[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
  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?

-------------- 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