<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I set up the matrix using MatCreateDense(), passing PETSC_DECIDE for the local dimensions.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
The same error appears with 8, 12, and 16 nodes (32 proc/node).</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I'll have to get back to you regarding a minimal example.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Peder</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Fra:</b> Jose E. Roman <jroman@dsic.upv.es><br>
<b>Sendt:</b> 18. januar 2024 19:28<br>
<b>Til:</b> Peder Jørgensgaard Olesen <pjool@dtu.dk><br>
<b>Cc:</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Emne:</b> Re: [petsc-users] ScaLAPACK EPS error</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">How are you setting up your input matrix? Are you giving the local sizes or setting them to PETSC_DECIDE?<br>
Do you get the same error for different number of MPI processes?<br>
Can you send a small code reproducing the error?<br>
<br>
Jose<br>
<br>
<br>
> El 18 ene 2024, a las 18:59, Peder Jørgensgaard Olesen via petsc-users <petsc-users@mcs.anl.gov> escribió:<br>
> <br>
> Hello,<br>
> <br>
> 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.<br>
> <br>
> 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?<br>
> <br>
> Thanks,<br>
> Peder<br>
> <br>
> Error message (excerpt):<br>
> <br>
> PETSC ERROR: #1 MatConvert_Dense_ScaLAPACK() at [...]/matscalapack.c:1032<br>
> PETSC ERROR: #2 MatConvert at [...]/matrix.c:4250<br>
> PETSC ERROR: #3 EPSSetUp_ScaLAPACK() at [...]/scalapack.c:47<br>
> PETSC ERROR: #4 EPSSetUp() at [...]/epssetup.c:323<br>
> PETSC ERROR: #5 EPSSolve at [...]/epssolve.c:134<br>
> PETSC ERROR: ------ Error message ------<br>
> PETSC ERROR: Error in external library<br>
> PETSC ERROR: Error in ScaLAPACK subroutine descinit: info=-9<br>
> (...)<br>
> <br>
> Log file (excerpt):<br>
> { 357, 0}: On entry to DESCINIT parameter number 9 had an illegal value<br>
> [and a few hundred lines similar to this]<br>
<br>
<br>
</div>
</span></font></div>
</body>
</html>