<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Oct 27, 2014 at 12:52 PM, Sharp Stone <span dir="ltr"><<a href="mailto:thronesf@gmail.com" target="_blank">thronesf@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>Dear All,<br><br></div>I'm using the TS theta solver for backward Euler method. But I found that the TSSetIJacobian() routine has never been called. I did the code just following what the example does: src/ts/examples/tutorials/ex17. I can run the example without errors. But when I run my code, it won't be as good as the example as stated above. I also attached the piece of my setups to this email. <br><br></div>How could I fix this problem? Thank you in advance!<br></div></blockquote><div><br></div><div>I cannot understand the paragraph above.</div><div><br></div><div>  1) TS ex17 does call TSSetIJacobian(), unless you choose a matrix-free version. You seem to imply above that it does not.</div><div><br></div><div>  2) Do you mean that you do not call it? You do not have too, since you can use the finite difference approximations offered by SNES</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>    TSCreate(PETSC_COMM_WORLD, &ts);<br>    TSSetProblemType(ts, TS_NONLINEAR);<br>    TSSetType(ts, TSTHETA);<br>    TSThetaSetTheta(ts, 1.0);<br>    TSSetIFunction(ts, NULL, FormIFunction, &appctx);<br>    DMSetMatType(da, MATAIJ);<br>    DMCreateMatrix(da, &J);<br>    TSSetDM(ts, da);<br>    TSSetIJacobian(ts, J, J, FormIJacobian, &appctx);<br>    TSSetDuration(ts, timeStepsMax, Tend); <br>    // Evaluate initial conditions<br>    FormInitialSolution(&appctx);<br>    TSSetSolution(ts, solution); <br>    dt = 1.0;<br>    TSSetInitialTimeStep(ts, 0.0, dt); <br>    TSSolve(ts, solution);<span class="HOEnZb"><font color="#888888"><br clear="all"><div><div><div><br>-- <br><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div>
</div></div></div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>