<div dir="ltr">ok, thank you very much!!<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Jan 19, 2014 at 10:10 PM, Karl Rupp <span dir="ltr"><<a href="mailto:rupp@mcs.anl.gov" target="_blank">rupp@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
ah, sorry, the matrix view was okay, I'm still confused from a recent debugging session. If you want to verify that the distribution is okay, verify the return values of MatGetOwnerShipRange() to get the rows 'owned' by the respective process. In your example you should get {0,1} for the first MPI rank and {2,3} for the second.<div class="im">
<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Example:<br>
<a href="http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex11.c.html" target="_blank">http://www.mcs.anl.gov/petsc/<u></u>petsc-current/src/ksp/ksp/<u></u>examples/tutorials/ex11.c.html</a><br>
Run it with<br>
  mpirun -np 2 ./Release/spin-scenarios-petsc  -n 2<br>
This is the results:<br>
1+1i<br>
1+1i<br>
Norm of error 1.19381e-16 iterations 4<br>
Vector Object: 2 MPI processes<br>
   type: mpi<br>
Process [0]<br>
4.16334e-17<br>
0 + 1.38778e-17 i<br>
Process [1]<br>
0<br>
1.11022e-16<br>
Matrix Object: 1 MPI processes<br>
   type: mpiaij<br>
row 0: (0, -7.11111 + 0.0800036 i) (1, -1)  (2, -1)<br>
row 1: (0, -1)  (1, -7.11111 + 0.00111359 i) (3, -1)<br>
row 2: (0, -1)  (2, -7.11111 + 0.0601271 i) (3, -1)<br>
row 3: (1, -1)  (2, -1)  (3, -7.11111 + 0.0943018 i)<br>
<br>
I am not sure about it, but the results is the same as I tried before<br>
with c++ class.<br>
</blockquote>
<br></div>
The output looks okay to me - sorry for the confusion.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
By the way, I install PETSc using<br>
--with-mpi-dir=/usr/local --with-scalar-type=complex<br>
--with-clanguage=cxx --with-fortran-kernels=1<br>
--with-blas-lapack-dir=/opt/<u></u>intel/composer_xe_2013_sp1.0.<u></u>080/mkl<br>
--with-debugging=no<br>
</blockquote>
<br></div>
Please perform *all* development using --with-debugging=yes and only use optimized builds once everything is working. Chances are high that this will save you a lot of headache :-)<br>
<br>
Best regards,<br>
Karli<br>
<br>
</blockquote></div><br></div>