overlaps_alignedplane_alignedbox.h

// Description: Overlap tests between an aligned plane and an aligned box

#ifndef PASTELGEOMETRY_OVERLAPS_ALIGNEDPLANE_ALIGNEDBOX_H
#define PASTELGEOMETRY_OVERLAPS_ALIGNEDPLANE_ALIGNEDBOX_H

#include "pastel/geometry/shape/alignedplane.h"
#include "pastel/geometry/shape/alignedbox.h"

namespace Pastel
{

    //! Tests if an aligned plane and an aligned box overlap.
    /*!
   Preconditions:
   plane.n() == box.n()

   Time complexity:
   O(1)
   */
    template <typename Real, int N>
    bool overlaps(
        const AlignedPlane<Real, N>& plane,
        const AlignedBox<Real, N>& box)
    {
        PENSURE_OP(plane.n(), ==, box.n());

        if (plane.position() >= box.max()[plane.axis()])
        {
            if (plane.position > box.max()[plane.axis()] ||
                box.maxTopology()[plane.axis()] == Topology::Open)
            {
                return false;
            }
        }

        if (plane.position() <= box.min()[plane.axis()])
        {
            if (plane.position() < box.min()[plane.axis()] ||
                box.minTopology()[plane.axis()] == Topology::Open)
            {
                return false;
            }
        }

        return true;
    }

    //! Tests if an aligned plane and an aligned box overlap.
    /*!
   boxOnPositiveSide:
   Filled with the information if the minimal point of the 
   aligned box is on the positive side of the plane.

   Returns:
   overlaps(plane, box)
   */
    template <typename Real, int N>
    bool overlaps(
        const AlignedPlane<Real, N>& plane,
        const AlignedBox<Real, N>& box,
        bool &boxOnPositiveSide)
    {
        boxOnPositiveSide =
            box.min()[plane.axis()] > plane.position();

        return overlaps(plane, box);
    }

}

#endif