代码
这里是我在过去的工作中曾经用到的代码。
-
求解形如
\[\partial_t f_i=\sum_{k=1}^N\nabla\cdot(h_{ik}\nabla f_k)+g_i(f_1,f_2,\cdots,f_N)\]的非线性反应-扩散方程。不过边界的发散问题尚未解决。
-
二维格点上的无相互作用run-and-tumble粒子，使用连续时间的Gillespie算法
另外还有高效地加入格点上的逻辑斯蒂群体动力学的版本，和使用CUDA的版本。
-
从对粒子拍摄的高速影片数据中提取粒子的运动信息，基于对模型的数值拉普拉斯逆变换。参见我的博士毕业论文。
-
-