libcruft-util/geom/aabb.hpp

/*
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
 * Copyright 2015-2017 Danny Robson <danny@nerdcruft.net>
 */


#ifndef CRUFT_UTIL_GEOM_AABB_HPP
#define CRUFT_UTIL_GEOM_AABB_HPP

#include "../point.hpp"
#include "../extent.hpp"

#include <cstdint>

namespace util::geom {
    ///////////////////////////////////////////////////////////////////////////
    /// represents an axis-aligned bounding-box through two opposing corners.
    ///
    /// p0 must be less-than-or-equal to p1. equality is allowed so that we
    /// can represent zero sized bounding-boxes.
    template <size_t S, typename T>
    struct aabb {
        aabb () = default;
        aabb (point<S,T>, point<S,T>);

        extent<S,T> magnitude (void) const;
        T diameter (void) const;

        /// tests whether a point lies within the region, inclusive of borders
        constexpr bool
        inclusive (point<S,T> p) const noexcept
        { return all (p0 <= p && p1 >= p); }

        point<S,T> closest (point<S,T>) const;

        void cover (point<S,T>);

        aabb<S,T> operator+ (vector<S,T>) const;
        aabb<S,T> operator- (vector<S,T>) const;

        ::util::point<S,T> p0;
        ::util::point<S,T> p1;
    };


    ///////////////////////////////////////////////////////////////////////////
    template <std::size_t S, typename T>
    constexpr bool
    operator== (const aabb<S,T> &a, const aabb<S,T> &b) noexcept
    {
        return a.p0 == b.p0 && a.p1 == b.p1;
    }


    ///////////////////////////////////////////////////////////////////////////
    typedef aabb<2,float> aabb2f;
    typedef aabb<2,unsigned> aabb2u;
    typedef aabb<2,int> aabb2i;

    typedef aabb<3,float> aabb3f;
    typedef aabb<3,unsigned> aabb3u;
    typedef aabb<3,int> aabb3i;
}


///////////////////////////////////////////////////////////////////////////////
#include "./sample.hpp"

#include <random>

namespace util::geom {
    template <size_t S, typename T, typename G>
    struct sampler<S,T,aabb,G> {
        static point<S,T>
        fn (aabb<S,T> b, G &g)
        {
            std::uniform_real_distribution<T> d;

            point<S,T> p;
            std::generate (p.begin (), p.end (), [&] (void) { return d (g); });

            return p * (b.p1 - b.p0) + b.p0.template as<util::vector> ();
        }
    };
}

#endif