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<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));