expoRkit (version 0.9.4)

orani: Australia Economic Model data, 1968-68.

Description

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.

References

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

See Also

expv

Examples

# NOT RUN {
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)
# }