// angle between p1 and p2 with respect to center point
// result is between -pi() and pi()
static double threePointAngle(IDoublePt p1, IDoublePt center, IDoublePt p2) {
  ret atan2(
      p1.x_double()-center.x_double(),
      p1.y_double()-center.y_double())
    -atan2(
      p2.x_double()-center.x_double(),
      p2.y_double()-center.y_double());
}