<div dir="ltr">Dear all<div><br></div><div>I am solving 3 coupled first order PDE in 2-D with SNES. The stencil is (3+1+3) in both directions due to weno scheme. I want to use matrix-free gmres. For preconditioner, I will use first order scheme with stencil (1+1+1) in both directions. I am using DMDA.</div><div><br></div><div>Two things that I am looking for are:</div><div><br></div><div>How to setup the matrix ?</div><div><br></div><div>Since I have 3 variables at each grid point, how to fill the jacobian matrix.</div><div><br></div><div>Is there an existing example that comes close to this situation. My search was not fruitful.<br></div><div><br></div><div>Thank you</div><div>praveen</div></div>