C++ 物理シミュレーション・テクニック - Qiita

for(int i = 0 ; i < Nbody ; ++ i){ for(int j = 0 ; j < Nbody ; ++ j){ vec3<double> dr_ij = r[i] - r[j]; double inv_r_ij = 1.0 / abs(dr_ij); a[i] += - m[j] * dr_ij * (inv_r_ij * inv_r_ij * inv_r_ij); } } なんと！たったの三行！ そう我々はバグを出さないためにこういう書き方をすべきなのだ。 各々次元にはどうせ同じ処理をしてしまうのだ、圧縮してしまえば良い。 そこをfor文で書くと、バグを生む可能性が出るし、for文の分計算コストがかかるかもしれない。 (もっとも、これはコンパイラの最適化で何とかなるかもしれないわけだが…) さて、では何故こんな事ができるか？ もちろんC++にはvec3とかいう

[kasumani]

C++ 物理シミュレーション・テクニック 物理演算をコンピュータにやらせようとなると、ベクトルの計算は頻出である。 たとえば天文学とか惑星科学の分野とかでは、よく出てくるのが重力である。 Tags: ifttt, kasumaniのスト