---->š petsc-3.1-p4 and petsc4py-1.1.1<br><br>I did only this change in <a href="http://code.google.com/p/petsc4py/source/browse/demo/ode/rober.py">http://code.google.com/p/petsc4py/source/browse/demo/ode/rober.py</a> ---<br>
--- (42) <br><pre class="prettyprint lang-py"><span class="pln">ts</span><span class="pun">.</span><span class="pln">setIFunction</span><span class="pun">(</span><span class="pln">ode</span><span class="pun">.</span><span class="pln">evalFunction</span><span class="pun">,</span><span class="pln"> f</span><span class="pun">) ---> </span><span class="pln">ts</span><span class="pun">.</span><span class="pln">setIFunction</span><span class="pun">(</span><span class="pln">ode</span><span class="pun">.</span><span class="pln">evalFunction</span><span class="pun">)</span><span class="pln"></span><span class="pun"></span>
</pre><br><div class="gmail_quote">2010/9/30 Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<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></p><blockquote type="cite"><div class="im">On Sep 30, 2010 4:25 AM, "Lisandro Dalcin" <<a href="mailto:dalcinl@gmail.com" target="_blank">dalcinl@gmail.com</a>> wrote:<br><br><p><font color="#500050">2010/9/29 îÉËĎĚÁĘ ëŐÔŐÚĎ× <<a href="mailto:kutuzovnp@gmail.com" target="_blank">kutuzovnp@gmail.com</a>>:</font></p>
</div><p><font color="#500050"><div class="im">> Where can i get it? I've tried this one --<br></div>> <a href="http://code.google.com/p/petsc4py/source/browse/demo/." target="_blank">http://code.google.com/p/petsc4py/source/browse/demo/.</a>..</font></p>
<div><div></div><div class="h5">
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></div></div></blockquote>
</blockquote></div><br>