Use approximate comparisons for floats

Author: GabrielMajeri

examples/stencil/src/lib.rs

@@ -133,22 +133,24 @@ fn assert_data_eq(a: &Data, b: &Data) {
             for x in 0..a.n.0 {
                 let idx = (x + y * a.n.1 + z * a.n.1 * a.n.0) as usize;
 
-                assert_eq!(
-                    a.vsq[idx], b.vsq[idx],
+                const EPSILON: f32 = 1E-4;
+
+                assert!(
+                    (a.vsq[idx] - b.vsq[idx]).abs() < EPSILON,
                     "vsq diff at idx = {} ({}, {}, {})",
-                    idx, x, y, z
+                    idx, x, y, z,
                 );
 
-                assert_eq!(
-                    a.a.0[idx], b.a.0[idx],
+                assert!(
+                    (a.a.0[idx] - b.a.0[idx]).abs() < EPSILON,
                     "a.0 diff at idx = {} ({}, {}, {})",
-                    idx, x, y, z
+                    idx, x, y, z,
                 );
 
-                assert_eq!(
-                    a.a.1[idx], b.a.1[idx],
+                assert!(
+                    (a.a.1[idx] - b.a.1[idx]).abs() < EPSILON,
                     "a.1 diff at idx = {} ({}, {}, {})",
-                    idx, x, y, z
+                    idx, x, y, z,
                 );
             }
         }

examples/stencil/src/scalar.rs

@@ -63,7 +63,8 @@ pub fn scalar(
 #[cfg(all(test, feature = "ispc"))]
 mod tests {
     use super::scalar;
-    use {ispc_loops::serial, Data};
+    use ispc_loops::serial;
+    use {assert_data_eq, Data};
 
     #[test]
 
@@ -74,6 +75,6 @@ mod tests {
         let mut data_ispc = Data::default();
         data_ispc.exec(serial);
 
-        assert_eq!(data_scalar, data_ispc);
+        assert_data_eq(&data_scalar, &data_ispc);
     }
 }