From 7ac4e3eb9bb3a9827a33e5cb3fd7bd95c723367c Mon Sep 17 00:00:00 2001 From: lisyarus Date: Sat, 23 Jan 2021 21:33:45 +0300 Subject: [PATCH] Fix normal_distribution bug --- libs/random/include/psemek/random/normal.hpp | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/libs/random/include/psemek/random/normal.hpp b/libs/random/include/psemek/random/normal.hpp index 54b94bb7..65d8626a 100644 --- a/libs/random/include/psemek/random/normal.hpp +++ b/libs/random/include/psemek/random/normal.hpp @@ -56,7 +56,9 @@ namespace psemek::random } uniform_real_distribution d; - T const u1 = d(rng); + T u1 = 0.f; + while (u1 == 0.f) + u1 = d(rng); T const u2 = d(rng); T const r = std::sqrt(- 2 * std::log(u1));