Que Cat <[email protected]> writes: > Hello Petsc-users, > > I have to solve a problem in which the jacobian is non-smooth ( incoporates > the heaviside step function behavior from the derivative of F(x) ). Could > you recommend any type of nonlinear solver in Petsc which handles this > issue ?
You can reformulate as a variational inequality and use SNESVI. Be warned that it's not for the faint of heart and it's less mature than the "smooth" solvers. If this nonsmoothness is absolutely fundamental to your physics, like contact or condensation, then you should expect to invest significant effort in finding the best formulation. If the nonsmoothness is not essential, you should put that effort into reformulating to avoid or lessen the nonsmoothness.
pgplr4UOfFmG5.pgp
Description: PGP signature
