<div class="gmail_quote">On Tue, Nov 8, 2011 at 04:14, Konstantinos Kontzialis <span dir="ltr">&lt;<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":18k">   ierr = TSSetSolution(sys.ts, sys.gsv);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = TSSetIFunction(sys.ts, PETSC_NULL, base_residual_implicit, &amp;sys);<br>
    CHKERRQ(ierr);<br></div></blockquote><div><br></div><div>Provide a residual vector here (instead of PETSC_NULL). Since you provided a state Vec in TSSetSolution(), it is possible to create a Vec for the residual. I&#39;ll add that logic and improve the error message, but for now, just pass in the residual Vec.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":18k">
<br>
    ierr = TSGetSNES(sys.ts, &amp;sys.snes);<br>
    CHKERRQ(ierr);<br>
<br>
    ierr = MatCreateSNESMF(sys.snes, &amp;sys.J);<br>
    CHKERRQ(ierr);</div></blockquote></div><br>