Expokit.jl
Expokit.jl copied to clipboard
expmv vs. expv
In Julia 1.0 the matrix exponential is done via exp
and not expm
as before. To have a consistent naming I would suggest to rename expmv
and phimv
.
But then what is the v for without the m?
Good question. To me, since exp
computes the exponential of something, expv
would give the exponential times a vector. But I don't know if that makes sense?
In the original Expokit package https://www.maths.uq.edu.au/expokit/download.html
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
||||||||||||||||||||||||||||||| MATLAB |||||||||||||||||||||||||||||||
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
chbv.m Chebyshev algorithm for w = exp(t*A)v
padm.m irreducible Pade algorithm for exp(t*A)
phiv.m Krylov algorithm for w = exp(t*A)*v + t*phi(t*A)u
expv.m Krylov algorithm for w = exp(t*A)*v
mexpv.m Markov algorithm for w = exp(t*A)*v
A better reason for renaming could be that we follow the original convention.