"Cholesky"
class is the class of Cholesky
decompositions of positive-semidefinite, real matrices. The
"BunchKaufman"
class is the class of Bunch-Kaufman
decompositions of symmetric, real matrices. The "pCholesky"
and "pBunchKaufman"
classes are their packed
storage versions.new("Cholesky",
...)
or new("BunchKaufman", ...)
, etc,
or rather by calls of the form chol(pm)
or
BunchKaufman(pm)
where pm
inherits from
the "dpoMatrix "
or "dsyMatrix "
class or as a side-effect of other functions
applied to "dpoMatrix"
objects (see dpoMatrix
)."MatrixFactorization"
and "dtrMatrix"
, directly.
Class "dgeMatrix"
, by class "dtrMatrix"
.
Class "Matrix"
, by class "dtrMatrix"
.dtrMatrix
, dpoMatrix
;
function chol
.(sm <- as(as(Matrix(diag(5) + 1), "dsyMatrix"), "dspMatrix"))
signif(csm <- chol(sm), 4)
(pm <- crossprod(Matrix(rnorm(18), nrow = 6, ncol = 3)))
(ch <- chol(pm))
if (toupper(ch@uplo) == "U") # which is TRUE
crossprod(ch)
stopifnot(all.equal(as(crossprod(ch), "matrix"),
as(pm, "matrix"), tol=1e-14))
Run the code above in your browser using DataLab