On Wed, Sep 5, 2012 at 12:24 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><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear All,<br>
<br>
I am trying to use the option '-snes_type test' to test my coded<br>
Jacobian. I tested with different snes options and it gives me<br>
different answers. I wonder if someone could give me a hint what is<br>
wrong with my settings.<br>
<br>
The command line looks like this:<br>
./my-code-opt -i test.i -snes_type test -snes_test_display<br>
<br>
1), when using 'petsc_option = -snes' in my input file, it says the<br>
Finite difference Jacobian is very different than the Hand-coded<br>
Jacobian<br></blockquote><div><br></div><div>This means your hand-coded routine is likely wrong.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2), when using 'petsc_option = -snes_fd' in my input file, it says the<br>
Finite difference Jacobian is idential to the Hand-coded Jacobian<br></blockquote><div><br></div><div>snes_fd replaces your hand-coded routine with our FD routine, so it</div><div>is of course the same as our FD routine.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
3), when using 'petsc_option = -snes_mf_operator', it gives error messages like:<br>
"Invalid argument! Cannot test with alternative preconditioner!"<br></blockquote><div><br></div><div>This is inappropriate for testing.</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">

Thanks in advance.<br>
<span class="HOEnZb"><font color="#888888"><br>
Ling<br>
</font></span></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<br>