<div class="gmail_quote">On Wed, Oct 26, 2011 at 12:45, Blaise Bourdin <span dir="ltr"><<a href="mailto:bourdin@lsu.edu">bourdin@lsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":1p5">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.<br>
<br>
Something else bugging me is that the Jacobian function does not seem to be ever called (even in ex22 from src/ts/examples/tutorial).</div></blockquote><div><br></div><div>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.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":1p5"> Is that expected? Is there a written description (or link to the literature) for each of the ODE solvers implemented in TS?</div>
</blockquote></div><br><div>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.</div>