moving variableWellStateIndices to StandardWells

Probably it is the time to introduce SeqQs and SeqP for the
multi-segment wells.

opm/autodiff/BlackoilModelBase.hpp

@@ -363,9 +363,6 @@ namespace Opm {
         std::vector<int>
         variableStateIndices() const;
 
-        std::vector<int>
-        variableWellStateIndices() const;
-
         SolutionState
         variableStateExtractVars(const ReservoirState& x,
                                  const std::vector<int>& indices,

opm/autodiff/BlackoilModelBase_impl.hpp

@@ -612,24 +612,6 @@ namespace detail {
 
 
 
-    template <class Grid, class WellModel, class Implementation>
-    std::vector<int>
-    BlackoilModelBase<Grid, WellModel, Implementation>::
-    variableWellStateIndices() const
-    {
-        // Black oil model standard is 5 equation.
-        // For the pure well solve, only the well equations are picked.
-        std::vector<int> indices(5, -1);
-        int next = 0;
-        indices[Qs] = next++;
-        indices[Bhp] = next++;
-        assert(next == 2);
-        return indices;
-    }
-
-
-
-
 
     template <class Grid, class WellModel, class Implementation>
     typename BlackoilModelBase<Grid, WellModel, Implementation>::SolutionState
@@ -1092,7 +1074,7 @@ namespace detail {
         V aliveWells;
         const int np = wells().number_of_phases;
         std::vector<ADB> cq_s(np, ADB::null());
-        std::vector<int> indices = variableWellStateIndices();
+        std::vector<int> indices = asImpl().stdWells().variableWellStateIndices();
         SolutionState state0 = state;
         WellState well_state0 = well_state;
         asImpl().makeConstantState(state0);

opm/autodiff/BlackoilMultiSegmentModel.hpp

@@ -223,7 +223,7 @@ namespace Opm {
         using Base::convergenceReduction;
         using Base::maxResidualAllowed;
         using Base::variableState;
-        using Base::variableWellStateIndices;
+        // using Base::variableWellStateIndices;
         using Base::asImpl;
 
         const std::vector<WellMultiSegmentConstPtr>& wellsMultiSegment() const { return wells_multisegment_; }

opm/autodiff/BlackoilMultiSegmentModel_impl.hpp

@@ -1559,7 +1559,7 @@ namespace Opm {
         V aliveWells;
         const int np = wells().number_of_phases;
         std::vector<ADB> cq_s(np, ADB::null());
-        std::vector<int> indices = variableWellStateIndices();
+        std::vector<int> indices = stdWells().variableWellStateIndices();
         SolutionState state0 = state;
         WellState well_state0 = well_state;
         makeConstantState(state0);

opm/autodiff/StandardWells.hpp

@@ -185,6 +185,8 @@ namespace Opm {
                                      int& next) const;
 
 
+            std::vector<int>
+            variableWellStateIndices() const;
 
         protected:
             bool wells_active_;

opm/autodiff/StandardWells_impl.hpp

@@ -1130,4 +1130,22 @@ namespace Opm
         indices[Bhp] = next++;
     }
 
+
+
+
+
+    std::vector<int>
+    StandardWells::variableWellStateIndices() const
+    {
+        // Black oil model standard is 5 equation.
+        // For the pure well solve, only the well equations are picked.
+        std::vector<int> indices(5, -1);
+        int next = 0;
+
+        variableStateWellIndices(indices, next);
+
+        assert(next == 2);
+        return indices;
+    }
+
 }