how should i use that patch?<br><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;">
<div class="im">On Thu, Sep 30, 2010 at 19:35, Jed Brown <<a href="mailto:jed@59a2.org">jed@59a2.org</a>> wrote:<br>
>> ts.setIFunction(ode.evalFunction, f) ---> ts.setIFunction(ode.evalFunction)<br>
><br>
> This won't work, you need to provide a residual vector, otherwise<br>
> PETSc does not know how large your vector is<br>
<br>
</div>Actually, you have to set the initial vector somewhere, and the<br>
residual vector can be obtained by duplicating it. PETSc works this<br>
way in C. Lisandro, why does TS.setIFunction require the residual<br>
vector to be provided? If you agree that it's not necessary, please<br>
apply the attached patch.<br>
<font color="#888888"><br>
Jed<br>
</font></blockquote></div><br>