+#define SCAL_EPSILON (S(0.000001))
+
+/*
+ * Check if two scalars are mostly equal, given a margin of error.
+ */
+INLINE_MAYBE
+bool scal_isequal2(scal_t a, scal_t b, scal_t e)
+{
+ return scal_abs(a -b) < e;
+}
+
+/*
+ * Check if two scalars are mostly equal.
+ */
+INLINE_MAYBE
+bool scal_isequal(scal_t a, scal_t b)
+{
+ return scal_isequal2(a, b, SCAL_EPSILON);
+}
+