From ffaa77af6be3294d7963eab06a58d56a89a64a25 Mon Sep 17 00:00:00 2001 From: Axel Kohlmeyer Date: Tue, 19 May 2020 09:02:17 -0400 Subject: [PATCH] reintroduce bugfix for bessel function distributed random numbers --- src/random_mars.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/random_mars.cpp b/src/random_mars.cpp index 66c12f8ff5..64d6808e43 100644 --- a/src/random_mars.cpp +++ b/src/random_mars.cpp @@ -151,7 +151,7 @@ double RanMars::besselexp(double theta, double alpha, double cp) { double first,v1,v2; - if (theta < 0.0 || alpha < 0.0 || alpha < 1.0) + if (theta < 0.0 || alpha < 0.0 || alpha > 1.0) error->all(FLERR,"Invalid Bessel exponential distribution parameters"); v1 = uniform();