WIP: 1次精度ラプラシアンを実装（粘性項のみ）

[LWisteria]

#30 を考えていたら、HLよりこっちの方が精度が良いのではと言う気がしてきたので

[LWisteria]

中心重力や水中崩壊のベンチマークを実行中です。

[LWisteria]

中心重力問題をやりました。#53 の問題があるのでSPPを外しています。あと半径を20にしました。

結果は全然駄目でした

1LじゃなくてHLの結果と比べると、

• 途中で圧力が消えたり高くなったりします
• そもそも圧力値高すぎ
• 表面もあんまり落ち着きません

単体テストは通っているんですけどね・・・。実装バグがありそうなので見直します。

[LWisteria]

調べてみたら単体試験通ってませんでした。

SPP	ラプラシアン	単体試験結果	備考
あり	HL	合格	今のmaster
あり	1L	合格	このPR最新
なし	1L	不合格	上記計算 #41 (comment) の設定
なし	HL	合格

不合格なのはNeighborDensityTestですね

[----------] 2 tests from NeighborDensityTest
[ RUN      ] NeighborDensityTest.ValueTestLinear
test_ComputerNeighborDensityVariation.cpp(159): error: The difference between std::abs((dndt - dndt_analy) / dndt_analy) and 0.0 is 0.75877817465439223, which exceeds testAccuracy, where
std::abs((dndt - dndt_analy) / dndt_analy) evaluates to 0.75877817465439223,
0.0 evaluates to 0, and
testAccuracy evaluates to 0.10000000000000001.
[  FAILED  ] NeighborDensityTest.ValueTestLinear (195 ms)
[ RUN      ] NeighborDensityTest.ValueTestPolynomial
test_ComputerNeighborDensityVariation.cpp(211): error: The difference between std::abs((dndt - dndt_analy) / dndt_analy) and 0.0 is 0.77584798548503164, which exceeds testAccuracy, where
std::abs((dndt - dndt_analy) / dndt_analy) evaluates to 0.77584798548503164,
0.0 evaluates to 0, and
testAccuracy evaluates to 0.10000000000000001.
[  FAILED  ] NeighborDensityTest.ValueTestPolynomial (196 ms)
[----------] 2 tests from NeighborDensityTest (394 ms total)

[LWisteria]

単体試験が落ちるやつ、SPPを導入していないとそもそもmasterでも落ちますね・・・（↑の「このPR最新」にはNeighborDensityVariationの試験は入ってなかった）。別問題っぽい。