<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jan 16, 2014 at 6:28 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
        I am solving a set of equations with SNES<br>
F1 (x1,x2,x3) = 0<br>
F2 (x1,x2,x3) = 0<br>
F3 (x1,x2,x3) = 0<br>
<br>
x1 and x2 are always > 0, how do i enforce that during the linesearch ?<br></blockquote><div><br></div><div>There are two options:</div><div><br></div><div>  a) Make your own line search that projects the step</div><div>
<br></div><div>      This is what we have traditionally done. It can destroy convergence, but its simple.</div><div><br></div><div>  b) Use SNESVI and set bounds on x1 and x2</div><div><br></div><div>      This is supposed to work, but has hardly been tested at all. When we get some time</div>
<div>      we plan on beating on this. If you feel like being a guinea pig, do this.</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">

The system of equations is shown on page 1 of pdf here<br>
<a href="http://dunham.ee.washington.edu/ee531/notes/SemiRev.pdf" target="_blank">http://dunham.ee.washington.edu/ee531/notes/SemiRev.pdf</a><br>
<br>
F1 = equation 1<br>
F2 = equation 2<br>
F3 = equation 5<br>
<br>
x1 = n, X2=p and X3 = psi<br>
<br>
Thanks<br>
<span class="HOEnZb"><font color="#888888">Reddy<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
</div></div>