Julia: isposdef()가 잘못되었습니다

isposdef()는 일부(전부는 아니지만) 양의 정부호 행렬이 다음과 같을 때 양의 정부호 행렬이라고 주장합니다.

julia> A = [1.69  2.21; 2.21  2.89]
2×2 Array{Float64,2}:
 1.69  2.21
 2.21  2.89
julia> det(A)
0.0
julia> isposdef(A)
true
julia> isposdef(ones(2,2))
false

이 동작의 원인이 무엇인지 확실하지 않습니다. isposdef()가 행렬의 Cholesky 분해를 시도하는 것처럼 보이며 모든 양의 준정부호 행렬에는 Cholesky 분해가 있습니다(IIRC에는 이를 찾는 일반적인 알고리즘이 없음) . 다음과 같이 Sylvestor의 기준을 구현하는 것이 더 간단할 것 같습니다.

sylvestor(A) = issymmetric(A) & all([det(A[1:j, 1:j]) > 0 for j = 1:size(A)[1]])

버전 정보:

julia> versioninfo()
Julia Version 0.5.0
Commit 3c9d753 (2016-09-19 18:14 UTC)
Platform Info:
  System: NT (x86_64-w64-mingw32)
  CPU: Intel(R) Xeon(R) CPU E3-1240 v3 @ 3.40GHz
  WORD_SIZE: 64
  BLAS: libopenblas (USE64BITINT DYNAMIC_ARCH NO_AFFINITY Haswell)
  LAPACK: libopenblas64_
  LIBM: libopenlibm
  LLVM: libLLVM-3.7.1 (ORCJIT, haswell)