----> petsc-3.1-p4 and petsc4py-1.1.1 I did only this change in http://code.google.com/p/petsc4py/source/browse/demo/ode/rober.py --- --- (42)
ts.setIFunction(ode.evalFunction, f) ---> ts.setIFunction(ode.evalFunction) 2010/9/30 Jed Brown <jed at 59a2.org> > And I get the same output as Lisandro if I comment the ts.setIJaconian line > (because that Jacobian is not being used). > > Jed > > On Sep 30, 2010 4:25 AM, "Lisandro Dalcin" <dalcinl at gmail.com> wrote: > > 2010/9/29 ??????? ??????? <kutuzovnp at gmail.com>: > > > Where can i get it? I've tried this one -- > > http://code.google.com/p/petsc4py/source/browse/demo/... > This is from branch release-1.1, a new release 1.1.2 is planned for > the weekend. As you can see from the output at the very end, SNES is > using matrix free. Note that I've not changed he code code, just > passed -snes_mf > > What petsc4py and PETSc versions are you using? > > [dalcinl at trantor petsc4py-release-1.1]$ python demo/ode/rober.py > -ts_view -snes_mf > TS Object: > type: theta > Theta=0.5 > Extrapolation=no > maximum steps=100 > maximum time=1e+30 > total number of nonlinear solver iterations=203 > total number of linear solver iterations=311 > SNES Object: > type: ls > line search variant: SNESLineSearchCubic > alpha=0.0001, maxstep=1e+08, minlambda=1e-12 > maximum iterations=50, maximum function evaluations=10000 > tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 > total number of linear solver iterations=3 > total number of function evaluations=8 > KSP Object: > type: gmres > GMRES: restart=30, using Classical (unmodified) Gram-Schmidt > Orthogonalization with no iterative refinement > GMRES: happy breakdown tolerance 1e-30 > maximum iterations=10000, initial guess is zero > tolerances: relative=1e-05, absolute=1e-50, divergence=10000 > left preconditioning > using PRECONDITIONED norm type for convergence test > PC Object: > type: none > linear system matrix = precond matrix: > Matrix Object: > type=mffd, rows=3, cols=3 > matrix-free approximation: > err=1e-07 (relative error in function evaluation) > Using wp compute h routine > Computes normA > Does not compute normU > > > -- > Lisandro Dalcin > --------------- > CIMEC (INTEC/CONICET-UNL) > Predio CONICET-Santa Fe > Colectora RN 168 Km 472, Paraje El Pozo > Tel: +54-342-4511594 (ext 1011) > Tel/Fax: +54-342-4511169 > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100930/70db7cab/attachment.htm>