Commit cb7e9ddb authored by Andreas Lauser's avatar Andreas Lauser
Browse files

box local jacobian: modify the base epsilon to 10^-10

and make sure that the base epsilon is still at least 4 orders of
magnitude larger than the machine precision...

git-svn-id: svn:// 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent a5f0c430
......@@ -267,9 +267,12 @@ public:
// resolution of the scalar type. E.g. for standard 64 bit
// floating point values, the resolution is about 10^-16 and
// the base epsilon is thus approximately 10^-8.
static const Scalar baseEps
= Dumux::geometricMean<Scalar>(std::numeric_limits<Scalar>::epsilon(), 1.0);
static const Scalar baseEps = 1e-10;
assert(std::numeric_limits<Scalar>::epsilon()*1e4 < baseEps);
// the epsilon value used for the numeric differentiation is
// now scaled by the absolute value of the primary variable...
Scalar pv = this->curVolVars_[scvIdx].primaryVar(pvIdx);
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment