46 S xDist = input[0] - this->_myMin[0];
47 S yDist = input[1] - this->_myMin[1];
48 S zDist = input[2] - this->_myMin[2];
50 int x = ((this->_myMin[0] + xDist)/_deltaR)+0.5;
51 int y = ((this->_myMin[1] + yDist)/_deltaR)+0.5;
52 int z = ((this->_myMin[2] + zDist)/_deltaR)+0.5;
54 if (x >= 0 && x < _blockData.getNx() && y >= 0 && y < _blockData.getNy() && z >= 0 && z < _blockData.getNz()) {
56 if (this->_blockData.get(input) > std::numeric_limits<S>::epsilon()) {
58 output[0] = S(this->_blockData.get(input));
70 output[0] = 1.-S(this->_blockData.get(input));
IndicatorBlockData3D(BlockData< 3, S, S > &blockData, Vector< S, 3 > extend, Vector< S, 3 > origin, S deltaR, bool invert=false)
constructor