You can also look into SLEPc to compute eigenvalues in parallel for the global system.
SLEPc is based on PETSc so it should not be too hard to use in your code.

Best,
Luc

On 01/21/2015 10:53 AM, Matthew Knepley wrote:
On Wed, Jan 21, 2015 at 9:51 AM, siddhesh godbole <[email protected] <mailto:[email protected]>> wrote:

    Hello,

    I want to execute Monte Carlo simulations for eigenvalue problem
    in parallel with PETSC. i am using LAPACKsygvx function to compute
    all the eigenvalues. Now the problem i am facing is that, i think
    LAPACKsygvx requires Matrices to be declared as MATSEQDENSE, and
    this sequential type is not allowing me to execute this program on
    multiple processors.

    How should i go about this problem? can i use different
    communicators or something? Please help!


Do you need all the eigenvalues? If so, you can try to use Elemental, with the development version of PETSc.

  Thanks,

    Matt


    *Siddhesh M Godbole*

    5th year Dual Degree,
    Civil Eng & Applied Mech.
    IIT Madras




--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener

Reply via email to