让是肯定的并且正对角线。让是未知的方阵。考虑方程组:
是否可以在不必完全重新计算 LU 分解的情况下对所有这些方程进行数值求解对于每个?
我一直在搞乱 Cholesky 和 Eigen 分解,但到目前为止还没有运气,我认为这可能是不可能的,但我想我会问。
如果有帮助,请注意解决方案都将是肯定的,因为
编辑:
众所周知,Brian Borchers 的这个答案似乎意味着这实际上是不可能的。
让是肯定的并且正对角线。让是未知的方阵。考虑方程组:
是否可以在不必完全重新计算 LU 分解的情况下对所有这些方程进行数值求解对于每个?
我一直在搞乱 Cholesky 和 Eigen 分解,但到目前为止还没有运气,我认为这可能是不可能的,但我想我会问。
如果有帮助,请注意解决方案都将是肯定的,因为
编辑:
众所周知,Brian Borchers 的这个答案似乎意味着这实际上是不可能的。
您可以随时尝试共轭渐变(https://en.wikipedia.org/wiki/Conjugate_gradient_method)。一旦你通过 LU 找到了解决方案,如果新的变化只影响少量的方程(甚至是一个很大的方程),那么收敛应该非常快。这里清楚地解释了如何做到这一点:ftp: //ftp.numerical.rl.ac.uk/pub/talks/isd_stanford50.pdf