Aliases: expm expm,Matrix-method expm,dMatrix-method expm,ddiMatrix-method expm,dgeMatrix-method expm,dspMatrix-method expm,dsparseMatrix-method expm,dsyMatrix-method expm,dtpMatrix-method expm,dtrMatrix-method expm,matrix-method
### ** Examples (m1 <- Matrix(c(1,0,1,1), ncol = 2))
2 x 2 Matrix of class "dtrMatrix"
[,1] [,2]
[1,] 1 1
[2,] . 1
(e1 <- expm(m1)) ; e <- exp(1)
2 x 2 Matrix of class "dtrMatrix"
[,1] [,2]
[1,] 2.718282 2.718282
[2,] . 2.718282
stopifnot(all.equal(e1@x, c(e,0,e,e), tolerance = 1e-15)) (m2 <- Matrix(c(-49, -64, 24, 31), ncol = 2))
2 x 2 Matrix of class "dgeMatrix"
[,1] [,2]
[1,] -49 24
[2,] -64 31
(e2 <- expm(m2))
2 x 2 Matrix of class "dgeMatrix"
[,1] [,2]
[1,] -0.7357588 0.5518191
[2,] -1.4715176 1.1036382
(m3 <- Matrix(cbind(0,rbind(6*diag(3),0))))# sparse!
4 x 4 sparse Matrix of class "dtCMatrix"
[1,] . 6 . .
[2,] . . 6 .
[3,] . . . 6
[4,] . . . .
(e3 <- expm(m3)) # upper triangular
4 x 4 Matrix of class "dtrMatrix"
[,1] [,2] [,3] [,4]
[1,] 1 6 18 36
[2,] . 1 6 18
[3,] . . 1 6
[4,] . . . 1