我在C++代码中使用Octave的EIG
类。事实上,它效果很好。我唯一的小问题是,它根据EIG
构造函数的参数列表自动确定算法。我想知道是否有一种方法可以确定EIG
类的算法。我注意到有一个名为info
的可选参数,我查看了C++源代码(在下面的链接上),看看它是否与算法有关,但我不明白。
http://octave.sourceforge.net/doxygen/html/d6/d42/EIG_8cc_source.html
我在Visual Studio 2008中使用Octave 3.6.1。
我研究了这个主题。显然Octave没有选择来确定eig
命令的算法。但对于条件较差的矩阵,MATLAB有'chol'
算法和'qz'
算法两种选择。
如果你想计算特征值&在C++中,可以实现LAPACK的DGGEV或DGGEVX例程。