n = sqrt(x*x + y*y); x /= n; y /= n;
// velocities along (x, y), or 0 if already moving away.
n = sqrt(x*x + y*y); x /= n; y /= n;
// velocities along (x, y), or 0 if already moving away.
- da = max(x*a->dx + y*a->dy, 0);
- db = min(x*b->dx + y*b->dy, 0);
+ va = max(x*a->dx + y*a->dy, 0);
+ vb = min(x*b->dx + y*b->dy, 0);
- a->dx += x*(db*mr - da); a->dy += y*(db*mr - da);
- b->dx += x*(da/mr - db); b->dy += y*(da/mr - db);
+ a->dx += x*(mb*vb - ma*va)/ma; a->dy += y*(mb*vb - ma*va)/ma;
+ b->dx += x*(ma*va - mb*vb)/mb; b->dy += y*(ma*va - mb*vb)/mb;