How are you setting up your input matrix? Are you giving the local sizes or
setting them to PETSC_DECIDE?
Do you get the same error for different number of MPI processes?
Can you send a small code reproducing the error?
Jose
> El 18 ene 2024, a las 18:59, Peder Jørgensgaard Olesen via petsc-users
> <petsc-users@mcs.anl.gov> escribió:
>
> 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]
