Refactors comparison functions to use math_utils implementations

Author: ndgrigorian

dpctl/tensor/libtensor/include/kernels/elementwise_functions/greater.hpp

@@ -30,6 +30,7 @@
 #include <cstdint>
 #include <type_traits>
 
+#include "utils/math_utils.hpp"
 #include "utils/offset_utils.hpp"
 #include "utils/type_dispatch.hpp"
 #include "utils/type_utils.hpp"
@@ -67,12 +68,8 @@ template <typename argT1, typename argT2, typename resT> struct GreaterFunctor
                       tu_ns::is_complex<argT2>::value)
         {
             static_assert(std::is_same_v<argT1, argT2>);
-            using realT = typename argT1::value_type;
-            realT real1 = std::real(in1);
-            realT real2 = std::real(in2);
-
-            return (real1 == real2) ? (std::imag(in1) > std::imag(in2))
-                                    : real1 > real2;
+            using dpctl::tensor::math_utils::greater_complex;
+            return greater_complex<argT1>(in1, in2);
         }
         else {
             return (in1 > in2);

dpctl/tensor/libtensor/include/kernels/elementwise_functions/greater_equal.hpp

@@ -30,6 +30,7 @@
 #include <cstdint>
 #include <type_traits>
 
+#include "utils/math_utils.hpp"
 #include "utils/offset_utils.hpp"
 #include "utils/type_dispatch.hpp"
 #include "utils/type_utils.hpp"
@@ -68,12 +69,8 @@ struct GreaterEqualFunctor
                       tu_ns::is_complex<argT2>::value)
         {
             static_assert(std::is_same_v<argT1, argT2>);
-            using realT = typename argT1::value_type;
-            realT real1 = std::real(in1);
-            realT real2 = std::real(in2);
-
-            return (real1 == real2) ? (std::imag(in1) >= std::imag(in2))
-                                    : real1 >= real2;
+            using dpctl::tensor::math_utils::greater_equal_complex;
+            return greater_equal_complex<argT1>(in1, in2);
         }
         else {
             return (in1 >= in2);

dpctl/tensor/libtensor/include/kernels/elementwise_functions/less.hpp

@@ -29,6 +29,7 @@
 #include <cstdint>
 #include <type_traits>
 
+#include "utils/math_utils.hpp"
 #include "utils/offset_utils.hpp"
 #include "utils/type_dispatch.hpp"
 #include "utils/type_utils.hpp"
@@ -66,12 +67,8 @@ template <typename argT1, typename argT2, typename resT> struct LessFunctor
                       tu_ns::is_complex<argT2>::value)
         {
             static_assert(std::is_same_v<argT1, argT2>);
-            using realT = typename argT1::value_type;
-            realT real1 = std::real(in1);
-            realT real2 = std::real(in2);
-
-            return (real1 == real2) ? (std::imag(in1) < std::imag(in2))
-                                    : real1 < real2;
+            using dpctl::tensor::math_utils::less_complex;
+            return less_complex<argT1>(in1, in2);
         }
         else {
             return (in1 < in2);

dpctl/tensor/libtensor/include/kernels/elementwise_functions/less_equal.hpp

@@ -30,6 +30,7 @@
 #include <cstdint>
 #include <type_traits>
 
+#include "utils/math_utils.hpp"
 #include "utils/offset_utils.hpp"
 #include "utils/type_dispatch.hpp"
 #include "utils/type_utils.hpp"
@@ -67,12 +68,8 @@ template <typename argT1, typename argT2, typename resT> struct LessEqualFunctor
                       tu_ns::is_complex<argT2>::value)
         {
             static_assert(std::is_same_v<argT1, argT2>);
-            using realT = typename argT1::value_type;
-            realT real1 = std::real(in1);
-            realT real2 = std::real(in2);
-
-            return (real1 == real2) ? (std::imag(in1) <= std::imag(in2))
-                                    : real1 <= real2;
+            using dpctl::tensor::math_utils::less_equal_complex;
+            return less_equal_complex<argT1>(in1, in2);
         }
         else {
             return (in1 <= in2);