[petsc-dev] controlling vector values while doing matrix free operations

Philip, Bobby philipb at ornl.gov
Wed Dec 14 16:07:03 CST 2011


When doing the following operation in PETSc using FD approximations to a  Jacobian vector product

J'(u)*a = [J(u+h*a) - J(u)]/h 

is there a way for the user to specify that the entries of u+h*a cannot be negative for example? 
Browsing the PETSc documentation it was not obvious whether I could do this.


More information about the petsc-dev mailing list