<br><br><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,<u></u>ierr);CHKERRQ(ierr)<br>
<br>
I get following crash:<br>
<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------<u></u>------<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><br></div><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: ------------------------------<u></u>------------------------------<u></u>------------<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: ------------------------------<u></u>------------------------------<u></u>------------<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/<u></u>openmpi-intel-complex-debug-f-<u></u>ds/lib<br>
[0]PETSC ERROR: Configure run at Wed Jul 18 12:14:30 2012<br>
[0]PETSC ERROR: ------------------------------<u></u>------------------------------<u></u>------------<br>
[0]PETSC ERROR: ISGetIndices() line 370 in /home/lib/petsc-3.3-p2/src/<u></u>vec/is/interface/index.c<br>
[0]PETSC ERROR: MatMatSolve_SeqAIJ() line 1109 in /home/lib/petsc-3.3-p2/src/<u></u>mat/impls/aij/seq/aijfact.c<br>
[0]PETSC ERROR: MatMatSolve() line 3122 in /home/lib/petsc-3.3-p2/src/<u></u>mat/interface/matrix.c<br>
<br>
The same code works fine on nproc > 1.<br>
<br>
Any ideas?<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Regards,<br>
Alexander<br>
</font></span></blockquote></div><br>