Reply To: MPI failed with 12 cores

[Adrian]

This is the patch that will be included in the next release:

diff –git a/src/functors/lattice/indicator/blockIndicatorF2D.hh b/src/functors/lattice/indicator/blockIndicatorF2D.hh
index 8d01a6b03c8cea37d1d9ec05033d1bda4a616efa..f52a48a3ec4b2b966e6341e186d253fe73bc453c 100644
— a/src/functors/lattice/indicator/blockIndicatorF2D.hh
+++ b/src/functors/lattice/indicator/blockIndicatorF2D.hh
@@ -124,7 +124,7 @@ bool BlockIndicatorMaterial2D<T>::operator() (bool output[], const int input[])
// of BlockGeometry2D<T>::get to avoid resetting geometry
// statistics:
const BlockGeometryStructure2D<T>& blockGeometry = this->_blockGeometryStructure;
– const int current = blockGeometry.get(input[0], input[1]);
+ const int current = blockGeometry.getMaterial(input[0], input[1]);
output[0] = std::any_of(_materials.cbegin(),
_materials.cend(),
[current](int material) { return current == material; });
diff –git a/src/functors/lattice/indicator/blockIndicatorF3D.hh b/src/functors/lattice/indicator/blockIndicatorF3D.hh
index dcf2f5fc0e7658c487d671a02322eea51669378f..60c00322c4f79dacf7cf88ab0687ee746f65bfb8 100644
— a/src/functors/lattice/indicator/blockIndicatorF3D.hh
+++ b/src/functors/lattice/indicator/blockIndicatorF3D.hh
@@ -136,7 +136,7 @@ bool BlockIndicatorMaterial3D<T>::operator() (bool output[], const int input[])
// of BlockGeometry3D<T>::get to avoid resetting geometry
// statistics:
const BlockGeometryStructure3D<T>& blockGeometry = this->_blockGeometryStructure;
– const int current = blockGeometry.get(input[0], input[1], input[2]);
+ const int current = blockGeometry.getMaterial(input[0], input[1], input[2]);
output[0] = std::any_of(_materials.cbegin(),
_materials.cend(),
[current](int material) { return current == material; });

i.e. change BlockGeometry(2,3)D::get calls to BlockGeometry(2,3)D::getMaterial for the block indicators.