Use RcppEigen matrices (no big effect). Derived variables calculated after integration of model using parallel programming. Some of the elements in the struct are const. The stats are stored in fixed-length arrays. Use pragma GCC ivdep (not much effect) in most for loops