Problem with MatMatMultTranspose

[Matthew Knepley]

On Thu, Sep 17, 2009 at 12:44 PM, Niall Moran <nmoran at thphys.nuim.ie> wrote:

> 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.
>

It is not curently in the todo list since we have only had one request. It
seems like it would just take being careful about the
complex case for SBAIJ if you would like the try the implementation. We can
answer questions.

  Thanks,

    Matt


> 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>
>>>>
>>>>
>>>
>>>
>


-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090917/88bdce57/attachment.htm>