40 #ifndef cellMDLimitedGrad_H
41 #define cellMDLimitedGrad_H
100 "cellMDLimitedGrad(const fvMesh&, Istream& schemeData)",
102 ) <<
"coefficient = " << k_
103 <<
" should be >= 0 and <= 1"
114 const Type& maxDelta,
115 const Type& minDelta,
136 const scalar& maxDelta,
137 const scalar& minDelta,
141 scalar extrapolate = dcf & g;
143 if (extrapolate > maxDelta)
145 g = g + dcf*(maxDelta - extrapolate)/
magSqr(dcf);
147 else if (extrapolate < minDelta)
149 g = g + dcf*(minDelta - extrapolate)/
magSqr(dcf);
158 const Type& maxDelta,
159 const Type& minDelta,
163 for(
direction cmpt=0; cmpt<Type::nComponents; cmpt++)
165 vector gi(g[cmpt], g[cmpt+3], g[cmpt+6]);
170 minDelta.component(cmpt),