Doolittle分解的算法实现问题

Yjc

各位老师好，我在写用“Doolittle分解法解线性方程组”这一数值算法的实现程序时遇到了问题，程序给出了错误的结果，想请教各位老师这个程序在算法上有什么错误，谢谢！
（附件中是实现doolittle分解法的子程序，这是我写的交互程序的一部分，A,X,B分别代表系数矩阵，解，常数矩阵；ndim代表方阵A的维度）
Doolittle_Decomposition.f90 (2.01 KB, 下载次数: 2)

2024-3-30 11:40 上传

点击文件名下载附件
源代码

Yjc

感谢老师解答，我的想法是在LU分解中结合列选主元，把每列中最大的元素放到主元所在的位置，以增加数值稳定性

Yjc

我发现问题所在了，在选列主元的时候，应该同时对矩阵A和矩阵L进行行变换，使得主元位置元素的绝对值最大。