<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Dec 19, 2013 at 4:26 PM, Que Cat <span dir="ltr"><<a href="mailto:quecat001@gmail.com" target="_blank">quecat001@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">Hello,<div><br></div><div>I use SNES to solve the non linear coupled equations. I have checked the hand-coded jacobian with -snes_type test and received the result:</div>
<div><div>Norm of matrix ratio 1.39606e-23 difference 1.53372e-14 (user-defined state)</div>
<div>Norm of matrix ratio 3.03157e-10 difference 0.33305 (constant state -1.0)</div><div>Norm of matrix ratio 3.03157e-10 difference 0.33305 (constant state 1.0)</div></div><div>It said that "  if the ratio is O(1.e-8), the hand-coded Jacobian is probably correct." What is the second paramenter "difference 1.53372e-14 ...."? Does it play any role in checking the accuracy of hand-coded jacobian?</div>
</div></blockquote><div><br></div><div>It appears that your Jacobian is correct, or at least has an error that is much smaller than most of the elements.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>The SNES solver for failed with the reason -6 "<span style>the line search failed". I tried with all other line search method stated in the PETSC manual, but they all failed. </span></div></div>
</blockquote><div><br></div><div>Do you know that the problem has a solution. Its possible that there is no real solution.</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"><div><span style>I tried to figure out what is wrong with my code? Could you please give me any recommendation to direct me how to check my code? I appreciate your help.</span></div><span class="HOEnZb"><font color="#888888">
<div><span style><br></span></div><div><span style>QUe</span></div><div><br></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>