add CUDSS test

Author: ChrisRackauckas

test/gpu/Project.toml

@@ -1,6 +1,7 @@
 [deps]
 BlockDiagonals = "0a1fb500-61f7-11e9-3c65-f5ef3456f9f0"
 CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
+CUDSS = "45b445bb-4962-46a0-9369-b4df9d0f772e"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"

test/gpu/cuda.jl

@@ -1,4 +1,5 @@
 using LinearSolve, CUDA, LinearAlgebra, SparseArrays, StableRNGs
+using CUDA.CUSPARSE, CUDSS
 using Test
 
 CUDA.allowscalar(false)
@@ -91,3 +92,17 @@ prob2 = LinearProblem(transpose(A), b)
     sol = solve(prob2, alg; alias = LinearAliasSpecifier(alias_A = false))
     @test norm(transpose(A) * sol.u .- b) < 1e-5
 end
+
+@testset "CUDSS" begin
+    T = Float32
+    n = 100
+    A_cpu = sprand(T, n, n, 0.05) + I
+    x_cpu = zeros(T, n)
+    b_cpu = rand(T, n)
+
+    A_gpu_csr = CuSparseMatrixCSR(A_cpu)
+    b_gpu = CuVector(b_cpu)
+
+    prob = LinearProblem(A_gpu_csr, b_gpu)
+    sol = solve(prob)
+end