Math.hpp

Go to the documentation of this file.

#pragma once
 
#include <cmath>
 
#include "utils/Types.hpp"
#include "utils/Macros.hpp"
 
#define APPROX_THRESHOLD 0.01

namespace ppforest2::math {
  template<typename A, typename B, typename T> inline bool is_approx(A a, B b, T threshold) {
    return fabs(a - b) < threshold;
  }

  template<typename A, typename B> inline bool is_approx(A a, B b) {
    return is_approx(a, b, APPROX_THRESHOLD);
  }

  template<typename A, typename B> inline bool is_module_approx(A a, B b) {
    return is_approx(fabs(a), fabs(b));
  }

  template<typename T> bool collinear(types::Vector<T> const& a, types::Vector<T> const& b) {
    return is_module_approx(a.dot(b) / (a.norm() * b.norm()), 1.0);
  }
}