33 Field<T> regionField(cellRegion_().nRegions(), pTraits<T>::zero);
35 forAll(cellRegion_(), cellI)
37 regionField[cellRegion_()[cellI]] += cellField[cellI];
51 const Field<T>& cellField,
56 const Field<T> summedField(regionSum(cellField));
68 label nlb2 = cellRegion_().nRegions()/2;
72 for (label j=0; j<nlb2; j++)
78 - regionField[cellRegion_().nRegions() - j - 1]
84 for (label j=0; j<nlb2; j++)
90 + regionField[cellRegion_().nRegions() - j - 1]
95 regionField.setSize(nlb2);