なるほど…ひとまず倍精度にしてみてどれくらい遅くなるか
見てみることにします。情報を小出しにしちゃって申し訳ないんですが、
一番問題なのは、上の処理をすると、領域ぎりぎりに粒子があったとき、
５０までしかセルがないのに５１セル目を見ようとしちゃうことですね。
このエラー自体が不可避ならとりあえずそこの例外処理だけは入れといたほうがよさそうですね。
いろいろありがとうございます。上記件について何かいいアイデアあればぜひ。