如何在Stan中得到非对称矩阵的特征向量？

Question

Stan提供的功能

vector eigenvalues_sym(matrix A)
matrix eigenvectors_sym(matrix A)

获取对称矩阵A的特征值和特征向量，但是如果你的矩阵是非对称的呢？你还能得到特征值和特征向量吗？还有，有没有办法测试我是否获得了正确的特征向量？如何测试Stan程序的各个部分？

Answer 1

Stan只提供对称案例，因为我们无法保证答案不会复杂。

虽然可以将复合体编码为真实对，并开始在Stan中实现其中一些，但这将是一个巨大的痛苦。 我担心我们只提供对称案例，因为我们无法保证答案是真实的，我们也不会支持。其他任何东西都必须在C ++或Stan中从头开始实现。

Ben Goodrich建议在C ++中使用一种解决方法，需要从源处重新编译Stan，因为非对称矩阵已知会产生真正的特征分解，但它几乎不是语言内的解决方案而不是我们真正推荐人们做的事情，因为它和＃39;每次更新Stan时都需要更新。