我正在使用 FORTRAN 来实现Gillespie 的随机模拟算法。我将并行运行许多模拟（具有不同种子和并行函数的并行实例）；如果我必须存储每个时间间隔，那么会消耗大量内存。所以我想我可以有一些固定数量的时间点来保存数据（预先动态或静态分配空间）。

如何确定时间间隔的大小？间隔应该具有固定大小还是应该基于局部动态（某些区域会有更多动态）？有没有办法在不实际运行确定性模拟的情况下基于确定性函数进行预测？

我有一个模糊的想法，即 tau-leap 公式可用于预测n发生重大反应次数的时间间隔。但这仍然不能帮助我预先预测步数。