Hello Matt
> No. You want 2 DAs, 1 for evaluating residuals with s = 3, and one for > making a Jacobian with s = 1. You give these > two DAs to > > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/ > DMDASNESSetFunctionLocal.html > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/ > DMDASNESSetJacobianLocal.html > > and I think it should work as you want. Its too bad we cannot handle this > case with SNESSetDM() since it only takes > a single DM and there is no way to specialize the Jacobian. > > ok, I was on wrong track trying to use SNESSetFunction and SNESSetJacobian which wouldnt have worked. 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 ? Thanks praveen
