50 const auto& blockGeometry = _indicatorF.getBlockGeometry();
51 const int nX = blockGeometry.getNx();
52 const int nY = blockGeometry.getNy();
53 const int nZ = blockGeometry.getNz();
54 const T weight =
util::pow(blockGeometry.getDeltaR(), 3);
57 W outputTmp[_f.getTargetDim()];
58 int inputTmp[_f.getSourceDim()];
60 for (inputTmp[0] = 0; inputTmp[0] < nX; ++inputTmp[0]) {
61 for (inputTmp[1] = 0; inputTmp[1] < nY; ++inputTmp[1]) {
62 for (inputTmp[2] = 0; inputTmp[2] < nZ; ++inputTmp[2]) {
63 if (_indicatorF(inputTmp)) {
64 _f(outputTmp, inputTmp);
65 for (
int iDim = 0; iDim < _f.getTargetDim(); ++iDim) {