Problem with MatMatMultTranspose

[zhifeng sheng]

Hi,

you mean the conjugate transpose for complex matrix is not supported?

then how can you implement the iterative solvers for complex matrices? 
because, some iterative solvers need it.

Thanks
Best regards



Hong Zhang wrote:
>
> 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> 
>>
>