let dist_from_point o v p =
  let v' = sub p o in
  let s = scal v' v in
    sqrt ((sqnorm v') -. s *. s)