[petsc-users] ScaLAPACK EPS error

[Peder Jørgensgaard Olesen]

It appears my setup doesn't allow me to use versions > 3.17.4, unfortunately (I believe I'll need to speak to admin for this).

Best,
Peder
________________________________
Fra: Barry Smith <bsmith at petsc.dev>
Sendt: 18. januar 2024 19:29
Til: Peder Jørgensgaard Olesen <pjool at dtu.dk>
Cc: petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
Emne: Re: [petsc-users] ScaLAPACK EPS error


   Looks like you are using an older version of PETSc. Could you please switch to the latest and try again and send same information if that also fails.

  Barry


On Jan 18, 2024, at 12:59 PM, Peder Jørgensgaard Olesen via petsc-users <petsc-users at mcs.anl.gov> wrote:

Hello,

I need to determine the full set of eigenpairs to a rather large (N=16,000) dense Hermitian matrix. I've managed to do this using SLEPc's standard Krylov-Schur EPS, but I think it could be done more efficiently using ScaLAPACK. I receive the following error when attempting this. As I understand it, descinit is used to initialize an array, and the variable in question designates the leading dimension of the array, for which it seems an illegal value is somehow passed.

I know ScaLAPACK is an external package, but it seems as if the error would be in the call from SLEPc. Any ideas as to what could cause this?

Thanks,
Peder

Error message (excerpt):

PETSC ERROR: #1 MatConvert_Dense_ScaLAPACK() at [...]/matscalapack.c:1032
PETSC ERROR: #2 MatConvert at [...]/matrix.c:4250
PETSC ERROR: #3 EPSSetUp_ScaLAPACK() at [...]/scalapack.c:47
PETSC ERROR: #4 EPSSetUp() at [...]/epssetup.c:323
PETSC ERROR: #5 EPSSolve at [...]/epssolve.c:134
PETSC ERROR: ------ Error message ------
PETSC ERROR: Error in external library
PETSC ERROR: Error in ScaLAPACK subroutine descinit: info=-9
(...)

Log file (excerpt):
{  357,    0}:  On entry to DESCINIT parameter number   9 had an illegal value
[and a few hundred lines similar to this]

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240118/78bd4331/attachment.html>