<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I would like to access the approximate Jacobian in quasi-Newton solvers (e.g., BFGS) to analyze its condition number. I guess for BFGS, where the approximate Jacobian is built from outer products of several vectors coming from current and previous gradient vectors, the approximate Jacobian might not be stored explicitly as a matrix and instead the Jacobian times a vector is performed in a matrix-free manner. My reason for believing so is the fact that when I try to retrieve the Jacobian matrix using SNESGetJacobian and feed it to SLEPc, I incur segmentation fault at EPSSetOperators. <br><br></div>If that's the case, is there a way to still get access to the approximate Jacobian and feed it to SLEPc for eigen analysis, like with the use of MatShell operations? <br><br></div>Thanks,<br></div>Bikash <br clear="all"><div><div><div><div><div><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><div><font color="#666666">Bikash S. Kanungo<br></font></div><font color="#666666">PhD Student<br></font></div><font color="#666666">Computational Materials Physics Group<br></font></div><font color="#666666">Mechanical Engineering <br></font></div><font color="#666666">University of Michigan<br><br></font></div></div>
</div></div></div></div></div></div>