Problem with MatMatMultTranspose

Jose E. Roman jroman at dsic.upv.es
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 mcs.anl.gov>  
> 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.

Jose



More information about the petsc-users mailing list