diff --git a/dumux/implicit/staggered/localjacobian.hh b/dumux/implicit/staggered/localjacobian.hh index 739b28f9c5dc7f23e1ab9f54a350f3ec39252515..12b3d34ee10c446932f638d4a884e7e082a32f5d 100644 --- a/dumux/implicit/staggered/localjacobian.hh +++ b/dumux/implicit/staggered/localjacobian.hh @@ -295,7 +295,6 @@ private: // update the global jacobian matrix with the current partial derivatives this->updateGlobalJacobian_(matrix[cellCenterIdx][cellCenterIdx], ccGlobalI_, globalJ, pvIdx, partialDeriv); -// printmatrix(std::cout,matrix[cellCenterIdx][cellCenterIdx], "part", ""); // restore the original volVars curVolVars = origVolVars; } @@ -328,6 +327,7 @@ private: for(auto pvIdx : PriVarIndices(faceIdx)) { PrimaryVariables priVars(CellCenterPrimaryVariables(0.0), FacePrimaryVariables(this->model_().curSol()[faceIdx][globalJ])); + const Scalar eps = numericEpsilon(priVars[pvIdx], cellCenterIdx, faceIdx); priVars[pvIdx] += eps; curFaceVars.update(priVars[faceIdx]);