[petsc-users] MatMultTranspose memory usage

[hong at aspiritech.org]

Kun:
How do you know 'MatMultTranpose creates an extra memory copy of matrix'?
Hong

Hi,
>
>
>
> I was using MatMultTranpose and MatMult to solver a linear system.
>
>
>
> However we found out, MatMultTranpose create an extra memory copy of
> matrix for its operation. This extra memory copy is not stated everywhere
> in petsc manual.
>
>
>
> This basically double my memory requirement to solve my system.
>
>
>
> I remember mkl’s routine can do inplace matrix transpose vector product,
> without transposing the matrix itself.
>
>
>
> Is this always the case? Or there is way to make petsc to do inplace
> matrix transpose vector product.
>
>
>
> Any help is greatly appreciated.
>
>
>
> Regards,
>
> Kun
>
>
>
> Schlumberger-Private
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190717/e6c16fa1/attachment.html>