Pivot on teh scaling function design

src/System/Random/Internal.hs

@@ -1469,18 +1469,14 @@ scaleFloating ::
   -- ^ Value in [0, 1] range that needs to be scaled.
   -> a
 scaleFloating l h x =
-  if abs l < abs h
-    -- Without clamping value furthest from zero we can get into a situation where scaled value can
-    -- get outside of it. This is not a problem for the value closest to zero, due to the
-    -- conditional above.
-    then let !y = l + x * (h - l)
-          in if l < h -- `l` is closer to zero
-               then min y h -- `l` is the true low, ensure `y` is not higher than `h`
-               else max y h -- `l` is the high, ensure `y` is not lower than `h`
-    else let !y = h + x * (l - h)
-          in if l < h -- `h` is closer to zero
-               then max y l -- `h` is the true high, ensure `y` is not lower than `l`
-               else min y l -- `h` is the low, ensure `y` is not higher than `l`
+  if isInfinite diff
+    then let !y = x * l + (1 - x) * h
+          in max (min y (max l h)) (min l h)
+    else if testBit w 31
+           then l + diff * x
+           else h + negate diff * x
+  where
+    !diff = h - l
 {-# INLINE scaleFloating #-}
 
 -- | Generates uniformly distributed 'Float' in the range \([0, 1]\).