On 12/28/23 12:22, Simon wrote:

I want to build the scalar product between rows of a
TrilinosWrappers::SparseMatrix and a TrilinosWrappers::MPI::Vector.

This sounds like you're computing a matrix-vector product -- the dot products you compute are the elements of the resulting vector.

As you have learned, it isn't entirely trivial to determine how information must be transmitted in parallel for these kinds of operations. But it is also not necessary: The parallel matrix classes all have matrix-vector ("vmult") operations implemented. Why not use those?

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/


--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/8e748b05-a6da-45d6-bd2a-23b552c818e6%40colostate.edu.

Reply via email to