Jump to content

  • Log In with Google      Sign In   
  • Create Account


taby

Member Since 10 Feb 2005
Offline Last Active May 11 2014 11:18 AM
***--

Profile Feed

Photo
Photo

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
PARTNERS