53 template<
class LimiterFunc>
69 <<
"coefficient = " << k_
70 <<
" should be >= 0 and <= 1"
76 k_ =
max(k_/2.0, 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 phict = LimiterFunc::phict
92 faceFlux, phiP, phiN, gradcP, gradcN, d
95 return min(
max(phict/k_, 0), 1);