|
12 | 12 |
|
13 | 13 | from functools import partial
|
14 | 14 |
|
15 |
| -import scipy.special |
16 |
| -import scipy.stats |
| 15 | +from scipy import special, stats |
17 | 16 |
|
18 | 17 | from pytensor import function, grad
|
19 | 18 | from pytensor import tensor as pt
|
|
45 | 44 |
|
46 | 45 |
|
47 | 46 | def scipy_special_gammau(k, x):
|
48 |
| - return scipy.special.gammaincc(k, x) * scipy.special.gamma(k) |
| 47 | + return special.gammaincc(k, x) * special.gamma(k) |
49 | 48 |
|
50 | 49 |
|
51 | 50 | def scipy_special_gammal(k, x):
|
52 |
| - return scipy.special.gammainc(k, x) * scipy.special.gamma(k) |
| 51 | + return special.gammainc(k, x) * special.gamma(k) |
53 | 52 |
|
54 | 53 |
|
55 | 54 | # Precomputing the result is brittle(it have been broken!)
|
56 | 55 | # As if we do any modification to random number here,
|
57 | 56 | # The input random number will change and the output!
|
58 |
| -expected_erf = scipy.special.erf |
59 |
| -expected_erfc = scipy.special.erfc |
60 |
| -expected_erfinv = scipy.special.erfinv |
61 |
| -expected_erfcinv = scipy.special.erfcinv |
62 |
| -expected_owenst = scipy.special.owens_t |
63 |
| -expected_gamma = scipy.special.gamma |
64 |
| -expected_gammaln = scipy.special.gammaln |
65 |
| -expected_psi = scipy.special.psi |
66 |
| -expected_tri_gamma = partial(scipy.special.polygamma, 1) |
67 |
| -expected_chi2sf = scipy.stats.chi2.sf |
68 |
| -expected_gammainc = scipy.special.gammainc |
69 |
| -expected_gammaincc = scipy.special.gammaincc |
| 57 | +expected_erf = special.erf |
| 58 | +expected_erfc = special.erfc |
| 59 | +expected_erfinv = special.erfinv |
| 60 | +expected_erfcinv = special.erfcinv |
| 61 | +expected_owenst = special.owens_t |
| 62 | +expected_gamma = special.gamma |
| 63 | +expected_gammaln = special.gammaln |
| 64 | +expected_psi = special.psi |
| 65 | +expected_tri_gamma = partial(special.polygamma, 1) |
| 66 | +expected_chi2sf = stats.chi2.sf |
| 67 | +expected_gammainc = special.gammainc |
| 68 | +expected_gammaincc = special.gammaincc |
70 | 69 | expected_gammau = scipy_special_gammau
|
71 | 70 | expected_gammal = scipy_special_gammal
|
72 |
| -expected_gammaincinv = scipy.special.gammaincinv |
73 |
| -expected_gammainccinv = scipy.special.gammainccinv |
74 |
| -expected_j0 = scipy.special.j0 |
75 |
| -expected_j1 = scipy.special.j1 |
76 |
| -expected_jv = scipy.special.jv |
77 |
| -expected_i0 = scipy.special.i0 |
78 |
| -expected_i1 = scipy.special.i1 |
79 |
| -expected_iv = scipy.special.iv |
80 |
| -expected_ive = scipy.special.ive |
81 |
| -expected_erfcx = scipy.special.erfcx |
82 |
| -expected_sigmoid = scipy.special.expit |
83 |
| -expected_hyp2f1 = scipy.special.hyp2f1 |
84 |
| -expected_betaincinv = scipy.special.betaincinv |
| 71 | +expected_gammaincinv = special.gammaincinv |
| 72 | +expected_gammainccinv = special.gammainccinv |
| 73 | +expected_j0 = special.j0 |
| 74 | +expected_j1 = special.j1 |
| 75 | +expected_jv = special.jv |
| 76 | +expected_i0 = special.i0 |
| 77 | +expected_i1 = special.i1 |
| 78 | +expected_iv = special.iv |
| 79 | +expected_ive = special.ive |
| 80 | +expected_erfcx = special.erfcx |
| 81 | +expected_sigmoid = special.expit |
| 82 | +expected_hyp2f1 = special.hyp2f1 |
| 83 | +expected_betaincinv = special.betaincinv |
85 | 84 |
|
86 | 85 | TestErfBroadcast = makeBroadcastTester(
|
87 | 86 | op=pt.erf,
|
@@ -837,14 +836,14 @@ def expected_log1mexp(x):
|
837 | 836 |
|
838 | 837 | TestBetaincBroadcast = makeBroadcastTester(
|
839 | 838 | op=pt.betainc,
|
840 |
| - expected=scipy.special.betainc, |
| 839 | + expected=special.betainc, |
841 | 840 | good=_good_broadcast_ternary_betainc,
|
842 | 841 | grad=_good_broadcast_ternary_betainc,
|
843 | 842 | )
|
844 | 843 |
|
845 | 844 | TestBetaincInplaceBroadcast = makeBroadcastTester(
|
846 | 845 | op=inplace.betainc_inplace,
|
847 |
| - expected=scipy.special.betainc, |
| 846 | + expected=special.betainc, |
848 | 847 | good=_good_broadcast_ternary_betainc,
|
849 | 848 | grad=_good_broadcast_ternary_betainc,
|
850 | 849 | inplace=True,
|
@@ -936,13 +935,13 @@ def test_beta_inc_stan_grad_combined(self):
|
936 | 935 |
|
937 | 936 | TestBetaincinvBroadcast = makeBroadcastTester(
|
938 | 937 | op=pt.betaincinv,
|
939 |
| - expected=scipy.special.betaincinv, |
| 938 | + expected=special.betaincinv, |
940 | 939 | good=_good_broadcast_ternary_betaincinv,
|
941 | 940 | )
|
942 | 941 |
|
943 | 942 | TestBetaincinvInplaceBroadcast = makeBroadcastTester(
|
944 | 943 | op=inplace.betaincinv_inplace,
|
945 |
| - expected=scipy.special.betaincinv, |
| 944 | + expected=special.betaincinv, |
946 | 945 | good=_good_broadcast_ternary_betaincinv,
|
947 | 946 | inplace=True,
|
948 | 947 | )
|
|
0 commit comments