Update Circle

Author: CarlosDaniel111

Geometry/Circle.cpp

@@ -8,6 +8,11 @@ Point circumCenter(Point a, Point b, Point c) {
   return a + (b * sq(c) - c * sq(b)).perp() / b.cross(c) / 2;
 }
 
+// Retorna el punto que se encuentra en el circulo dado el angulo
+Point circlePoint(Point c, double r, double ang) {
+  return Point{c.x + cos(ang) * r, c.y + sin(ang) * r};
+}
+
 // Retorna el numero de intersecciones de la linea l con el circulo (o,r)
 // y los pone en out. Si solo hay una interseccion el par de out es igual
 int circleLine(Point o, double r, Line l, pair<Point, Point> &out) {
@@ -24,7 +29,7 @@ int circleLine(Point o, double r, Line l, pair<Point, Point> &out) {
 // la interseccion con una linea
 int circleCircle(Point o1, double r1, Point o2, double r2, pair<Point, Point> &out) {
   Point d = o2 - o1;
-  double d2 = sq(d);
+  double d2 = d.sq();
   if (d2 == 0) {
     assert(r1 != r2);  // los circulos son iguales
     return 0;

Geometry/Point.cpp

@@ -37,12 +37,13 @@ struct Point {
   Point translate(Point v) { return *this + v; }
   Point scale(Point c, double factor) { return c + (*this - c) * factor; }
   Point rotate(double ang) { return {x * cos(ang) - y * sin(ang), x * sin(ang) + y * cos(ang)}; }
-  Point rot_around(double ang, Point c) { return p + (*this - p).rotate(ang); }
+  Point rot_around(double ang, Point c) { return c + (*this - c).rotate(ang); }
   Point perp() { return {-y, x}; }
 
   T dot(Point p) { return x * p.x + y * p.y; }
   T cross(Point p) const { return x * p.y - y * p.x; }
   T cross(Point a, Point b) const { return (a - *this).cross(b - *this); }
+  double angle() const { return atan2(y, x); }
 
   friend ostream& operator<<(ostream& os, Point p) {
     return os << "(" << p.x << "," << p.y << ")";