added maximum based tensor method

The addition of the experimental maximum based tensor based load-balancing method.

include/ALL_Tensor.hpp

@@ -78,7 +78,7 @@ public:
 
   /// method to execute a load-balancing step
   /// @param[in] step number of the load-balancing step
-  virtual void balance(int step) override;
+  virtual void balance(int step) override {balance(step, MPI_SUM);}
 
   // getter for variables (passed by reference to avoid
   // direct access to private members of the object)
@@ -92,6 +92,9 @@ public:
   /// @param[in] data pointer to the data structure
   virtual void setAdditionalData(known_unused const void *data) override {}
 
+  /// method to execute a load-balancing step
+  /// @param[in] step number of the load-balancing step
+  virtual void balance(int step, MPI_Op reductionMode);
 private:
   // type for MPI communication
   MPI_Datatype MPIDataTypeT;
@@ -104,9 +107,10 @@ private:
   // list of neighbors
   std::vector<int> neighbors;
   std::vector<int> nNeighbors;
+  
+protected:
+
 
-  // used MPI reduction method
-  const MPI_Op reductionMode = MPI_SUM;
 };
 
 // setup routine for the tensor-based load-balancing scheme
@@ -176,17 +180,20 @@ template <class T, class W> void Tensor_LB<T, W>::setup() {
   }
 }
 
-template <class T, class W> void Tensor_LB<T, W>::balance(int) {
+template <class T, class W> void Tensor_LB<T, W>::balance(int step, MPI_Op reductionMode) {
   int dim = this->getDimension();
   std::vector<Point<T>> newVertices = this->vertices;
   this->prevVertices = this->vertices;
 
+  bool localIsSum = reductionMode == MPI_SUM;
+  bool localIsMax = reductionMode == MPI_MAX;
+
   // loop over all available dimensions
   for (int i = 0; i < dim; ++i) {
     W work_local_plane;
     // collect work from all processes in the same plane
     MPI_Allreduce(this->getWork().data(), &work_local_plane, 1, MPIDataTypeW,
-                  this->reductionMode, communicators.at(i));
+                  reductionMode, communicators.at(i));
 
     // correct right border:
 
@@ -277,10 +284,7 @@ template <class T, class W> class TensorMax_LB : public Tensor_LB<T, W> {
 public:
   TensorMax_LB() {}
   TensorMax_LB(int d, W w, T g) : Tensor_LB<T, W>(d, w, g) {}
-
-private:
-  // used MPI reduction method
-  const MPI_Op reductionMode = MPI_MAX;
+  virtual void balance(int step) override {Tensor_LB<T,W>::balance(step, MPI_MAX);}
 };
 
 } // namespace ALL