> > > You can use SNESSetJacobian, but DMSNESSetJacobianLocal is more > convenient in my opinion. > > > I still have some doubt. If > > > > da0 --> for residual > > da1 --> for jacobian > > > > I still need to call SNESetDM. So here I pass da0 ? But then how will > snes > > know about da1 ? > > This should work. > > DMCreateMatrix(da1, &Jpre); > SNESSetDM(snes, da0); > SNESSetJacobian(snes, NULL, Jpre, NULL, NULL); > DMDASNESSetJacobianLocal(da0, func, ctx); > > > In the last line above, should it not be da1 ?
Thanks a lot for the examples. Best praveen
