On Mon, Aug 8, 2011 at 2:17 AM, Gabriele Jost <span dir="ltr"><<a href="mailto:gjost@tacc.utexas.edu">gjost@tacc.utexas.edu</a>></span> wrote:<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="direction: ltr; color: rgb(0, 0, 0); "><div style="color: rgb(0, 0, 0); ">
<div style="font-family: 'Times New Roman'; font-size: 16px; "></div>
<div style="font-family: 'Times New Roman'; font-size: 16px; ">
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hello!
<div><br>
</div>
<div>I am working with a user of XSEDE (forerly know as TeraGrid) on converting his code from Scalapack to using PETSc. His physical problem yields matrices which are sparse, but the structure is rather general.</div>
<div>We have a version going that works very well for one type of problems. I can use MUMPS or SuperLU and I get results which are pretty close to what he gets when using Scalapack. </div>
<div>However, I am running into problems with the matrices that correspond to a different problem type. I get results when using MUMPS, but they differ widely from the Scalapack results. With MUMPS I get an error like:</div>

<div><br>
</div>
<div>Error reported by MUMPS in numerical factorization phase: INFO(1)=-1, INFO(2)=359</div>
<div><br>
</div>
<div>I noticed that for the matrix types where I get good results, the nonzeros reported for the PC Object type exactly the number of entries in my input data for the matrix (MPIAIJ format). For the type that does not work, petsc reports more nonzeros than
 what I enter. The message looks like this:</div>
<div><br>
</div>
<div>
<div>  linear system matrix = precond matrix:</div>
<div>  Matrix Object:</div>
<div>    type=mpiaij, rows=123120, cols=123120</div>
<div>    total: nonzeros=59094408, allocated nonzeros=59173613</div>
</div>
<div><br>
</div>
<div>My questions are:</div>
<div><br>
</div>
<div>1) Does anybody know what the error INFO(1)=-1 means in MUMPS?</div></div></div></div></div></div></blockquote><div><br></div><div>-1 means the error code is in INFO(2). I cannot find any online list of the error codes. I suggest mailing the MUMPS list.</div>
<div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="direction: ltr; color: rgb(0, 0, 0); ">
<div style="color: rgb(0, 0, 0); "><div style="font-family: 'Times New Roman'; font-size: 16px; "><div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<div>2) How does PETSc calculate the number of nonzeros for the PC matrix?</div>
<div><br>
</div>
<div>Have banging my head against this problem for a while. Any advice is highly appreciated.</div>
<div><br>
</div>
<div>Many thanks in advance!</div>
<div><br>
</div><font color="#888888">
<div>Gabriele Jost</div>
<div>Texas Advanced Computing Center</div>
<div><br>
</div>
<div><br>
</div>
</font></div>
</div>
</div>
</div>
</div>

</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>