25 OGS_FATAL(
"An invalid norm type has been passed");
41 if (norm_x > std::numeric_limits<double>::epsilon())
43 return norm_dx / norm_x;
47 if (norm_dx < std::numeric_limits<double>::epsilon())
53 return norm_dx / std::numeric_limits<double>::epsilon();
Global vector based on Eigen vector.
double computeRelativeChangeFromPreviousTimestep(GlobalVector const &x, GlobalVector const &x_old, MathLib::VecNormType norm_type)
std::unique_ptr< GlobalVector > _dx
Used to store .
double norm(MatrixOrVector const &x, MathLib::VecNormType type)
void axpy(PETScVector &y, double const a, PETScVector const &x)
void copy(PETScVector const &x, PETScVector &y)
VecNormType
Norm type. Not declared as class type in order to use the members as integers.