110 label sizeby2 = this->size()/2;
112 if (this->doTransform())
114 for (label facei=0; facei<sizeby2; facei++)
118 this->forwardT()[0], iField[faceCells[facei + sizeby2]]
123 this->reverseT()[0], iField[faceCells[facei]] + jf[facei]
129 for (label facei=0; facei<sizeby2; facei++)
131 pnf[facei] = iField[faceCells[facei + sizeby2]] - jf[facei];
132 pnf[facei + sizeby2] = iField[faceCells[facei]] + jf[facei];
153 label sizeby2 = this->size()/2;
156 if (&psiInternal == &this->internalField())
161 for (label facei=0; facei<sizeby2; facei++)
163 pnf[facei] = psiInternal[faceCells[facei + sizeby2]] - jf[facei];
164 pnf[facei + sizeby2] = psiInternal[faceCells[facei]] + jf[facei];
169 for (label facei=0; facei<sizeby2; facei++)
171 pnf[facei] = psiInternal[faceCells[facei + sizeby2]];
172 pnf[facei + sizeby2] = psiInternal[faceCells[facei]];
177 this->transformCoupleField(pnf, cmpt);
182 result[faceCells[elemI]] -= coeffs[elemI]*pnf[elemI];