§ `ƒi0ãó—dZddlmZejdddd¬¦«Zejdd d d ¬¦«Zejd d dd¬¦«Zejdddd¬¦«Zejdddd¬¦«Zejdddd¬¦«Z dZ ejddde d ¬!¦«Z d"Z ejd#d$d%e d&¬!¦«Z ejd'd(d)e d*¬!¦«Ze d+zZejd,dd-ed.¬!¦«Zejd/dd0ed1¬!¦«Zd2Zejd3d4d5ed6¬!¦«Zd7S)8aIbasic special functions cotdg and tandg implementations are adapted from the following SciPy code: https://github.com/scipy/scipy/blob/master/scipy/special/cephes/tandg.c radian is from https://github.com/scipy/scipy/blob/master/scipy/special/cephes/sindg.c cosm1 is from https://github.com/scipy/scipy/blob/main/scipy/special/cephes/unity.c polevl is from https://github.com/scipy/scipy/blob/main/scipy/special/cephes/polevl.h Cephes Math Library Release 2.0: April, 1987 Copyright 1984, 1987 by Stephen L. Moshier Direct inquiries to 30 Frost Street, Cambridge, MA 02140 é)Ú_coreÚcupyx_scipy_log1p))úf->fzout0 = log1pf(in0))úF->Fz-out0 = out0_type(log1p((complex)in0))úd->dúD->Dzout0 = log1p(in0);zöElementwise function for scipy.special.log1p Calculates log(1 + x) for use when `x` is near zero. Notes ----- This implementation currently does not support complex-valued `x`. .. seealso:: :meth:`scipy.special.log1p` )ÚdocÚcupyx_scipy_special_cbrt))rzout0 = cbrtf(in0)rzout0 = cbrt(in0)z=Cube root. .. seealso:: :meth:`scipy.special.cbrt` Úcupyx_scipy_special_exp2))rzout0 = exp2f(in0)rzout0 = exp2(in0)zEComputes ``2**x``. .. seealso:: :meth:`scipy.special.exp2` Úcupyx_scipy_special_exp10))rzout0 = exp10f(in0)rzout0 = exp10(in0)zGComputes ``10**x``. .. seealso:: :meth:`scipy.special.exp10` Úcupyx_scipy_special_expm1))rzout0 = expm1f(in0)rrrzout0 = expm1(in0)zLComputes ``exp(x) - 1``. .. seealso:: :meth:`scipy.special.expm1` Úcupyx_scipy_special_exprel)rrrrz/out0 = abs(in0) >= 1e-16 ? expm1(in0) / in0 : 1zSComputes ``(exp(x) - 1) / x``. .. seealso:: :meth:`scipy.special.exprel` aŽ //Define from npy_math.h //https://github.com/numpy/numpy/blob/main/numpy/core/include/numpy/npy_math.h #define NPY_PI_4 0.785398163397448309615660845819875721 /* pi/4 */ __constant__ double coscof[] = { 4.7377507964246204691685E-14, -1.1470284843425359765671E-11, 2.0876754287081521758361E-9, -2.7557319214999787979814E-7, 2.4801587301570552304991E-5, -1.3888888888888872993737E-3, 4.1666666666666666609054E-2, }; __device__ static double polevl(double x, const double coef[], int N) { double ans; int i; const double *p; p = coef; ans = *p++; i = N; do ans = ans * x + *p++; while (--i); return (ans); } __device__ static double cosm1(double x) { double xx; if ((x < -NPY_PI_4) || (x > NPY_PI_4)) return (cos(x) - 1.0); xx = x * x; xx = -0.5 * xx + xx * xx * polevl(xx, coscof, 6); return xx; } Úcupyx_scipy_special_cosm1)rrzout0 = cosm1(in0)zLComputes ``cos(x) - 1``. .. seealso:: :meth:`scipy.special.cosm1` )Úpreambler zG __constant__ double PI180 = 1.74532925199432957692E-2; // pi/180 Úcupyx_scipy_special_cosdg))rzout0 = cosf(PI180 * in0)rzout0 = cos(PI180 * in0)zRCosine of x with x in degrees. .. seealso:: :meth:`scipy.special.cosdg` Úcupyx_scipy_special_sindg))rzout0 = sinf(PI180 * in0)rzout0 = sin(PI180 * in0)zPSine of x with x in degrees. .. seealso:: :meth:`scipy.special.sindg` aü // include for CUDART_INF #include __constant__ double lossth = 1.0e14; __device__ static double tancot(double, int); __device__ static double tandg(double x) { return tancot(x, 0); } __device__ static double cotdg(double x) { return tancot(x, 1); } __device__ static double tancot(double xx, int cotflg) { double x; int sign; /* make argument positive but save the sign */ if (xx < 0) { x = -xx; sign = -1; } else { x = xx; sign = 1; } if (x > lossth) { // sf_error("tandg", SF_ERROR_NO_RESULT, NULL); return 0.0; } /* modulo 180 */ x = x - 180.0 * floor(x / 180.0); if (cotflg) { if (x <= 90.0) { x = 90.0 - x; } else { x = x - 90.0; sign *= -1; } } else { if (x > 90.0) { x = 180.0 - x; sign *= -1; } } if (x == 0.0) { return 0.0; } else if (x == 45.0) { return sign * 1.0; } else if (x == 90.0) { // sf_error((cotflg ? "cotdg" : "tandg"), SF_ERROR_SINGULAR, NULL); return CUDART_INF; } /* x is now transformed into [0, 90) */ return sign * tan(x * PI180); } Úcupyx_scipy_special_tandgzout0 = tandg(in0)zSTangent of x with x in degrees. .. seealso:: :meth:`scipy.special.tandg` Úcupyx_scipy_special_cotdgzout0 = cotdg(in0)zUCotangent of x with x in degrees. .. seealso:: :meth:`scipy.special.cotdg` zè /* 1 arc second, in radians*/ __constant__ double P64800 = 4.848136811095359935899141023579479759563533023727e-6; template __device__ T radian(T d, T m, T s) { return (((d * 60.0 + m) * 60.0 + s) * P64800); } Úcupyx_scipy_special_radian)zfff->fzddd->dzout0 = radian(in0, in1, in2)zZDegrees, minutes, seconds to radians: .. seealso:: :meth:`scipy.special.radian` N)Ú__doc__ÚcupyrÚ create_ufuncÚlog1pÚcbrtÚexp2Úexp10Úexpm1ÚexprelÚcosm1_implementationÚcosm1Úpi180_preambleÚcosdgÚsindgÚtancot_implementationÚtandgÚcotdgÚradian_implementationÚradian©óún/home/jaya/work/projects/VOICE-AGENT/VIET/agent-env/lib/python3.11/site-packages/cupyx/scipy/special/_basic.pyúr,sVðððð2ÐÐÐÐÐð ˆÔØððð ð  ñ ô €ð(€uÔØØ+Øð ð  ñ ô €ð€uÔØØ+Øð ð  ñ ô €ð ˆÔØØ,Øð ð  ñ ô €ð ˆÔØØ<Øð ð  ñ ô €ð ˆÔ Ø Ø&Ø5ð ð  ñ ô €ð)ÐðV ˆÔØÐ!1ØØ !ð ð  ñ ô €ð€ð ˆÔØØ2ØØ ð ð ñ ô €ð ˆÔØØ2ØØ ð ð ñ ô €ð'ðF*ñFÐðP ˆÔØÐ!1ØØ "ð ð  ñ ô €ð ˆÔØÐ!1ØØ "ð ð  ñ ô €ð Ðð ˆÔ Ø Ð"6Ø"Ø "ð ð  ñ ô €€€r*