diff --git a/libs/random/include/psemek/random/uniform_real.hpp b/libs/random/include/psemek/random/uniform_real.hpp
index 377197dc..531d7524 100644
--- a/libs/random/include/psemek/random/uniform_real.hpp
+++ b/libs/random/include/psemek/random/uniform_real.hpp
@@ -52,25 +52,12 @@ namespace psemek::random
 	template <typename RNG>
 	T uniform_real_distribution<T>::generate(RNG & rng)
 	{
-		// see https://en.cppreference.com/w/cpp/numeric/random/generate_canonical
-
-		static constexpr std::size_t digits = std::numeric_limits<T>::digits;
-		static T const limit = std::pow<T>(2, digits);
+		// Previously used an algorithm like in https://en.cppreference.com/w/cpp/numeric/random/generate_canonical
+		// But the simple algorithm here is way faster, so...
 
 		T const R = static_cast<T>(rng.max()) - static_cast<T>(rng.min()) + 1;
 
-		T mult = T(1);
-
-		T sum = T(0);
-
-		while (mult < limit)
-		{
-			T inc = static_cast<T>(rng() - rng.min());
-			sum += inc * mult;
-			mult *= R;
-		}
-
-		return sum / mult;
+		return static_cast<T>(rng() - rng.min()) / R;
 	}
 
 }