Problem with MatMatMultTranspose

[zhifeng sheng]

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>