On Wed, Dec 14, 2011 at 4:07 PM, Philip, Bobby <span dir="ltr"><<a href="mailto:philipb@ornl.gov">philipb@ornl.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">
Hi:<br>
<br>
When doing the following operation in PETSc using FD approximations to a  Jacobian vector product<br>
<br>
J'(u)*a = [J(u+h*a) - J(u)]/h<br></blockquote><div><br></div><div>Usually, we would put in F where you have J.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

is there a way for the user to specify that the entries of u+h*a cannot be negative for example?<br>
Browsing the PETSc documentation it was not obvious whether I could do this.<br></blockquote><div><br></div><div>You can do whatever you want in your F function, but you should really investigate SNESVI if</div><div>you have bounds on your variables.</div>
<div><br></div><div>  Thanks,</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,<br>
Bobby</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>