Skip to content

Commit d735f8d

Browse files
committed
More fixes for Julia 1
1 parent 9ef2bc3 commit d735f8d

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

src/eigenSelfAdjoint.jl

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -21,7 +21,7 @@ function Base.getproperty(S::SymmetricTridiagonalFactorization, s::Symbol)
2121
if s == :Q
2222
return EigenQ(S.uplo, S.factors, S.τ)
2323
else
24-
throw(ArgumentError("no property $s"))
24+
return getfield(S, s)
2525
end
2626
end
2727

@@ -492,8 +492,8 @@ function symtriUpper!(AS::StridedMatrix{T},
492492
n = LinearAlgebra.checksquare(AS)
493493

494494
@inbounds for k = 1:(n - 2 + !(T<:Real))
495-
# FixMe! Avoid the deprecated ConjArray
496-
τk = LinearAlgebra.reflector!(LinearAlgebra._conj(view(AS, k, (k + 1):n)))
495+
# This is a bit convoluted method to get the conjugated vector but conjugation is required for correctness of arrays of quaternions. Eventually, it should be sufficient to write vec(x') but it currently (July 10, 2018) hits a bug in LinearAlgebra
496+
τk = LinearAlgebra.reflector!(vec(transpose(view(AS, k, (k + 1):n)')))
497497
τ[k] = τk'
498498

499499
for j in (k + 1):n
@@ -536,7 +536,7 @@ LinearAlgebra.eigvals!(A::SymTridiagonal ) = _eigvals!(A)
536536
LinearAlgebra.eigvals!(A::Hermitian ) = _eigvals!(A)
537537

538538
_eigen!(A::SymmetricTridiagonalFactorization) =
539-
LinearAlgebra.Eigen(eigQL!(A.diagonals, Array(getq(A)), eps(eltype(A.diagonals)), false)...)
539+
LinearAlgebra.Eigen(eigQL!(A.diagonals, Array(A.Q), eps(eltype(A.diagonals)), false)...)
540540
_eigen!(A::SymTridiagonal) =
541541
LinearAlgebra.Eigen(eigQL!(A, Matrix{eltype(A)}(I, size(A, 1), size(A, 1)), eps(eltype(A)), false)...)
542542
_eigen!(A::Hermitian) = _eigen!(symtri!(A))

0 commit comments

Comments
 (0)