marcusps / ExpmV.jl

Julia package to compute the result of expm(t*A)*v when A is a sparse matrix, without computing expm(t*A).
Other
22 stars 12 forks source link

Avoid globals in benchmark #3

Closed acroy closed 9 years ago

acroy commented 9 years ago

Just noticed that you use global variables in the benchmark. Currently, Julia isn't very good optimizing code which involves globals (basically since their type can change any time). Maybe you can wrap everything in a function?

marcusps commented 9 years ago

I fixed it and will push it soon. It does help Expokit quite a bit, and it makes me want to properly implement normAm ...