Update cusolverrf.jl

Author: ChrisRackauckas

test/gpu/cusolverrf.jl

@@ -54,20 +54,6 @@ using Test
         @test norm(A2 * sol2.u - b2) / norm(b2) < 1e-10
     end
 
-    # Test multiple right-hand sides
-    @testset "Multiple RHS" begin
-        nrhs = 5
-        B = rand(n, nrhs)
-        
-        prob_multi = LinearProblem(A, B)
-        sol_multi = solve(prob_multi, CUSOLVERRFFactorization())
-        
-        # Check each solution
-        for i in 1:nrhs
-            @test norm(A * sol_multi.u[:, i] - B[:, i]) / norm(B[:, i]) < 1e-10
-        end
-    end
-    
     # Test adjoint solve
     @testset "Adjoint Solve" begin
         prob_adj = LinearProblem(A', b)
@@ -85,4 +71,4 @@ using Test
         # This should error since CUSOLVERRF only supports Float64
         @test_throws Exception solve(prob_f32, CUSOLVERRFFactorization())
     end
-end
+end