More tests

Author: Katharine Hyatt

src/MatrixAlgebraKit.jl

@@ -19,6 +19,7 @@ export eigh_full, eigh_vals, eigh_trunc
 export eigh_full!, eigh_vals!, eigh_trunc!
 export eig_full, eig_vals, eig_trunc
 export eig_full!, eig_vals!, eig_trunc!
+export gen_eig, gen_eig!
 export gen_eig_full, gen_eig_vals
 export gen_eig_full!, gen_eig_vals!
 export schur_full, schur_vals

test/gen_eig.jl

@@ -14,7 +14,7 @@ using LinearAlgebra: Diagonal
         B_init = copy(B)
         Tc = complex(T)
 
-        W, V = @constinferred gen_eig_full(A, B; alg=($alg))
+        W, V = @constinferred gen_eig(A, B; alg=($alg))
         @test eltype(W) == eltype(V) == Tc
         @test A_init == A
         @test B_init == B
@@ -24,7 +24,7 @@ using LinearAlgebra: Diagonal
 
         Ac = similar(A)
         Bc = similar(B)
-        W2, V2 = @constinferred gen_eig_full!(copy!(Ac, A), copy!(Bc, B), (W, V), alg′)
+        W2, V2 = @constinferred gen_eig!(copy!(Ac, A), copy!(Bc, B), (W, V), alg′)
         @test W2 === W
         @test V2 === V
         @test A * V ≈ B * V * Diagonal(W)
@@ -51,4 +51,12 @@ using LinearAlgebra: Diagonal
     # a tuple of the input types is passed to `default_algorithm`
     @test_throws MethodError MatrixAlgebraKit.default_algorithm(gen_eig_full, A, B)
     @test_throws MethodError MatrixAlgebraKit.default_algorithm(gen_eig_vals, A, B)
+    if T <: Real
+        # Float16 isn't supported
+        Afp16 = Float16.(A)
+        Bfp16 = Float16.(B)
+        @test_throws MethodError MatrixAlgebraKit.default_algorithm(gen_eig_full, (Afp16, Bfp16))
+        @test_throws MethodError MatrixAlgebraKit.default_gen_eig_algorithm(Afp16, Bfp16)
+    end
 end
+