pageexpm

Paged matrix exponential function

pageexpm extends the expm (matrix exponential) function for efficient operation on multicore and multiprocessor computers and GPUs by using paged matrix operations pagemtimes and pagemldivide.

For large-size, large-norm, matrices (2-D), pageexpm typically runs about 30% faster than expm and can exhibit order-of-magnitude accuracy improvement. The included demo_pageexpm.mlx illustrates the comparative performance of pageexpm and expm. Documentation of the underlying algorithms is cited in the code comments (https://vixra.org/abs/2203.0072).

pageexpm supports automatic differentiation via class mpoly (https://www.mathworks.com/matlabcentral/fileexchange/182185-class-mpoly-multivariate-polynomial).