diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimetic.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimetic.hh index 76e7c90d00671d284d80b3997b954307a57cfb0d..7c98f76ab66ea22aa42721e24d01bcdd881d56cb 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimetic.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimetic.hh @@ -361,6 +361,8 @@ private: case Sn: sat = cellData.saturation(nPhaseIdx); break; + default: + DUNE_THROW(Dune::NotImplemented, "Only saturation formulation Sw and Sn are implemented!"); } Scalar error = (sat > 1.0) ? sat - 1.0 : 0.0; @@ -641,6 +643,8 @@ void MimeticTwoPLocalStiffness<TypeTag>::assembleElementMatrices(const Element& fracFlow = -cellData.fracFlowFunc(wPhaseIdx); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } @@ -669,6 +673,8 @@ void MimeticTwoPLocalStiffness<TypeTag>::assembleElementMatrices(const Element& fracFlow = -cellData.fracFlowFunc(wPhaseIdx); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } } else if (flux < 0. && bctype.isDirichlet(satEqIdx)) @@ -695,6 +701,8 @@ void MimeticTwoPLocalStiffness<TypeTag>::assembleElementMatrices(const Element& / (krw / viscosity_[wPhaseIdx] + krn / viscosity_[nPhaseIdx]); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } } diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimeticadaptive.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimeticadaptive.hh index 896605001e23cba58013eb5495d2c609f3e90fbb..5bf9b3bcc7a74f3d6db7f3f55a7f1ff94e172835 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimeticadaptive.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/mimeticadaptive.hh @@ -377,6 +377,8 @@ private: case Sn: sat = cellData.saturation(nPhaseIdx); break; + default: + DUNE_THROW(Dune::NotImplemented, "Only saturation formulation Sw and Sn are implemented!"); } Scalar error = (sat > 1.0) ? sat - 1.0 : 0.0; @@ -661,6 +663,8 @@ void MimeticTwoPLocalStiffnessAdaptive<TypeTag>::assembleElementMatrices(const E fracFlow = -cellData.fracFlowFunc(wPhaseIdx); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } rhs_[eIdxGlobal] -= (faceVol[idx] * fracFlow * flux); @@ -688,6 +692,8 @@ void MimeticTwoPLocalStiffnessAdaptive<TypeTag>::assembleElementMatrices(const E fracFlow = -cellData.fracFlowFunc(wPhaseIdx); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } } else if (flux < 0. && bctype.isDirichlet(satEqIdx)) @@ -714,6 +720,8 @@ void MimeticTwoPLocalStiffnessAdaptive<TypeTag>::assembleElementMatrices(const E / (krw / viscosity_[wPhaseIdx] + krn / viscosity_[nPhaseIdx]); break; } + default: + DUNE_THROW(Dune::NotImplemented, "Only pressure formulation pw and pn are implemented!"); } } diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh index 4160138ba7ceb1a62381727359771ea4f90d4716..d23db57673491dbc5bb674b5d3acf4ad5b9655d2 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh @@ -24,6 +24,8 @@ #ifndef DUMUX_MIMETICPRESSURE2P_HH #define DUMUX_MIMETICPRESSURE2P_HH +#include <dune/common/exceptions.hh> + // dumux environment #include <dumux/porousmediumflow/sequential/mimetic/properties.hh> #include <dumux/porousmediumflow/sequential/cellcentered/pressure.hh> @@ -136,6 +138,8 @@ template<class TypeTag> class MimeticPressure2P case sn: sat = problem_.variables().cellData(i).saturation(nPhaseIdx); break; + default: + DUNE_THROW(Dune::NotImplemented, "Only saturation formulation sw and sn are implemented!"); } if (sat > 1.0) { diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh index a774ed7b1019b4cfc1a516a35195eca3b6a0f887..948bc3dc83df5970d420757467fb931838695cd7 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh @@ -24,6 +24,8 @@ #ifndef DUMUX_MIMETICPRESSURE2PADAPTIVE_HH #define DUMUX_MIMETICPRESSURE2PADAPTIVE_HH +#include <dune/common/exceptions.hh> + // dumux environment #include <dumux/porousmediumflow/sequential/mimetic/properties.hh> #include <dumux/porousmediumflow/sequential/cellcentered/pressure.hh> @@ -137,6 +139,8 @@ template<class TypeTag> class MimeticPressure2PAdaptive case Sn: sat = problem_.variables().cellData(i).saturation(nPhaseIdx); break; + default: + DUNE_THROW(Dune::NotImplemented, "Only saturation formulation Sw and Sn are implemented!"); } if (sat > 1.0) {