let kernel_to_distance k =
  fun x y ->
    k x x -. 2.0 *. k x y +. k y y