OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-25 18:50:19 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #864 from lisajulia/feature/deterministic-indicesSyncer

Browse Source 
Construct the matrices in the AMG hierarchy with deterministic indices if possible.
This commit is contained in:
Markus Blatt 2024-01-29 13:59:44 +01:00 committed by GitHub
commit 3aa2e0969f
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/simulators/linalg/parallelamgbackend.hh
View File

@ -255,6 +255,18 @@ protected:
}
#endif
// trailing return type with decltype used for detecting existence of setUseFixedOrder member function by overloading the setUseFixedOrder function
template <typename C>
auto setUseFixedOrder(C criterion, bool booleanValue) -> decltype(criterion.setUseFixedOrder(booleanValue))
{
return criterion.setUseFixedOrder(booleanValue); // Set flag to ensure that the matrices in the AMG hierarchy are constructed with deterministic indices.
}
template <typename C>
void setUseFixedOrder(C, ...)
{
// do nothing, since the function setUseFixedOrder does not exist yet
}
void setupAmg_()
{
if (amg_)
@ -291,6 +303,7 @@ protected:
// coarsenCriterion.setAccumulate(Dune::Amg::noAccu);
coarsenCriterion.setAccumulate(Dune::Amg::atOnceAccu);
coarsenCriterion.setSkipIsolated(false);
setUseFixedOrder(coarsenCriterion, true); // If possible, set flag to ensure that the matrices in the AMG hierarchy are constructed with deterministic indices.
// instantiate the AMG preconditioner
#if HAVE_MPI