Move all logic with cyl_bessel_i0 support in i0.hpp

Author: antonwolfy

dpnp/backend/extensions/window/kaiser.cpp

@@ -32,22 +32,6 @@
 
 #include <sycl/sycl.hpp>
 
-/**
- * Version of SYCL DPC++ 2025.1 compiler where an issue with
- * sycl::ext::intel::math::cyl_bessel_i0(x) is fully resolved.
- */
-#ifndef __SYCL_COMPILER_BESSEL_I0_SUPPORT
-#define __SYCL_COMPILER_BESSEL_I0_SUPPORT 20241208L
-#endif
-
-// The <sycl/ext/intel/math.hpp> header and related types like
-// _iml_half_internal are not compatible with NVIDIA GPUs and cause compilation
-// errors when building with -fsycl-targets=nvptx64-nvidia-cuda.
-#if !defined(__NVPTX__) &&                                                     \
-    (__SYCL_COMPILER_VERSION >= __SYCL_COMPILER_BESSEL_I0_SUPPORT)
-#include <sycl/ext/intel/math.hpp>
-#endif
-
 #include "../kernels/elementwise_functions/i0.hpp"
 
 namespace dpnp::extensions::window
@@ -78,12 +62,7 @@ class KaiserFunctor
 
     void operator()(sycl::id<1> id) const
     {
-#if !defined(__NVPTX__) &&                                                     \
-    (__SYCL_COMPILER_VERSION >= __SYCL_COMPILER_BESSEL_I0_SUPPORT)
-        using sycl::ext::intel::math::cyl_bessel_i0;
-#else
-        using dpnp::kernels::i0::impl::cyl_bessel_i0;
-#endif
+        using dpnp::kernels::i0::cyl_bessel_i0;
 
         const auto i = id.get(0);
         const T alpha = (N - 1) / T(2);

dpnp/backend/kernels/elementwise_functions/i0.hpp

@@ -35,16 +35,28 @@
 #define __SYCL_COMPILER_BESSEL_I0_SUPPORT 20241208L
 #endif
 
-// The <sycl/ext/intel/math.hpp> header and related types like
-// _iml_half_internal are not compatible with NVIDIA GPUs and cause compilation
-// errors when building with -fsycl-targets=nvptx64-nvidia-cuda.
-#if !defined(__NVPTX__) &&                                                     \
+/**
+ * The <sycl/ext/intel/math.hpp> header and related types like
+ * _iml_half_internal are not compatible with NVIDIA GPUs and cause compilation
+ * errors when building with -fsycl-targets=nvptx64-nvidia-cuda.
+ */
+#if !defined(__NVPTX__)
+#define __SYCL_EXT_INTEL_MATH_SUPPORT
+#endif
+
+#if defined(__SYCL_EXT_INTEL_MATH_SUPPORT) &&                                  \
     (__SYCL_COMPILER_VERSION >= __SYCL_COMPILER_BESSEL_I0_SUPPORT)
 #include <sycl/ext/intel/math.hpp>
 #endif
 
 namespace dpnp::kernels::i0
 {
+#if defined(__SYCL_EXT_INTEL_MATH_SUPPORT) &&                                  \
+    (__SYCL_COMPILER_VERSION >= __SYCL_COMPILER_BESSEL_I0_SUPPORT)
+using sycl::ext::intel::math::cyl_bessel_i0;
+
+#else
+
 /**
  * The below implementation of Bessel function of order 0
  * is based on the source code from https://github.com/gcc-mirror/gcc
@@ -243,6 +255,10 @@ inline Tp cyl_bessel_i0(Tp x)
 }
 } // namespace impl
 
+using impl::cyl_bessel_i0;
+
+#endif
+
 template <typename argT, typename resT>
 struct I0Functor
 {
@@ -257,13 +273,6 @@ struct I0Functor
 
     resT operator()(const argT &x) const
     {
-#if !defined(__NVPTX__) &&                                                     \
-    (__SYCL_COMPILER_VERSION >= __SYCL_COMPILER_BESSEL_I0_SUPPORT)
-        using sycl::ext::intel::math::cyl_bessel_i0;
-#else
-        using impl::cyl_bessel_i0;
-#endif
-
         if constexpr (std::is_same_v<resT, sycl::half>) {
             return static_cast<resT>(cyl_bessel_i0<float>(float(x)));
         }