[petsc-dev] Problem mixing sieve and TS

[Jed Brown]

On Wed, Oct 26, 2011 at 12:45, Blaise Bourdin <bourdin at lsu.edu> wrote:

> Any suggestion? FWIW, in this test, both the IJacobian matrix and RHS
> Vector are constant, but this will not be the case once in the future.
>
> Something else bugging me is that the Jacobian function does not seem to be
> ever called (even in ex22 from src/ts/examples/tutorial).
>

How are you coming to that conclusion? The default method for ex22 is a
method with an implicit part. If you don't set any type at all, then the
default method is forward Euler (we should probably change this, maybe to a
Rosenbrock-W), which does not evaluate a Jacobian, but that code sets the
default type to ARKIMEX which does evaluate the IJacobian.


> Is that expected? Is there a written description (or link to the
> literature) for each of the ODE solvers implemented in TS?
>

The man page for each method should have a citation to a paper that
introduced it or discussed it. Some of the methods in there were designed by
Emil and have not been published yet.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111026/8df6b7cf/attachment.html>