下一节:稀疏矩阵的直接求解器
上一级:直接线性求解器简述
上一节:稠密矩阵的直接求解器
带状矩阵的直接求解器
本节类似于稠密矩阵的§10.3.1部分:
- A-\sigma I 是厄米正定矩阵。
- 在这种情况下,Cholesky分解是首选算法。它在LAPACK计算例程中实现,
xPBTRF用于计算分解,xPBTRS用于使用分解求解(两者在LAPACK驱动例程xPBSVX中结合)。
Cholesky分解在类似的ScaLAPACK例程中实现,包括PxPBTRF、PxPBTRS和PxPBSV。
- A-\sigma I 是非厄米矩阵。
- 在这种情况下,高斯消去法是首选算法。它在LAPACK计算例程中实现,
xGBTRF用于计算分解,xGBTRS用于使用分解求解(两者在LAPACK驱动例程xGBSVX中结合)。
它在类似的ScaLAPACK例程中实现,包括带部分主元消去的PxGBTRF、PxGBTRS和PxGBSV,
以及不带主元消去的PxDBTRF、PxDBTRS和PxDBSV(风险较大但速度比使用主元消去快)。
没有利用厄米不定结构的例程可用(因为没有关于带宽因主元消去而增长多少的简单界限)。
MATLAB中没有带状矩阵的例程。如果A-\sigma I是厄米矩阵且带宽足够窄(如三对角矩阵),
应使用§4.2中的直接特征值求解方法。
下一节:稀疏矩阵的直接求解器
上一级:直接线性求解器简述
上一节:稠密矩阵的直接求解器
Susan Blackford
2000-11-20