Better explanation

Author: awf

docs/source/05-stochastic-rounding.ipynb

@@ -91,9 +91,10 @@ def round_ndarray(
             ## RTNE delta to srbits
             d = delta * 2.0**srnumbits
             floord = np.floor(d).astype(np.int64)
-            d = floord + ((d - floord > 0.5) | ((d - floord == 0.5) & _isodd(floord)))
+            dd = d - floord
+            drnd = floord + (dd > 0.5) + ((dd == 0.5) & ~_isodd(floord))
 
-            should_round_away = d > srbits
+            should_round_away = drnd > srbits
         case RoundMode.StochasticFast:
             assert srbits is not None
             should_round_away = delta > (2 * srbits + 1) * 2.0 ** -(1 + srnumbits)