Problem with MatMatMultTranspose

Hong Zhang hzhang at mcs.anl.gov
Mon Oct 20 14:06:16 CDT 2008


Zhifeng,

We do not have support for matrix operations on Hermitian matrix yet.
Hong

On Mon, 20 Oct 2008, zhifeng sheng wrote:

> Dear all
>
> I am using this MatMatMultTranspose function for complex matrices, but it 
> seems to be doing something weird.
>
> for instance, if I have complex matrix A, and I compute A^T*A with this 
> function, it does not generate a Hermitian 
> <http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html> 
> matrix.
>
> I am thinking that maybe the function take the transpose of A instead of the 
> conjugate transpose ....
>
> Do you know how I can get an A^H*A instead of A^T*A for complex matrices?
>
> Thanks a lot
> Best regards
> Zhifeng
> <http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html> 
>




More information about the petsc-users mailing list