On Thu, Sep 17, 2009 at 12:44 PM, Niall Moran <span dir="ltr"><<a href="mailto:nmoran@thphys.nuim.ie">nmoran@thphys.nuim.ie</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
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.<br></blockquote>
<div><br>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<br>complex case for SBAIJ if you would like the try the implementation. We can answer questions.<br>
<br> Thanks,<br><br> Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Regards,<br>
<br>
Niall.<br>
<br>
zhifeng sheng wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
you mean the conjugate transpose for complex matrix is not supported?<br>
<br>
then how can you implement the iterative solvers for complex matrices? because, some iterative solvers need it.<br>
<br>
Thanks<br>
Best regards<br>
<br>
<br>
<br>
Hong Zhang wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Zhifeng,<br>
<br>
We do not have support for matrix operations on Hermitian matrix yet.<br>
Hong<br>
<br>
On Mon, 20 Oct 2008, zhifeng sheng wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear all<br>
<br>
I am using this MatMatMultTranspose function for complex matrices, but it seems to be doing something weird.<br>
<br>
for instance, if I have complex matrix A, and I compute A^T*A with this function, it does not generate a Hermitian <<a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html" target="_blank">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html</a>> matrix.<br>
<br>
I am thinking that maybe the function take the transpose of A instead of the conjugate transpose ....<br>
<br>
Do you know how I can get an A^H*A instead of A^T*A for complex matrices?<br>
<br>
Thanks a lot<br>
Best regards<br>
Zhifeng<br>
<<a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html" target="_blank">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html</a>> <br>
</blockquote>
<br>
</blockquote></blockquote>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>