いまいち安定して計算できなかったMPSのプログラムだが、 粒子速度に合わせてCFL条件をきちんと守る可変タイムステップ導入（CFL Divergence-free形式を用いる圧力ソース項にする振動対策（＋粒子衝突対策など） 自由表面粒子の選択アルゴリズム（オリジナル…