MueLu: filtered matrix lumping may be buggy
Created by: aprokop
A potential bug in MueLu's
FilteredAFactory lumping procedure affecting both non-kokkos and kokkos branches.
MueLu_FilteredAFactory_def.hpp the lumping procedure seem to only use local indices for determining which element is diagonal. This is a problem for situations when row map and column map do not "match". In that situation, the dropped entries may be added to non-diagonal element thus exacerbating situation.
Who is affected
Potentially anybody using non-standard column map in conjuction with filtering. Panzer?
How to address
As a first step, we need to identify if anybody is getting hit by this. I propose adding a check
isSameMap similarly to
TentativePFactory. If we do get somebody affected, we need to decide what to do. Going to global indices is going to be very expensive.