improve numerical stability of inverse sigmoid

Author: LarsKue

bayesflow/utils/numpy_utils.py

@@ -4,7 +4,7 @@
 
 def inverse_sigmoid(x: np.ndarray) -> np.ndarray:
     """Inverse of the sigmoid function."""
-    return np.log(x / (1 - x))
+    return np.log(x) - np.log1p(-x)
 
 
 def inverse_shifted_softplus(