BBox.h

Go to the documentation of this file.

/*
 * Copyright (C) <2012-2024> <EDF-DTG> <FRANCE>
 * This file is part of Code_TYMPAN (R).
 * Code_TYMPAN (R) is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * Code_TYMPAN (R) is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 * See the GNU General Public License for more details.
 * You should have received a copy of the GNU General Public License along
 * with Code_TYMPAN (R). If not, see <https://www.gnu.org/licenses/>.
 */
 
#ifndef BBOX_H
#define BBOX_H
 
#include <iostream>
#include "Geometry/mathlib.h"
#include <algorithm>
 
// using namespace std;
 
class BBox
{
public:
    // public data
    vec3 pMin;     
    vec3 pMax;     
    vec3 centroid; 
    bool isNull;   
 
    // public methods
    BBox() : isNull(true) {}
 
    BBox(const vec3& _pMin, const vec3& _pMax) : isNull(false)
    {
        pMin = _pMin;
        pMax = _pMax;
        centroid = (pMin + pMax) / 2.0f;
    }
 
    BBox(const BBox& other)
    {
        pMin = vec3(other.pMin);
        pMax = vec3(other.pMax);
        centroid = (pMin + pMax) / 2.0f;
        isNull = other.isNull;
    }
 
    vec3 getBBMin()
    {
        return pMin;
    }
 
    vec3 getBBMax()
    {
        return pMax;
    }
 
    vec3 getCentroid()
    {
        return centroid;
    }
 
    void setBBMin(vec3& _pMin)
    {
        pMin = _pMin;
        centroid = (pMin + pMax) / 2.0f;
    }
 
    void setBBMax(vec3& _pMax)
    {
        pMax = _pMax;
        centroid = (pMin + pMax) / 2.0f;
    }
 
    double SurfaceArea() const
    {
        vec3 d = pMax - pMin;
        return 2 * (d.x * d.y + d.x * d.z + d.y * d.z);
    }
 
    BBox Union(const BBox& b, const vec3& p)
    {
        BBox ret = b;
        ret.pMin.x = std::min(b.pMin.x, p.x);
        ret.pMin.y = std::min(b.pMin.y, p.y);
        ret.pMin.z = std::min(b.pMin.z, p.z);
        ret.pMax.x = std::max(b.pMax.x, p.x);
        ret.pMax.y = std::max(b.pMax.y, p.y);
        ret.pMax.z = std::max(b.pMax.z, p.z);
        ret.centroid = (ret.pMin + ret.pMax) * 0.5;
        return ret;
    }
 
    BBox Union(const vec3& p)
    {
        BBox ret;
        ret.isNull = false;
        ret.pMin.x = std::min(pMin.x, p.x);
        ret.pMin.y = std::min(pMin.y, p.y);
        ret.pMin.z = std::min(pMin.z, p.z);
        ret.pMax.x = std::max(pMax.x, p.x);
        ret.pMax.y = std::max(pMax.y, p.y);
        ret.pMax.z = std::max(pMax.z, p.z);
        ret.centroid = (ret.pMin + ret.pMax) * 0.5;
        return ret;
    }
 
    BBox Union(const BBox& b2)
    {
        BBox ret;
        if (isNull)
        {
            ret = b2;
        }
        else if (b2.isNull)
        {
            ret = *this;
        }
        else
        {
            ret.isNull = false;
            ret.pMin.x = std::min(pMin.x, b2.pMin.x);
            ret.pMin.y = std::min(pMin.y, b2.pMin.y);
            ret.pMin.z = std::min(pMin.z, b2.pMin.z);
            ret.pMax.x = std::max(pMax.x, b2.pMax.x);
            ret.pMax.y = std::max(pMax.y, b2.pMax.y);
            ret.pMax.z = std::max(pMax.z, b2.pMax.z);
            ret.centroid = (ret.pMin + ret.pMax) * 0.5;
        }
        return ret;
    }
 
    BBox Union(const BBox& b1, const BBox& b2)
    {
        BBox ret = b1;
        ret.pMin.x = std::min(b1.pMin.x, b2.pMin.x);
        ret.pMin.y = std::min(b1.pMin.y, b2.pMin.y);
        ret.pMin.z = std::min(b1.pMin.z, b2.pMin.z);
        ret.pMax.x = std::max(b1.pMax.x, b2.pMax.x);
        ret.pMax.y = std::max(b1.pMax.y, b2.pMax.y);
        ret.pMax.z = std::max(b1.pMax.z, b2.pMax.z);
        ret.centroid = (ret.pMin + ret.pMax) * 0.5;
        return ret;
    }
 
