<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>Hello Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class=""><div></div></span><div>No. You want 2 DAs, 1 for evaluating residuals with s = 3, and one for making a Jacobian with s = 1. You give these</div><div>two DAs to</div><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/DMDASNESSetFunctionLocal.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>petsc-current/docs/<wbr>manualpages/SNES/<wbr>DMDASNESSetFunctionLocal.html</a><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/DMDASNESSetJacobianLocal.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>petsc-current/docs/<wbr>manualpages/SNES/<wbr>DMDASNESSetJacobianLocal.html</a><br></div><div><br></div><div>and I think it should work as you want. Its too bad we cannot handle this case with SNESSetDM() since it only takes</div><div>a single DM and there is no way to specialize the Jacobian.</div><div><br></div></div></div></div></blockquote><div><br></div><div>ok, I was on wrong track trying to use SNESSetFunction and SNESSetJacobian which wouldnt have worked.</div><div><br></div><div>I still have some doubt. If</div><div><br></div><div>da0 --> for residual</div><div>da1 --> for jacobian </div><div><br></div><div>I still need to call SNESetDM. So here I pass da0 ? But then how will snes know about da1 ?</div><div><br></div><div>Thanks</div><div>praveen</div></div></div></div>