about KSP based on parallel dense matrix

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


Yujie,

>
> 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.

Hong

> 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