3#ifndef DUNE_PDELAB_LOCALOPERATOR_DGINTERIORPENALTYPARAMETER_HH
4#define DUNE_PDELAB_LOCALOPERATOR_DGINTERIORPENALTYPARAMETER_HH
6#include <dune/common/parametertreeparser.hh>
29 sigma = config.get<RF>(
"sigma");
32 template<
typename GEO,
typename IGEO,
typename OGEO>
33 RF
getFaceIP(
const GEO& geo,
const IGEO& igeo,
const OGEO& ogeo)
const
37 return sigma * geo.volume()/min(igeo.volume(), ogeo.volume());
40 template<
typename GEO,
typename IGEO>
41 RF
getFaceIP(
const GEO& geo,
const IGEO& igeo)
const
44 return sigma * geo.volume()/igeo.volume();
65 sigma = config.get<RF>(
"sigma");
66 beta = config.get<RF>(
"beta");
69 template<
typename GEO,
typename IGEO,
typename OGEO>
70 RF
getFaceIP(
const GEO& geo,
const IGEO& igeo,
const OGEO& ogeo)
const
75 return sigma * pow(geo.volume()/min(igeo.volume(),ogeo.volume()), beta);
78 template<
typename GEO,
typename IGEO>
79 RF
getFaceIP(
const GEO& geo,
const IGEO& igeo)
const
83 return sigma * pow(geo.volume()/igeo.volume(), beta);
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
Default implementation of the interior penalty factor.
Definition: dginteriorpenaltyparameter.hh:23
RF getFaceIP(const GEO &geo, const IGEO &igeo) const
Definition: dginteriorpenaltyparameter.hh:41
DefaultInteriorPenalty(const Dune::ParameterTree &config)
Definition: dginteriorpenaltyparameter.hh:27
RF getFaceIP(const GEO &geo, const IGEO &igeo, const OGEO &ogeo) const
Definition: dginteriorpenaltyparameter.hh:33
Implementation of overpenalized interior penalty.
Definition: dginteriorpenaltyparameter.hh:58
RF getFaceIP(const GEO &geo, const IGEO &igeo) const
Definition: dginteriorpenaltyparameter.hh:79
RF getFaceIP(const GEO &geo, const IGEO &igeo, const OGEO &ogeo) const
Definition: dginteriorpenaltyparameter.hh:70
OverPenalizedInteriorPenalty(const Dune::ParameterTree &config)
Definition: dginteriorpenaltyparameter.hh:63