Prevent SciPy error by using float64 point in test_dirichlet_with_batch_shapes

Author: brandonwillard

pymc3/tests/test_distributions.py

@@ -1920,8 +1920,11 @@ def test_dirichlet_with_batch_shapes(self, dist_shape):
         with pm.Model() as model:
             d = pm.Dirichlet("d", a=a)
 
+        # Generate sample points to test
         d_value = d.tag.value_var
-        d_point = d.eval()
+        d_point = d.eval().astype("float64")
+        d_point /= d_point.sum(axis=-1)[..., None]
+
         if hasattr(d_value.tag, "transform"):
             d_point_trans = d_value.tag.transform.forward(d, aet.as_tensor(d_point)).eval()
         else: