Problem overflow when computing the power of higher values

[salmuz]

Hello,

Thanks for your library. I catch a bug in calculating the power of higher values.

When we have higher mu, scale and x values, for instance:

x = [2881890.0, 1520335.0, 338717.0, 1842502.0, 325209.0, 282761.0, 847733.0, 1183272.0, 16533530.0]
mu = [1540092.433139984, 1157811.842, 775235.501, 660113.3310, 876908.8905, 1098942.967, 1275848.455, 1136826.756, 17490555.9386]
tweedie(mu=mu, p=1.0275417, phi=216098.00079).logpdf(x).sum()

We obtain an error: RuntimeWarning: overflow encountered in multiply tweedie
I tested that with Python 3.10

[thequackdaddy]

Thanks for the PR!

So I haven't been around recently ... sorry. There used to be a travis test suite that would run, but I know GitHub has been changing things around. I guess we should get the test suite working first (and made up to date) before merging. Once that is done, I'm happy to make some wheels and publish to pypi.