关于mpm88行的stress，以及C矩阵的计算

[jlxy11]

大家好，我在蒋老师的mpmcourse中没有找到这两处代码所对应的公式，大家能为我解答一下吗？
stress = -dt * 4 * E * p_vol * (J[p] - 1) / dx**2
affine = ti.Matrix([[stress, 0], [0, stress]]) + p_mass * C[p]

new_C += 4 * weight * g_v.outer_product(dpos) / dx**2

[jlxy11]

为什么前面的系数是4呢？

[BillXu2000]

列一些可能相关的链接：

1. APIC https://www.math.ucla.edu/~cffjiang/research/apic/paper.pdf
2. MLS-MPM A Moving Least Squares Material Point Method with Displacement Discontinuity and Two-Way Rigid Body Coupling | Yuanming Hu
3. GAMES 201 lec7, 8 GAMES201：高级物理引擎实战指南2020_哔哩哔哩_bilibili

[jlxy11]

谢谢