R updown-methods 更新和降級稀疏 Cholesky 分解

說明

計算稀疏 Cholesky 分解的排名 - k 更新或降級

對於某些 k 列矩陣 C 來說是因式分解

此處，s = 1 表示更新，s = -1 表示降級。

用法

updown(update, C, L)

參數

值

維度與 L 匹配的稀疏 Cholesky 分解，通常屬於 dCHMsimpl 類。

例子

m <- sparseMatrix(i = c(3, 1, 3:2, 2:1), p = c(0:2, 4, 4, 6), x = 1:6,
                  dimnames = list(LETTERS[1:3], letters[1:5]))
uc0 <- Cholesky(A <- crossprod(m) + Diagonal(5))
uc1 <- updown("+", Diagonal(5, 1), uc0)
uc2 <- updown("-", Diagonal(5, 1), uc1)
stopifnot(all.equal(uc0, uc2))

作者

Initial implementation by Nicholas Nagle, University of Tennessee.

參考

Davis, T. A., Hager, W. W. (2001). Multiple-rank modifications of a sparse Cholesky factorization. SIAM Journal on Matrix Analysis and Applications, 22(4), 997-1013. doi:10.1137/S0895479899357346

也可以看看

類 dCHMsimpl 和 dCHMsuper 及其方法，特別是通用函數 update ，它不等於 updown(update = TRUE) 。

通用函數Cholesky。