<p>And I get the same output as Lisandro if I comment the ts.setIJaconian line (because that Jacobian is not being used).</p>
<p>Jed</p>
<p><blockquote type="cite">On Sep 30, 2010 4:25 AM, "Lisandro Dalcin" <<a href="mailto:dalcinl@gmail.com">dalcinl@gmail.com</a>> wrote:<br><br><p><font color="#500050">2010/9/29 Николай Кутузов <<a href="mailto:kutuzovnp@gmail.com">kutuzovnp@gmail.com</a>>:</font></p>
<p><font color="#500050">> Where can i get it? I've tried this one --<br>> <a href="http://code.google.com/p/petsc4py/source/browse/demo/.">http://code.google.com/p/petsc4py/source/browse/demo/.</a>..</font></p>
This is from branch release-1.1, a new release 1.1.2 is planned for<br>
the weekend. As you can see from the output at the very end, SNES is<br>
using matrix free. Note that I've not changed he code code, just<br>
passed -snes_mf<br>
<br>
What petsc4py and PETSc versions are you using?<br>
<br>
[dalcinl@trantor petsc4py-release-1.1]$ python demo/ode/rober.py<br>
-ts_view -snes_mf<br>
TS Object:<br>
type: theta<br>
Theta=0.5<br>
Extrapolation=no<br>
maximum steps=100<br>
maximum time=1e+30<br>
total number of nonlinear solver iterations=203<br>
total number of linear solver iterations=311<br>
SNES Object:<br>
type: ls<br>
line search variant: SNESLineSearchCubic<br>
alpha=0.0001, maxstep=1e+08, minlambda=1e-12<br>
maximum iterations=50, maximum function evaluations=10000<br>
tolerances: relative=1e-08, absolute=1e-50, solution=1e-08<br>
total number of linear solver iterations=3<br>
total number of function evaluations=8<br>
KSP Object:<br>
type: gmres<br>
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt<br>
Orthogonalization with no iterative refinement<br>
GMRES: happy breakdown tolerance 1e-30<br>
maximum iterations=10000, initial guess is zero<br>
tolerances: relative=1e-05, absolute=1e-50, divergence=10000<br>
left preconditioning<br>
using PRECONDITIONED norm type for convergence test<br>
PC Object:<br>
type: none<br>
linear system matrix = precond matrix:<br>
Matrix Object:<br>
type=mffd, rows=3, cols=3<br>
matrix-free approximation:<br>
err=1e-07 (relative error in function evaluation)<br>
Using wp compute h routine<br>
Computes normA<br>
Does not compute normU<br>
<font color="#888888"><br>
<br>
--<br>
Lisandro Dalcin<br>
---------------<br>
CIMEC (INTEC/CONICET-UNL)<br>
Predio CONICET-Santa Fe<br>
Colectora RN 168 Km 472, Paraje El Pozo<br>
Tel: +54-342-4511594 (ext 1011)<br>
Tel/Fax: +54-342-4511169<br>
</font></blockquote></p>