scipy 1.15.3__cp313-cp313-macosx_12_0_arm64.whl → 1.16.0__cp313-cp313-macosx_12_0_arm64.whl

scipy/special/xsf/cephes/shichi.h

@@ -1,248 +0,0 @@
-/* Translated into C++ by SciPy developers in 2024.
- * Original header with Copyright information appears below.
- */
-
-/*                                                     shichi.c
- *
- *     Hyperbolic sine and cosine integrals
- *
- *
- *
- * SYNOPSIS:
- *
- * double x, Chi, Shi, shichi();
- *
- * shichi( x, &Chi, &Shi );
- *
- *
- * DESCRIPTION:
- *
- * Approximates the integrals
- *
- *                            x
- *                            -
- *                           | |   cosh t - 1
- *   Chi(x) = eul + ln x +   |    -----------  dt,
- *                         | |          t
- *                          -
- *                          0
- *
- *               x
- *               -
- *              | |  sinh t
- *   Shi(x) =   |    ------  dt
- *            | |       t
- *             -
- *             0
- *
- * where eul = 0.57721566490153286061 is Euler's constant.
- * The integrals are evaluated by power series for x < 8
- * and by Chebyshev expansions for x between 8 and 88.
- * For large x, both functions approach exp(x)/2x.
- * Arguments greater than 88 in magnitude return INFINITY.
- *
- *
- * ACCURACY:
- *
- * Test interval 0 to 88.
- *                      Relative error:
- * arithmetic   function  # trials      peak         rms
- *    IEEE         Shi      30000       6.9e-16     1.6e-16
- *        Absolute error, except relative when |Chi| > 1:
- *    IEEE         Chi      30000       8.4e-16     1.4e-16
- */
-
-/*
- * Cephes Math Library Release 2.0:  April, 1987
- * Copyright 1984, 1987 by Stephen L. Moshier
- * Direct inquiries to 30 Frost Street, Cambridge, MA 02140
- */
-#pragma once
-
-#include "../config.h"
-
-#include "chbevl.h"
-#include "const.h"
-
-namespace xsf {
-namespace cephes {
-
-    namespace detail {
-
-        /* x exp(-x) shi(x), inverted interval 8 to 18 */
-        constexpr double shichi_S1[] = {
-            1.83889230173399459482E-17,  -9.55485532279655569575E-17, 2.04326105980879882648E-16,
-            1.09896949074905343022E-15,  -1.31313534344092599234E-14, 5.93976226264314278932E-14,
-            -3.47197010497749154755E-14, -1.40059764613117131000E-12, 9.49044626224223543299E-12,
-            -1.61596181145435454033E-11, -1.77899784436430310321E-10, 1.35455469767246947469E-9,
-            -1.03257121792819495123E-9,  -3.56699611114982536845E-8,  1.44818877384267342057E-7,
-            7.82018215184051295296E-7,   -5.39919118403805073710E-6,  -3.12458202168959833422E-5,
-            8.90136741950727517826E-5,   2.02558474743846862168E-3,   2.96064440855633256972E-2,
-            1.11847751047257036625E0};
-
-        /* x exp(-x) shi(x), inverted interval 18 to 88 */
-        constexpr double shichi_S2[] = {
-            -1.05311574154850938805E-17, 2.62446095596355225821E-17,  8.82090135625368160657E-17,
-            -3.38459811878103047136E-16, -8.30608026366935789136E-16, 3.93397875437050071776E-15,
-            1.01765565969729044505E-14,  -4.21128170307640802703E-14, -1.60818204519802480035E-13,
-            3.34714954175994481761E-13,  2.72600352129153073807E-12,  1.66894954752839083608E-12,
-            -3.49278141024730899554E-11, -1.58580661666482709598E-10, -1.79289437183355633342E-10,
-            1.76281629144264523277E-9,   1.69050228879421288846E-8,   1.25391771228487041649E-7,
-            1.16229947068677338732E-6,   1.61038260117376323993E-5,   3.49810375601053973070E-4,
-            1.28478065259647610779E-2,   1.03665722588798326712E0};
-
-        /* x exp(-x) chin(x), inverted interval 8 to 18 */
-        constexpr double shichi_C1[] = {
-            -8.12435385225864036372E-18, 2.17586413290339214377E-17, 5.22624394924072204667E-17,
-            -9.48812110591690559363E-16, 5.35546311647465209166E-15, -1.21009970113732918701E-14,
-            -6.00865178553447437951E-14, 7.16339649156028587775E-13, -2.93496072607599856104E-12,
-            -1.40359438136491256904E-12, 8.76302288609054966081E-11, -4.40092476213282340617E-10,
-            -1.87992075640569295479E-10, 1.31458150989474594064E-8,  -4.75513930924765465590E-8,
-            -2.21775018801848880741E-7,  1.94635531373272490962E-6,  4.33505889257316408893E-6,
-            -6.13387001076494349496E-5,  -3.13085477492997465138E-4, 4.97164789823116062801E-4,
-            2.64347496031374526641E-2,   1.11446150876699213025E0};
-
-        /* x exp(-x) chin(x), inverted interval 18 to 88 */
-        constexpr double shichi_C2[] = {
-            8.06913408255155572081E-18,  -2.08074168180148170312E-17, -5.98111329658272336816E-17,
-            2.68533951085945765591E-16,  4.52313941698904694774E-16,  -3.10734917335299464535E-15,
-            -4.42823207332531972288E-15, 3.49639695410806959872E-14,  6.63406731718911586609E-14,
-            -3.71902448093119218395E-13, -1.27135418132338309016E-12, 2.74851141935315395333E-12,
-            2.33781843985453438400E-11,  2.71436006377612442764E-11,  -2.56600180000355990529E-10,
-            -1.61021375163803438552E-9,  -4.72543064876271773512E-9,  -3.00095178028681682282E-9,
-            7.79387474390914922337E-8,   1.06942765566401507066E-6,   1.59503164802313196374E-5,
-            3.49592575153777996871E-4,   1.28475387530065247392E-2,   1.03665693917934275131E0};
-
-        /*
-         * Evaluate 3F0(a1, a2, a3; z)
-         *
-         * The series is only asymptotic, so this requires z large enough.
-         */
-        XSF_HOST_DEVICE inline double hyp3f0(double a1, double a2, double a3, double z) {
-            int n, maxiter;
-            double err, sum, term, m;
-
-            m = std::pow(z, -1.0 / 3);
-            if (m < 50) {
-                maxiter = m;
-            } else {
-                maxiter = 50;
-            }
-
-            term = 1.0;
-            sum = term;
-            for (n = 0; n < maxiter; ++n) {
-                term *= (a1 + n) * (a2 + n) * (a3 + n) * z / (n + 1);
-                sum += term;
-                if (std::abs(term) < 1e-13 * std::abs(sum) || term == 0) {
-                    break;
-                }
-            }
-
-            err = std::abs(term);
-
-            if (err > 1e-13 * std::abs(sum)) {
-                return std::numeric_limits<double>::quiet_NaN();
-            }
-
-            return sum;
-        }
-
-    } // namespace detail
-
-    /* Sine and cosine integrals */
-    XSF_HOST_DEVICE inline int shichi(double x, double *si, double *ci) {
-        double k, z, c, s, a, b;
-        short sign;
-
-        if (x < 0.0) {
-            sign = -1;
-            x = -x;
-        } else {
-            sign = 0;
-        }
-
-        if (x == 0.0) {
-            *si = 0.0;
-            *ci = -std::numeric_limits<double>::infinity();
-            return (0);
-        }
-
-        if (x >= 8.0) {
-            goto chb;
-        }
-
-        if (x >= 88.0) {
-            goto asymp;
-        }
-
-        z = x * x;
-
-        /*     Direct power series expansion   */
-        a = 1.0;
-        s = 1.0;
-        c = 0.0;
-        k = 2.0;
-
-        do {
-            a *= z / k;
-            c += a / k;
-            k += 1.0;
-            a /= k;
-            s += a / k;
-            k += 1.0;
-        } while (std::abs(a / s) > detail::MACHEP);
-
-        s *= x;
-        goto done;
-
-    chb:
-        /* Chebyshev series expansions */
-        if (x < 18.0) {
-            a = (576.0 / x - 52.0) / 10.0;
-            k = std::exp(x) / x;
-            s = k * chbevl(a, detail::shichi_S1, 22);
-            c = k * chbevl(a, detail::shichi_C1, 23);
-            goto done;
-        }
-
-        if (x <= 88.0) {
-            a = (6336.0 / x - 212.0) / 70.0;
-            k = std::exp(x) / x;
-            s = k * chbevl(a, detail::shichi_S2, 23);
-            c = k * chbevl(a, detail::shichi_C2, 24);
-            goto done;
-        }
-
-    asymp:
-        if (x > 1000) {
-            *si = std::numeric_limits<double>::infinity();
-            *ci = std::numeric_limits<double>::infinity();
-        } else {
-            /* Asymptotic expansions
-             * http://functions.wolfram.com/GammaBetaErf/CoshIntegral/06/02/
-             * http://functions.wolfram.com/GammaBetaErf/SinhIntegral/06/02/0001/
-             */
-            a = detail::hyp3f0(0.5, 1, 1, 4.0 / (x * x));
-            b = detail::hyp3f0(1, 1, 1.5, 4.0 / (x * x));
-            *si = std::cosh(x) / x * a + std::sinh(x) / (x * x) * b;
-            *ci = std::sinh(x) / x * a + std::cosh(x) / (x * x) * b;
-        }
-        if (sign) {
-            *si = -*si;
-        }
-        return 0;
-
-    done:
-        if (sign) {
-            s = -s;
-        }
-
-        *si = s;
-
-        *ci = detail::SCIPY_EULER + std::log(x) + c;
-        return (0);
-    }
-
-} // namespace cephes
-} // namespace xsf

