2009-11-05 123 views
2

根据我的理解,需要分解/因式分解(LU,QR,Cholesky等),接下来是基于因式分解的矩阵反演计算。有没有其他的方法来解决它(我试图弄清楚我是否可以坚持在CULAtools的试用版中免费赠送的6个功能)?预先感谢您的意见。Lapack调用反转矩阵

+0

看看MAGMA项目 - 它为混合CPU/GPU系统提供LAPACK例程 - 是开源的 – prettygully 2013-02-19 22:38:54

回答

3

当然,找到Adjugate矩阵;这是一个反转小基础的简单方法。辅助矩阵只是辅因子矩阵的转置矩阵,而矩阵矩阵的逆矩阵就是由(标量)行列式分割的偶数矩阵。如果他们不熟悉,请在Wikipedia上查找这些条款。

如果你正在工作瓦特/大矩阵,我会买包。

保罗

5

的LAPACK例程计算矩阵逆是xyyTRI,其中x指示数据类型(“S”为单精度实,“d”的双精度实,“C”对于单精度复杂,'Z'表示双精度复数),yy表示矩阵的类型(对于非对称矩阵的一般情况,'GE';对于其他矩阵类型,有20多个其他双字母代码)。对于实数矩阵,通常使用DGETRI,对于复数矩阵,通常使用ZGETRI