Deal with empty sparse matrices

Author: michel2323

lib/mkl/wrappers_sparse.jl

@@ -23,7 +23,10 @@ for (fname, elty, intty) in ((:onemklSsparse_set_csr_data   , :Float32   , :Int3
             m, n = dims
             nnzA = length(nzVal)
             queue = global_queue(context(nzVal), device(nzVal))
-            $fname(sycl_queue(queue), handle_ptr[], m, n, 'O', rowPtr, colVal, nzVal)
+            # Don't update handle if matrix is empty
+            if m != 0 && n != 0
+                $fname(sycl_queue(queue), handle_ptr[], m, n, 'O', rowPtr, colVal, nzVal)
+            end
             dA = oneSparseMatrixCSR{$elty,$intty}(handle_ptr[], rowPtr, colVal, nzVal, (m,n), nnzA)
             finalizer(sparse_release_matrix_handle, dA)
             return dA
@@ -38,8 +41,10 @@ for (fname, elty, intty) in ((:onemklSsparse_set_csr_data   , :Float32   , :Int3
             onemklXsparse_init_matrix_handle(handle_ptr)
             m, n = dims
             nnzA = length(nzVal)
-            @show dims
-            $fname(sycl_queue(queue), handle_ptr[], n, m, 'O', colPtr, rowVal, nzVal)  # CSC of A is CSR of Aᵀ
+            # Don't update handle if matrix is empty
+            if m != 0 && n != 0
+                $fname(sycl_queue(queue), handle_ptr[], n, m, 'O', colPtr, rowVal, nzVal)  # CSC of A is CSR of Aᵀ
+            end
             dA = oneSparseMatrixCSC{$elty,$intty}(handle_ptr[], colPtr, rowVal, nzVal, dims, nnzA)
             finalizer(sparse_release_matrix_handle, dA)
             return dA