<div dir="ltr"><div><div><div>Hello everyone,<br><br></div>From what I read online, one could use a matrix-free method if using the conjugate gradient method for a symmetric and positive definite problem. Some of the PETSc SNES examples like ex12 compute the Jacobian A explicitly, but I was wondering if it's possible to convert this problem into matrix free form. That is, to directly assemble the product of A*x into vector form without first forming A and then invoking MatMult().<br><br></div>Thanks,<br></div>Justin<br></div>