By looking at src/mat/examples/tests/output/ex88_1.out, it seems that the output of MatGetDiagonal is wrong for MATSHELL.
it should be 2.96782e+08 1.41731e+09 and not (as computed) 1.10918e+08 2.06459e+08 -- Stefano
By looking at src/mat/examples/tests/output/ex88_1.out, it seems that the output of MatGetDiagonal is wrong for MATSHELL.
it should be 2.96782e+08 1.41731e+09 and not (as computed) 1.10918e+08 2.06459e+08 -- Stefano