<div dir="ltr">Does this look good:<div><br></div><div><div>Norm of matrix ratio 2.19793e-09, difference 2.14739e-07 (user-defined state)</div><div>Norm of matrix ratio 1.10356e-08, difference 1.07917e-06 (constant state -1.0)</div><div>Norm of matrix ratio 1.15084e-08, difference 1.12541e-06 (constant state 1.0)</div></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 17, 2017 at 12:44 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.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">I see, I had a tiny time step that make the Jacobian big and the errors small.</div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 17, 2017 at 12:30 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.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"><br><div class="gmail_extra"><br><div class="gmail_quote"><span>On Thu, Aug 17, 2017 at 12:22 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="m_8155288774901714469m_6859141497539562307gmail-">On Wed, Aug 16, 2017 at 10:14 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I just want to check the interpretation of snes_type test. Is this pretty conclusively good?</div></blockquote><div><br></div></span><div>It looks right, but the entires in the Jacobian are crazy. You should non-dimensionalize first.</div></div></div></div></blockquote><div><br></div></span><div>This is non-dimensionalize and it is 2D FE 2nd order equation so everything should be O(1), that is, there is no 1/h like term. And the domain is (1) and the grid is small so h is O(1) anyway.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br></div><div>  Matt</div><span class="m_8155288774901714469m_6859141497539562307gmail-"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Testing hand-coded Jacobian, if the ratio is</div><div>O(1.e-8), the hand-coded Jacobian is probably correct.</div><div>Run with -snes_test_display to show difference</div><div>of hand-coded and finite difference Jacobian.</div><div>Norm of matrix ratio 1.38299e-16, difference 1.11738 (user-defined state)</div><div>Norm of matrix ratio 9.50529e-09, difference 7.67972e+07 (constant state -1.0)</div><div>Norm of matrix ratio 1.07438e-08, difference 8.68036e+07 (constant state 1.0)</div></div><div><br></div></div>
</blockquote></span></div><span class="m_8155288774901714469m_6859141497539562307gmail-HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_8155288774901714469m_6859141497539562307gmail-m_-717686420364861377gmail_signature"><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">http://www.caam.rice.edu/~mk51<wbr>/</a><br></div></div></div>
</font></span></div></div>
</blockquote></span></div><br></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>