| 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 |