tests/inverse.R

 set.seed(42)

 library(aster)

 a <- matrix(rnorm(1000), 5)
 a <- apply(a, 1, cumsum)
 a <- var(a)

 b <- .Call(aster:::C_pos_def_mat_inv, a)

 all.equal(b %*% a, diag(nrow(a)))

Try the aster package in your browser

Any scripts or data that you put into this service are public.

aster documentation built on June 8, 2025, 1:19 p.m.