<div dir="ltr">Thank you Peter.<div><br></div><div style>Can I simply ignore this step:</div><div style><p style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,38)">TSSetRHSJacobian(ts, J, J, SNESDefaultComputeJacobian, PETSC_NULL);</p>
<p style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,38)"><br></p><p style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,38)"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">And get the SNES context from TS, then set the Jacobian option this way:</span><br>
</p><p style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,38)">SNES snes;<br></p><p style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,38)">TSGetSNES(ts, &snes);</p><p style="margin:0px">
</p><p style="color:rgb(0,132,38);font-family:Menlo;font-size:11px;margin:0px">SNESSetJacobian(snes, J, J, SNESDefaultComputeJacobian, PETSC_NULL);</p><p style="color:rgb(0,132,38);font-family:Menlo;font-size:11px;margin:0px">
<br></p><p style="margin:0px"><font face="arial, sans-serif">(Guess my last email didn't work well)</font></p><p style="margin:0px"><font face="arial, sans-serif"><br></font></p><p style="margin:0px"><font face="arial, sans-serif">Ling</font></p>
<p></p></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 22, 2013 at 3:17 PM, Peter Brune <span dir="ltr"><<a href="mailto:prbrune@gmail.com" target="_blank">prbrune@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>This may be done with the command-line option:<br><br></div>-snes_fd<span class="HOEnZb"><font color="#888888"><br>
<br></font></span></div><span class="HOEnZb"><font color="#888888">- Peter<br></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 22, 2013 at 4:03 PM, Zou (Non-US), Ling <span dir="ltr"><<a href="mailto:ling.zou@inl.gov" target="_blank">ling.zou@inl.gov</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi All,<br><br>In 'SNESSetJacobian' I could pass 'SNESDefaultComputeJacobian' to let PETSc do the finite differencing for me to calculate Jacobian, i.e.,<br>

<br>SNESSetJacobian(snes, J, J, SNESDefaultComputeJacobian, PETSC_NULL);<br>
<br><br>In TSSetRHSJacobian, this options seems to be invalid, i.e., it requires a function to evaluate the Jacobian, i.e.,<br><br><br>TSSetRHSJacobian(ts, J, J, FormRHSJacobian, PETSC_NULL);<br><br><br>Is it possible that TS could also use PETSc finite differencing Jacobian as SNES does?<br>


<br><br>Best,<br><br><br>Ling<br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>