Most of the matrix decompositions in Julia can be applied to matrices with element types that are subtypes of AbstractFloat. For example
A = BigFloat.(rand(3,3))
qr(A)
lu(A)
svd(A)
schur(A)
all work fine. However,
gives a method error because lq() seems to require eltype(A) to be either Float16, Float32, or Float64.
Can the lq() function be generalized so that it can be applied to matrices whose elements are a subtype of AbstractFloat?
Most of the matrix decompositions in Julia can be applied to matrices with element types that are subtypes of AbstractFloat. For example
all work fine. However,
gives a method error because lq() seems to require eltype(A) to be either Float16, Float32, or Float64.
Can the lq() function be generalized so that it can be applied to matrices whose elements are a subtype of AbstractFloat?