<a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatShellSetOperation.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatShellSetOperation.html</a><br><div><br></div><div>Note that without an assembled matrix, you'll generally need to write your own preconditioner, which you can do using PCShell.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Nov 30, 2012 at 11:24 PM, Hesameddin Ilatikhameneh <span dir="ltr"><<a href="mailto:hilatikh@purdue.edu" target="_blank">hilatikh@purdue.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
We have a very large system (100M variables) which we want to minimize its energy. We are using Petsc extensively for matrix operations.<br>
I know that it is theoretically possible to use conjugate gradient method without storing Jacobian matrix (which is very huge in this case) since the conjugate gradient needs Jacobian times a vector not Jacobian itself.<br>

<br>
I read your manual about matrix free methods but I couldn't find where I can determine that function to Jacobian times a vector.<br>
<br>
I appreciate your help.<br>
<br>
Thanks<br>
Hesam<br>
</blockquote></div><br></div>