[petsc-users] petsc4py

[Jed Brown]

On Thu, Sep 30, 2010 at 19:35, Jed Brown <jed at 59a2.org> wrote:
>> ts.setIFunction(ode.evalFunction, f) ---> ts.setIFunction(ode.evalFunction)
>
> This won't work, you need to provide a residual vector, otherwise
> PETSc does not know how large your vector is

Actually, you have to set the initial vector somewhere, and the
residual vector can be obtained by duplicating it.  PETSc works this
way in C.  Lisandro, why does TS.setIFunction require the residual
vector to be provided?  If you agree that it's not necessary, please
apply the attached patch.

Jed
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jed-setifunction.patch
Type: text/x-patch
Size: 1172 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100930/efa97b15/attachment.bin>