about KSP based on parallel dense matrix

Hong Zhang hzhang at mcs.anl.gov
Tue Sep 23 10:26:08 CDT 2008


> Barry has tried to establish an interface for Plapack. However, there are
> some bugs in Plapack. Therefore, it doesn't work. I am wondering if CG in

The Plapack interface in the realeased petsc-2.3.3 should work fine.
Additional work is needed to make it work under reorganization
in petsc-dev, which I'm working on.

Note, Plapack implements LU precontioner.


> Petsc can work with parallel dense matrix. When using the same matrix, which
> one is faster, sequential or parallel? thanks.
> Regards,
> Yujie
> On Mon, Sep 22, 2008 at 6:51 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
>> Well, any iterative solver will actually work, but expect a really
>> poor scalability :-). I believe (never used dense matrices) that you
>> could use a direct method (PLAPACK?), but again, be prepared for long
>> running times if your problem is (even moderately) large.
>> On Mon, Sep 22, 2008 at 10:35 PM, Yujie <recrusader at gmail.com> wrote:
>>> To my knowledge, PETsc doesn't provide parallel dense matrix-based
>> solvers,
>>> such as for CG, GMRES and so on. If it is, how to deal with this problem?
>>> Thanks.
>>> Regards,
>>> Yujie
>> --
>> Lisandro Dalcín
>> ---------------
>> Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
>> Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
>> Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
>> PTLC - Güemes 3450, (3000) Santa Fe, Argentina
>> Tel/Fax: +54-(0)342-451.1594

More information about the petsc-users mailing list