We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent b431b32 commit bfc594aCopy full SHA for bfc594a
src/factorizations/svd.jl
@@ -11,9 +11,11 @@ struct BlockPermutedDiagonalAlgorithm{A<:MatrixAlgebraKit.AbstractAlgorithm} <:
11
alg::A
12
end
13
14
+# TODO: this is a hardcoded for now to get around this function not being defined in the
15
+# type domain
16
function MatrixAlgebraKit.default_svd_algorithm(A::AbstractBlockSparseMatrix; kwargs...)
- @assert blocktype(A) <: StridedMatrix{<:LinearAlgebra.BLAS.BlasFloat} "unsupported type:
- $(blocktype(A))"
17
+ blocktype(A) <: StridedMatrix{<:LinearAlgebra.BLAS.BlasFloat} ||
18
+ error("unsupported type: $(blocktype(A))")
19
alg = MatrixAlgebraKit.LAPACK_DivideAndConquer(; kwargs...)
20
return BlockPermutedDiagonalAlgorithm(alg)
21
0 commit comments