// angle between p2-p1 and the x axis (?)
static double twoPointAngle(IDoublePt p1, IDoublePt p2) {
  ret atan2(
    p2.x_double()-p1.x_double(),
    p2.y_double()-p1.y_double());
}