Function circumcircle

code/geometry/circles.cpp:39–42

  assert(tangent(ip, A, rA, P.first, Q.first) == 2);
  assert(tangent(ip, B, rB, P.second, Q.second) == 2); }
pair<point,double> circumcircle(point a, point b, point c) {
  b -= a, c -= a;
  point p = perp(b*norm(c)-c*norm(b))/2.0/cross(b, c);
  return make_pair(a+p,abs(p)); }
// vim: cc=60 ts=2 sts=2 sw=2: