[petsc-users] MatMultTranspose memory usage

[Zhang, Hong]

Karl:
What matrix format do you use? Run it in parallel or sequential?
Hong

We used /proc/self/stat to track the resident set size during program run, and we saw the resident set size jumped by the size of the matrix right after we did matmulttranspose.

On Wed, Jul 17, 2019 at 12:04 PM hong--- via petsc-users <petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov>> wrote:
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/94e64e69/attachment-0001.html>