intersect2d
pub fn ulps_eq<T>(a: &T, b: &T) -> bool where T: Float + CoordNum + PartialOrd + AbsDiffEq + UlpsEq, T::Epsilon: Copy,