Jose :
> sbaij is for symmetric matrix, not Hermitian, as you noticed. > > Petsc does not have support for Hermitian operations for Cholesky solve yet. > We need MatSolve_SeqSBAIJ_1_Hermitian().
Then why not give a PETSC_ERR_SUP error?
Jose