計算復矩陣的行列式fortran90

我需要計算對稱的復矩陣的行列式。矩陣的大小範圍從500 * 500到2000 * 2000。有我打電話的子程序嗎？順便說一句，我用ifort編譯。計算復矩陣的行列式fortran90

2013-09-26 someone

你確定你確實需要計算行列式嗎？很可能你想要一個相關的數量，如條件編號。在LAPACK中存在幾個相關例程（這包括在英特爾的MKL庫中，如果您有ifort，則可能有這些例程）。根據LAPACK FAQ，使用決定因素存在穩定性問題http://www.netlib.org/lapack/faq.html#_are_there_routines_in_lapack_to_compute_determinants – DaveP

最簡單的方法是做一個LU分解as described here。我會建議使用LAPACK這個任務......

這article在C一些代碼做，對於一個實數對稱矩陣，所以你需要通過zspsv交換dspsv處理雙精度複雜基質。

來源

2013-09-26 05:22:35

計算復矩陣的行列式fortran90

回答

相關問題