I do research in the field of mathematical ecology and use R for numerical computations and simulations. Now to deal with my new research problem, I have to first reproduce a numerical simulation from a published paper ( https://royalsocietypublishing.org/doi/pdf/10.1098/rstb.2015.0281). I`m having trouble doing this. I am telling the process briefly: The simulation starts with a system of 2 odes, which I have to integrate for 10000 time steps and need the steady state. After that, the number of variables in the system of ode may increase to 3 (depending on some conditions) and again I have to integrate it for 10000 time steps. After that, the number of variables may be 4 and so on. Each step of increasing the number of variables is called the mutation steps in the language of ecology.
Now I have to do 50000 mutation steps (each step contains ode integration for 10000 time points) to get a complete data set. I call the functions from C++ (using Rcpp), then integrate it using the RK4 method of numerical integration of ode but it is taking a huge amount of time (near about 50-60 hrs) to complete the process. Can one suggest to me what I should do to reduce the time of the simulation process.
_______________________________________________ Rcpp-devel mailing list Rcpp-devel@lists.r-forge.r-project.org https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel