Add some tests

Author: mtfishman

test/Project.toml

@@ -1,5 +1,6 @@
 [deps]
 Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
+FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
 KroneckerArrays = "05d0b138-81bc-4ff7-84be-08becefb1ccc"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 MatrixAlgebraKit = "6c742aac-3347-4629-af66-fc926824e5e4"
@@ -9,6 +10,7 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
 Aqua = "0.8"
+FillArrays = "1"
 KroneckerArrays = "0.1"
 LinearAlgebra = "1.10"
 MatrixAlgebraKit = "0.2"

test/test_matrixalgebrakit.jl

@@ -1,5 +1,6 @@
+using FillArrays: Eye
 using KroneckerArrays: ⊗
-using LinearAlgebra: Hermitian, I, diag, norm
+using LinearAlgebra: Hermitian, I, diag, hermitianpart, norm
 using MatrixAlgebraKit:
   eig_full,
   eig_trunc,
@@ -23,28 +24,35 @@ using MatrixAlgebraKit:
   svd_vals
 using Test: @test, @test_throws, @testset
 
+herm(a) = hermitianpart(a).data
+
 @testset "MatrixAlgebraKit" begin
   elt = Float32
 
   a = randn(elt, 2, 2) ⊗ randn(elt, 3, 3)
   d, v = eig_full(a)
   @test a * v ≈ v * d
 
+  a = randn(elt, 2, 2) ⊗ Eye(3)
+  d, v = eig_full(a)
+  @test d.b == v.b == Eye(3)
+  @test a * v ≈ v * d
+
   a = randn(elt, 2, 2) ⊗ randn(elt, 3, 3)
   @test_throws MethodError eig_trunc(a)
 
   a = randn(elt, 2, 2) ⊗ randn(elt, 3, 3)
   d = eig_vals(a)
   @test d ≈ diag(eig_full(a)[1])
 
-  a = Hermitian(randn(elt, 2, 2)) ⊗ Hermitian(randn(elt, 3, 3))
+  a = herm(randn(elt, 2, 2)) ⊗ herm(randn(elt, 3, 3))
   d, v = eigh_full(a)
   @test a * v ≈ v * d
 
-  a = Hermitian(randn(elt, 2, 2)) ⊗ Hermitian(randn(elt, 3, 3))
+  a = herm(randn(elt, 2, 2)) ⊗ herm(randn(elt, 3, 3))
   @test_throws MethodError eigh_trunc(a)
 
-  a = Hermitian(randn(elt, 2, 2)) ⊗ Hermitian(randn(elt, 3, 3))
+  a = herm(randn(elt, 2, 2)) ⊗ herm(randn(elt, 3, 3))
   d = eigh_vals(a)
   @test d ≈ diag(eigh_full(a)[1])
 
@@ -102,6 +110,11 @@ using Test: @test, @test_throws, @testset
   @test collect(u'u) ≈ I
   @test collect(v * v') ≈ I
 
+  a = randn(elt, 2, 2) ⊗ Eye(3)
+  u, s, v = svd_compact(a)
+  @test u * s * v ≈ a
+  @test u.b == s.b == v.b == Eye(3)
+
   a = randn(elt, 2, 2) ⊗ randn(elt, 3, 3)
   u, s, v = svd_full(a)
   @test u * s * v ≈ a