FIX: random.gumbel segfault on GPU (#560)

* FIX: random.gumbel segfault on GPU

Author: samir-nasibli

dpnp/backend/kernels/dpnp_krnl_random.cpp

@@ -298,6 +298,9 @@ void dpnp_rng_geometric_c(void* result, const float p, const size_t size)
     event_out.wait();
 }
 
+template <typename _KernelNameSpecialization>
+class dpnp_blas_scal_c_kernel;
+
 template <typename _DataType>
 void dpnp_rng_gumbel_c(void* result, const double loc, const double scale, const size_t size)
 {
@@ -308,15 +311,33 @@ void dpnp_rng_gumbel_c(void* result, const double loc, const double scale, const
     }
 
     const _DataType alpha = (_DataType(-1.0));
-    const _DataType stride = (_DataType(1.0));
+    std::int64_t incx = 1;
     _DataType* result1 = reinterpret_cast<_DataType*>(result);
     double negloc = loc * (double(-1.0));
 
     mkl_rng::gumbel<_DataType> distribution(negloc, scale);
-    // perform generation
     event = mkl_rng::generate(distribution, DPNP_RNG_ENGINE, size, result1);
     event.wait();
-    event = mkl_blas::scal(DPNP_QUEUE, size, alpha, result1, stride);
+
+    // OK for CPU and segfault for GPU device
+    // event = mkl_blas::scal(DPNP_QUEUE, size, alpha, result1, incx);
+    if (dpnp_queue_is_cpu_c())
+    {
+        event = mkl_blas::scal(DPNP_QUEUE, size, alpha, result1, incx);
+    }
+    else
+    {
+        // for (size_t i = 0; i < size; i++) result1[i] *= alpha;
+        cl::sycl::range<1> gws(size);
+        auto kernel_parallel_for_func = [=](cl::sycl::id<1> global_id) {
+            size_t i = global_id[0];
+            result1[i] *= alpha;
+        };
+        auto kernel_func = [&](cl::sycl::handler& cgh) {
+            cgh.parallel_for<class dpnp_blas_scal_c_kernel<_DataType>>(gws, kernel_parallel_for_func);
+        };
+        event = DPNP_QUEUE.submit(kernel_func);
+    }
     event.wait();
 }
 

tests/skipped_tests_gpu.tbl

@@ -1,8 +1,3 @@
-// gumbel distribution issue https://github.com/IntelPython/dpnp/issues/569
-tests/test_random.py::TestDistributionsGumbel::test_extreme_value
-tests/test_random.py::TestDistributionsGumbel::test_invalid_args
-tests/test_random.py::TestDistributionsGumbel::test_moments
-tests/test_random.py::TestDistributionsGumbel::test_seed
 //----------------------------------------------------------------------
 // eig/eigvals/svd issue https://github.com/IntelPython/dpnp/issues/567
 tests/test_linalg.py::test_eig_arange[2-float64]