Julia中LowRankApprox的问题

Question

我正在尝试使用Julia v0.6.0中LowRankApprox.jl包提供的pheigfact函数进行Hermitian特征分解。 基本上，这只是一行代码，例如：

(E, F) = pheigfact(A);

其中A是实对称正定矩阵。 但是，出现以下错误：

MethodError: no method matching
start(::LowRankApprox.PartialHermitianEigen{Float64,Float64})
Closest candidates are:
  start(!Matched::SimpleVector) at essential.jl:258
  start(!Matched::Base.MethodList) at reflection.jl:560
  start(!Matched::ExponentialBackOff) at error.jl:107 

感谢任何帮助！

Answer 1

TL; DR

使用函数pheig not pheigfact返回值和向量的元组

---

完整答案

我没有这个包，但是从文档看来， pheigfact返回一个元素，您可以使用getindex(x,ind::Symbol)从中访问值/向量。

例如

F = pheigfact(A)
values=F[:values]
vectors=F[:vectors]

如果您尝试为管分配单个元素，它将尝试遍历不支持该元素的类型，从而给您带来错误（即，该类型没有方法start ）。 我可能在执行x,y = :onetwo或start(:onetwo)遇到类似的错误

解

使用pheig函数返回一个元组。

E, F = pheig(A)