<div dir="ltr">On Thu, May 2, 2013 at 7:04 AM, Abdul-Lateef Haji-Ali <span dir="ltr"><<a href="mailto:abdullateef.hajiali@kaust.edu.sa" target="_blank">abdullateef.hajiali@kaust.edu.sa</a>></span> wrote:<br><div class="gmail_extra">
<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-family:arial,helvetica,sans-serif">Greetings,<br><br></div><div style="font-family:arial,helvetica,sans-serif">
I am trying to solve a system of 4 non-linear equations to find 4 unknowns.<br>

The equations are slightly complicated  and unstable for slightly large values of the unknowns.<br>In most of the cases the correct solution is found. However sometimes the values that are fed to my residual function are big and nan/inf values are produced.<br>
</div></div></blockquote><div><br></div><div style>It could be overflow, so scaling the equations may help. However, it could also be the</div><div style>the math is poorly organized so that you create huge intermediate values (this seems</div>
<div style>likely since all your other values are quite small).</div><div style><br></div><div style>   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 style="font-family:arial,helvetica,sans-serif">This is the output of my residual function (x are the parameters and f is the computed residual)<br>x= [1, 1, 2, 4]; f= [-1.32993, 6.46116, 0.490443, -4.47854];<br>
x= [2.32993, -5.46116, 1.50956, 8.47854]; f= [-1.77078, -133.872, 4.26217, 71.9728];<br>

x= [1.06488, 0.684771, 1.97607, 4.2185]; f= [-0.431978, 2.79625, 0.190123, 0.0503181];<br>x= [4.10071, 128.411, -2.75261, -63.4942]; f= [62.0141, 2548.28, -95.0523, -1354.74];<br>x= [1.74294, 0.960827, 1.63558, 4.54888]; f= [19.6641, 4.52508, -9.48464, 7.44219];<br>


x= [-57.9134, -2419.86, 92.2997, 1291.24]; f= [inf, -inf, inf, -inf];<br></div><div style="font-family:arial,helvetica,sans-serif">This was using ngmres, I noticed similar behaviour with other types.<br>

</div><div style="font-family:arial,helvetica,sans-serif"><br>What can I do to improve the solver?<br></div><div style="font-family:arial,helvetica,sans-serif"></div>N<div style="font-family:arial,helvetica,sans-serif;display:inline">


ote that my experience with petsc in new and limited.<br><br>Thank you,</div><span class="HOEnZb"><font color="#888888"><br><div>-- <br><div dir="ltr">Abdul Lateef Haji Ali<br></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>