add defaults and more tests

Author: Jutho

src/interface/lq.jl

@@ -72,6 +72,9 @@ default_lq_algorithm(A; kwargs...) = default_lq_algorithm(typeof(A); kwargs...)
 function default_lq_algorithm(T::Type; kwargs...)
     throw(MethodError(default_lq_algorithm, (T,)))
 end
+function default_lq_algorithm(::Type{T}; kwargs...) where {T <: AbstractMatrix}
+    return Native_HouseholderLQ(; kwargs...)
+end
 function default_lq_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.BlasMat}
     return LAPACK_HouseholderLQ(; kwargs...)
 end

src/interface/qr.jl

@@ -72,6 +72,9 @@ default_qr_algorithm(A; kwargs...) = default_qr_algorithm(typeof(A); kwargs...)
 function default_qr_algorithm(T::Type; kwargs...)
     throw(MethodError(default_qr_algorithm, (T,)))
 end
+function default_qr_algorithm(::Type{T}; kwargs...) where {T <: AbstractMatrix}
+    return Native_HouseholderQR(; kwargs...)
+end
 function default_qr_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.BlasMat}
     return LAPACK_HouseholderQR(; kwargs...)
 end

test/lq.jl

@@ -6,7 +6,7 @@ using LinearAlgebra: diag, I, Diagonal
 using MatrixAlgebraKit: LQViaTransposedQR, LAPACK_HouseholderQR
 
 lapack_eltypes = (Float32, Float64, ComplexF32, ComplexF64)
-native_eltypes = (lapack_eltypes..., BigFloat, Complex{BigFloat})
+native_eltypes = (lapack_eltypes..., Float16, BigFloat, Complex{BigFloat})
 
 @testset "lq_compact! for T = $T" for T in lapack_eltypes
     rng = StableRNG(123)
@@ -299,5 +299,10 @@ end
         @test Q2 === Q
         @test L * Q ≈ A
         @test isunitary(Q)
+
+        if !(T ∈ lapack_eltypes)
+            alg = MatrixAlgebraKit.select_algorithm(lq_compact!, A)
+            @test alg == Native_HouseholderLQ()
+        end
     end
 end

test/qr.jl

@@ -5,7 +5,7 @@ using StableRNGs
 using LinearAlgebra: diag, I, Diagonal
 
 lapack_eltypes = (Float32, Float64, ComplexF32, ComplexF64)
-native_eltypes = (lapack_eltypes..., BigFloat, Complex{BigFloat})
+native_eltypes = (lapack_eltypes..., Float16, BigFloat, Complex{BigFloat})
 
 @testset "qr_compact! and qr_null! for T = $T" for T in lapack_eltypes
     rng = StableRNG(123)
@@ -269,5 +269,10 @@ end
         @test R2 === R
         @test Q * R ≈ A
         @test isunitary(Q)
+
+        if !(T ∈ lapack_eltypes)
+            alg = MatrixAlgebraKit.select_algorithm(qr_compact!, A)
+            @test alg == Native_HouseholderQR()
+        end
     end
 end