MueLu: identify duplicates at region interfaces on coarse levels
Created by: mayrmt
We want to identify pairs of GIDs on coarse levels that refer to the same point on region interfaces.
We do not identify pairs of GIDs, but rather rely on knowing the coarsening of both regions attached to an interface, so we can figure out the mapping by integer arithmetic.
Motivation and Context
To enable transferring vectors from regional to composite layout (and vice versa), we need to identify at region interfaces which GIDs in the region layout refer to which GID in the composite layout. This identification must not rely on a priori knowledge and assumptions about the coarsening.
Definition of Done
recursively identify interface GIDs on coarse levels
form coarse level maps in composite, quasiRegional and regional layout
allow for multi-level case with recursive V-cycle implementation
- Inject coarse GIDs from region layout into fine level.
- Move to composite layout with a combineMode Minimum.
- Identify mapping by looking at aggregates
I'll work on this in branch
mayrmt-regions-multigrid in my fork and will share via a pull request.