pdf nan outside domain

Author: dweindl

petab/v1/distributions.py

@@ -177,10 +177,15 @@ def _pdf_transformed_untruncated(self, x) -> np.ndarray | float:
             )
 
             return np.where(
-                x > 0,
-                self._pdf_untransformed_untruncated(self._log(x))
-                * chain_rule_factor,
-                0,
+                x >= 0,
+                np.where(
+                    x > 0,
+                    self._pdf_untransformed_untruncated(self._log(x))
+                    * chain_rule_factor,
+                    0,
+                ),
+                # NaN outside its domain
+                np.nan,
             )
 
     @property

tests/v1/test_priors.py

@@ -272,7 +272,7 @@ def cdf(x):
     # check that the integral of the PDF is 1 for the unscaled parameters
     integral, abserr = quad(
         lambda x: prior.pdf(x, x_scaled=False),
-        -np.inf,
+        -np.inf if prior.distribution.logbase is False else 0,
         np.inf,
         limit=100,
         epsabs=1e-10,