[petsc-users] MatGetDiagonal

[Uwe Schlifkowitz]

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