Problem with MatMatMultTranspose

[Barry Smith]

    There is only support for CG with Hermitian transpose, ksp_cg_type  
symmetric or hermitian, KSPCGSetType()
the others only support complex, no Hermitian transpose.

    Barry

On Oct 22, 2008, at 2:40 AM, 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 
>>> >
>>
>