Skip to content

Commit

Permalink
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Fixed symmetric
Browse files Browse the repository at this point in the history
kshyatt committed Aug 13, 2015
1 parent c94fa7b commit df4c774
Showing 2 changed files with 10 additions and 6 deletions.
4 changes: 2 additions & 2 deletions base/linalg/symmetric.jl
Original file line number Diff line number Diff line change
@@ -57,9 +57,9 @@ ctranspose(A::Hermitian) = A
trace(A::Hermitian) = real(trace(A.data))

#tril/triu
tril(A::Hermitian,k::Integer=0) = tril(A.data,k)
tril(A::Hermitian,k::Integer=0) = triu(A.data,-k)'
triu(A::Hermitian,k::Integer=0) = triu(A.data,k)
tril(A::Symmetric,k::Integer=0) = tril(A.data,k)
tril(A::Symmetric,k::Integer=0) = triu(A.data,-k)'
triu(A::Symmetric,k::Integer=0) = triu(A.data,k)

## Matvec
12 changes: 8 additions & 4 deletions test/linalg/symmetric.jl
Original file line number Diff line number Diff line change
@@ -60,10 +60,14 @@ let n=10
@test ctranspose(Hermitian(asym)) == asym

#tril/triu
@test triu(Symmetric(asym),1) == triu(asym,1)
@test tril(Symmetric(asym),1) == tril(asym,1)
@test triu(Hermitian(asym),1) == triu(asym,1)
@test tril(Hermitian(asym),1) == tril(asym,1)
@test triu(Symmetric(asym),1) == triu(asym,1)
@test triu(Symmetric(asym),-1) == triu(asym,-1)
@test tril(Symmetric(asym),1) == tril(asym,1)
@test tril(Symmetric(asym),-1) == tril(asym,-1)
@test triu(Hermitian(asym),1) == triu(asym,1)
@test triu(Hermitian(asym),-1) == triu(asym,-1)
@test tril(Hermitian(asym),1) == tril(asym,1)
@test tril(Hermitian(asym),-1) == tril(asym,-1)

eltya == BigFloat && continue # Revisit when implemented in julia
d, v = eig(asym)

0 comments on commit df4c774

Please sign in to comment.