Problem with MatMatMultTranspose

Jose E. Roman jroman at
Fri Sep 18 09:19:31 CDT 2009

On 18/09/2009, Matthew Knepley wrote:

> On Fri, Sep 18, 2009 at 8:48 AM, Hong Zhang <hzhang at>  
> wrote:
> Should we create a new matrix type, say hbaij,
> for operations of Hermitian matrices?
> No, it should be a minor perturbation.
>   Matt

I think that the complex SBAIJ can handle both cases by taking into  
account the MAT_SYMMETRIC and MAT_HERMITIAN flags (set with  
MatSetOption), then decide whether elements should be conjugated or not.


More information about the petsc-users mailing list