如果我有一个 psd,对称矩阵我需要做 LU decomps on(在哪里是一个对角 psd 矩阵,其中随机变化)。有没有一种简单的方法可以从 LU 出发+ 的一些功能得到 LU?
文献中的任何参考资料将不胜感激。
如果我有一个 psd,对称矩阵我需要做 LU decomps on(在哪里是一个对角 psd 矩阵,其中随机变化)。有没有一种简单的方法可以从 LU 出发+ 的一些功能得到 LU?
文献中的任何参考资料将不胜感激。
使用 Cholesky 分解或 LDL 分解代替 LU。
从以下答案来看:
没有好的方法可以更新任何这些分解。除非对角矩阵的秩非常低,否则每次都从头开始重新计算分解不会更糟,在这种情况下,也许您可以执行类似于 Sherman-Morrison(-Woodbury) 类型的更新。