inet::SpatialGrid::LineSegmentIterator Class Reference

#include <SpatialGrid.h>

Public Member Functions
	LineSegmentIterator (const SpatialGrid *spatialGrid, const LineSegment &lineSegment, const Triplet< double > &voxelSizes, const Triplet< int > &numVoxels)
LineSegmentIterator &	operator++ ()
const Triplet< int > &	getMatrixIndices () const
bool	end () const

Protected Attributes
Triplet< int >	index
Triplet< int >	endPoint
Triplet< int >	step
Triplet< double >	tDelta
Triplet< double >	tExit
bool	reachedEnd

Constructor & Destructor Documentation

inet::SpatialGrid::LineSegmentIterator::LineSegmentIterator	(	const SpatialGrid *	spatialGrid,
		const LineSegment &	lineSegment,
		const Triplet< double > &	voxelSizes,
		const Triplet< int > &	numVoxels
	)

{
    reachedEnd = false;
    Coord p0 = lineSegment.getPoint1();
    Coord p1 = lineSegment.getPoint2();
    Coord segmentDirection = p1 - p0;
    Triplet<double> point0 = Triplet<double>(p0.x - spatialGrid->constraintAreaMin.x, p0.y - spatialGrid->constraintAreaMin.y, p0.z - spatialGrid->constraintAreaMin.z);
    Triplet<double> direction = Triplet<double>(segmentDirection.x, segmentDirection.y, segmentDirection.z);
    index = spatialGrid->coordToMatrixIndices(p0);
    endPoint = spatialGrid->coordToMatrixIndices(p1);
    for (int i = 0; i < 3; i++)
    {
        tDelta[i] = voxelSizes[i] / std::abs(direction[i]);
        double ithdirection = direction[i];
        if (ithdirection > 0)
            step[i] = 1;
        else if (ithdirection < 0)
            step[i] = -1;
        else
            step[i] = 0;
        if (ithdirection != 0) {
            double d = point0[i] - index[i] * voxelSizes[i];
            if (step[i] > 0)
                d = voxelSizes[i] - d;
            ASSERT(d >= 0 && d <= voxelSizes[i]);
            tExit[i] = d / std::abs(ithdirection);
        }
        else
            tExit[i] = std::numeric_limits<double>::max();
    }
}

Member Function Documentation

bool inet::SpatialGrid::LineSegmentIterator::end ( ) const

inline

Referenced by inet::SpatialGrid::lineSegmentQuery().

{ return reachedEnd; }

const Triplet<int>& inet::SpatialGrid::LineSegmentIterator::getMatrixIndices ( ) const

inline

{ return index; }

SpatialGrid::LineSegmentIterator & inet::SpatialGrid::LineSegmentIterator::operator++

(

)

{
    if (index.x != endPoint.x || index.y != endPoint.y || index.z != endPoint.z)
    {
        int axis = 0;
        if (tExit.x < tExit.y)
        {
            if (tExit.x < tExit.z)
                axis = 0;
            else
                axis = 2;
        }
        else if (tExit.y < tExit.z)
            axis = 1;
        else
            axis = 2;
        index[axis] += step[axis];
        tExit[axis] += tDelta[axis];
    }
    else
        reachedEnd = true;
    return *this;
}

Member Data Documentation

Triplet<int> inet::SpatialGrid::LineSegmentIterator::endPoint

protected

Triplet<int> inet::SpatialGrid::LineSegmentIterator::index

protected

bool inet::SpatialGrid::LineSegmentIterator::reachedEnd

protected

Triplet<int> inet::SpatialGrid::LineSegmentIterator::step

protected

Triplet<double> inet::SpatialGrid::LineSegmentIterator::tDelta

protected

Triplet<double> inet::SpatialGrid::LineSegmentIterator::tExit

protected

---

The documentation for this class was generated from the following files:

• SpatialGrid.h
• SpatialGrid.cc