El 26/07/2012, a las 21:34, Hong Zhang escribió:

Jose :

Then why not give a PETSC_ERR_SUP error?

Patched 3.3 and petsc-dev. http://petsc.cs.iit.edu/petsc/petsc-dev/rev/ff096aff2a39

When I find time, I'll implement MatSolve() for seqsbaij_hermitian :-)

Hong

Thanks. I think the same error should go to MatGetFactor_seqaij_petsc otherwise my second sample run gives wrong results: $ ./ex10 -f $PETSC_DIR/share/petsc/datafiles/matrices/hpd-complex-int32-float64 -ksp_type preonly -pc_type cholesky Number of iterations = 1 Residual norm 1.64361 Jose