[petsc-dev] [petsc-maint #42811] TSPSEUDO and snes ex27

Jed Brown jed at 59A2.org
Tue Mar 9 04:05:20 CST 2010

Apparently I'm being excessively dense.  As long as the predictor for
the next step is just the old state, then we have

  \dot{x} = (x^{n+1}_{\text{predicted}}-x^n)/dt = (x^n - x^n)/dt = 0

so that


is in fact the steady state residual, in which case doing one Newton
step (with no line search, and no need to even check residual) is
equivalent to the algorithm described in the paper.  With respect to
ex27, this means that no code change is necessary to perform the
"correct" algorithm, and indeed the differences are minor.  Compared to
Figure 3.1, the DAE formulation converges in the same number of steps,
albeit with different residuals along the way, and the fully parabolized
form takes two extra steps (but both of these take many fewer linear
solves than when Newton is converged on the transient term at each
step).  I've put a note to this effect in the ex27 and TSPSEUDO


More information about the petsc-dev mailing list