scipy/special/xsf/cephes/sici.h

@@ -1,224 +0,0 @@
-/* Translated into C++ by SciPy developers in 2024.
- * Original header with Copyright information appears below.
- */
-
-/*                                                     sici.c
- *
- *     Sine and cosine integrals
- *
- *
- *
- * SYNOPSIS:
- *
- * double x, Ci, Si, sici();
- *
- * sici( x, &Si, &Ci );
- *
- *
- * DESCRIPTION:
- *
- * Evaluates the integrals
- *
- *                          x
- *                          -
- *                         |  cos t - 1
- *   Ci(x) = eul + ln x +  |  --------- dt,
- *                         |      t
- *                        -
- *                         0
- *             x
- *             -
- *            |  sin t
- *   Si(x) =  |  ----- dt
- *            |    t
- *           -
- *            0
- *
- * where eul = 0.57721566490153286061 is Euler's constant.
- * The integrals are approximated by rational functions.
- * For x > 8 auxiliary functions f(x) and g(x) are employed
- * such that
- *
- * Ci(x) = f(x) sin(x) - g(x) cos(x)
- * Si(x) = pi/2 - f(x) cos(x) - g(x) sin(x)
- *
- *
- * ACCURACY:
- *    Test interval = [0,50].
- * Absolute error, except relative when > 1:
- * arithmetic   function   # trials      peak         rms
- *    IEEE        Si        30000       4.4e-16     7.3e-17
- *    IEEE        Ci        30000       6.9e-16     5.1e-17
- */
-
-/*
- * Cephes Math Library Release 2.1:  January, 1989
- * Copyright 1984, 1987, 1989 by Stephen L. Moshier
- * Direct inquiries to 30 Frost Street, Cambridge, MA 02140
- */
-#pragma once
-
-#include "../config.h"
-
-#include "const.h"
-#include "polevl.h"
-
-namespace xsf {
-namespace cephes {
-
-    namespace detail {
-
-        constexpr double sici_SN[] = {
-            -8.39167827910303881427E-11, 4.62591714427012837309E-8,  -9.75759303843632795789E-6,
-            9.76945438170435310816E-4,   -4.13470316229406538752E-2, 1.00000000000000000302E0,
-        };
-
-        constexpr double sici_SD[] = {
-            2.03269266195951942049E-12, 1.27997891179943299903E-9, 4.41827842801218905784E-7,
-            9.96412122043875552487E-5,  1.42085239326149893930E-2, 9.99999999999999996984E-1,
-        };
-
-        constexpr double sici_CN[] = {
-            2.02524002389102268789E-11, -1.35249504915790756375E-8, 3.59325051419993077021E-6,
-            -4.74007206873407909465E-4, 2.89159652607555242092E-2,  -1.00000000000000000080E0,
-        };
-
-        constexpr double sici_CD[] = {
-            4.07746040061880559506E-12, 3.06780997581887812692E-9, 1.23210355685883423679E-6,
-            3.17442024775032769882E-4,  5.10028056236446052392E-2, 4.00000000000000000080E0,
-        };
-
-        constexpr double sici_FN4[] = {
-            4.23612862892216586994E0,  5.45937717161812843388E0,  1.62083287701538329132E0,  1.67006611831323023771E-1,
-            6.81020132472518137426E-3, 1.08936580650328664411E-4, 5.48900223421373614008E-7,
-        };
-
-        constexpr double sici_FD4[] = {
-            /*  1.00000000000000000000E0, */
-            8.16496634205391016773E0,  7.30828822505564552187E0,  1.86792257950184183883E0,  1.78792052963149907262E-1,
-            7.01710668322789753610E-3, 1.10034357153915731354E-4, 5.48900252756255700982E-7,
-        };
-
-        constexpr double sici_FN8[] = {
-            4.55880873470465315206E-1, 7.13715274100146711374E-1,  1.60300158222319456320E-1,
-            1.16064229408124407915E-2, 3.49556442447859055605E-4,  4.86215430826454749482E-6,
-            3.20092790091004902806E-8, 9.41779576128512936592E-11, 9.70507110881952024631E-14,
-        };
-
-        constexpr double sici_FD8[] = {
-            /*  1.00000000000000000000E0, */
-            9.17463611873684053703E-1,  1.78685545332074536321E-1,  1.22253594771971293032E-2,
-            3.58696481881851580297E-4,  4.92435064317881464393E-6,  3.21956939101046018377E-8,
-            9.43720590350276732376E-11, 9.70507110881952025725E-14,
-        };
-
-        constexpr double sici_GN4[] = {
-            8.71001698973114191777E-2, 6.11379109952219284151E-1, 3.97180296392337498885E-1, 7.48527737628469092119E-2,
-            5.38868681462177273157E-3, 1.61999794598934024525E-4, 1.97963874140963632189E-6, 7.82579040744090311069E-9,
-        };
-
-        constexpr double sici_GD4[] = {
-            /*  1.00000000000000000000E0, */
-            1.64402202413355338886E0,  6.66296701268987968381E-1, 9.88771761277688796203E-2, 6.22396345441768420760E-3,
-            1.73221081474177119497E-4, 2.02659182086343991969E-6, 7.82579218933534490868E-9,
-        };
-
-        constexpr double sici_GN8[] = {
-            6.97359953443276214934E-1, 3.30410979305632063225E-1,  3.84878767649974295920E-2,
-            1.71718239052347903558E-3, 3.48941165502279436777E-5,  3.47131167084116673800E-7,
-            1.70404452782044526189E-9, 3.85945925430276600453E-12, 3.14040098946363334640E-15,
-        };
-
-        constexpr double sici_GD8[] = {
-            /*  1.00000000000000000000E0, */
-            1.68548898811011640017E0,  4.87852258695304967486E-1,  4.67913194259625806320E-2,
-            1.90284426674399523638E-3, 3.68475504442561108162E-5,  3.57043223443740838771E-7,
-            1.72693748966316146736E-9, 3.87830166023954706752E-12, 3.14040098946363335242E-15,
-        };
-
-    } // namespace detail
-
-    XSF_HOST_DEVICE inline int sici(double x, double *si, double *ci) {
-        double z, c, s, f, g;
-        short sign;
-
-        if (x < 0.0) {
-            sign = -1;
-            x = -x;
-        } else {
-            sign = 0;
-        }
-
-        if (x == 0.0) {
-            *si = 0.0;
-            *ci = -std::numeric_limits<double>::infinity();
-            return (0);
-        }
-
-        if (x > 1.0e9) {
-            if (std::isinf(x)) {
-                if (sign == -1) {
-                    *si = -M_PI_2;
-                    *ci = std::numeric_limits<double>::quiet_NaN();
-                } else {
-                    *si = M_PI_2;
-                    *ci = 0;
-                }
-                return 0;
-            }
-            *si = M_PI_2 - std::cos(x) / x;
-            *ci = std::sin(x) / x;
-        }
-
-        if (x > 4.0) {
-            goto asympt;
-        }
-
-        z = x * x;
-        s = x * polevl(z, detail::sici_SN, 5) / polevl(z, detail::sici_SD, 5);
-        c = z * polevl(z, detail::sici_CN, 5) / polevl(z, detail::sici_CD, 5);
-
-        if (sign) {
-            s = -s;
-        }
-        *si = s;
-        *ci = detail::SCIPY_EULER + std::log(x) + c; /* real part if x < 0 */
-        return (0);
-
-        /* The auxiliary functions are:
-         *
-         *
-         * *si = *si - M_PI_2;
-         * c = cos(x);
-         * s = sin(x);
-         *
-         * t = *ci * s - *si * c;
-         * a = *ci * c + *si * s;
-         *
-         * *si = t;
-         * *ci = -a;
-         */
-
-    asympt:
-
-        s = std::sin(x);
-        c = std::cos(x);
-        z = 1.0 / (x * x);
-        if (x < 8.0) {
-            f = polevl(z, detail::sici_FN4, 6) / (x * p1evl(z, detail::sici_FD4, 7));
-            g = z * polevl(z, detail::sici_GN4, 7) / p1evl(z, detail::sici_GD4, 7);
-        } else {
-            f = polevl(z, detail::sici_FN8, 8) / (x * p1evl(z, detail::sici_FD8, 8));
-            g = z * polevl(z, detail::sici_GN8, 8) / p1evl(z, detail::sici_GD8, 9);
-        }
-        *si = M_PI_2 - f * c - g * s;
-        if (sign) {
-            *si = -(*si);
-        }
-        *ci = f * s - g * c;
-
-        return (0);
-    }
-
-} // namespace cephes
-} // namespace xsf

