在不完全重新计算每个方程的 LU 分解的情况下求解相关方程组

计算科学 线性代数 线性系统 矩阵方程
2021-12-18 07:30:17

Σ是肯定的并且Di正对角线。Xi是未知的方阵。考虑方程组:

(I+ΣDi)Xi=Σfori=1,...n.

是否可以在不必完全重新计算 LU 分解的情况下对所有这些方程进行数值求解(I+ΣDi)对于每个i?

我一直在搞乱 Cholesky 和 ​​Eigen 分解Σ,但到目前为止还没有运气,我认为这可能是不可能的,但我想我会问。

如果有帮助,请注意解决方案Xi都将是肯定的,因为

Xi1=Σ1+Di.

编辑:

众所周知,Brian Borchers 的这个答案似乎意味着这实际上是不可能的。

1个回答

您可以随时尝试共轭渐变(https://en.wikipedia.org/wiki/Conjugate_gradient_method)。一旦你通过 LU 找到了解决方案,如果新的变化只影响少量的方程(甚至是一个很大的方程),那么收敛应该非常快。这里清楚地解释了如何做到这一点:ftp: //ftp.numerical.rl.ac.uk/pub/talks/isd_stanford50.pdf