Add a few more Diagonal tests for coverage

Author: kshyatt

src/TensorKit.jl

@@ -76,6 +76,8 @@ export leftorth, rightorth, leftnull, rightnull,
        qr_full, qr_compact, qr_null, lq_full, lq_compact, lq_null,
        qr_full!, qr_compact!, qr_null!, lq_full!, lq_compact!, lq_null!,
        tsvd!, tsvd, eigen, eigen!, eig, eig!, eigh, eigh!, exp, exp!,
+       eigh_full!, eigh_full, eig_full!, eig_full, eigh_vals!, eigh_vals,
+       eig_vals!, eig_vals,
        isposdef, isposdef!, ishermitian, isisometry, isunitary, sylvester, rank, cond
 export braid, braid!, permute, permute!, transpose, transpose!, twist, twist!, repartition,
        repartition!

test/diagonal.jl

@@ -188,14 +188,19 @@ diagspacelist = ((ℂ^4)', ℂ[Z2Irrep](0 => 2, 1 => 3),
     @timedtestset "Factorization" begin
         for T in (Float32, ComplexF64)
             t = DiagonalTensorMap(rand(T, reduceddim(V)), V)
-            @testset "eig" begin
+            @testset "eig/eigh" begin
+                D, W = @constinferred eig_full(t)
+                @test t * W ≈ W * D
                 D, W = @constinferred eig(t)
                 @test t * W ≈ W * D
                 t2 = t + t'
                 D2, V2 = @constinferred eigh(t2)
                 VdV2 = V2' * V2
                 @test VdV2 ≈ one(VdV2)
                 @test t2 * V2 ≈ V2 * D2
+                
+                D3 = @constinferred eigh_vals(t2)
+                @test D2 ≈ D3
 
                 @test rank(D) ≈ rank(t)
                 @test cond(D) ≈ cond(t)
@@ -231,6 +236,26 @@ diagspacelist = ((ℂ^4)', ℂ[Z2Irrep](0 => 2, 1 => 3),
                     @test isposdef(L)
                 end
             end
+            @testset "qr_full with $alg" for alg in (TensorKit.DiagonalAlgorithm(),)
+                Q, R = @constinferred qr_full(t; alg=alg)
+                @test isisometry(Q)
+                @test Q * R ≈ t
+            end
+            @testset "qr_compact with $alg" for alg in (TensorKit.DiagonalAlgorithm(),)
+                Q, R = @constinferred qr_compact(t; alg=alg)
+                @test isisometry(Q)
+                @test Q * R ≈ t
+            end
+            @testset "lq_full with $alg" for alg in (TensorKit.DiagonalAlgorithm(),)
+                L, Q = @constinferred lq_full(t; alg=alg)
+                @test isisometry(Q; side=:right)
+                @test L * Q ≈ t
+            end
+            @testset "lq_compact with $alg" for alg in (TensorKit.DiagonalAlgorithm(),)
+                L, Q = @constinferred lq_compact(t; alg=alg)
+                @test isisometry(Q; side=:right)
+                @test L * Q ≈ t
+            end
             @testset "tsvd with $alg" for alg in (TensorKit.DiagonalAlgorithm(),)
                 U, S, Vᴴ = @constinferred tsvd(t; alg=alg)
                 UdU = U' * U