throw when calling addWellContributions() for MSW

opm/simulators/wells/MultisegmentWell.hpp

@@ -152,6 +152,8 @@ namespace Opm
                                                  const WellState& well_state,
                                                  Opm::DeferredLogger& deferred_logger) override; // should be const?
 
+        virtual void  addWellContributions(Mat& mat) const override;
+
         /// number of segments for this well
         /// int number_of_segments_;
         int numberOfSegments() const;

opm/simulators/wells/MultisegmentWell_impl.hpp

@@ -875,6 +875,18 @@ namespace Opm
 
 
 
+    template<typename TypeTag>
+    void
+    MultisegmentWell<TypeTag>::
+    addWellContributions(Mat& /* mat */) const
+    {
+        OPM_THROW(std::runtime_error, "addWellContributions is not supported by multisegment well yet");
+    }
+
+
+
+
+
     template <typename TypeTag>
     const WellSegments&
     MultisegmentWell<TypeTag>::

opm/simulators/wells/WellInterface.hpp

@@ -208,8 +208,7 @@ namespace Opm
         void calculateReservoirRates(WellState& well_state) const;
 
         // Add well contributions to matrix
-        virtual void addWellContributions(Mat&) const
-        {}
+        virtual void addWellContributions(Mat&) const = 0;
 
         void addCellRates(RateVector& rates, int cellIdx) const;