Problem with MatMatMultTranspose

[Niall Moran]

Hi,

I am just wondering if anything has changed on the status of this 
feature. Would be great to be able to perform matrix vector 
multiplications on complex Hermitian matrices by only providing one half 
of the matrix.

Regards,

Niall.

zhifeng sheng wrote:
> 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> 
>>>
>>