src/integer-geometry/unit-test/unit-test-det.cpp

View this file on GitHub
Last update: 2022-05-14 03:05:02+09:00

Depends on

Code

#include <cassert>

#include "../det.hpp"

int main() {
  using intgeometry2d::i64;
  using intgeometry2d::lattice_point;
  using intgeometry2d::det;

  using LP64 = lattice_point<i64>;
  assert(det(LP64(12345, 67890), LP64(3141592, 653589)) == -205'214'124'675ll);
}

#line 1 "src/integer-geometry/unit-test/unit-test-det.cpp"
#include <cassert>

#line 2 "src/integer-geometry/det.hpp"

#line 2 "src/integer-geometry/point.hpp"

#line 2 "src/integer-geometry/alias.hpp"

#include <cstddef>
#include <cstdint>

namespace intgeometry2d {
  using isize = std::ptrdiff_t;
  using usize = std::size_t;

  using i32 = std::int_fast32_t;
  using i64 = std::int_fast64_t;
  using u32 = std::uint_fast32_t;
  using u64 = std::uint_fast64_t;
} // intgeometry2d
#line 4 "src/integer-geometry/point.hpp"

namespace intgeometry2d {

  template< typename Z >
  class lattice_point {
    Z x_, y_;

  public:
    lattice_point() {}
    lattice_point(Z x_, Z y_) : x_(x_), y_(y_) {}
    
    Z x() const { return x_; }
    Z y() const { return y_; }

    bool operator==(const lattice_point &p) const { return x_ == p.x_ and y_ == p.y_; }
    bool operator!=(const lattice_point &p) const { return x_ != p.x_ or  y_ != p.y_; }

    lattice_point operator+(lattice_point p) { return lattice_point(x_ + p.x_, y_ + p.y_); }
    lattice_point operator-(lattice_point p) { return lattice_point(x_ - p.x_, y_ - p.y_); }

    Z norm() const { return x_ * x_ + y_ * y_; }
  };

} // intgeometry2d
#line 4 "src/integer-geometry/det.hpp"

namespace intgeometry2d {

  template< typename Z >
  Z det(const lattice_point<Z> &a, const lattice_point<Z> &b) {
    return a.x() * b.y() - a.y() * b.x();
  }

} // intgeometry2d
#line 4 "src/integer-geometry/unit-test/unit-test-det.cpp"

int main() {
  using intgeometry2d::i64;
  using intgeometry2d::lattice_point;
  using intgeometry2d::det;

  using LP64 = lattice_point<i64>;
  assert(det(LP64(12345, 67890), LP64(3141592, 653589)) == -205'214'124'675ll);
}