对于公式：
C_p^{n+1} = \frac 4 {{\Delta x}^2} \sum_i w_{ip} \mathbf{v}_i^{n+1}{\left(x_i - x_p^n\right)}^T
第一个问题：4/dx^2
GAMES201似乎并没有讲这个 \frac 4 {{\Delta x}^2} 如何推出来的，只说了MPM中这一项是从APIC里用Quadratic得来的，但是在上一讲APIC时也没有推导这一项
在APIC论文原文中，C矩阵是两个矩阵计算得出的：
C_p^n=B_p^n(D_p^n)^{-1} \\
D_p^n=\sum_i w_{ip}^n{\left(x_i - x_p^n\right)} {\left(x_i - x_p^n\right)}^T \\
B_p^{n+1}=\sum_i w_{ip}^n \tilde{\mathbf{v}}_i^{n+1} {\left(x_i - x_p^n\right)}^T
论文在后面又说在Quadratic差值形式下很简单可以化为 D_p^n = \frac 1 4 \Delta x^2 I ， 从而得出GAMES201中计算 C_p^{n+1} 的形式。但是论文原文也没有说如何推出这个 \frac 1 4 \Delta x^2 I
我试了两整天也没办法把这个 \frac 1 4 \Delta x^2 I 推出来
问题2：张量积
APIC论文原文中提到使用 v(x) = v_p + C_p(x - x_p) 来近似估计粒子p邻域的速度场，并且GAMES201中也使用了一张直观的图片来解释 C_p 矩阵。
照这样看 C_p 矩阵应该类似于速度关于位置的导数关系， 在应用中应该使用解线性方程 v(x) = v_p + C_p(x - x_p) 的方式来求解 C_p
但是实际上，应用中使用了张量积来近似这个导数关系，并且通过除以 \Delta x^2 来解决了量纲对不上的问题。这非常的奇怪
以上两个问题，希望有高人解答