35 Field<T> regionField(cellRegion_().nRegions(),
Zero);
37 forAll(cellRegion_(), celli)
39 regionField[cellRegion_()[celli]] += cellField[celli];
53 const Field<T>& cellField,
58 const Field<T> summedField(regionSum(cellField));
70 label nlb2 = cellRegion_().nRegions()/2;
74 for (label j=0; j<nlb2; j++)
80 - regionField[cellRegion_().nRegions() - j - 1]
86 for (label j=0; j<nlb2; j++)
92 + regionField[cellRegion_().nRegions() - j - 1]
97 regionField.setSize(nlb2);