§ `ƒiïãóJ—ddlmZdZdZejdddeezd¬¦«Zd S) é)Ú_corezP #include __constant__ double EUL = 0.5772156649015328; a template static __device__ T exp1(T x) { if (x == 0) { return CUDART_INF; } else if (x <= 1) { T e1 = 1.0; T R = 1.0; for (int k = 1; k <= 25; k++) { int den = (k + 1) * (k + 1); R = -R * k * x / den; e1 += R; } return -EUL - log(x) + x * e1; } int M = 20 + 80.0 / x; T t0 = 0; for (int k = M; k != 0; k--) { t0 = k / (1.0 + k / (x + t0)); } T t = 1.0 / (x + t0); return exp(-x) * t; } Úcupyx_scipy_special_exp1)zf->fzd->dzout0 = exp1(in0)aExponential integral E1. Parameters ---------- x : cupy.ndarray Real argument Returns ------- y : scalar or cupy.ndarray Values of the exponential integral E1 See Also -------- :func:`scipy.special.exp1` )ÚpreambleÚdocN)ÚcupyrÚmath_constants_and_eulÚexp1_definitionÚ create_ufuncÚexp1©óúm/home/jaya/work/projects/VOICE-AGENT/VIET/agent-env/lib/python3.11/site-packages/cupyx/scipy/special/_exp1.pyúrs]ððÐÐÐÐÐðÐð €ð8€uÔØØØØ # oÑ 5ð ð ñô€€€r