static double ptLength(Pt p) {
  ret sqrt(p.x*p.x+p.y*p.y);
}