pub fn plane_side_3d_epsilon<S>( a: Point3<S>, b: Point3<S>, c: Point3<S>, d: Point3<S>, epsilon: S, ) -> Swhere S: OrderedField + AbsDiffEq<Epsilon = S>,