Change the type of multiplier for well connectivities to double

Author: michal-toth

opm/grid/common/ZoltanGraphFunctions.cpp

@@ -189,7 +189,7 @@ EdgeWeightType calculateWellEdgeWeight(const Dune::CpGrid& grid,
     }
 
     // when multipltWellConnectivities is provided, set the well weight to the average of grid weight times that coefficient
-    EdgeWeightType mWC = graph.getMultiplyWellConnectivities();
+    double mWC = graph.getMultiplyWellConnectivities();
     if (mWC >= 0) {
         if (total != std::numeric_limits<double>::infinity()) {
             total /= grid.numFaces();

opm/grid/common/ZoltanGraphFunctions.hpp

@@ -194,12 +194,12 @@ class CombinedGridWellGraph
             return 1.0;
     }
 
-    void setMultiplyWellConnectivities(const float& mWC)
+    void setMultiplyWellConnectivities(const double& mWC)
     {
         multiplyWellConnectivities = mWC;
     }
 
-    float getMultiplyWellConnectivities() const
+    double getMultiplyWellConnectivities() const
     {
         return multiplyWellConnectivities;
     }
@@ -250,16 +250,15 @@ class CombinedGridWellGraph
     int edgeWeightsMethod_;
     WellConnections well_indices_;
     double log_min_;
-    float multiplyWellConnectivities = -1;
+    double multiplyWellConnectivities = -1;
 };
 
 /// \brief Get the number of edges of the graph of the grid and the wells for one cell
 int getNumberOfEdgesForSpecificCellForGridWithWells(const CombinedGridWellGraph& graph, int localCellId);
 
-/// \brief Iterate over the grid and to calculate the edge weights
+/// \brief Iterate over the grid to calculate the edge weights between well connections
 ///
-/// Used as a weight for edges between cells of a well
-/// The result is the sum of all grid edges, or multiplyWellConnectivities*grid_average if it is positive
+/// The result is the sum of all grid edges, or multiplyWellConnectivities*grid_average if the multiplier is positive
 template <typename EdgeWeightType>
 EdgeWeightType calculateWellEdgeWeight(const Dune::CpGrid& grid,
                                        const CombinedGridWellGraph& graph);

opm/grid/common/ZoltanPartition.cpp

@@ -326,10 +326,10 @@ zoltanGraphPartitionGridOnRoot(const CpGrid& cpgrid,
     }
     setDefaultZoltanParameters(zz);
     Zoltan_Set_Param(zz, "IMBALANCE_TOL", std::to_string(zoltanImbalanceTol).c_str());
-    float mWC = -1;
+    double mWC = -1;
     for (const auto& [key, value] : params) {
         if (key=="MultiplyWellConnectivities") {
-            mWC = std::stof(value);
+            mWC = std::stod(value);
         } else {
             Zoltan_Set_Param(zz, key.c_str(), value.c_str());
         }
@@ -532,7 +532,7 @@ class ZoltanSerialPartitioner
 
         for (const auto& [key, value] : params) {
             if (key=="MultiplyWellConnectivities") {
-                gridAndWells->setMultiplyWellConnectivities(std::stof(value));
+                gridAndWells->setMultiplyWellConnectivities(std::stod(value));
             } else {
                  Zoltan_Set_Param(zz, key.c_str(), value.c_str());
             }