/* stepSystem: advance N-body system one time-step. */ NBodyStatus nbStepSystemPlain(const NBodyCtx* ctx, NBodyState* st) { NBodyStatus rc; const real dt = ctx->timestep; advancePosVel(st, st->nbody, dt); rc = nbGravMap(ctx, st); advanceVelocities(st, st->nbody, dt); st->step++; return rc; }
/* stepSystem: advance N-body system one time-step. */ NBodyStatus nbStepSystemPlain(const NBodyCtx* ctx, NBodyState* st) { NBodyStatus rc; const real dt = ctx->timestep; advancePosVel(st, st->nbody, dt); rc = nbGravMap(ctx, st); apply_dynamical_friction(st,ctx); advanceVelocities(st, st->nbody, dt); st->step++; #ifdef NBODY_BLENDER_OUTPUT blenderPrintBodies(st, ctx); // printf("Frame: %d (%f%%)\n", (int)(st->step),100.0*st->step/ctx->nStep); #endif return rc; }