#pragma once

#include <math.h>

struct vec {
	double x, y;

	vec(double xx, double yy) : x(xx), y(yy) {}

	double norm() const {
		return sqrt(x*x + y*y);
	}

	double angle() const {
		double xx = x / norm();
		double a = acos(x);
		return (y > 0 ? a : -a);
	}

	vec normalize() const {
		double n = norm();
		return vec(x / n, y / n);
	}

	bool is_nil() const {
		return x == 0 && y == 0;
	}

	vec operator+ (const vec& o) const {
		return vec(x + o.x, y + o.y);
	}
	vec operator- (const vec& o) const {
		return vec(x - o.x, y - o.y);
	}
	vec operator* (double m) const {
		return vec(m * x, m * y);
	}
	vec operator/ (double m) const {
		return vec(m / x, m / y);
	}

	static vec from_polar(double r, double theta) {
		return vec(r * cos(theta), r * sin(theta));
	}

	static vec dot(vec a, vec b) {		// dot product (produit scalaire)
		return a.x * b.x + a.y * b.y;
	}
	static vec cross(vec a, vec b) {	// cross product (déterminant 2x2)
		return a.x * b.y - a.y * b.x;
	}
	static vec angle(vec a, vec b) {	// oriented angle between two vectors
		if (a.is_nil() || b.is_nil()) return 0;
		float cos = dot(a.normalize(), b.normalize());
		if (cos <= -1) return M_PI;
		float uangle = acos(cos);
		if (cross(a, b) > 0) {
			return uangle;
		} else {
			return -uangle;
		}
	}
};

struct line {
	// Line defined by ax + by + c = 0
	double a, b, c;

	line(double aa, double bb, double cc) : a(aa), b(bb), c(cc) {}
	line(vec p1, vec p2) {
		a = p1.x-p2.x ;
		b = p1.y-p2.y ;
		c = p1.x*(p2.x-p1.x)+p1.y*(p2.y-p1.y);
	}

	bool on_line(vec p) const {
		return a * p.x + b * p.y + c == 0;
	}

	double dist(vec p) const {
		// calculate distance from p to the line
	    return abs(a*p.x+b*p.y+c)/sqrt(a*a+b*b);
	}

	vec proj(vec p) const {
		// calculate orthogonal projection of point p on the line
		// TODO
		return vec(0, 0);
	}

	double angle() const {
		return vec(-b, a).angle();
	}

};

struct segment {
	vec a, b;

	segment(vec pa, vec pb) : a(pa), b(pb) {}

	bool on_segment(vec p) const {
		// TODO
		// does point intersect segment?
		return false;
	}

	double dist(vec p) const {
		
		return 1;
	}
};

struct circle {
	vec c;
	double r;

	circle(double x, double y, double rr) : c(x, y), r(rr) {}
	circle(vec cc, double rr) : c(cc), r(rr) {}

	bool on_circle(vec.p) const {
		return ((p - c).norm() == r);
	}

	double dist(vec p) const {
		// TODO
		return 1;
	}
};

struct circpoint {
	circle c;
	double theta;

	circpoint(circle cc, double th) : c(cc), theta(th) {}
	circpoint(vec cc, double rr, double th : c(cc, rr), theta(th) {}

	vec pos() const {
		return c.c + vec(c.r * cos(theta), c.r * sin(theta));
	}
};

struct circarc {
	circle c;
	double theta1, theta2;

	circarc(circle cc, double tha, double thb) : c(cc), theta1(tha), theta2(thb) {}

	double dist(vec p) const {
		// TODO
		return 1;
	}
};

/* vim: set ts=4 sw=4 tw=0 noet :*/