58template<
class LimiterFunc>
74 <<
"coefficient = " << k_
75 <<
" should be >= 0 and <= 1"
81 k_ =
max(k_/2.0, SMALL);
86 const scalar cdWeight,
87 const scalar faceFlux,
88 const typename LimiterFunc::phiType& phiP,
89 const typename LimiterFunc::phiType& phiN,
90 const typename LimiterFunc::gradPhiType& gradcP,
91 const typename LimiterFunc::gradPhiType& gradcN,
95 scalar phict = LimiterFunc::phict
97 faceFlux, phiP, phiN, gradcP, gradcN, d
100 return min(
max(phict/k_, 0), 1);
Class with limiter function which returns the limiter for the Gamma differencing scheme based on phic...
GammaLimiter(Istream &is)
scalar limiter(const scalar cdWeight, const scalar faceFlux, const typename LimiterFunc::phiType &phiP, const typename LimiterFunc::phiType &phiN, const typename LimiterFunc::gradPhiType &gradcP, const typename LimiterFunc::gradPhiType &gradcN, const vector &d) const
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
errorManipArg< error, int > exit(error &err, const int errNo=1)