AABB.hpp

Go to the documentation of this file.

/**************************************************************************************/
/*                                                                                    */
/*  Visualization Library                                                             */
/*  http://visualizationlibrary.org                                                   */
/*                                                                                    */
/*  Copyright (c) 2005-2020, Michele Bosi                                             */
/*  All rights reserved.                                                              */
/*                                                                                    */
/*  Redistribution and use in source and binary forms, with or without modification,  */
/*  are permitted provided that the following conditions are met:                     */
/*                                                                                    */
/*  - Redistributions of source code must retain the above copyright notice, this     */
/*  list of conditions and the following disclaimer.                                  */
/*                                                                                    */
/*  - Redistributions in binary form must reproduce the above copyright notice, this  */
/*  list of conditions and the following disclaimer in the documentation and/or       */
/*  other materials provided with the distribution.                                   */
/*                                                                                    */
/*  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND   */
/*  ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED     */
/*  WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE            */
/*  DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR  */
/*  ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES    */
/*  (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;      */
/*  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON    */
/*  ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT           */
/*  (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS     */
/*  SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.                      */
/*                                                                                    */
/**************************************************************************************/

#ifndef AABB_INCLUDE_ONCE
#define AABB_INCLUDE_ONCE

#include <vlCore/Vector3.hpp>
#include <vlCore/Matrix4.hpp>

namespace vl
{
  //-----------------------------------------------------------------------------
  // AABB
  //-----------------------------------------------------------------------------
  class VLCORE_EXPORT AABB
  {
  public:
    AABB();

    AABB( const vec3& center, real radius );

    AABB( const vec3& pt1, const vec3& pt2, real displace=0);

    void setNull() { mMin = 1; mMax = -1; }

    bool isNull()  const { return mMin.x() > mMax.x() || mMin.y() > mMax.y() || mMin.z() > mMax.z(); }

    bool isPoint() const { return mMin == mMax; }

    void enlarge(real displace);

    bool intersects(const AABB & bb) const;

    vec3 clip(const vec3& p, bool clipx=true, bool clipy=true, bool clipz=true) const;

    bool isInside(const vec3& p, bool clipx, bool clipy, bool clipz) const;

    bool isInside(const vec3& p) const;

    real width() const;

    real height() const;

    real depth() const;

    bool operator==(const AABB& aabb) const
    {
      return mMin == aabb.mMin && mMax == aabb.mMax;
    }

    bool operator!=(const AABB& aabb) const
    {
      return !operator==(aabb);
    }

    AABB operator+(const AABB& aabb) const;

    AABB& operator+=(const AABB& other)
    {
      *this = *this + other;
      return *this;
    }

    AABB operator+(const vec3& p)
    {
      AABB aabb = *this;
      aabb += p;
      return aabb;
    }

    const AABB& operator+=(const vec3& p)
    {
      addPoint(p);
      return *this;
    }

    vec3 center() const;

    real longestSideLength() const
    {
      real side = width();
      if (height() > side)
        side = height();
      if (depth() > side)
        side = depth();
      return side;
    }

    void addPoint(const vec3& p, real radius);

    void addPoint(const vec3& p)
    {
      if (isNull())
      {
        mMax = p;
        mMin = p;
        return;
      }

      if ( mMax.x() < p.x() ) mMax.x() = p.x();
      if ( mMax.y() < p.y() ) mMax.y() = p.y();
      if ( mMax.z() < p.z() ) mMax.z() = p.z();
      if ( mMin.x() > p.x() ) mMin.x() = p.x();
      if ( mMin.y() > p.y() ) mMin.y() = p.y();
      if ( mMin.z() > p.z() ) mMin.z() = p.z();
    }

    void transformed(AABB& out, const mat4& mat) const
    {
      out.setNull();
      if ( !isNull() )
      {
        out.addPoint( mat * vec3(minCorner().x(), minCorner().y(), minCorner().z()) );
        out.addPoint( mat * vec3(minCorner().x(), maxCorner().y(), minCorner().z()) );
        out.addPoint( mat * vec3(maxCorner().x(), maxCorner().y(), minCorner().z()) );
        out.addPoint( mat * vec3(maxCorner().x(), minCorner().y(), minCorner().z()) );
        out.addPoint( mat * vec3(minCorner().x(), minCorner().y(), maxCorner().z()) );
        out.addPoint( mat * vec3(minCorner().x(), maxCorner().y(), maxCorner().z()) );
        out.addPoint( mat * vec3(maxCorner().x(), maxCorner().y(), maxCorner().z()) );
        out.addPoint( mat * vec3(maxCorner().x(), minCorner().y(), maxCorner().z()) );
      }
    }

    AABB transformed(const mat4& mat) const
    {
      AABB aabb;
      transformed(aabb, mat);
      return aabb;
    }

    const vec3& minCorner() const { return mMin; }

    const vec3& maxCorner() const { return mMax; }

    void setMinCorner(real x, real y, real z) { mMin.x() = x; mMin.y() = y; mMin.z() = z; }

    void setMinCorner(const vec3& v) { mMin = v; }

    void setMaxCorner(real x, real y, real z) { mMax.x() = x; mMax.y() = y; mMax.z() = z; }

    void setMaxCorner(const vec3& v) { mMax = v; }

    real volume() const { return width() * height() * depth(); }

  protected:
    vec3 mMin;
    vec3 mMax;
  };
}

#endif