R/expm.r

In nmatzke/rexpokit: R Wrappers for EXPOKIT; Other Matrix Functions

Defines functions expm

Documented in expm

expm <- function(x, t=1)
{
  if (nrow(x) != ncol(x))
    stop("Matrix exponentiation is only defined for square matrices.")
  
  ret <- expokit_dgpadm_Qmat(Qmat=x, t=t, transpose_needed=TRUE)
  
  return( ret )
}