Merge pull request #36617 from marload:master
PiperOrigin-RevId: 294916821 Change-Id: I1ace34ff994058b84a714d1002d3d4637127effc
This commit is contained in:
commit
5bad62a486
@ -191,56 +191,44 @@ def _ndtri(p):
|
||||
# Constants used in piece-wise rational approximations. Taken from the cephes
|
||||
# library:
|
||||
# https://root.cern.ch/doc/v608/SpecFuncCephesInv_8cxx_source.html
|
||||
p0 = list(reversed([-5.99633501014107895267E1,
|
||||
9.80010754185999661536E1,
|
||||
-5.66762857469070293439E1,
|
||||
1.39312609387279679503E1,
|
||||
-1.23916583867381258016E0]))
|
||||
q0 = list(reversed([1.0,
|
||||
1.95448858338141759834E0,
|
||||
4.67627912898881538453E0,
|
||||
8.63602421390890590575E1,
|
||||
-2.25462687854119370527E2,
|
||||
2.00260212380060660359E2,
|
||||
-8.20372256168333339912E1,
|
||||
1.59056225126211695515E1,
|
||||
-1.18331621121330003142E0]))
|
||||
p1 = list(reversed([4.05544892305962419923E0,
|
||||
3.15251094599893866154E1,
|
||||
5.71628192246421288162E1,
|
||||
4.40805073893200834700E1,
|
||||
1.46849561928858024014E1,
|
||||
2.18663306850790267539E0,
|
||||
-1.40256079171354495875E-1,
|
||||
-3.50424626827848203418E-2,
|
||||
-8.57456785154685413611E-4]))
|
||||
q1 = list(reversed([1.0,
|
||||
1.57799883256466749731E1,
|
||||
4.53907635128879210584E1,
|
||||
4.13172038254672030440E1,
|
||||
1.50425385692907503408E1,
|
||||
2.50464946208309415979E0,
|
||||
-1.42182922854787788574E-1,
|
||||
-3.80806407691578277194E-2,
|
||||
-9.33259480895457427372E-4]))
|
||||
p2 = list(reversed([3.23774891776946035970E0,
|
||||
6.91522889068984211695E0,
|
||||
3.93881025292474443415E0,
|
||||
1.33303460815807542389E0,
|
||||
2.01485389549179081538E-1,
|
||||
1.23716634817820021358E-2,
|
||||
3.01581553508235416007E-4,
|
||||
2.65806974686737550832E-6,
|
||||
6.23974539184983293730E-9]))
|
||||
q2 = list(reversed([1.0,
|
||||
6.02427039364742014255E0,
|
||||
3.67983563856160859403E0,
|
||||
1.37702099489081330271E0,
|
||||
2.16236993594496635890E-1,
|
||||
1.34204006088543189037E-2,
|
||||
3.28014464682127739104E-4,
|
||||
2.89247864745380683936E-6,
|
||||
6.79019408009981274425E-9]))
|
||||
|
||||
p0 = [
|
||||
-1.23916583867381258016E0, 1.39312609387279679503E1,
|
||||
-5.66762857469070293439E1, 9.80010754185999661536E1,
|
||||
-5.99633501014107895267E1
|
||||
]
|
||||
q0 = [
|
||||
-1.18331621121330003142E0, 1.59056225126211695515E1,
|
||||
-8.20372256168333339912E1, 2.00260212380060660359E2,
|
||||
-2.25462687854119370527E2, 8.63602421390890590575E1,
|
||||
4.67627912898881538453E0, 1.95448858338141759834E0, 1.0
|
||||
]
|
||||
p1 = [
|
||||
-8.57456785154685413611E-4, -3.50424626827848203418E-2,
|
||||
-1.40256079171354495875E-1, 2.18663306850790267539E0,
|
||||
1.46849561928858024014E1, 4.40805073893200834700E1,
|
||||
5.71628192246421288162E1, 3.15251094599893866154E1,
|
||||
4.05544892305962419923E0
|
||||
]
|
||||
q1 = [
|
||||
-9.33259480895457427372E-4, -3.80806407691578277194E-2,
|
||||
-1.42182922854787788574E-1, 2.50464946208309415979E0,
|
||||
1.50425385692907503408E1, 4.13172038254672030440E1,
|
||||
4.53907635128879210584E1, 1.57799883256466749731E1, 1.0
|
||||
]
|
||||
p2 = [
|
||||
6.23974539184983293730E-9, 2.65806974686737550832E-6,
|
||||
3.01581553508235416007E-4, 1.23716634817820021358E-2,
|
||||
2.01485389549179081538E-1, 1.33303460815807542389E0,
|
||||
3.93881025292474443415E0, 6.91522889068984211695E0,
|
||||
3.23774891776946035970E0
|
||||
]
|
||||
q2 = [
|
||||
6.79019408009981274425E-9, 2.89247864745380683936E-6,
|
||||
3.28014464682127739104E-4, 1.34204006088543189037E-2,
|
||||
2.16236993594496635890E-1, 1.37702099489081330271E0,
|
||||
3.67983563856160859403E0, 6.02427039364742014255E0, 1.0
|
||||
]
|
||||
|
||||
def _create_polynomial(var, coeffs):
|
||||
"""Compute n_th order polynomial via Horner's method."""
|
||||
|
Loading…
Reference in New Issue
Block a user