[petsc-users] A^TA

[Jed Brown]

On Sat, Feb 4, 2012 at 01:47, Kuhlemann, Verena <vkuhlem at emory.edu> wrote:

>  Hi,
>
>  one quick question: I have a matrix P and need to compute P^t*P.
> Is there a function available for this? I couldn't find anything.
> Right now I am computing P^t*P by using MatPtAP with A=identity,
> but I was hoping there is another way where I don't need to set up
> the identity matrix.
>

Use petsc-dev and MatTransposeMatMult(P,P,...,&PtP)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120204/716f74c5/attachment.htm>