Results with seldom updates of matrix exponential about consim HOT 2 CLOSED

Here it is a more complete version of the plot above:

The take home message is always the same: we don't need to update the matrix exponential (roughly) more than once every ms because doing so doesn't improve accuracy (in the range of time steps tested). Therefore, if the integration accuracy you get with a time step of 1 ms is not enough, you can simply update more frequently the other terms of the dynamics (M, h, J), while keeping updating the matrix exponential only every 1 ms. That's gonna look very good in a plot accuracy VS computation time, where we would see a fast decrease of the error for a small increase of computation!

from consim.

I've added the possibility to set the frequency of the updates of the matrix exponential in the c++ implementation of the Exponential simulator. However, at the moment I haven't optimized the code to exploit this feature in terms of computation time because it would mainly benefit the high accuracy simulations, whereas I am currently more interested in the low-accuracy but high-speed simulations.

from consim.