plane.hpp

Go to the documentation of this file.

 
/* COPYRIGHT (c) 2024 John Reppy (https://cs.uchicago.edu/~jhr)
 * All rights reserved.
 */
 
#ifndef _CS237_PLANE_HPP_
#define _CS237_PLANE_HPP_
 
#ifndef _CS237_HPP_
#error "cs237/plane.hpp should not be included directly"
#endif
 
namespace cs237 {
 
namespace __detail {
 
template <typename T>
struct Plane {
    using vec3 = glm::vec<3, T, glm::defaultp>;
    using vec4 = glm::vec<4, T, glm::defaultp>;
 
    vec4 _nd;     
 
    Plane () { }
 
    Plane (vec3 n, T d) : _nd(vec4(n, T(d))) { }
 
    Plane (vec3 n, vec3 p)
    {
        vec3 norm = glm::normalize(n);
        double d = -glm::dot(norm, p);
        this->_nd = vec4(norm, d);
    }
 
    vec3 norm () const { return vec3(this->_nd); }
 
    T dist () const { return this->_nd.w; }
 
    T distanceToPt (vec3 const &p) const
    {
    return glm::dot(this->_nd, vec4(p, 1.0));
    }
 
    vec3 project (vec3 const &p) const
    {
        return p - this->distanceToPt(p) * this->norm();
    }
 
    std::string toString () const
    {
        return "Plane(" + glm::to_string(this->norm()) + ", "
            + to_string(this->dist()) + ")";
    }
 
};
 
/***** Output *****/
 
template <typename T>
inline std::ostream& operator<< (std::ostream& s, Plane<T> const &plane)
{
    return (s << plane.toString());
}
 
} // namespace __detail
 
using Planef_t = __detail::Plane<float>;
using Planed_t = __detail::Plane<double>;
 
} // namespace cs237
 
#endif // !_CS237_PLANE_HPP_