50 const auto& geometry = _indicatorF.getBlockGeometry();
52 for (
int i = 0; i < this->getTargetDim(); ++i) {
56 if (!_indicatorF(input)) {
61 geometry.getPhysR(centerOfCircle, input);
65 _indicatorF.getBlockGeometry());
67 std::size_t voxels(0);
70 for (inputTmp[0] = 0; inputTmp[0] < geometry.getNx(); ++inputTmp[0]) {
71 for (inputTmp[1] = 0; inputTmp[1] < geometry.getNy(); ++inputTmp[1]) {
72 if (latticeCircle(inputTmp) && _indicatorF(inputTmp)) {
73 T outputTmp[_f.getTargetDim()];
74 _f(outputTmp, inputTmp);
75 for (
int i = 0; i < this->getTargetDim(); ++i) {
76 output[i] += outputTmp[i];
84 for (
int i = 0; i < this->getTargetDim(); ++i) {