let div_check x y =
  classify_float (x /. y) = FP_infinite