Skip to content
Snippets Groups Projects

Fix scaling of the Jacobian with linear solver PaStiX and nonlinear solver NOX

Merged Reinhard Resch requested to merge fix-pastix-scale into develop
4 files
+ 18
10
Compare changes
  • Side-by-side
  • Inline
Files
4
@@ -614,6 +614,7 @@ bool LinSol::SetTolerance(doublereal dToleranceRes)
{
switch (currSolver) {
case LinSol::AZTECOO_SOLVER:
case LinSol::PASTIX_SOLVER:
dTolRes = dToleranceRes;
break;
default:
@@ -941,12 +942,12 @@ LinSol::GetSolutionManager(integer iNLD,
case LinSol::SOLVER_FLAGS_ALLOWS_GRAD: {
SAFENEWWITHCONSTRUCTOR(pCurrSM,
PastixSolutionManager<SpGradientSparseMatrixHandler>,
PastixSolutionManager<SpGradientSparseMatrixHandler>(iNLD, nThreads, iMaxIter, scale, solverFlags, dLowRankCompressTol, dLowRankCompressMinRatio, iVerbose));
PastixSolutionManager<SpGradientSparseMatrixHandler>(iNLD, nThreads, iMaxIter, dTolRes, scale, solverFlags, dLowRankCompressTol, dLowRankCompressMinRatio, iVerbose));
} break;
default:
SAFENEWWITHCONSTRUCTOR(pCurrSM,
PastixSolutionManager<SpMapMatrixHandler>,
PastixSolutionManager<SpMapMatrixHandler>(iNLD, nThreads, iMaxIter, scale, solverFlags, dLowRankCompressTol, dLowRankCompressMinRatio, iVerbose));
PastixSolutionManager<SpMapMatrixHandler>(iNLD, nThreads, iMaxIter, dTolRes, scale, solverFlags, dLowRankCompressTol, dLowRankCompressMinRatio, iVerbose));
}
} break;
#else /* !USE_PASTIX */
Loading