From 77684898dd8d06d5ad58a6a2275e65521d0730e9 Mon Sep 17 00:00:00 2001 From: Peter Teuben Date: Wed, 25 Dec 2024 19:42:07 -0500 Subject: [PATCH] fix corner case --- src/kernel/misc/xrandom.c | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/kernel/misc/xrandom.c b/src/kernel/misc/xrandom.c index 3e69aa3e3..0a1ecaef7 100644 --- a/src/kernel/misc/xrandom.c +++ b/src/kernel/misc/xrandom.c @@ -190,6 +190,8 @@ double xrandom(double xl, double xh) { double retval; + if (xl==xh) return xl; + for(;;) { #if defined(HAVE_GSL) if (my_r == NULL) error("GSL init_xrandom was never called"); @@ -237,6 +239,8 @@ double grandom(double mean, double sdev) static double v1, v2, s; static int gcount = 0; + if (sdev <= 0) return mean; + if (gcount) { gcount = 0; return mean + sdev * v2 * s;