point2d.h

#pragma once
 
#include <Eigen/Dense>
#include "cevalm.hpp"

class Point2d {
   private:
    int x_, y_;
 
   public:
    constexpr Point2d() : x_(0), y_(0) {}

    constexpr Point2d(int x, int y) : x_(x), y_(y) {}

    constexpr Point2d(Eigen::Vector2i vector)
     : x_(vector(0)), y_(vector(1)) {}

    constexpr int x() const { return x_; }

    constexpr void set_x(int val) { x_ = val; }

    constexpr int y() const { return y_; }

    constexpr void set_y(int val) { y_ = val; }

    constexpr Eigen::Vector2i as_vector() const { return Eigen::Vector2i(x_, y_); }

    constexpr Eigen::Vector2d as_vector(Eigen::Vector2d X, Eigen::Vector2d Y) const {
        return Eigen::Vector2d(
                x_ * X(0) + y_ * Y(0),
                x_ * X(1) + y_ * Y(1)
        );
    }

    constexpr int manhattan_distance(Point2d other) const {
        return cevalm::abs(x_ - other.x_) + cevalm::abs(y_ - other.y_);
    }

    constexpr int manhattan_norm() const { return cevalm::abs(x_) + cevalm::abs(y_); }

    constexpr double distance(Point2d other) const { return cevalm::hypot(x_ - other.x_, y_ - other.y_);  }

    constexpr double norm() const { return cevalm::hypot(x_, y_); }

    constexpr Point2d inverse() const { return Point2d{-x_, -y_};; }

    constexpr int dot(Point2d other) const { return (x_ * other.x_) + (y_ * other.y_); }

    constexpr Point2d operator+(Point2d other) const { return Point2d{x_ + other.x_, y_ + other.y_}; }

    constexpr Point2d& operator+=(Point2d other) { return *this = *this + other; }

    constexpr Point2d operator-(Point2d other) const { return Point2d{x_ - other.x_, y_ - other.y_}; }

    constexpr Point2d& operator-=(Point2d other) { return *this = *this - other; }

    constexpr Point2d operator-() const { return inverse(); }

    constexpr Point2d operator*(int scalar) const { return Point2d{x_ * scalar, y_ * scalar}; }

    friend constexpr Point2d operator*(int scalar, Point2d point) { return point * scalar; }

    constexpr int operator*(Point2d other) const { return dot(other); }

    constexpr Point2d& operator*=(int scalar) { return *this = *this * scalar; }

    constexpr bool operator==(Point2d other) const { return x_ == other.x_ && y_ == other.y_; }
};