39 #ifndef limitedLinear_H
40 #define limitedLinear_H
53 template<
class LimiterFunc>
70 <<
"coefficient = " << k_
71 <<
" should be >= 0 and <= 1"
76 twoByk_ = 2.0/
max(k_, SMALL);
81 const scalar cdWeight,
82 const scalar faceFlux,
83 const typename LimiterFunc::phiType& phiP,
84 const typename LimiterFunc::phiType& phiN,
85 const typename LimiterFunc::gradPhiType& gradcP,
86 const typename LimiterFunc::gradPhiType& gradcN,
90 scalar r = LimiterFunc::r
92 faceFlux, phiP, phiN, gradcP, gradcN, d
95 return max(
min(twoByk_*r, 1), 0);