    inline bool isInBox(const vec3& p) const
    {
        return p.x >= pMin.x && p.y >= pMin.y && p.z >= pMin.z && p.x <= pMax.x && p.y <= pMax.y &&
               p.z <= pMax.z;
    }
 
    inline bool intersectBox(const BBox& box)
    {
        if (isNull || box.isNull)
        {
            return false;
        }
        bool x = (pMax.x >= box.pMin.x) && (pMin.x <= box.pMax.x);
        bool y = (pMax.y >= box.pMin.y) && (pMin.y <= box.pMax.y);
        bool z = (pMax.z >= box.pMin.z) && (pMin.z <= box.pMax.z);
        return (x && y && z);
    }
 
    int MaximumExtend() const
    {
        vec3 diag = pMax - pMin;
        if (diag.x > diag.y && diag.x > diag.z)
        {
            return 0;
        }
        else if (diag.y > diag.z)
        {
            return 1;
        }
        else
        {
            return 2;
        }
    }
 
    bool IntersectP(vec3 rayPos, vec3 rayDir, decimal* hitt0 = NULL, decimal* hitt1 = NULL) const
    {
        if (isNull)
        {
            return false;
        }
 
        // float t0 = ray->mint; float t1 = ray->maxt;
        decimal t0 = 0, t1 = 1000000;
        decimal invRayDir, tNear, tFar;
 
        // Composante sur X
        // TODO : If denominator is Zero
        invRayDir = 1 / rayDir.x;
        tNear = (pMin.x - rayPos.x) * invRayDir;
        tFar = (pMax.x - rayPos.x) * invRayDir;
 
        if (tNear > tFar)
        {
            std::swap(tNear, tFar);
        }
        t0 = tNear > t0 ? tNear : t0;
        t1 = tFar < t1 ? tFar : t1;
        if (t0 > t1)
        {
            return false;
        }
 
        // Composante sur Y
        // TODO : If denominator is Zero
        invRayDir = 1 / rayDir.y;
        tNear = (pMin.y - rayPos.y) * invRayDir;
        tFar = (pMax.y - rayPos.y) * invRayDir;
 
        if (tNear > tFar)
        {
            std::swap(tNear, tFar);
        }
        t0 = tNear > t0 ? tNear : t0;
        t1 = tFar < t1 ? tFar : t1;
        if (t0 > t1)
        {
            return false;
        }
 
        // Composante sur Z
        // TODO : If denominator is Zero
        invRayDir = 1 / rayDir.z;
        tNear = (pMin.z - rayPos.z) * invRayDir;
        tFar = (pMax.z - rayPos.z) * invRayDir;
 
        if (tNear > tFar)
        {
            std::swap(tNear, tFar);
        }
        t0 = tNear > t0 ? tNear : t0;
        t1 = tFar < t1 ? tFar : t1;
        if (t0 > t1)
        {
            return false;
        }
 
        if (hitt0)
        {
            *hitt0 = t0;
        }
        if (hitt1)
        {
            *hitt1 = t1;
        }
        return true;
    }
 
    /* Not used (and seems not coded)
    void print()
    {
        if (isNull)   //std::cout << "Non initialisee." << std::endl;
        {
        }
        else   //std::cout << "Initialisee." << std::endl;
        {
        }
    } */
 
    vec3& operator[](int i)
    {
        if (i == 0)
        {
            return pMin;
        }
        else
        {
            return pMax;
        }
    }
    const vec3& operator[](int i) const
    {
        if (i == 0)
        {
            return pMin;
        }
        else
        {
            return pMax;
        }
    }
 
    decimal diag()
    {
        return sqrt((pMax.x - pMin.x) * (pMax.x - pMin.x) + (pMax.y - pMin.y) * (pMax.y - pMin.y) +
                    (pMax.z - pMin.z) * (pMax.z - pMin.z));
    }
 
    bool Inside(const vec3& pt) const
    {
        return (pt.x >= pMin.x && pt.x <= pMax.x && pt.y >= pMin.y && pt.y <= pMax.y && pt.z >= pMin.z &&
                pt.z <= pMax.z);
    }
};
 
#endif // BBOX_H