Dear list, as explained in the last thread (see http://lists.mcs.anl.gov/pipermail/petsc-users/2011-August/009706.html ) about MatGetDiagonal, the fix for a bug was pushed to petsc-dev.
My results from running example 52 from ~petsc-dev/src/ksp/ksp/examples/tutorials/ex52.c are attached below. Are these the correct results? I ran the example on Max OS X 10.6.8 as well as on CentOS with both petsc-3.2-p5 and petsc-dev (pulled from hg yesterday) and obtained the same results. Please let me know if you need the full output. (I recently started working on Clemens' project, i think he left in August.) Regards, Uwe $ ./ex52 -use_petsc_ch Matrix Object: 1 MPI processes type: seqaij row 0: (0, 4) (5, -1) (28, -1) (29, -1) row 1: (1, 4) (5, -1) (12, -1) (13, -1) ... Vector Object: 1 MPI processes type: seq 4 4 4 4 4 3.25 4 3.5 4 4 3.75 3.46548 3.1138 2.70534 ... 3.10993 2.55974 2.38793 2.32336 2.34087 2.5002 3.30916 Norm of error < 1.e-12 iterations 1
