From 1fed2813e5f20e92f16f79df64209e5c50c32705 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Adrien=20Gallou=C3=ABt?= <adrien@gallouet.fr>
Date: Thu, 3 Oct 2019 15:10:37 +0000
Subject: [PATCH] Fix the unprobable s.v==1 case
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: Adrien Gallouët <adrien@gallouet.fr>
---
 src/bench.c | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/bench.c b/src/bench.c
index a12312b..394ef7b 100644
--- a/src/bench.c
+++ b/src/bench.c
@@ -91,13 +91,13 @@ gt_bench(int argc, char **argv)
                 s.m = mbps;
                 s.d = 0;
                 continue;
-            } else {
-                int64_t d1 = mbps - s.m; s.m += d1 / s.n;
-                int64_t d2 = mbps - s.m; s.d += d1 * d2;
             }
 
+            int64_t d1 = mbps - s.m; s.m += d1 / s.n;
+            int64_t d2 = mbps - s.m; s.d += d1 * d2;
+
             s.v = s.d / (s.n - 1);
-            s.sigma = s.v / 2;
+            s.sigma = 1 + (s.v - 1) / 2;
 
             while (s.sigma && s.sigma * s.sigma > s.v)
                 s.sigma--;
@@ -114,7 +114,7 @@ gt_bench(int argc, char **argv)
             printf("bench %"PRIi64" %"PRIi64" %"PRIi64"\n", size, s.m, s.sigma);
         }
 
-        size += 2*11*13;
+        size += 2 * 11 * 13;
     }
 
     return 0;