[petsc-users] petsc4py

[Николай Кутузов]

---->  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>