Fix normal_distribution bug

2021-01-23 21:33:45 +03:00 · 2021-01-23 21:33:45 +03:00 · 7ac4e3eb9b
commit 7ac4e3eb9b
parent 5e7aa7601d
1 changed files with 3 additions and 1 deletions
--- 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<T> 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));