49 const auto& geometry = _indicatorF.getBlockGeometry();
51 for (
int i = 0; i < this->getTargetDim(); ++i) {
55 if (!_indicatorF(input)) {
60 geometry.getPhysR(centerOfSphere, input);
64 _indicatorF.getBlockGeometry());
66 std::size_t voxels(0);
69 for (inputTmp[0] = 0; inputTmp[0] < geometry.getNx(); ++inputTmp[0]) {
70 for (inputTmp[1] = 0; inputTmp[1] < geometry.getNy(); ++inputTmp[1]) {
71 for (inputTmp[2] = 0; inputTmp[2] < geometry.getNz(); ++inputTmp[2]) {
72 if (latticeSphere(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];
85 for (
int i = 0; i < this->getTargetDim(); ++i) {