Find content
Not Telling
Profile Feed
This status update has reached the maximum number of replies allowed
taby
inline void proceed_rk2(const double &dt) { vector_3 k1_velocity = mercury.velocity; vector_3 k1_acceleration = gravitation(mercury.position, k1_velocity); vector_3 k2_velocity = mercury.velocity + k1_acceleration*dt*0.5; vector_3 k2_acceleration = gravitation(mercury.position + k1_velocity*dt*0.5, k2_velocity); mercury.velocity += k2_acceleration*dt; mercury.position += k2_velocity*dt; } inline void proceed_RK4(const double &dt) { // Performance is increased on some archite...
Jun 14 2012 05:26 PM