Add WeibullRV

aesara/tensor/random/basic.py

@@ -179,6 +179,17 @@ class ExponentialRV(RandomVariable):
 exponential = ExponentialRV()
 
 
+class WeibullRV(RandomVariable):
+    name = "weibull"
+    ndim_supp = 0
+    ndims_params = [0]
+    dtype = "floatX"
+    _print_name = ("Weibull", "\\operatorname{Weibull}")
+
+
+weibull = WeibullRV()
+
+
 def safe_multivariate_normal(mean, cov, size=None, rng=None):
     """A shape consistent multivariate normal sampler.
 

tests/tensor/random/test_basic.py

@@ -44,6 +44,7 @@ from aesara.tensor.random.basic import (
     truncexpon,
     uniform,
     wald,
+    weibull,
 )
 from aesara.tensor.type import iscalar, scalar, tensor
 
@@ -259,6 +260,13 @@ def test_exponential_samples():
     rv_numpy_tester(exponential, test_lambda, size=[2, 3])
 
 
+def test_weibull_samples():
+    test_alpha = np.array(3, dtype=config.floatX)
+
+    rv_numpy_tester(weibull, test_alpha)
+    rv_numpy_tester(weibull, test_alpha, size=[2, 3])
+
+
 def test_pareto_samples():
     test_alpha = np.array(0.5, dtype=config.floatX)