CFD/hypre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed a problem with setting absolute tolerances.

Browse Source
This commit is contained in:
falgout 2013-08-01 15:22:21 +00:00
parent a239552b48
commit 0b0076ca35
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
FEI_mv/fei-hypre/HYPRE_LinSysCore.cxx
View File

@ -4786,13 +4786,12 @@ int HYPRE_LinSysCore::launchSolver(int& solveStatus, int &iterations)
HYPRE_PCGSetRecomputeResidual(HYSolver_, pcgRecomputeRes_); HYPRE_PCGSetRecomputeResidual(HYSolver_, pcgRecomputeRes_);
if ( normAbsRel_ == 0 ) if ( normAbsRel_ == 0 )
{ {
HYPRE_PCGSetStopCrit(HYSolver_,0);
HYPRE_PCGSetTol(HYSolver_, tolerance_); HYPRE_PCGSetTol(HYSolver_, tolerance_);
} }
else else
{ {
HYPRE_PCGSetStopCrit(HYSolver_,1);
HYPRE_PCGSetAbsoluteTol(HYSolver_, tolerance_); HYPRE_PCGSetAbsoluteTol(HYSolver_, tolerance_);
HYPRE_PCGSetTol(HYSolver_, 0.0);
} }
if ( (HYOutputLevel_ & HYFEI_SPECIALMASK) >= 1 ) if ( (HYOutputLevel_ & HYFEI_SPECIALMASK) >= 1 )
{ {
@ -5029,13 +5028,12 @@ int HYPRE_LinSysCore::launchSolver(int& solveStatus, int &iterations)
HYPRE_ParCSRGMRESSetMaxIter(HYSolver_, maxIterations_); HYPRE_ParCSRGMRESSetMaxIter(HYSolver_, maxIterations_);
if ( normAbsRel_ == 0 ) if ( normAbsRel_ == 0 )
{ {
HYPRE_GMRESSetStopCrit(HYSolver_,0);
HYPRE_GMRESSetTol(HYSolver_, tolerance_); HYPRE_GMRESSetTol(HYSolver_, tolerance_);
} }
else else
{ {
HYPRE_GMRESSetStopCrit(HYSolver_,1);
HYPRE_GMRESSetAbsoluteTol(HYSolver_, tolerance_); HYPRE_GMRESSetAbsoluteTol(HYSolver_, tolerance_);
HYPRE_GMRESSetTol(HYSolver_, 0.0);
} }
if ( (HYOutputLevel_ & HYFEI_SPECIALMASK) >= 1 ) if ( (HYOutputLevel_ & HYFEI_SPECIALMASK) >= 1 )
{ {