// NOTE(Brendan): This represents the piecewise function:
//
// F(d) = {
-// 0 < d < 1: (1 / d) - 1
+// 0 < d < 1: r * (1 - d) / d
// 1 <= d < distance_range + 1: (-4 * max_force / (distance_range ^ 2)) * (d - (distance_range / 2) - 1) ^ 2 + max_force,
// otherwise: 0
// }
persist const f64 step = 0.01;
foreach (i, low, high) body_accumulate_move(&state->bodies[i], &state->qt_bodies, TEMP_dt);
+ pthread_barrier_wait(&state->thread_sync_barrier);
foreach (i, low, high) body_apply_move(&state->bodies[i], TEMP_dt);
}