SNESGetKSP( KSPSetType(ksp,KSPPREONLY) KSPGetPC( PCSetType(pc,PCSHELL); PCShellSetApply(pc,YourSolver);
On Mon, 29 Oct 2007, Lisandro Dalcin wrote: > Some time ago I made a request for adding SNESSetLinearSolve(), in > order to let users set a custom linear solve routine to SNES. I was > not clear to me if this idea was finally accepted, so I want to ask > again. > > I really needs this feature in the near future, I I would like this to > go in the next petsc release. My use case is rather simple: I'm > helping a coworker to use SNES for solving a nonlinear optimization > problem related to mesh movement (that is, move nodes retaining > element connectivities). This problem uses a regularized functional, > wich introduces a global degree of freedom coupled with all others > dofs. To avoid assembling a sparse parallel AIJ matrix with a dense > column and row, we endup needing to solve two linear systems with the > same matrix for actually solving the Hessian of the regularized > functional. > > > >
