Skip to content
GitLab

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