On Mon, Jun 20, 2011 at 9:40 PM, Tian(ICT) <span dir="ltr"><<a href="mailto:rongtian@ncic.ac.cn">rongtian@ncic.ac.cn</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
Good morning/evening!<br>
I have this issue for quite a while. Using MSVC with Petsc, when debuging the code, the FormJacobian was never called, instead calling "MatMFFDComputeJacobian(SNES snes,Vec x,Mat *jac,Mat *B,MatStructure *flag,void *dummy)", though I donot use matrix free option.<br>
<br>
However, when runing the same thing in command line (still a debug compilation), I can see FormJacobian was called correctly.<br>
<br>
Any hint?<br></blockquote><div><br></div><div>My guess is that you are getting the matrix-free option from somewhere (maybe a .petscrc file?). You can use -options_dump to see</div><div>everything passed in. If not, look at SNESSetFromOptions() in the debugger, and see why it is choosing MF.</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;">
Rong<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <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<br>