Reenable truncated CUDA tests

Author: kshyatt

test/cuda/factorizations.jl

@@ -229,17 +229,17 @@ for V in spacelist
                 @test isisometric(N)
                 @test norm(N' * t) ≈ 0 atol = 100 * eps(norm(t))
 
-                #N = @constinferred left_null(t; trunc = (; atol = 100 * eps(norm(t))))
-                #@test isisometric(N)
-                #@test norm(N' * t) ≈ 0 atol = 100 * eps(norm(t))
+                N = @constinferred left_null(t; trunc = (; atol = 100 * eps(norm(t))))
+                @test isisometric(N)
+                @test norm(N' * t) ≈ 0 atol = 100 * eps(norm(t))
 
                 Nᴴ = @constinferred right_null(t; alg = :svd)
                 @test isisometric(Nᴴ; side = :right)
                 @test norm(t * Nᴴ') ≈ 0 atol = 100 * eps(norm(t))
 
-                #Nᴴ = @constinferred right_null(t; trunc = (; atol = 100 * eps(norm(t))))
-                #@test isisometric(Nᴴ; side = :right)
-                #@test norm(t * Nᴴ') ≈ 0 atol = 100 * eps(norm(t))
+                Nᴴ = @constinferred right_null(t; trunc = (; atol = 100 * eps(norm(t))))
+                @test isisometric(Nᴴ; side = :right)
+                @test norm(t * Nᴴ') ≈ 0 atol = 100 * eps(norm(t))
             end
 
             # empty tensor
@@ -258,15 +258,15 @@ for V in spacelist
             end
         end
 
-        #=@testset "truncated SVD" begin
+        @testset "truncated SVD" begin
             for T in eltypes,
                     t in (
                         CUDA.randn(T, W, W),
-                        #CUDA.randn(T, W, W)',
+                        CUDA.randn(T, W, W)',
                         CUDA.randn(T, W, V4),
                         CUDA.randn(T, V4, W),
-                        #CUDA.randn(T, W, V4)',
-                        #CUDA.randn(T, V4, W)',
+                        CUDA.randn(T, W, V4)',
+                        CUDA.randn(T, V4, W)',
                         DiagonalTensorMap(CUDA.randn(T, reduceddim(V1)), V1),
                     )
 
@@ -327,7 +327,7 @@ for V in spacelist
                 @test minimum(diagview(S5)) >= λ
                 @test dim(domain(S5)) ≤ nvals
             end
-        end=# # TODO
+        end
 
         @testset "Eigenvalue decomposition" begin
             for T in eltypes,
@@ -349,10 +349,10 @@ for V in spacelist
                 @test @constinferred isposdef(vdv)
                 t isa DiagonalTensorMap || @test !isposdef(t) # unlikely for non-hermitian map
 
-                #=nvals = round(Int, dim(domain(t)) / 2)
+                nvals = round(Int, dim(domain(t)) / 2)
                 d, v = @constinferred eig_trunc(t; trunc = truncrank(nvals))
                 @test t * v ≈ v * d
-                @test dim(domain(d)) ≤ nvals=#
+                @test dim(domain(d)) ≤ nvals
 
                 t2 = @constinferred project_hermitian(t)
                 D, V = eigen(t2)
@@ -380,10 +380,9 @@ for V in spacelist
                 @test isposdef(t2 - λ * one(t) + 0.1 * one(t2))
                 @test !isposdef(t2 - λ * one(t) - 0.1 * one(t2))
 
-                # TODO
-                #=d, v = @constinferred eigh_trunc(t2; trunc = truncrank(nvals))
+                d, v = @constinferred eigh_trunc(t2; trunc = truncrank(nvals))
                 @test t2 * v ≈ v * d
-                @test dim(domain(d)) ≤ nvals=#
+                @test dim(domain(d)) ≤ nvals
             end
         end
 

test/cuda/tensors.jl

@@ -258,14 +258,14 @@ for V in spacelist
                 end
             end
         end
-        @timedtestset "Tensor conversion" begin # TODO adjoint conversion methods don't work yet
+        @timedtestset "Tensor conversion" begin
             W = V1 ⊗ V2
             t = @constinferred CUDA.randn(W ← W)
-            #@test typeof(convert(TensorMap, t')) == typeof(t) # TODO Adjoint not supported yet
+            @test typeof(convert(TensorMap, t')) == typeof(t)
             tc = complex(t)
             @test convert(typeof(tc), t) == tc
             @test typeof(convert(typeof(tc), t)) == typeof(tc)
-            # @test typeof(convert(typeof(tc), t')) == typeof(tc) # TODO Adjoint not supported yet
+            @test typeof(convert(typeof(tc), t')) == typeof(tc)
             @test Base.promote_typeof(t, tc) == typeof(tc)
             @test Base.promote_typeof(tc, t) == typeof(tc + t)
         end