/// This method rescales each statistic and combines all the results
void MultiGrid2D::reduceStatistics(){
std::vector<BlockStatistics*> levelStatistics(this->getNumLevels());
std::vector<int> dimensionsX, dimensionsT;
dimensionsX = statsSubscriber.getDimensionsX();
dimensionsT = statsSubscriber.getDimensionsT();
for (plint iLevel=0; iLevel<this->getNumLevels(); ++iLevel) {
int dxScale = this->getReferenceLevel() - iLevel;
int dtScale = dxScale; // TODO: here, we assume convective scaling; general case should be considered.
std::vector<double> scales(dimensionsX.size());
for (pluint iScale=0; iScale<scales.size(); ++iScale) {
scales[iScale] = scaleToReference(dxScale, dimensionsX[iScale], dtScale, dimensionsT[iScale]);
}
// copy the statistics
levelStatistics[iLevel] = new BlockStatistics(getComponent(iLevel).getInternalStatistics());
// rescale the statistics to the reference level
levelStatistics[iLevel]->rescale(scales);
}
combine(levelStatistics, getInternalStatistics() );
internalStatistics.incrementStats();
for (plint iLevel=0; iLevel<this->getNumLevels(); ++iLevel){
delete levelStatistics[iLevel];
}
}
void executeDataProcessor( ReductiveDataProcessorGenerator3D& generator,
std::vector<MultiGrid3D*> multiGrids,
plint referenceLevel )
{
if(multiGrids.empty()) return;
plint numLevels = (plint)multiGrids[0]->getNumLevels();
std::vector<int> dimensionsX, dimensionsT;
generator.getDimensionsX(dimensionsX);
generator.getDimensionsT(dimensionsT);
std::vector<ReductiveDataProcessorGenerator3D*> localGenerators(numLevels);
std::vector<BlockStatistics*> localStatistics(numLevels);
for (plint iLevel=0; iLevel<numLevels; ++iLevel) {
int dxScale = (int)referenceLevel - (int)iLevel;
int dtScale = dxScale; // TODO: here, we assume convective scaling; general case could be considered.
localGenerators[iLevel] = generator.clone();
plint boxRescaleFactor = util::roundToInt(util::twoToThePowerPlint(std::abs(referenceLevel-iLevel)));
if (dxScale < 0)
generator.divide(boxRescaleFactor);
else
generator.multiply(boxRescaleFactor);
std::vector<MultiBlock3D*> localBlocks(multiGrids.size());
for (plint iBlock=0; iBlock<(plint)localBlocks.size(); ++iBlock) {
localBlocks[iBlock] = &multiGrids[iBlock]->getComponent(iLevel);
}
executeDataProcessor(*localGenerators[iLevel], localBlocks);
std::vector<double> scales(dimensionsX.size());
for (pluint iScale=0; iScale<scales.size(); ++iScale) {
scales[iScale] = scaleToReference(dxScale, dimensionsX[iScale], dtScale, dimensionsT[iScale]);
}
localGenerators[iLevel]->getStatistics().rescale(scales);
localStatistics[iLevel] = &(localGenerators[iLevel]->getStatistics());
}
combine(localStatistics, generator.getStatistics());
}
