[問題] 科學計算函式庫GSL中eigen問題請教

作者: profyang (prof)   2017-07-22 16:35:51
開發平台(Platform): (Ex: Win10, Linux, ...)
Win7
編譯器(Ex: GCC, clang, VC++...)+目標環境(跟開發平台不同的話需列出)
MinGW
額外使用到的函數庫(Library Used): (Ex: OpenGL, ...)
GNU Scientific Library(GSL)
問題(Question):
主要是最近Matlab用膩了,想換換比較有深度(?)的,就回來用用C++
雖然以前C++學得很爛,但基本的迴圈之類的不太可能忘
而我主要會需要用到程式的地方其實就如標題所說:科學計算,所以也夠用了
只要有個科學計算的函式庫應該可以很快上手?
所以看到有GSL這種東西當然就很高興地想給他下載
但是看到他的Documentation就卡住了
https://www.gnu.org/software/gsl/doc/html/index.html
我有一個一定要用到的功能就是解eigen value & eigen vector
但是doc上面:
Eigensystems
Real Symmetric Matrices
Complex Hermitian Matrices
Real Nonsymmetric Matrices
Real Generalized Symmetric-Definite Eigensystems
Complex Generalized Hermitian-Definite Eigensystems
Real Generalized Nonsymmetric Eigensystems
Sorting Eigenvalues and Eigenvectors
Examples
References and Further Reading
似乎是說他提供解Eigensystems只有上面幾種特殊矩陣
但我的矩陣不會是以上幾種有這麼好性質的特殊矩陣
查了wiki似乎也是說對於任意矩陣求eigen沒有什麼好的加速演算法
等於只能硬是去解他的eigen polynomial的root (det(A-xI)=0的根)
但是我找了半天好像也沒找到GSL有可以自動幫我列出eigen polynomial的係數的...
也就是我得自己寫個程式把任意矩陣的eigen polynomial的係數算出來嗎...?
不知道有沒有大大有遇過類似情形的?還是是我哪裡搞錯了請指證!謝謝!
作者: sunev (Veritas)   2017-07-22 16:50:00
Real Nonsymmetric Matrices ?
作者: profyang (prof)   2017-07-22 18:05:00
很不幸的我的矩陣是complex...
作者: freeunixer (御劍客)   2017-07-22 18:48:00
直接用 python 吧,計算模組就是仿 matlab 的.
作者: longlongint (華哥爾)   2017-07-22 19:45:00
解 eigen value 當然要用 EIGEN 啊(誤https://goo.gl/pp1kqf
作者: soheadsome (師大狗鼻哥)   2017-07-23 08:55:00
用python吧 一堆用gpu最佳化的函式庫都有人幫你包好
作者: longlongint (華哥爾)   2017-07-23 12:45:00
推個python C語言是拿來寫別的語言的編譯器用的(誤
作者: gozule (好冷啊~~)   2017-07-26 19:12:00
C++封裝blas和lpack的函式庫Armadillo也可試試,很好上手http://arma.sourceforge.net/
作者: Neisseria (Neisseria)   2017-07-26 21:51:00
先用 Python 或 MATLAB 試模型會比較簡單忘記原 po 先前就是用 MATLAB,歹勢,請忽略小弟留言
作者: Killercat (殺人貓™)   2017-07-26 23:37:00
intel那個被告翻了 現在不敢再這樣搞了而且主要影響是TBL比較多 而非MKL抱歉 是TBB XD 不是TBL

Links booklink

Contact Us: admin [ a t ] ucptt.com