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 c0d06d3 commit 46cf572Copy full SHA for 46cf572
stdlib/LinearAlgebra/src/bidiag.jl
@@ -886,7 +886,7 @@ eigvals(M::Bidiagonal) = M.dv
886
function eigvecs(M::Bidiagonal{T}) where T
887
n = length(M.dv)
888
Q = Matrix{T}(undef, n,n)
889
- blks = [0; findall(x -> x == 0, M.ev); n]
+ blks = [0; findall(iszero, M.ev); n]
890
v = zeros(T, n)
891
if M.uplo == 'U'
892
for idx_block = 1:length(blks) - 1, i = blks[idx_block] + 1:blks[idx_block + 1] #index of eigenvector
0 commit comments