[petsc-users] SNES realted

[Barry Smith]

On Jan 16, 2014, at 6:28 PM, Dharmendar Reddy <dharmareddy84 at gmail.com> wrote:

> Hello,
>        I am solving a set of equations with SNES
> F1 (x1,x2,x3) = 0
> F2 (x1,x2,x3) = 0
> F3 (x1,x2,x3) = 0
> 
> x1 and x2 are always > 0, how do i enforce that during the line search ?

  You don’t 

   This is an example of a variational inequality.   You use the SNESType of vinewtonrsls  to solve it and call SNESVISetVariableBounds() to set the bounds on the variables.

    Good luck,

    Barry



> 
> The system of equations is shown on page 1 of pdf here
> http://dunham.ee.washington.edu/ee531/notes/SemiRev.pdf
> 
> F1 = equation 1
> F2 = equation 2
> F3 = equation 5
> 
> x1 = n, X2=p and X3 = psi
> 
> Thanks
> Reddy