external solvers - scalapack

[Nichols A. Romero]

Barry,

PLAPACK has some dense eigensolvers. The public release version has QR
(which is probably too slow for us), but a developmental version (available
by contacting the authors) implements MR^3.

Does petsc have support for PLAPACK QR? If so, would it be possible to put support
for MR^3?

Nichols A. Romero, Ph.D.
Argonne Leadership Computing Facility
Argonne National Laboratory
Building 360 Room L-146
9700 South Cass Avenue
Argonne, IL 60490
(630) 252-3441


----- Original Message -----
From: "Barry Smith" <bsmith at mcs.anl.gov>
To: petsc-users at mcs.anl.gov
Sent: Wednesday, December 17, 2008 10:09:42 PM GMT -06:00 US/Canada Central
Subject: Re: external solvers - scalapack


    Whoever wrote a PETSc TACC2008 tutorial that says PETSc interfaces  
to scalapack
is a fool. It does not interface to scalapack in any way (since we  
support building
mumps which needs scalapack we support building scalapack but PETSc  
itself makes
no calls to scalapack.)*

    We do have support for using Cholesky in PLAPACK with our dense  
matrices.
You should use petsc-dev or wait a couple of days for the petsc-3.0  
release
before using it.

    Barry

* One reason we don't screw with scalapack is the complicated messy  
way the dense
matrix has to be spread across processes which is painful.


On Dec 17, 2008, at 9:36 PM, Nichols A. Romero wrote:

> Hi,
>
> I looked through the mailing lists to see if this question
> was already asked but could not find anything.
>
> In the PETSc TACC2008 tutorial, it states that PETSc interfaces
> with Scalapack?
>
> I would like to know what is supported from Scalapack. I am  
> particularly
> interested in:
> 1. Cholesky factorization (I know there are a couple of alternative  
> options here, e.g. PLAPACK).
> 2. Dense eigenvalue solver
>
>
>
> Nichols A. Romero, Ph.D.
> Argonne Leadership Computing Facility
> Argonne National Laboratory
> Building 360 Room L-146
> 9700 South Cass Avenue
> Argonne, IL 60490
> (630) 252-3441
>