From 54bf292307e0297301aadbb115d4e7878b4a45ed Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jakub=20Klinkovsk=C3=BD?= <klinkjak@fjfi.cvut.cz>
Date: Sat, 14 Nov 2015 20:10:12 +0100
Subject: [PATCH] Fixed bug in GMRES
---
src/solvers/linear/krylov/tnlGMRESSolver_impl.h | 13 ++++++-------
1 file changed, 6 insertions(+), 7 deletions(-)
diff --git a/src/solvers/linear/krylov/tnlGMRESSolver_impl.h b/src/solvers/linear/krylov/tnlGMRESSolver_impl.h
index 8d738b1b74..770cb26b31 100644
--- a/src/solvers/linear/krylov/tnlGMRESSolver_impl.h
+++ b/src/solvers/linear/krylov/tnlGMRESSolver_impl.h
@@ -251,15 +251,14 @@ bool tnlGMRESSolver< Matrix, Preconditioner > :: solve( const Vector& b, Vector&
sn[ i ] );
this->setResidue( fabs( s[ i + 1 ] ) / normb );
- this->refreshSolverMonitor();
-
- /*if( this->getResidue() < this->getConvergenceResidue() )
- {
+ if( this->nextIteration() ) {
+ this->refreshSolverMonitor();
+ }
+ else {
update( i, m, _H, _s, _v, x );
- return true;
+ this->refreshSolverMonitor( true );
+ return this->checkConvergence();
}
- if( ! this->nextIteration() )
- return false;*/
}
update( m - 1, m, _H, _s, _v, x );
--
GitLab