33 const word& derivedType
37 geomDecomDict_(decompositionDict.
subDict(derivedType +
"Coeffs")),
38 n_(geomDecomDict_.lookup(
"n")),
39 delta_(
readScalar(geomDecomDict_.lookup(
"delta"))),
44 if (nProcessors_ != n_.x()*n_.y()*n_.z())
48 "geomDecomp::geomDecomp"
49 "(const dictionary& decompositionDict)"
50 ) <<
"Wrong number of processor divisions in geomDecomp:" <<
nl
51 <<
"Number of domains : " << nProcessors_ <<
nl
52 <<
"Wanted decomposition : " << n_
56 scalar
d = 1 - 0.5*delta_*delta_;
65 a*d - a2*d, a*a2 + d2, -2*a*d,
66 a*d2 + a2, a*d - a2*d, d2 - a2