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

[initindicator] use new neumann interface

parent c4b29ffe
...@@ -28,6 +28,7 @@ ...@@ -28,6 +28,7 @@
#include <dune/geometry/type.hh> #include <dune/geometry/type.hh>
#include <dumux/common/properties.hh> #include <dumux/common/properties.hh>
#include <dumux/common/parameters.hh> #include <dumux/common/parameters.hh>
#include <dumux/common/deprecated.hh>
#include <dumux/discretization/method.hh> #include <dumux/discretization/method.hh>
namespace Dumux { namespace Dumux {
...@@ -163,10 +164,14 @@ public: ...@@ -163,10 +164,14 @@ public:
// get the fvGeometry and elementVolVars needed for the bc and source interfaces // get the fvGeometry and elementVolVars needed for the bc and source interfaces
auto fvGeometry = localView(*fvGridGeometry_); auto fvGeometry = localView(*fvGridGeometry_);
fvGeometry.bindElement(element); fvGeometry.bind(element);
auto elemVolVars = localView(gridVariables_->curGridVolVars()); auto elemVolVars = localView(gridVariables_->curGridVolVars());
elemVolVars.bindElement(element, fvGeometry, sol); elemVolVars.bind(element, fvGeometry, sol);
// elemFluxVarsCache for neumann interface
auto elemFluxVarsCache = localView(gridVariables_->gridFluxVarsCache());
elemFluxVarsCache.bind(element, fvGeometry, elemVolVars);
//! Check if we have to refine around a source term //! Check if we have to refine around a source term
if (refineAtSource_) if (refineAtSource_)
...@@ -208,7 +213,7 @@ public: ...@@ -208,7 +213,7 @@ public:
// we are on a pure Neumann boundary // we are on a pure Neumann boundary
else if(refineAtFluxBC_) else if(refineAtFluxBC_)
{ {
const auto fluxes = problem_->neumann(element, fvGeometry, elemVolVars, scvf); const auto fluxes = Deprecated::neumann(*problem_, element, fvGeometry, elemVolVars, elemFluxVarsCache, scvf);
if (fluxes.infinity_norm() > eps_) if (fluxes.infinity_norm() > eps_)
{ {
indicatorVector_[eIdx] = true; indicatorVector_[eIdx] = true;
...@@ -242,7 +247,7 @@ public: ...@@ -242,7 +247,7 @@ public:
//! check if scvf is on Neumann boundary //! check if scvf is on Neumann boundary
if (scvf.boundary() && bcTypes[scvf.insideScvIdx()].hasNeumann()) if (scvf.boundary() && bcTypes[scvf.insideScvIdx()].hasNeumann())
{ {
const auto fluxes = problem_->neumann(element, fvGeometry, elemVolVars, scvf); const auto fluxes = Deprecated::neumann(*problem_, element, fvGeometry, elemVolVars, elemFluxVarsCache, scvf);
if (fluxes.infinity_norm() > eps_) if (fluxes.infinity_norm() > eps_)
{ {
indicatorVector_[eIdx] = true; indicatorVector_[eIdx] = true;
......
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