[Trilinos-Users] Tpetra::CrsMatrix caching related to C= A^T B and D = A*C
Denis Davydov
davydden at gmail.com
Mon Sep 10 10:07:58 EDT 2018
Hi all,
I was looking at the documentation of TPetra and could not immediately see whether or not
something is cached for matrix-matrix products.
Say I have tall and skinny sparse matrices A, B, D and need to calculate the products
C=A^T B
and
D = A*C
on each iteration of a non-linear solver.
Both the intersection of sparsities of A and B as well as MPI communication pattern does not change from iteration to iteration.
But it looks like Tpetra::MatrixMatrix::Multiply() will be calculating those things every time anew.
Is it indeed the case or is there some mechanism behind-the-scenes which caches this information?
Regards,
Denis.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://trilinos.org/pipermail/trilinos-users/attachments/20180910/fd005fce/attachment.html>
More information about the Trilinos-Users
mailing list