Fix test_weibull_logp flakyness

michaelosthege · michaelosthege · commit b213cf3278b9 · 2023-02-03T14:29:22.000+01:00
diff --git a/pymc/tests/distributions/test_continuous.py b/pymc/tests/distributions/test_continuous.py
@@ -618,10 +618,13 @@ def test_pareto(self):
         reason="Fails on float32 due to numerical issues",
     )
     def test_weibull_logp(self):
+        # SciPy has new (?) precision issues at {alpha=20, beta=2, x=100}
+        # We circumvent it by skipping alpha=20:
+        rplusbig = Domain([0, 0.5, 0.9, 0.99, 1, 1.5, 2, np.inf])
         check_logp(
             pm.Weibull,
             Rplus,
-            {"alpha": Rplusbig, "beta": Rplusbig},
+            {"alpha": rplusbig, "beta": Rplusbig},
             lambda value, alpha, beta: st.exponweib.logpdf(value, 1, alpha, scale=beta),
         )
 

Original file line number	Diff line number	Diff line change
`@@ -618,10 +618,13 @@ def test_pareto(self):`
`618`	`618`	`reason="Fails on float32 due to numerical issues",`
`619`	`619`	`)`
`620`	`620`	`def test_weibull_logp(self):`
	`621`	`+ # SciPy has new (?) precision issues at {alpha=20, beta=2, x=100}`
	`622`	`+ # We circumvent it by skipping alpha=20:`
	`623`	`+ rplusbig = Domain([0, 0.5, 0.9, 0.99, 1, 1.5, 2, np.inf])`
`621`	`624`	`check_logp(`
`622`	`625`	`pm.Weibull,`
`623`	`626`	`Rplus,`
`624`		`- {"alpha": Rplusbig, "beta": Rplusbig},`
	`627`	`+ {"alpha": rplusbig, "beta": Rplusbig},`
`625`	`628`	`lambda value, alpha, beta: st.exponweib.logpdf(value, 1, alpha, scale=beta),`
`626`	`629`	`)`
`627`	`630`