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? Thanks very much! Cheers, Rebecca