线性方程组 · 把消元记录留下来
换一行,也要换它以前的记录
每一列都从当前尚未处理的行里选主元。看工作矩阵如何变化,再检查排列和 L 中已经保存的乘子有没有一起跟上。
| 被更新的量 | 原值 | 乘法结果 | 相减后的存储值 |
|---|
L:已经存下的乘子
当前 U 工作数组
增长因子 ρ
max |(PA−LU)ᵢⱼ|
已执行动作数
增长因子记录所有消元阶段出现过的最大系数绝对值 / 原矩阵最大系数绝对值,含原始阶段,排除右端。PA−LU 和最终残差都用原始数据与普通浮点重新核对,不拿已舍入的因子自证正确。
回代得到的解,放回原方程
下一次该把哪一行换到主元位置?
三位十进制模型使用 Number(value.toPrecision(3)),每次乘、减、加、除后单独舍入。它是低精度教学模型,不是 IEEE 二进制格式,也不承诺十进制平局按偶数舍入;本章关键算例不靠平局规则。