49 label sizeby2 = deltas.
size()/2;
51 for (label facei = 0; facei < sizeby2; facei++)
53 scalar avFa = (magFa[facei] + magFa[facei + sizeby2])/2.0;
55 if (
mag(magFa[facei] - magFa[facei + sizeby2])/avFa > 1
e-4)
57 FatalErrorIn(
"cyclicFvPatch::makeWeights(scalarField& w) const")
58 <<
"face " << facei <<
" and " << facei + sizeby2
59 <<
" areas do not match by "
60 << 100*
mag(magFa[facei] - magFa[facei + sizeby2])/avFa
61 <<
"% -- possible face ordering problem"
65 scalar di = deltas[facei];
66 scalar dni = deltas[facei + sizeby2];
68 w[facei] = dni/(di + dni);
69 w[facei + sizeby2] = 1 - w[facei];
78 label sizeby2 = deltas.
size()/2;
80 for (label facei = 0; facei < sizeby2; facei++)
82 scalar di = deltas[facei];
83 scalar dni = deltas[facei + sizeby2];
85 dc[facei] = 1.0/(di + dni);
86 dc[facei + sizeby2] = dc[facei];
95 label sizeby2 = patchD.
size()/2;
103 for (label facei = 0; facei < sizeby2; facei++)
105 vector ddi = patchD[facei];
106 vector dni = patchD[facei + sizeby2];
108 pdv[facei] = ddi - dni;
109 pdv[facei + sizeby2] = -pdv[facei];
114 for (label facei = 0; facei < sizeby2; facei++)
116 vector ddi = patchD[facei];
117 vector dni = patchD[facei + sizeby2];
133 return patchInternalField(internalData);
146 label sizeby2 = this->size()/2;
148 for (label facei=0; facei<sizeby2; facei++)
150 pnf[facei] = interfaceData[facei + sizeby2];
151 pnf[facei + sizeby2] = interfaceData[facei];
169 label sizeby2 = this->size()/2;
171 for (label facei=0; facei<sizeby2; facei++)
173 pnf[facei] = iF[faceCells[facei + sizeby2]];
174 pnf[facei + sizeby2] = iF[faceCells[facei]];