From 89c6a75fe3e3d2d328f41dcb252cea28c33bdf2d Mon Sep 17 00:00:00 2001
From: Antonio Carzaniga <antonio.carzaniga@usi.ch>
Date: Fri, 5 Jul 2019 15:58:10 +0200
Subject: [PATCH] added init function for truly uniform random velocity.

---
 balls.c | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/balls.c b/balls.c
index 8b154d3..6f01c84 100644
--- a/balls.c
+++ b/balls.c
@@ -39,6 +39,15 @@ unsigned int n_balls = 50;
 static double g_y = 20;
 static double g_x = 0;
 
+void random_velocity(struct ball * p) {
+    double r2;
+    do {
+	p->v_x = v_min + rand() % (v_max + 1 - v_min);
+	p->v_y = v_min + rand() % (v_max + 1 - v_min);
+	r2 = p->v_x*p->v_x + p->v_y*p->v_y;
+    } while (r2 > v_max*v_max || r2 < v_min*v_min);
+}
+
 void balls_init_state () {
     srand(time(NULL));
     static const unsigned int border = 10;
@@ -48,8 +57,7 @@ void balls_init_state () {
     for (unsigned int i = 0; i < n_balls; ++i) {
 	balls[i].x = border + rand() % w;
 	balls[i].y = border + rand() % h;
-	balls[i].v_x = v_min + rand() % (v_max + 1 - v_min);
-	balls[i].v_y = v_min + rand() % (v_max + 1 - v_min);
+	random_velocity(balls + i);
 	if (rand() % 2)
 	    balls[i].v_x = -balls[i].v_x;
 	if (rand() % 2)