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.
1 2 3 4 5 6 7 8 9 10 11 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.