Commit a78a2590 authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[facet][box]darcyslaw] throw error if xi != 1.0 wants to be used

parent 1dab0494
......@@ -26,6 +26,7 @@
#include <dune/common/exceptions.hh>
#include <dune/common/fvector.hh>
#include <dune/common/float_cmp.hh>
#include <dumux/common/math.hh>
#include <dumux/common/parameters.hh>
......@@ -72,6 +73,10 @@ public:
if (!scvf.interiorBoundary())
return DefaultBoxDarcysLaw::flux(problem, element, fvGeometry, elemVolVars, scvf, phaseIdx, elemFluxVarCache);
static const Scalar xi = getParamFromGroup<Scalar>(problem.paramGroup(), "FacetCoupling.Xi", 1.0);
if ( !Dune::FloatCmp::eq(xi, 1.0, 1e-6) )
DUNE_THROW(Dune::NotImplemented, "Xi != 1.0 cannot be used with the Box-Facet-Coupling scheme");
// get some references for convenience
const auto& fluxVarCache = elemFluxVarCache[scvf];
const auto& shapeValues = fluxVarCache.shapeValues();
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment