Re: Eigensystems in Fortran with MPI

From: Ian Bush (I.J.Bush_at_nospam.dl.ac.uk)
Date: 02/22/04 

Date: Sun, 22 Feb 2004 13:16:51 +0000

Maxime Boissonneault wrote:

> Hello,
> We want to diagonalize a matrix (so find its eigenvectors and eigenvalues)
> on multiple processors. We want to use ScaLAPACK library to do so.
> However, we have problems doing it. Does anyone have some examples of code
> that do this ? Our matrix are hemitians.
>

Could you be more specific about "problems" ? Most of the errors I have
seen are to do with not getting the data distribution correct - it's
a little tricky to get your head around it ( at least form me, it was )
but once you do it's quite convenient,

Ian

Relevant Pages

  • Eigensystems in Fortran with MPI
    ... We want to diagonalize a matrix (so find its eigenvectors and eigenvalues) ... on multiple processors. ... Maxime Boissonneault ...
    (comp.lang.fortran)
  • Re: solving this impasse could result in solving the Navier-Stokes
    ... where u is an incompressible flow field in wavenumber space and A ... where A_n are matrices that are a function of wavenumber alone and c_n ... eigenvectors of Ain can be a function of time. ... diagonalize nicely are completely correct, ...
    (sci.math)
  • Re: Matrix Diagonalization
    ... LAPACK and virtually everything else get that one wrong! ... Run it, extracting all eigenvectors, and see :-) ... then the routine will ignore the lower triangle of the ... input and diagonalize a unit matrix of dimension 2x2. ...
    (comp.lang.fortran)