demo/expm_check.r

library(rexpokit, quietly=TRUE)
require(Matrix)

set.seed(1234)

n <- 5
x <- matrix(rnorm(n*n), nrow=n, ncol=n)

A <- rexpokit::expm(x, t=1)
B <- Matrix::expm(x)@x
dim(B) <- c(n, n)

all.equal(A, B)
nmatzke/rexpokit documentation built on Nov. 28, 2023, 9:35 p.m.