[petsc-users] How to get diagonal entries ...

[Clemens Domanig]

Hi everyone,

I did a Cholesky-decomposition using MUMPS and now I want to get the 
diagonal entries of the Cholesky-decomposition but MatGetDiagonal 
doesnt't seem to work.

Thanks for your help - Clemens Domanig

...
ierr=PCSetType( prec, PCCHOLESKY); CHKERRQ( ierr);
...
ierr=PCFactorGetMatrix( prec, &M); CHKERRQ( ierr);
ierr=MatGetDiagonal( M, z); CHKERRQ( ierr);
printf("\nPivots:\n\n");
ierr=VecView( z, PETSC_VIEWER_STDOUT_SELF); CHKERRQ( ierr);

output
Process [0]
inf
0
0