R/TDtT.R

Defines functions TDtT

TDtT <- function(A, ...){
    ch <- chol(A)
    dd <- diag(ch)
 return(list(T = t(drop0(zapsmall(ch / dd, ...))), D = Diagonal(nrow(A), dd^2)))
}
matthewwolak/nadiv documentation built on July 7, 2023, 1:24 p.m.