Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

RFC: sort eigenvalues in a canonical order #21598

Merged
merged 15 commits into from
Feb 5, 2019
Merged

RFC: sort eigenvalues in a canonical order #21598

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
726b309
sort eigenvalues in a canonical order
stevengj Jan 10, 2019
c493100
no sortby for special matrix types
stevengj Jan 10, 2019
75cc29d
simplification
stevengj Jan 10, 2019
e639171
fixes
stevengj Jan 24, 2019
c0576ea
pass through keyword args from eigen to eigen! and similar
stevengj Jan 31, 2019
8c6072e
more fixes
stevengj Feb 1, 2019
bc1ce98
whoops
stevengj Feb 1, 2019
167424e
fix bidiag test
stevengj Feb 1, 2019
c3fc98d
simplify lapack test fixes
stevengj Feb 1, 2019
79fb517
don't discard sortby if the matrix happens to be hermitian
stevengj Feb 1, 2019
cf08bda
missing check in permutecols
stevengj Feb 1, 2019
e37586e
Merge branch 'master' into eigsort
stevengj Feb 1, 2019
d6acbb0
test fixes
stevengj Feb 4, 2019
118bcf5
Merge branch 'eigsort' of git://github.com/stevengj/julia into eigsort
stevengj Feb 4, 2019
1b913d5
NEWS [ci skip]
stevengj Feb 5, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
fixes
  • Loading branch information
@stevengj
stevengj committed Jan 24, 2019
commit e63917133a3c29018195be10ddb054d59ad57e3a
9 changes: 5 additions & 4 deletions base/combinatorics.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,9 +66,10 @@ isperm(p::Tuple{Int,Int}) = ((p[1] == 1) & (p[2] == 2)) | ((p[1] == 2) & (p[2] =
# swap columns i and j of a, in-place
function swapcols!(a::AbstractMatrix, i, j)
i == j && return
cols = indices(a,2)
(i in cols && j in cols) || throw(BoundsError())
for k in indices(a,1)
cols = axes(a,2)
@boundscheck i in cols || throw(BoundsError(a, (:,i)))
@boundscheck j in cols || throw(BoundsError(a, (:,j)))
for k in axes(a,1)
@inbounds a[k,i],a[k,j] = a[k,j],a[k,i]
end
end
Expand All @@ -77,7 +78,7 @@ function permutecols!!(a::AbstractMatrix, p::AbstractVector{<:Integer})
count = 0
start = 0
while count < length(p)
ptr = start = findnext(p, start+1)
ptr = start = findnext(!iszero, p, start+1)::Int
next = p[start]
count += 1
while next != start
Expand Down
4 changes: 2 additions & 2 deletions stdlib/LinearAlgebra/docs/src/index.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,13 +38,13 @@ julia> A = [-4. -17.; 2. 2.]

julia> eigvals(A)
2-element Array{Complex{Float64},1}:
-1.0 + 5.0im
-1.0 - 5.0im
-1.0 + 5.0im

julia> eigvecs(A)
2×2 Array{Complex{Float64},2}:
0.945905+0.0im 0.945905-0.0im
-0.166924-0.278207im -0.166924+0.278207im
-0.166924+0.278207im -0.166924-0.278207im
```

In addition, Julia provides many [factorizations](@ref man-linalg-factorizations) which can be used to
Expand Down