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`.