diff --git a/dumux/implicit/problem.hh b/dumux/implicit/problem.hh index ec5788175ec08f06b723acc97370d81ef28c9a36..460968ce5d7dbed8d3b1cbe84dc4921fd12312d5 100644 --- a/dumux/implicit/problem.hh +++ b/dumux/implicit/problem.hh @@ -536,6 +536,32 @@ public: "a initialAtPos() method."); } + /*! + * \brief Evaluate the initial phase state inside a control volume. + * + * \param vertex The vertex + * \param vIdxGlobal The global index of the vertex + * \param globalPos The global position + */ + int initialPhasePresence(const Vertex &vertex, + int &vIdxGlobal, + const GlobalPosition &globalPos) const + { + // forward to generic interface + return asImp_().initialPhasePresenceAtPos(globalPos); + } + + /*! + * \brief Evaluate the initial value for a control volume. + * + * \param globalPos The global position + */ + void initialPhasePresenceAtPos(const GlobalPosition &globalPos) const + { + DUNE_THROW(Dune::InvalidStateException, + "The problem does not provide a initialPhasePresenceAtPos() method."); + } + /*! * \brief Return how much the domain is extruded at a given sub-control volume. *