西の方のがくせーくんが困ってる問題の原因が大分しぼれてきたので, りっつせんせーの timestep の limitter とPrice の人工熱伝導項を入れてみた. 後者については, がじぇっとのいろんな物理量のスケールファクター依存性が独特でちょっと苦労した.
で, 走らせてみてるわけだが, ワタシのみたいに強烈なフィードバックをいれてると limitter は結構結果に影響する. そして計算もすげー遅くなる. 何しろ一番 timestep が長い粒子の次に数が多いのが一番 timestep の短い粒子とかいう狂った hierarchy.
人工熱伝導はワタシのみたいな銀河形成のシミュレーションの結果に影響するのかどうか今のところよーわからん. まーでもぴよぴよと動いてはいる.