#1202 Termination criterion in RestartedGMRes wrong
Metadata
Property | Value |
---|---|
Reported by | Uli Sack (usack@math.fu-berlin.de) |
Reported at | Oct 24, 2012 14:54 |
Type | Bug Report |
Version | Git (pre2.4) [autotools] |
Operating System | Unspecified / All |
Last edited by | Markus Blatt (markus@dr-blatt.de) |
Last edited at | Aug 12, 2013 13:28 |
Closed by | Markus Blatt (markus@dr-blatt.de) |
Closed at | Aug 12, 2013 13:28 |
Closed in version | Unknown |
Resolution | Fixed |
Comment | fixed with jurgis' patch in rev d3fe65f |
Description
AFAIK the residual-norm available during GMRes-iteration is that of the preconditioned residual, which for the convergence check should be compared to the preconditioned initial residual if anything. In the istl implementation however it is compared to the UNpreconditioned initial residual.
Therefore I propose the attached patch. It sets norm_0 to norm after that initial convergence check right after the initial preconditioning (couldn't we just kick that one out altogether?). To avoid confusing output we must move the output around a little.
Greets, Uli