From 752593a5893f27078570bc4dddd0ea8d31ec6f43 Mon Sep 17 00:00:00 2001 From: lisyarus Date: Mon, 5 Feb 2024 13:32:35 +0300 Subject: [PATCH] Fix random::normal_distribution: actually use mean & stddev instead of always generating N(0,1) --- libs/random/include/psemek/random/normal.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/libs/random/include/psemek/random/normal.hpp b/libs/random/include/psemek/random/normal.hpp index 65d8626a..cb761fb1 100644 --- a/libs/random/include/psemek/random/normal.hpp +++ b/libs/random/include/psemek/random/normal.hpp @@ -64,8 +64,8 @@ namespace psemek::random T const r = std::sqrt(- 2 * std::log(u1)); T const th = 2 * geom::pi * u2; - T const z1 = r * std::cos(th); - T const z2 = r * std::sin(th); + T const z1 = mean_ + stddev_ * r * std::cos(th); + T const z2 = mean_ + stddev_ * r * std::sin(th); cached_ = z2; return z1;