This computation of C as following is not consistent with [this line](https://github.com/taichi-dev/taichi/blob/215d6d13c1733cb297c4ffebc84f896c5c8a5576/examples/mpm88.py#L68). https://github.com/yuanming-hu/taichi_mpm/blob/3bb90fbe4c901aafc048dbb2d8d8aa388226d011/mls-mpm88-explained.cpp#L171 The above line is lack of one `1/dx` according to the last page of [this slide](https://github.com/yuanming-hu/taichi_mpm/releases/download/SIGGRAPH2018/mls-mpm-cpic-slides.pdf). 
This computation of C as following is not consistent with this line.
taichi_mpm/mls-mpm88-explained.cpp
Line 171 in 3bb90fb
The above line is lack of one
1/dxaccording to the last page of this slide.