[petsc-users] MatMultTranspose memory usage

[Kun Jiao]

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/bf4dab10/attachment-0001.html>