Skip to content
Snippets Groups Projects
Commit 6cae5a88 authored by Robert Kloefkorn's avatar Robert Kloefkorn
Browse files

NoFlux class for elliptic problems.

parent 8453c255
No related branches found
No related tags found
No related merge requests found
fluxesdir = $(includedir)/dune/fem-dg/operator/fluxes
fluxes_HEADERS = averageflux.hh diffusionflux.hh ldgflux.hh \
mhd_fluxes.hh rotator.hh dgprimalfluxes.hh eulerfluxes.hh \
mhd_eqns.hh mhdfluxes.hh mhd_eqns.cc mhd_fluxes.cc
mhd_eqns.hh mhdfluxes.hh mhd_eqns.cc mhd_fluxes.cc noflux.hh
include $(top_srcdir)/am/global-rules
#ifndef DUNE_FEMDG_NOFLUX_HH
#define DUNE_FEMDG_NOFLUX_HH
#include <string>
namespace Dune
{
template <class ModelType>
class NoFlux {
public:
typedef ModelType Model;
typedef typename Model::Traits Traits;
enum { dimRange = Model::dimRange };
typedef typename Model :: DomainType DomainType;
typedef typename Model :: RangeType RangeType;
typedef typename Model :: FluxRangeType FluxRangeType;
typedef typename Model :: FaceDomainType FaceDomainType;
typedef typename Model :: EntityType EntityType;
typedef typename Model :: IntersectionType IntersectionType;
public:
/**
* @brief constructor
*/
NoFlux(const Model& mod) : model_(mod) {}
static std::string name () { return "UpwindFlux"; }
const Model& model() const {return model_;}
/**
* @brief evaluates the flux \f$g(u,v)\f$
*
* @return maximum wavespeed * normal
*/
template <class QuadratureImp>
inline double numericalFlux( const IntersectionType& it,
const EntityType& inside,
const EntityType& outside,
const double time,
const QuadratureImp& faceQuadInner,
const QuadratureImp& faceQuadOuter,
const int quadPoint,
const RangeType& uLeft,
const RangeType& uRight,
RangeType& gLeft,
RangeType& gRight ) const
{
gLeft = 0;
gRight = 0;
return 0;
}
protected:
const Model& model_;
};
}
#endif
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment