Commit 83337960 authored by Dennis Gläser's avatar Dennis Gläser Committed by Kilian Weishaupt
Browse files

[veloutput] use new neumann interface

parent 561ed6e6
...@@ -29,6 +29,7 @@ ...@@ -29,6 +29,7 @@
#include <dune/geometry/referenceelements.hh> #include <dune/geometry/referenceelements.hh>
#include <dumux/common/parameters.hh> #include <dumux/common/parameters.hh>
#include <dumux/common/deprecated.hh>
#include <dumux/io/velocityoutput.hh> #include <dumux/io/velocityoutput.hh>
#include <dumux/discretization/method.hh> #include <dumux/discretization/method.hh>
#include <dumux/discretization/elementsolution.hh> #include <dumux/discretization/elementsolution.hh>
...@@ -310,7 +311,7 @@ public: ...@@ -310,7 +311,7 @@ public:
if (bcTypes.hasNeumann()) if (bcTypes.hasNeumann())
{ {
// check if we have Neumann no flow, we can just use 0 // check if we have Neumann no flow, we can just use 0
const auto neumannFlux = problem_.neumann(element, fvGeometry, elemVolVars, scvf); const auto neumannFlux = Deprecated::neumann(problem_, element, fvGeometry, elemVolVars, elemFluxVarsCache, scvf);
using NumEqVector = std::decay_t<decltype(neumannFlux)>; using NumEqVector = std::decay_t<decltype(neumannFlux)>;
if (Dune::FloatCmp::eq<NumEqVector, Dune::FloatCmp::CmpStyle::absolute>(neumannFlux, NumEqVector(0.0), 1e-30)) if (Dune::FloatCmp::eq<NumEqVector, Dune::FloatCmp::CmpStyle::absolute>(neumannFlux, NumEqVector(0.0), 1e-30))
scvfFluxes[scvfIndexInInside[localScvfIdx]] = 0; scvfFluxes[scvfIndexInInside[localScvfIdx]] = 0;
......
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