MueLu: add local values in regionalToComposite manually
Created by: mayrmt
@trilinos/muelu
Expectations
Replace current use of Epetra_CombineMode Epetra_AddLocalAlso
in the MueLu HHG driver by an implementation that performs the local adds manually when transforming a region matrix back to a composite matrix.
Current Behavior
In the HHG driver's routine regionToComposite()
, we rely on the Epetra_CombineMode Epetra_AddLocalAlso
to not only add values accross processes, but also to add process-local values.
Motivation and Context
Tpetra/Xpetra
does not support a AddLocalAlso
-type CombineMode.
Definition of Done
-
Add manual add operation for processor-local values. -
Remove use of Epetra_CombineMode Epetra_AddLocalAlso
.
Possible Solution
This is already implemented for vector operations. We might follow a similar idea.
There's already some code that's trying to do that by perform inter-process add operations based on a Xpetra::CombineMode Add
, but performs on-process operations separately. We just have to get it right.
Related Issues
- Blocks #4084 (closed)