On Sun, Apr 10, 2011 at 2:13 AM, TAY wee-beng <span dir="ltr">&lt;<a href="mailto:zonexo@gmail.com">zonexo@gmail.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,<br>
<br>
I am already using PETSc to solve my momentum and poisson equations. However in some parts of my code, I need to solve a dense (usually) or sparse matrix, which arises from the radial basis function interpolation. Depending on the problem, it can be a big or small matrix.<br>

<br>
I am thinking whether to use PETSc or just a simple solver.<br>
<br>
Can you recommend the minimum size of sparse or dense matrix to use PETSc?<br></blockquote><div><br></div><div>For a dense matrix, we just call LAPACK in serial. And you can can change from dense to sparse if it is</div>
<div>sparse enough by changing the matrix type. The only regime to worry about is very large, dense matrices,</div><div>but I do not think you have those.</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;">

Thank you.<br>
<br>
-- <br>
Yours sincerely,<br><font color="#888888">
<br>
TAY wee-beng<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>