piLaboratory / GillesCom

Gillespie Simulation of Ecological Communities Dynamics
2 stars 0 forks source link

Rcpp trajectories [ Depends on #21 ] #22

Closed andrechalom closed 8 years ago

andrechalom commented 8 years ago

Implemented the saving of time and cycles at each snapshot, this now is implemented as an Rcpp::List and returned as a list with three elements to R. This might impact all older code that referenced to x$trajectories, which should now read x$trajectories[[1]]. My initial performance testings suggest that this modification does not impact significantly the performance of the code - at least for reasonably large saving intervals.

The plotting functions have been redesigned to work with the new trajectories, and now use the exact time in their plots. This will probably have minimal effects, but can be seen when the community is very small or the saving interval is very short. Notice how the lines are not aligned with the "index" 0, 1, 2, 3, etc, but with the actual time 0, 0.22, 1.99, 2.38, etc.

newp