orani: Australia Economic Model data, 1968-68.
In expoRkit: Expokit in R

Description Author(s) References See Also Examples

This sparse matrix of order 2,529 and with 90,158 non-zero elements was used in Sidje (1998) to illustrate the application of Expokit.

Niels Richard Hansen Niels.R.Hansen@math.ku.dk

Sidje, R. B. (1998) Expokit. Software Package for Computing Matrix Exponentials. ACM Trans. Math. Softw. 24(1), 130-156.

expv

data(orani)  ## Load the data as a 'dgCMatrix' (CCS format)
v <- rep(1, 2529)
### Solving a system of 2529 coupled linear differential equations
system.time(wCCS <- expv(orani, v = v, t = 10))
oraniCOO <- as(orani, "TsparseMatrix")  ## Coerce to COO format
### Sidje (1998) reported that the COO format gives a slight increase
### in computational time. Wether this is the case for this example
### depends on the machine used.
system.time(wCOO <- expv(oraniCOO, v = v, t = 10))

print(cbind(wCCS[1:5], wCOO[1:5]), digits = 14)