src/integer-geometry/dot.hpp

View this file on GitHub
Last update: 2022-05-14 03:05:02+09:00
Include: #include "src/integer-geometry/dot.hpp"

Depends on

Required by

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

Code

#pragma once

#include "point.hpp"

namespace intgeometry2d {

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

} // intgeometry2d

#line 2 "src/integer-geometry/dot.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/dot.hpp"

namespace intgeometry2d {

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

} // intgeometry2d