define default algorithm for base wrapper types

Author: lkdvos

src/interface/eig.jl

@@ -161,12 +161,18 @@ See also [`eig_full(!)`](@ref eig_full) and [`eig_trunc(!)`](@ref eig_trunc).
 # -------------------
 default_eig_algorithm(A; kwargs...) = default_eig_algorithm(typeof(A); kwargs...)
 default_eig_algorithm(T::Type; kwargs...) = throw(MethodError(default_eig_algorithm, (T,)))
-function default_eig_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasMat}
+function default_eig_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasVecOrMat}
     return LAPACK_Expert(; kwargs...)
 end
 function default_eig_algorithm(::Type{T}; kwargs...) where {T <: Diagonal}
     return DiagonalAlgorithm(; kwargs...)
 end
+function default_eig_algorithm(::Type{<:Base.ReshapedArray{T, N, A}}) where {T, N, A}
+    return default_eig_algorithm(A)
+end
+function default_eig_algorithm(::Type{SubArray{T, N, A}}) where {T, N, A}
+    return default_eig_algorithm(A)
+end
 
 for f in (:eig_full!, :eig_vals!)
     @eval function default_algorithm(::typeof($f), ::Type{A}; kwargs...) where {A}

src/interface/eigh.jl

@@ -167,12 +167,18 @@ default_eigh_algorithm(A; kwargs...) = default_eigh_algorithm(typeof(A); kwargs.
 function default_eigh_algorithm(T::Type; kwargs...)
     throw(MethodError(default_eigh_algorithm, (T,)))
 end
-function default_eigh_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasMat}
+function default_eigh_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasVecOrMat}
     return LAPACK_MultipleRelativelyRobustRepresentations(; kwargs...)
 end
 function default_eigh_algorithm(::Type{T}; kwargs...) where {T <: Diagonal}
     return DiagonalAlgorithm(; kwargs...)
 end
+function default_eigh_algorithm(::Type{<:Base.ReshapedArray{T, N, A}}) where {T, N, A}
+    return default_eigh_algorithm(A)
+end
+function default_eigh_algorithm(::Type{SubArray{T, N, A}}) where {T, N, A}
+    return default_eigh_algorithm(A)
+end
 
 for f in (:eigh_full!, :eigh_vals!)
     @eval function default_algorithm(::typeof($f), ::Type{A}; kwargs...) where {A}

src/interface/lq.jl

@@ -75,12 +75,18 @@ 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.MaybeBlasMat}
+function default_lq_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasVecOrMat}
     return LAPACK_HouseholderLQ(; kwargs...)
 end
 function default_lq_algorithm(::Type{T}; kwargs...) where {T <: Diagonal}
     return DiagonalAlgorithm(; kwargs...)
 end
+function default_lq_algorithm(::Type{<:Base.ReshapedArray{T, N, A}}) where {T, N, A}
+    return default_lq_algorithm(A)
+end
+function default_lq_algorithm(::Type{SubArray{T, N, A}}) where {T, N, A}
+    return default_lq_algorithm(A)
+end
 
 for f in (:lq_full!, :lq_compact!, :lq_null!)
     @eval function default_algorithm(::typeof($f), ::Type{A}; kwargs...) where {A}

src/interface/qr.jl

@@ -75,12 +75,18 @@ 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.MaybeBlasMat}
+function default_qr_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasVecOrMat}
     return LAPACK_HouseholderQR(; kwargs...)
 end
 function default_qr_algorithm(::Type{T}; kwargs...) where {T <: Diagonal}
     return DiagonalAlgorithm(; kwargs...)
 end
+function default_qr_algorithm(::Type{<:Base.ReshapedArray{T, N, A}}) where {T, N, A}
+    return default_qr_algorithm(A)
+end
+function default_qr_algorithm(::Type{SubArray{T, N, A}}) where {T, N, A}
+    return default_qr_algorithm(A)
+end
 
 for f in (:qr_full!, :qr_compact!, :qr_null!)
     @eval function default_algorithm(::typeof($f), ::Type{A}; kwargs...) where {A}

src/interface/svd.jl

@@ -161,12 +161,18 @@ default_svd_algorithm(A; kwargs...) = default_svd_algorithm(typeof(A); kwargs...
 function default_svd_algorithm(T::Type; kwargs...)
     throw(MethodError(default_svd_algorithm, (T,)))
 end
-function default_svd_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasMat}
+function default_svd_algorithm(::Type{T}; kwargs...) where {T <: YALAPACK.MaybeBlasVecOrMat}
     return LAPACK_DivideAndConquer(; kwargs...)
 end
 function default_svd_algorithm(::Type{T}; kwargs...) where {T <: Diagonal}
     return DiagonalAlgorithm(; kwargs...)
 end
+function default_svd_algorithm(::Type{<:Base.ReshapedArray{T, N, A}}) where {T, N, A}
+    return default_svd_algorithm(A)
+end
+function default_svd_algorithm(::Type{SubArray{T, N, A}}) where {T, N, A}
+    return default_svd_algorithm(A)
+end
 
 for f in (:svd_full!, :svd_compact!, :svd_vals!)
     @eval function default_algorithm(::typeof($f), ::Type{A}; kwargs...) where {A}