09/11/29 13:14:09
みなさん、レスありがとうございます。まだ全部追えて無いんですが…
>>730
実は無次元化してません…
>>731
733が仰ってますが、単に変位です。負値は大丈夫です。
>>732
もちろんそうなんですが、どこで間違えてるのか分からなくて。
>>734
計算の一行前にwrite文で表示させてるので一回前の値と言うことはないはずなんですが、
この値で止まるのも同程度におかしいので何か間違えてるんだと思います。
>>735
ver.11です。デフォだとNaNで続行されますけど、NaNだと粒子が行方不明になってしまって困るので
トラップして原因を突き止めようとしてるのです。開放系ならごく稀に粒子が居なくなるくらいは問題ないんですが、
閉鎖系だと粒子がどんどん減ってしまうので…
>>737
こういう問題が発生しないアルゴリズムってどういうものでしょうか？原理的に発生しないんですか？
今自分で思いつくところだとNaNを検出して適当な座標を振りなおしてやるくらいしか思いつかないんですが…
ちなみにゼロ割した場合、floating invalidじゃなくてfloating divide by zeroとでました。分子も0だと
floating invalidになるのでこっちなんでしょうか？まだ不可解なところがあるのでもう少し見直してみます。