43 #ifndef filteredLinear3_H
44 #define filteredLinear3_H
57 template<
class LimiterFunc>
78 <<
"coefficient = " << k_
79 <<
" should be >= 0 and <= 1"
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,
96 scalar df = phiN - phiP;
99 scalar dP = 2*(d & gradcP);
100 scalar dN = 2*(d & gradcN);
103 scalar
limiter = 1 - k_*(dN - df)*(dP - df)/
max(
sqr(dN + dP), SMALL);
106 return max(
min(limiter, 1), 0);