GPvecchia (version 0.1.7)

ichol: Wrapper for incomplete Cholesky decomposition

Description

Wrapper for incomplete Cholesky decomposition

Usage

ichol(M, S = NULL)

Value

the incomplete Cholesky factor in the sparse format

Arguments

M

the matrix to be decomposed

S

sparsity pattern matrix given

Examples

A = matrix(runif(25), ncol = 5)
A = t(A) * A + 2 * Matrix::Diagonal(5)
S = Matrix::Matrix(c(rep(1, 5), c(0, 1, 1, 0, 0), c(0, 0, 1, 0, 1),
c(0, 0, 0, 1, 0), c(0, 0, 0, 0, 1)), ncol = 5, byrow = TRUE)
I1 = ichol(A, S)
I2 = ichol(A * S)