tests/expm.r

In wrathematics/libexpm: Matrix Exponentiation

library(libexpm)

n <- 250
n <- 10
x <- matrix(rnorm(n*n), n, n)

t1 <- system.time(exp1 <- expm(x))[3]

library(Matrix)
t2 <- system.time(exp2 <- Matrix::expm(x))[3]

t1
t2

all.equal(exp1, as.matrix(exp2), check.attributes=FALSE)