line.h

/*
 * Copyright (c) 2017 FilipeCN
 *
 * The MIT License (MIT)
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to deal
 * in the Software without restriction, including without limitation the rights
 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 * copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 * THE SOFTWARE.
 *
 */
 
#ifndef HERMES_GEOMETRY_LINE_H
#define HERMES_GEOMETRY_LINE_H
 
#include <hermes/geometry/point.h>
#include <hermes/geometry/vector.h>
 
namespace hermes {
 
// *********************************************************************************************************************
//                                                                                                               Line
// *********************************************************************************************************************
class Line2 {
public:
  // *******************************************************************************************************************
  //                                                                                                     CONSTRUCTORS
  // *******************************************************************************************************************
  Line2() = default;
  Line2(point2 _a, vec2 _d) {
    a = _a;
    d = normalize(_d);
  }
  // *******************************************************************************************************************
  //                                                                                                        OPERATORS
  // *******************************************************************************************************************
  point2 operator()(float t) const { return a + d * t; }
  // *******************************************************************************************************************
  //                                                                                                          METHODS
  // *******************************************************************************************************************
  //                                                                                                           access
  [[nodiscard]] vec2 direction() const { return normalize(d); }
  //                                                                                                          queries
  [[nodiscard]] float projection(point2 p) const { return dot((p - a), d); }
  [[nodiscard]] point2 closestPoint(point2 p) const { return (*this)(projection(p)); }
  // *******************************************************************************************************************
  //                                                                                                            DEBUG
  // *******************************************************************************************************************
  friend std::ostream &operator<<(std::ostream &os, const Line2 &p) {
    os << "[Line]\n";
    os << p.a << p.d;
    return os;
  }
  // *******************************************************************************************************************
  //                                                                                                    PUBLIC FIELDS
  // *******************************************************************************************************************
  point2 a; 
  vec2 d;   
};
 
// *********************************************************************************************************************
//                                                                                                               Line
// *********************************************************************************************************************
class Line {
public:
  // *******************************************************************************************************************
  //                                                                                                     CONSTRUCTORS
  // *******************************************************************************************************************
  Line() = default;
  Line(point3 _a, vec3 _d) {
    a = _a;
    d = normalize(_d);
  }
  // *******************************************************************************************************************
  //                                                                                                        OPERATORS
  // *******************************************************************************************************************
  point3 operator()(float t) const { return a + d * t; }
  // *******************************************************************************************************************
  //                                                                                                          METHODS
  // *******************************************************************************************************************
  //                                                                                                           access
  [[nodiscard]] vec3 direction() const { return normalize(d); }
  //                                                                                                          queries
  [[nodiscard]] float projection(point3 p) const { return dot((p - a), d); }
  [[nodiscard]] point3 closestPoint(point3 p) const { return (*this)(projection(p)); }
  // *******************************************************************************************************************
  //                                                                                                            DEBUG
  // *******************************************************************************************************************
  friend std::ostream &operator<<(std::ostream &os, const Line &p) {
    os << "[Line]\n";
    os << p.a << p.d;
    return os;
  }
  // *******************************************************************************************************************
  //                                                                                                    PUBLIC FIELDS
  // *******************************************************************************************************************
  point3 a; 
  vec3 d;   
};
 
} // hermes namespace
 
#endif