/// 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()); }