<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Constantin:</div><div class="gmail_quote">I can reproduce your report, and am trying to debug it.</div><div class="gmail_quote">One thing I found is memory problem. W<span style="color:rgb(80,0,80)">ith one processor on o-build,</span></div><div class="gmail_quote"><font color="#500050">I get </font></div><div class="gmail_quote"><font color="#500050"><div class="gmail_quote">$ valgrind --tool=memcheck ./ex51f</div><div class="gmail_quote">==28659== Memcheck, a memory error detector</div><div class="gmail_quote">==28659== Copyright (C) 2002-2011, and GNU GPL'd, by Julian Seward et al.</div><div class="gmail_quote">==28659== Using Valgrind-3.7.0 and LibVEX; rerun with -h for copyright info</div><div class="gmail_quote">==28659== Command: ./ex51f</div><div class="gmail_quote">==28659==</div><div class="gmail_quote"> BEGIN PROC 0</div><div class="gmail_quote"> ITERATION 1</div><div class="gmail_quote">==28659== Invalid write of size 4</div><div class="gmail_quote">==28659== at 0x509B519: matload_ (in /scratch/hzhang/petsc/arch-linux-gcc-gfortran-o/lib/libpetsc.so.3.07.1)</div><div class="gmail_quote">==28659== by 0x402204: MAIN__ (in /scratch/hzhang/petsc/src/ksp/ksp/examples/tests/ex51f)</div><div class="gmail_quote">==28659== by 0x4030DC: main (in /scratch/hzhang/petsc/src/ksp/ksp/examples/tests/ex51f)</div><div class="gmail_quote">==28659== Address 0xe is not stack'd, malloc'd or (recently) free'd</div><div class="gmail_quote"><br></div><div class="gmail_quote">g-build does not show this error, plenty memory leak though.</div><div class="gmail_quote"><br></div><div class="gmail_quote">I'm further investing it...</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div><div class="gmail_quote"><br></div></font><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div class=""><div class="h5">Hi,<br>
<br>
I'm a new user of PETSc and I try to use it with MUMPS functionalities to compute a nullbasis.<br>
I wrote a code where I compute 4 times the same nullbasis. It does work well when I run it with several procs but with only one processor I get an error on the 2nd iteration when KSPSetUp is called. Furthermore when it is run with a debugger ( --with-debugging=yes), it works fine with one or several processors.<br>
Have you got any idea about why it doesn't work with one processor and no debugger?<br>
<br>
Thanks.<br>
Constantin.<br>
<br>
PS: You can find the code and the files required to run it enclosed.</div></div></blockquote></div><br></div></div>