计算复数矩阵fortran90的行列式

Question

我需要计算对称的复矩阵的行列式。矩阵的大小范围从500 * 500到2000 * 2000。我有什么子程序可以打电话吗？顺便说一句，我使用ifort进行编译。

Answer 1

最简单的方法是进行LU分解as described here。我建议使用LAPACK来完成这项任务......

此article在C中有一些代码用于实值对称矩阵，因此您需要通过dspsv交换zspsv来处理双精度复数矩阵。