Alexander :<div>Hmm, isrow is a matrix ordering, likely a bug :-(</div><div>Can you send us a stand-alone code that reproduces this error?</div><div>Do you use an external package for MatMatSolve()?</div><div><br></div><div>
Hong</div><div><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"><div class="im">
<div>On 13.08.2012 16:54, Hong Zhang wrote:<br>
</div>
<blockquote type="cite">
<div class="gmail_quote">On Mon, Aug 13, 2012 at 8:17 AM,
Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de" target="_blank">agrayver@gfz-potsdam.de</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
Trying to call this routine on 1 proc:<br>
<br>
call MatMatSolve(F,B,solutions,ierr);CHKERRQ(ierr)<br>
<br>
I get following crash:<br>
<br>
[0]PETSC ERROR: --------------------- Error Message
------------------------------------<br>
[0]PETSC ERROR: Null argument, when expecting valid pointer!<br>
[0]PETSC ERROR: Null Object: Parameter # 1!<br>
</blockquote>
<div> </div>
<div>This indicates matrix F (factored matrix) is not provided.</div>
<div><br>
</div>
<div>I tested petsc-3.3/src/mat/examples/tests/ex129.c with
MatMatSovle():</div>
<div>./ex129 -mat_type aij -dof 2</div>
<div>which works fine.</div>
</div>
</blockquote>
<br></div>
Hong,<br>
<br>
I get matrix F as:<br>
call PCFactorGetMatrix(pc,F,ierr);CHKERRQ(ierr)<br>
<br>
Please note that all checks in MatMatSolve are done before this
crash. <br>
So matrix F is provided correctly and it is not zero. <br>
<br>
The reason is that iscol and isrow == 0 here:<br>
<br>
PetscErrorCode MatMatSolve_SeqAIJ(Mat A,Mat B,Mat X)<br>
{<br>
Mat_SeqAIJ *a = (Mat_SeqAIJ*)A->data;<br>
IS iscol = a->col,isrow = a->row;<br>
<br>
...<br>
ierr = ISGetIndices(isrow,&rout);CHKERRQ(ierr); r = rout;
<<<< crash<br>
<br>
Again, same code runs fine when nproc > 1<div class="im"><br>
<br>
<blockquote type="cite">
<div class="gmail_quote">
<div>Hong</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 2, Fri Jul
13 15:42:00 CDT 2012<br>
[0]PETSC ERROR: See docs/changes/index.html for recent
updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble
shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: /app on a openmpi-i named glic1 by ag Mon Aug
13 15:06:36 2012<br>
[0]PETSC ERROR: Libraries linked from /home/lib/petsc-3.3-p2/openmpi-intel-complex-debug-f-ds/lib<br>
[0]PETSC ERROR: Configure run at Wed Jul 18 12:14:30 2012<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: ISGetIndices() line 370 in
/home/lib/petsc-3.3-p2/src/vec/is/interface/index.c<br>
[0]PETSC ERROR: MatMatSolve_SeqAIJ() line 1109 in
/home/lib/petsc-3.3-p2/src/mat/impls/aij/seq/aijfact.c<br>
[0]PETSC ERROR: MatMatSolve() line 3122 in
/home/lib/petsc-3.3-p2/src/mat/interface/matrix.c<br>
<br>
The same code works fine on nproc > 1.<br>
<br>
Any ideas?<span><font color="#888888"><br>
<br>
-- <br>
Regards,<br>
Alexander<br>
</font></span></blockquote>
</div>
<br>
</blockquote>
<br>
<br>
</div><span class="HOEnZb"><font color="#888888"><pre cols="72">--
Regards,
Alexander</pre>
</font></span></div>
</blockquote></div><br></div>