Non-allocating versions of expm and phim

acroy / Expokit.jl

Julia implementation of EXPOKIT routines

Other

26 stars 11 forks source link

Non-allocating versions of expm and phim #19

Open ChrisRackauckas opened 6 years ago

ChrisRackauckas commented 6 years ago

It would be nice if there could be an optional argument for the Krylov vectors so that way if this is called in a loop the caches don't need to be created each time.

mforets commented 6 years ago

It looks like this feature has been implemented recently in OrdinaryDiffEq.

ChrisRackauckas commented 6 years ago

Yup it has been. We may want to refactor that part out later, but for now it's just easier to work on them for the specific purpose (exponential integrator) all together.