<div dir="ltr">thank you very much, Hong.<br><br><div class="gmail_quote">On Tue, Sep 23, 2008 at 8:26 AM, Hong Zhang <span dir="ltr">&lt;<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
Yujie,<div class="Ih2E3d"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Barry has tried to establish an interface for Plapack. However, there are<br>
some bugs in Plapack. Therefore, it doesn&#39;t work. I am wondering if CG in<br>
</blockquote>
<br></div>
The Plapack interface in the realeased petsc-2.3.3 should work fine.<br>
Additional work is needed to make it work under reorganization<br>
in petsc-dev, which I&#39;m working on.<br>
<br>
Note, Plapack implements LU precontioner.<br><font color="#888888">
<br>
Hong</font><div><div class="Wj3C7c"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Petsc can work with parallel dense matrix. When using the same matrix, which<br>
one is faster, sequential or parallel? thanks.<br>
<br>
Regards,<br>
<br>
Yujie<br>
<br>
On Mon, Sep 22, 2008 at 6:51 PM, Lisandro Dalcin &lt;<a href="mailto:dalcinl@gmail.com" target="_blank">dalcinl@gmail.com</a>&gt; wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Well, any iterative solver will actually work, but expect a really<br>
poor scalability :-). I believe (never used dense matrices) that you<br>
could use a direct method (PLAPACK?), but again, be prepared for long<br>
running times if your problem is (even moderately) large.<br>
<br>
On Mon, Sep 22, 2008 at 10:35 PM, Yujie &lt;<a href="mailto:recrusader@gmail.com" target="_blank">recrusader@gmail.com</a>&gt; wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
To my knowledge, PETsc doesn&#39;t provide parallel dense matrix-based<br>
</blockquote>
solvers,<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
such as for CG, GMRES and so on. If it is, how to deal with this problem?<br>
Thanks.<br>
<br>
Regards,<br>
<br>
Yujie<br>
<br>
</blockquote>
<br>
<br>
<br>
--<br>
Lisandro Dalcín<br>
---------------<br>
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)<br>
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)<br>
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)<br>
PTLC - Güemes 3450, (3000) Santa Fe, Argentina<br>
Tel/Fax: +54-(0)342-451.1594<br>
<br>
<br>
</blockquote>
</blockquote>
</div></div></blockquote></div><br></div>