So this turned out to be a bug in my own code after I made the original change in trilinos_sparse_matrix.cc. When I was constructing the sparsity patterns for the matrices using DoFTools::make_sparsity_pattern, there is a flag I had set keep_constrained_dofs = false. Setting this flag to true has fixed my matrix multiplication problem in parallel (though it's not really clear to me why this would only cause a problem in parallel).
-- 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. For more options, visit https://groups.google.com/d/optout.