On Sat, Sep 4, 2010 at 10:26 AM, jinshan wu <span dir="ltr">&lt;<a href="mailto:wujinshan@yahoo.com">wujinshan@yahoo.com</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi all, I am new to this list. But it seems to have a lot of resources and<br>
friendly people here.<br>
<br>
I need to solve a huge linear system. In order to save some memory, I am<br>
thinking to store the matrices in its own data structure other than in real<br>
matrices. I understand that petsc can do that via GMRES (and others)<br>
without accessing matrix elements . So my question is where can we find some<br>
example? I am perfectly fine with supplying the matrix-vector product part by my<br>
own. I need more infor on how to connect it to the solver in petsc.<br></blockquote><div><br></div><div>You can use</div><div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatCreateShell.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatCreateShell.html</a></div>
<div><br></div><div>There are links to examples on that page.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Also another technical question: do we have to perform the preconditioner?<br></blockquote><div><br></div><div>Usually, if you want it to be efficient.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

thanks,<br>
<font color="#888888">Jinshan<br>
<br>
<br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>