scipy/special/xsf/cephes/sindg.h

@@ -1,221 +0,0 @@
-/* Translated into C++ by SciPy developers in 2024.
- * Original header with Copyright information appears below.
- */
-
-/*                                                     sindg.c
- *
- *     Circular sine of angle in degrees
- *
- *
- *
- * SYNOPSIS:
- *
- * double x, y, sindg();
- *
- * y = sindg( x );
- *
- *
- *
- * DESCRIPTION:
- *
- * Range reduction is into intervals of 45 degrees.
- *
- * Two polynomial approximating functions are employed.
- * Between 0 and pi/4 the sine is approximated by
- *      x  +  x**3 P(x**2).
- * Between pi/4 and pi/2 the cosine is represented as
- *      1  -  x**2 P(x**2).
- *
- *
- *
- * ACCURACY:
- *
- *                      Relative error:
- * arithmetic   domain      # trials      peak         rms
- *    IEEE      +-1000       30000      2.3e-16      5.6e-17
- *
- * ERROR MESSAGES:
- *
- *   message           condition        value returned
- * sindg total loss   x > 1.0e14 (IEEE)     0.0
- *
- */
-/*							cosdg.c
- *
- *	Circular cosine of angle in degrees
- *
- *
- *
- * SYNOPSIS:
- *
- * double x, y, cosdg();
- *
- * y = cosdg( x );
- *
- *
- *
- * DESCRIPTION:
- *
- * Range reduction is into intervals of 45 degrees.
- *
- * Two polynomial approximating functions are employed.
- * Between 0 and pi/4 the cosine is approximated by
- *      1  -  x**2 P(x**2).
- * Between pi/4 and pi/2 the sine is represented as
- *      x  +  x**3 P(x**2).
- *
- *
- * ACCURACY:
- *
- *                      Relative error:
- * arithmetic   domain      # trials      peak         rms
- *    IEEE     +-1000        30000       2.1e-16     5.7e-17
- *  See also sin().
- *
- */
-
-/* Cephes Math Library Release 2.0:  April, 1987
- * Copyright 1985, 1987 by Stephen L. Moshier
- * Direct inquiries to 30 Frost Street, Cambridge, MA 02140 */
-#pragma once
-
-#include "../config.h"
-#include "../error.h"
-
-#include "const.h"
-#include "polevl.h"
-
-namespace xsf {
-namespace cephes {
-
-    namespace detail {
-
-        constexpr double sincof[] = {1.58962301572218447952E-10, -2.50507477628503540135E-8,
-                                     2.75573136213856773549E-6,  -1.98412698295895384658E-4,
-                                     8.33333333332211858862E-3,  -1.66666666666666307295E-1};
-
-        constexpr double coscof[] = {1.13678171382044553091E-11, -2.08758833757683644217E-9, 2.75573155429816611547E-7,
-                                     -2.48015872936186303776E-5, 1.38888888888806666760E-3,  -4.16666666666666348141E-2,
-                                     4.99999999999999999798E-1};
-
-        constexpr double sindg_lossth = 1.0e14;
-
-    } // namespace detail
-
-    XSF_HOST_DEVICE inline double sindg(double x) {
-        double y, z, zz;
-        int j, sign;
-
-        /* make argument positive but save the sign */
-        sign = 1;
-        if (x < 0) {
-            x = -x;
-            sign = -1;
-        }
-
-        if (x > detail::sindg_lossth) {
-            set_error("sindg", SF_ERROR_NO_RESULT, NULL);
-            return (0.0);
-        }
-
-        y = std::floor(x / 45.0); /* integer part of x/M_PI_4 */
-
-        /* strip high bits of integer part to prevent integer overflow */
-        z = std::ldexp(y, -4);
-        z = std::floor(z);        /* integer part of y/8 */
-        z = y - std::ldexp(z, 4); /* y - 16 * (y/16) */
-
-        j = z; /* convert to integer for tests on the phase angle */
-        /* map zeros to origin */
-        if (j & 1) {
-            j += 1;
-            y += 1.0;
-        }
-        j = j & 07; /* octant modulo 360 degrees */
-        /* reflect in x axis */
-        if (j > 3) {
-            sign = -sign;
-            j -= 4;
-        }
-
-        z = x - y * 45.0;   /* x mod 45 degrees */
-        z *= detail::PI180; /* multiply by pi/180 to convert to radians */
-        zz = z * z;
-
-        if ((j == 1) || (j == 2)) {
-            y = 1.0 - zz * polevl(zz, detail::coscof, 6);
-        } else {
-            y = z + z * (zz * polevl(zz, detail::sincof, 5));
-        }
-
-        if (sign < 0)
-            y = -y;
-
-        return (y);
-    }
-
-    XSF_HOST_DEVICE inline double cosdg(double x) {
-        double y, z, zz;
-        int j, sign;
-
-        /* make argument positive */
-        sign = 1;
-        if (x < 0)
-            x = -x;
-
-        if (x > detail::sindg_lossth) {
-            set_error("cosdg", SF_ERROR_NO_RESULT, NULL);
-            return (0.0);
-        }
-
-        y = std::floor(x / 45.0);
-        z = std::ldexp(y, -4);
-        z = std::floor(z);        /* integer part of y/8 */
-        z = y - std::ldexp(z, 4); /* y - 16 * (y/16) */
-
-        /* integer and fractional part modulo one octant */
-        j = z;
-        if (j & 1) { /* map zeros to origin */
-            j += 1;
-            y += 1.0;
-        }
-        j = j & 07;
-        if (j > 3) {
-            j -= 4;
-            sign = -sign;
-        }
-
-        if (j > 1)
-            sign = -sign;
-
-        z = x - y * 45.0;   /* x mod 45 degrees */
-        z *= detail::PI180; /* multiply by pi/180 to convert to radians */
-
-        zz = z * z;
-
-        if ((j == 1) || (j == 2)) {
-            y = z + z * (zz * polevl(zz, detail::sincof, 5));
-        } else {
-            y = 1.0 - zz * polevl(zz, detail::coscof, 6);
-        }
-
-        if (sign < 0)
-            y = -y;
-
-        return (y);
-    }
-
-    /* Degrees, minutes, seconds to radians: */
-
-    /* 1 arc second, in radians = 4.848136811095359935899141023579479759563533023727e-6 */
-
-    namespace detail {
-        constexpr double sindg_P64800 = 4.848136811095359935899141023579479759563533023727e-6;
-    }
-
-    XSF_HOST_DEVICE inline double radian(double d, double m, double s) {
-        return (((d * 60.0 + m) * 60.0 + s) * detail::sindg_P64800);
-    }
-
-} // namespace cephes
-} // namespace xsf