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
