2024/08/10/DEM%E7%B2%97%E7%B2%92%E5%8C%96-CODE/

[giscus[bot]]

2024/08/10/DEM%E7%B2%97%E7%B2%92%E5%8C%96-CODE/

基于这篇的内容，我用 Python 写了一份对单种颗粒系统进行粗粒化的代码。 其具体使用流程为： obj = CG_mono_3D(...): 初始化。 obj.set_domain(...): 设置矩形的计算域。 obj.get_CG_data(R, CG_width): 对 R 点进行粗粒化，其粗粒化的长度为 CG_width。 DEM_CG_mono.py 1234567891011

https://ives0721.github.io/2024/08/10/DEM%E7%B2%97%E7%B2%92%E5%8C%96-CODE/

[Ives0721]

这里需要吐槽一下。在这份代码里，我的 self._get_distance_vec 函数算的是 $\boldsymbol{r}_{i} - \boldsymbol{r}$，和公式里要的东西（ $\boldsymbol{r} - \boldsymbol{r}_{i}$ ）是反过来的，导致代码看着有点怪😓

[Ives0721]

代码已经修改为正确的版本。