解线性方程组的直接方法 --- 选主元的三角分解法（C语言实现）

  从直接三角分解法公式能够看出，当U矩阵中的Urr=0时，计算将中断。或者当Urr绝对值很小时，按分解公式计算可能引发舍入偏差的累积。code 当A满秩的时候（非奇异）， 咱们能够经过交换A的行实现矩阵PA的LU分解，所以能够采用与列主元消去法相似的方法， 将直接三角分解法修改成（部分）选主元的三角分解法。对象 void GaussTri(double (*Arr)[3], int N){

>>阅读原文<<