254 |
IF ( krylovIter.EQ.SEAICEkrylovIterMax ) THEN |
IF ( krylovIter.EQ.SEAICEkrylovIterMax ) THEN |
255 |
krylovFails = krylovFails + 1 |
krylovFails = krylovFails + 1 |
256 |
ENDIF |
ENDIF |
257 |
C Set the stopping criterion for the Newton iteration |
C Set the stopping criterion for the Newton iteration and the |
258 |
IF ( newtonIter .EQ. 1 ) JFNKtol=JFNKgamma_nonlin*JFNKresidual |
C criterion for the transition from accurate to approximate FGMRES |
259 |
|
IF ( newtonIter .EQ. 1 ) THEN |
260 |
|
JFNKtol=JFNKgamma_nonlin*JFNKresidual |
261 |
|
IF ( JFNKres_tFac .NE. UNSET_RL ) |
262 |
|
& JFNKres_t = JFNKresidual * JFNKres_tFac |
263 |
|
ENDIF |
264 |
C Update linear solution vector and return to Newton iteration |
C Update linear solution vector and return to Newton iteration |
265 |
C Do a linesearch if necessary, and compute a new residual. |
C Do a linesearch if necessary, and compute a new residual. |
266 |
C Note that it should be possible to do the following operations |
C Note that it should be possible to do the following operations |