On Tue, Jan 17, 2012 at 2:51 PM, Xuefei (Rebecca) Yuan <xyuan at lbl.gov>wrote:
> Hello all, > > I have a piece of code that need to store the RHS vector b of the ksp > object via vecview. > > However, the routine is like: > > > ierr = DMMGSetSNESLocal(dmmg, FormFunctionLocal, > FormJacobianLocal,0,0);CHKERRQ(ierr); > ierr = DMMGSolve(dmmg);CHKERRQ(ierr); > > The rhs evaluation is in FormFunctionLocal(), how could store the rhs > vector b? > 1) In your nonlinear problem, what do you mean by b? 2) Don't use DMMG, use petsc-dev with regular SNES and SNESSetDM() Matt > Thanks very much! > > Cheers, > > Rebecca > > > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120117/abbd3946/attachment.htm>
