<p>Note that in the code in ts/impls/implicit/sundials it says:</p>
<p>This uses its own nonlinear solver and krylov method so PETSc SNES and KSP options do not apply...</p>
<p>- Peter</p>
<div class="gmail_quote">On Jun 21, 2012 7:59 AM, "Geoff Oxberry" <<a href="mailto:goxberry@gmail.com">goxberry@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Running the following example from PETSC 3.3.0-dev (changeset: 23631:0e86ac5e4170)<br>
<br>
/path/to/petsc-dev/src/ts/examples/tutorials/ex8 -problem_type rober -snes_fd -ts_type sundials<br>
<br>
gives the following output<br>
<br>
steps 1000 (0 rejected, 0 SNES fails), ftime 744.845, nonlinits 3739, linits 3739<br>
WARNING! There are options you set that were not used!<br>
WARNING! could be spelling mistake, etc!<br>
Option left: name:-snes_fd no value<br>
<br>
Just to confirm, is it currently impossible to use a finite difference Jacobian matrix in concert with CVODE? If so, could this feature be implemented in a future release? I currently rely on Sundials to integrate stiff systems of ODEs, and for my application, it is impractical to derive an analytical Jacobian matrix. (It is an issue I've discussed both with Jed and Matt on another forum.)<br>

<br>
Cheers,<br>
<br>
Geoff</blockquote></div>