Closed mforets closed 6 years ago
unless you had sth else mind, i'd be happy to continue with this issue.
Please go ahead 🙂 The only thing I thought was that it might be better to start from the Julia implementation of expmv
instead of the original phimv
in EXPOKIT?
So i started trying the minimal changes for expmv
to get phimv
, but failed.
Then i stepped back and wrote phimv
using your notation, but still using Matlab's k1
and mb
variables. Do we want to get rid of k1
? I guess we can continue the discussion in PR #18
The
phimv
function, phiv in Expokit, computes the solution of a nonhomogeneous linear ODE without direct evaluation of a matrix exponential.