<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 id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;">

<p>Hi, <br>

</p>

<p><br>

</p>

<p>I want to diagonalize a huge sparse matrix and I'm using the Kryov-Schur method with harmonic extraction (command line option -eps_harmonic ) implemented in SLEPc.

<br>

</p>

<p><br>

</p>

<p>I manually distribute a sparse matrix across several CPUs and everything works fine when:</p>

<p>- I do _not_ use harmonic extraction<br>

</p>

<p>- I use harmonic extraction on only a single CPU <br>

</p>

<p><br>

</p>

<p>If I try do use harmonic extraction on multiple CPUs, I get a memory corruption.</p>

<p>I'm not quite sure where to look at, but somewhere in the output, I find:</p>

<p><br>

</p>

<p><br>

</p>

<p></p>

<div>[1]PETSC ERROR: PetscMallocValidate: error detected at PetscSignalHandlerDefault() line 145 in /home/applications/sources/libraries/petsc-3.9.3/src/sys/error/signal.c<br>

[1]PETSC ERROR: Memory [id=0(9072)] at address 0x145bcd0 is corrupted (probably write past end of array)<br>

[1]PETSC ERROR: Memory originally allocated in DSAllocateWork_Private() line 74 in /home/applications/sources/libraries/slepc-3.9.2/src/sys/classes/ds/interface/dspriv.c<br>

<br>

</div>

<p></p>

<p><br>

</p>

<p>Now, I have the feeling that this might be a bug in SLEPc because, if I had messed up the matrix initialization and distribution, I should also get a memory corruption when I don't use harmonic extraction, right?</p>

<p><br>

Any suggestions what might be going on?</p>

<p><br>

</p>

<p>Regards, <br>

</p>

<p>Moritz<br>

</p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

</div>

</body>

</html>