equation.cpp

Go to the documentation of this file.

#include "../../include/monolish_blas.hpp"
#include "../../include/monolish_equation.hpp"
#include "../internal/monolish_internal.hpp"
 
namespace monolish {
 
template <typename MATRIX, typename T>
T solver::solver<MATRIX, T>::get_residual(vector<T> &x) {
  switch (resid_method) {
  case 0:
    return blas::nrm2(x);
    break;
  default:
    throw std::runtime_error("error vector size is not same");
    break;
  }
}
 
template double
solver::solver<matrix::Dense<double>, double>::get_residual(vector<double> &x);
template float
solver::solver<matrix::Dense<float>, float>::get_residual(vector<float> &x);
 
template double
solver::solver<matrix::CRS<double>, double>::get_residual(vector<double> &x);
template float
 
solver::solver<matrix::CRS<float>, float>::get_residual(vector<float> &x);
template double
solver::solver<matrix::LinearOperator<double>, double>::get_residual(
    vector<double> &x);
template float
solver::solver<matrix::LinearOperator<float>, float>::get_residual(
    vector<float> &x);
} // namespace monolish