1 files changed, 46 insertions, 0 deletions
diff --git a/compiler/nike/src/RelativeEpsilonEqual.cpp b/compiler/nike/src/RelativeEpsilonEqual.cpp
new file mode 100644
index 000000000..ac8d2ad82
--- /dev/null
+++ b/compiler/nike/src/RelativeEpsilonEqual.cpp
@@ -0,0 +1,46 @@
+/*
+ * Copyright (c) 2018 Samsung Electronics Co., Ltd. All Rights Reserved
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *    http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+#include "nike/RelativeEpsilonEqual.h"
+
+#include <cmath>
+#include <cfloat>
+
+#include <algorithm>
+
+namespace nike
+{
+
+bool RelativeEpsilonEqualFunctor::operator()(float lhs, float rhs) const
+{
+  if (std::isnan(lhs) && std::isnan(rhs))
+  {
+    return true;
+  }
+
+  // TODO How to handle sign difference?
+  auto const delta = std::fabs(lhs - rhs);
+  auto const max = std::max(std::fabs(lhs), std::fabs(rhs));
+
+  return delta <= (max * FLT_EPSILON * _tolerance);
+}
+
+RelativeEpsilonEqualFunctor relative_epsilon_equal(unsigned tolerance)
+{
+  return RelativeEpsilonEqualFunctor{tolerance};
+}
+
+} // namespace nike