ʽXc@`sdZddlmZmZmZdZdZdZdZdZ dZ dd l Z dd l Z dd l Z dd lZdd lZdd lZed Zi d d 6d d6dd6dd6dd6dd6dd6dd6dd6Zi dd 6dd6dd6d d6d!d6d"d6d#d6d$d6dd6Zid%d 6d&d6d'd6d(d6d)d6d*d6d+d6d,d6Zi d-d.gd/6d0gd16ddgd26dgd36ddgd46d0d5gd66dgd76dgd86ddgd96Zd:Zd;Zeddfddfdd fddfddfddfddfdd fddfddfddfddfdd fddfddfd dfd dfddfddfddfddfddfgZid<d 6d<d6d<d6d<d6d<d6d<d6d=d6d=d6Zd>Zd?Zd@ZdAZdBZ dCZ!dDZ"dEe#fdFYZ$dGe$fdHYZ%dIe$fdJYZ&dKZ'dLZ(dMZ)dNZ*dOZ+e,dPkre+nd S(Qsz generate_ufuncs.py Generate Ufunc definition source files for scipy.special. Produces files '_ufuncs.c' and '_ufuncs_cxx.c' by first producing Cython. This will generate both calls to PyUFunc_FromFuncAndData and the required ufunc inner loops. The syntax in the ufunc signature list is : '--' '--' : [',' ]* : ':' '*' '->' '*' : * : * : ? : ? : [',' ]* The input parameter types are denoted by single character type codes, according to 'f': 'float' 'd': 'double' 'g': 'long double' 'F': 'float complex' 'D': 'double complex' 'G': 'long double complex' 'i': 'int' 'l': 'long' 'v': 'void' If multiple kernel functions are given for a single ufunc, the one which is used is determined by the standard ufunc mechanism. Kernel functions that are listed first are also matched first against the ufunc input types, so functions listed earlier take precedence. In addition, versions with casted variables, such as d->f,D->F and i->d are automatically generated. There should be either a single header that contains all of the kernel functions listed, or there should be one header for each kernel function. Cython pxd files are allowed in addition to .h files. Cython functions may use fused types, but the names in the list should be the specialized ones, such as 'somefunc[float]'. Function coming from C++ should have ``++`` appended to the name of the header. Floating-point exceptions inside these Ufuncs are converted to special function errors --- which are separately controlled by the user, and off by default, as they are usually not especially useful for the user. The C++ module -------------- In addition to ``_ufuncs`` module, a second module ``_ufuncs_cxx`` is generated. This module only exports function pointers that are to be used when constructing some of the ufuncs in ``_ufuncs``. The function pointers are exported via Cython's standard mechanism. This mainly avoids build issues --- Python distutils has no way to figure out what to do if you want to link both C++ and Fortran code in the same shared library. i(tdivisiontprint_functiontabsolute_importsjL _sf_error_test_function -- _sf_error_test_function: i->i -- sf_error.pxd sph_harm -- sph_harmonic: iidd->D, sph_harmonic_unsafe: dddd->D -- sph_harm.pxd, _legacy.pxd _lambertw -- lambertw_scalar: Dld->D -- lambertw.pxd wrightomega -- wrightomega : D->D -- _wright.h++ _ellip_harm -- ellip_harmonic: ddiiddd->d, ellip_harmonic_unsafe: ddddddd->d --_ellip_harm.pxd, _legacy.pxd logit -- logitf: f->f, logit: d->d, logitl: g->g -- _logit.h expit -- expitf: f->f, expit: d->d, expitl: g->g -- _logit.h bdtrc -- bdtrc: iid->d, bdtrc_unsafe: ddd->d -- cephes.h, _legacy.pxd bdtr -- bdtr: iid->d, bdtr_unsafe: ddd->d -- cephes.h, _legacy.pxd bdtri -- bdtri: iid->d, bdtri_unsafe: ddd->d -- cephes.h, _legacy.pxd binom -- binom: dd->d -- orthogonal_eval.pxd btdtr -- btdtr: ddd->d -- cephes.h btdtri -- incbi: ddd->d -- cephes.h fdtrc -- fdtrc: ddd->d -- cephes.h fdtr -- fdtr: ddd->d -- cephes.h fdtri -- fdtri: ddd->d -- cephes.h gdtrc -- gdtrc: ddd->d -- cephes.h gdtr -- gdtr: ddd->d -- cephes.h hyp0f1 -- _hyp0f1_real: dd->d, _hyp0f1_cmplx: dD->D -- _hyp0f1.pxd hyp2f1 -- hyp2f1: dddd->d, chyp2f1_wrap: dddD->D -- cephes.h, specfun_wrappers.h hyp1f1 -- hyp1f1_wrap: ddd->d, chyp1f1_wrap: ddD->D -- specfun_wrappers.h hyperu -- hypU_wrap: ddd->d -- specfun_wrappers.h hyp2f0 -- hyp2f0: dddi*d->d, hyp2f0_unsafe: dddd*d->d -- cephes.h, _legacy.pxd hyp1f2 -- onef2: dddd*d->d -- cephes.h hyp3f0 -- threef0: dddd*d->d -- cephes.h betainc -- incbet: ddd->d -- cephes.h betaincinv -- incbi: ddd->d -- cephes.h nbdtrc -- nbdtrc: iid->d, nbdtrc_unsafe: ddd->d -- cephes.h, _legacy.pxd nbdtr -- nbdtr: iid->d, nbdtr_unsafe: ddd->d -- cephes.h, _legacy.pxd nbdtri -- nbdtri: iid->d, nbdtri_unsafe: ddd->d -- cephes.h, _legacy.pxd beta -- beta: dd->d -- cephes.h betaln -- lbeta: dd->d -- cephes.h cbrt -- cbrt: d->d -- cephes.h chdtrc -- chdtrc: dd->d -- cephes.h chdtr -- chdtr: dd->d -- cephes.h chdtri -- chdtri: dd->d -- cephes.h ellipeinc -- ellie: dd->d -- cephes.h ellipkinc -- ellik: dd->d -- cephes.h ellipe -- ellpe: d->d -- cephes.h ellipkm1 -- ellpk: d->d -- cephes.h eval_jacobi -- eval_jacobi[double]: dddd->d, eval_jacobi[double complex]: dddD->D, eval_jacobi_l: lddd->d -- orthogonal_eval.pxd eval_sh_jacobi -- eval_sh_jacobi[double]: dddd->d, eval_sh_jacobi[double complex]: dddD->D, eval_sh_jacobi_l: lddd->d -- orthogonal_eval.pxd eval_gegenbauer -- eval_gegenbauer[double]: ddd->d, eval_gegenbauer[double complex]: ddD->D, eval_gegenbauer_l: ldd->d -- orthogonal_eval.pxd eval_chebyt -- eval_chebyt[double]: dd->d, eval_chebyt[double complex]: dD->D, eval_chebyt_l: ld->d -- orthogonal_eval.pxd eval_chebyu -- eval_chebyu[double]: dd->d, eval_chebyu[double complex]: dD->D, eval_chebyu_l: ld->d -- orthogonal_eval.pxd eval_chebyc -- eval_chebyc[double]: dd->d, eval_chebyc[double complex]: dD->D, eval_chebyc_l: ld->d -- orthogonal_eval.pxd eval_chebys -- eval_chebys[double]: dd->d, eval_chebys[double complex]: dD->D, eval_chebys_l: ld->d -- orthogonal_eval.pxd eval_sh_chebyt -- eval_sh_chebyt[double]: dd->d, eval_sh_chebyt[double complex]: dD->D, eval_sh_chebyt_l:ld->d -- orthogonal_eval.pxd eval_sh_chebyu -- eval_sh_chebyu[double]: dd->d, eval_sh_chebyu[double complex]: dD->D, eval_sh_chebyu_l:ld->d -- orthogonal_eval.pxd eval_legendre -- eval_legendre[double]: dd->d, eval_legendre[double complex]: dD->D, eval_legendre_l: ld->d -- orthogonal_eval.pxd eval_sh_legendre -- eval_sh_legendre[double]: dd->d, eval_sh_legendre[double complex]: dD->D, eval_sh_legendre_l:ld->d -- orthogonal_eval.pxd eval_genlaguerre -- eval_genlaguerre[double]: ddd->d, eval_genlaguerre[double complex]: ddD->D, eval_genlaguerre_l:ldd->d -- orthogonal_eval.pxd eval_laguerre -- eval_laguerre[double]: dd->d, eval_laguerre[double complex]: dD->D, eval_laguerre_l:ld->d -- orthogonal_eval.pxd eval_hermite -- eval_hermite: ld->d -- orthogonal_eval.pxd eval_hermitenorm -- eval_hermitenorm: ld->d -- orthogonal_eval.pxd exp10 -- exp10: d->d -- cephes.h exp2 -- exp2: d->d -- cephes.h gamma -- Gamma: d->d, cgamma: D->D -- cephes.h, _loggamma.pxd _gammaln -- lgam: d->d, clngamma_wrap: D->D -- cephes.h, specfun_wrappers.h gammasgn -- gammasgn: d->d -- c_misc/misc.h i0 -- i0: d->d -- cephes.h i0e -- i0e: d->d -- cephes.h i1 -- i1: d->d -- cephes.h i1e -- i1e: d->d -- cephes.h gammaincc -- igamc: dd->d -- cephes.h gammainc -- igam: dd->d -- cephes.h gammaincinv -- gammaincinv: dd->d -- cephes.h gammainccinv -- igami: dd->d -- cephes.h iv -- iv: dd->d, cbesi_wrap: dD->D -- cephes.h, amos_wrappers.h ive -- cbesi_wrap_e_real: dd->d, cbesi_wrap_e: dD->D -- amos_wrappers.h ellipj -- ellpj: dd*dddd->*i -- cephes.h expn -- expn: id->d, expn_unsafe: dd->d -- cephes.h, _legacy.pxd exp1 -- exp1_wrap: d->d, cexp1_wrap: D->D -- specfun_wrappers.h expi -- expi_wrap: d->d, cexpi_wrap: D->D -- specfun_wrappers.h kn -- cbesk_wrap_real_int: id->d, kn_unsafe: dd->d -- cephes.h, _legacy.pxd pdtrc -- pdtrc: id->d, pdtrc_unsafe: dd->d -- cephes.h, _legacy.pxd pdtr -- pdtr: id->d, pdtr_unsafe: dd->d -- cephes.h, _legacy.pxd pdtri -- pdtri: id->d, pdtri_unsafe: dd->d -- cephes.h, _legacy.pxd yn -- yn: id->d, yn_unsafe: dd->d -- cephes.h, _legacy.pxd smirnov -- smirnov: id->d, smirnov_unsafe: dd->d -- cephes.h, _legacy.pxd smirnovi -- smirnovi: id->d, smirnovi_unsafe: dd->d -- cephes.h, _legacy.pxd airy -- airy_wrap: d*dddd->*i, cairy_wrap: D*DDDD->*i -- amos_wrappers.h itairy -- itairy_wrap: d*dddd->*i -- specfun_wrappers.h airye -- cairy_wrap_e_real: d*dddd->*i, cairy_wrap_e: D*DDDD->*i -- amos_wrappers.h fresnel -- fresnl: d*dd->*i, cfresnl_wrap: D*DD->*i -- cephes.h, specfun_wrappers.h shichi -- shichi: d*dd->*i, cshichi: D*DD->*i -- cephes.h, _sici.pxd sici -- sici: d*dd->*i, csici: D*DD->*i -- cephes.h, _sici.pxd itj0y0 -- it1j0y0_wrap: d*dd->*i -- specfun_wrappers.h it2j0y0 -- it2j0y0_wrap: d*dd->*i -- specfun_wrappers.h iti0k0 -- it1i0k0_wrap: d*dd->*i -- specfun_wrappers.h it2i0k0 -- it2i0k0_wrap: d*dd->*i -- specfun_wrappers.h j0 -- j0: d->d -- cephes.h y0 -- y0: d->d -- cephes.h j1 -- j1: d->d -- cephes.h y1 -- y1: d->d -- cephes.h jv -- cbesj_wrap_real: dd->d, cbesj_wrap: dD->D -- amos_wrappers.h jve -- cbesj_wrap_e_real: dd->d, cbesj_wrap_e: dD->D -- amos_wrappers.h yv -- cbesy_wrap_real: dd->d, cbesy_wrap: dD->D -- amos_wrappers.h yve -- cbesy_wrap_e_real: dd->d, cbesy_wrap_e: dD->D -- amos_wrappers.h k0 -- k0: d->d -- cephes.h k0e -- k0e: d->d -- cephes.h k1 -- k1: d->d -- cephes.h k1e -- k1e: d->d -- cephes.h kv -- cbesk_wrap_real: dd->d, cbesk_wrap: dD->D -- amos_wrappers.h kve -- cbesk_wrap_e_real: dd->d, cbesk_wrap_e: dD->D -- amos_wrappers.h hankel1 -- cbesh_wrap1: dD->D -- amos_wrappers.h hankel1e -- cbesh_wrap1_e: dD->D -- amos_wrappers.h hankel2 -- cbesh_wrap2: dD->D -- amos_wrappers.h hankel2e -- cbesh_wrap2_e: dD->D -- amos_wrappers.h ndtr -- ndtr: d->d, faddeeva_ndtr: D->D -- cephes.h, _faddeeva.h++ log_ndtr -- log_ndtr: d->d, faddeeva_log_ndtr: D->D -- cephes.h, _faddeeva.h++ ndtri -- ndtri: d->d -- cephes.h psi -- digamma: d->d, cdigamma: D->D -- _digamma.pxd, _digamma.pxd rgamma -- rgamma: d->d, crgamma: D->D -- cephes.h, _loggamma.pxd round -- round: d->d -- cephes.h sindg -- sindg: d->d -- cephes.h cosdg -- cosdg: d->d -- cephes.h radian -- radian: ddd->d -- cephes.h tandg -- tandg: d->d -- cephes.h cotdg -- cotdg: d->d -- cephes.h log1p -- log1p: d->d, clog1p: D->D -- cephes.h, _cunity.pxd expm1 -- expm1: d->d, cexpm1: D->D -- cephes.h, _cunity.pxd cosm1 -- cosm1: d->d -- cephes.h spence -- spence: d->d, cspence: D-> D -- cephes.h, _spence.pxd zetac -- zetac: d->d -- cephes.h struve -- struve_h: dd->d -- misc.h modstruve -- struve_l: dd->d -- misc.h _struve_power_series -- struve_power_series: ddi*d->d -- misc.h _struve_asymp_large_z -- struve_asymp_large_z: ddi*d->d -- misc.h _struve_bessel_series -- struve_bessel_series: ddi*d->d -- misc.h itstruve0 -- itstruve0_wrap: d->d -- specfun_wrappers.h it2struve0 -- it2struve0_wrap: d->d -- specfun_wrappers.h itmodstruve0 -- itmodstruve0_wrap: d->d -- specfun_wrappers.h kelvin -- kelvin_wrap: d*DDDD->*i -- specfun_wrappers.h ber -- ber_wrap: d->d -- specfun_wrappers.h bei -- bei_wrap: d->d -- specfun_wrappers.h ker -- ker_wrap: d->d -- specfun_wrappers.h kei -- kei_wrap: d->d -- specfun_wrappers.h berp -- berp_wrap: d->d -- specfun_wrappers.h beip -- beip_wrap: d->d -- specfun_wrappers.h kerp -- kerp_wrap: d->d -- specfun_wrappers.h keip -- keip_wrap: d->d -- specfun_wrappers.h _zeta -- zeta: dd->d -- cephes.h kolmogorov -- kolmogorov: d->d -- cephes.h kolmogi -- kolmogi: d->d -- cephes.h besselpoly -- besselpoly: ddd->d -- c_misc/misc.h btdtria -- cdfbet3_wrap: ddd->d -- cdf_wrappers.h btdtrib -- cdfbet4_wrap: ddd->d -- cdf_wrappers.h bdtrik -- cdfbin2_wrap: ddd->d -- cdf_wrappers.h bdtrin -- cdfbin3_wrap: ddd->d -- cdf_wrappers.h chdtriv -- cdfchi3_wrap: dd->d -- cdf_wrappers.h chndtr -- cdfchn1_wrap: ddd->d -- cdf_wrappers.h chndtrix -- cdfchn2_wrap: ddd->d -- cdf_wrappers.h chndtridf -- cdfchn3_wrap: ddd->d -- cdf_wrappers.h chndtrinc -- cdfchn4_wrap: ddd->d -- cdf_wrappers.h fdtridfd -- cdff4_wrap: ddd->d -- cdf_wrappers.h ncfdtr -- cdffnc1_wrap: dddd->d -- cdf_wrappers.h ncfdtri -- cdffnc2_wrap: dddd->d -- cdf_wrappers.h ncfdtridfn -- cdffnc3_wrap: dddd->d -- cdf_wrappers.h ncfdtridfd -- cdffnc4_wrap: dddd->d -- cdf_wrappers.h ncfdtrinc -- cdffnc5_wrap: dddd->d -- cdf_wrappers.h gdtrix -- cdfgam2_wrap: ddd->d -- cdf_wrappers.h gdtrib -- cdfgam3_wrap: ddd->d -- cdf_wrappers.h gdtria -- cdfgam4_wrap: ddd->d -- cdf_wrappers.h nbdtrik -- cdfnbn2_wrap: ddd->d -- cdf_wrappers.h nbdtrin -- cdfnbn3_wrap: ddd->d -- cdf_wrappers.h nrdtrimn -- cdfnor3_wrap: ddd->d -- cdf_wrappers.h nrdtrisd -- cdfnor4_wrap: ddd->d -- cdf_wrappers.h pdtrik -- cdfpoi2_wrap: dd->d -- cdf_wrappers.h stdtr -- cdft1_wrap: dd->d -- cdf_wrappers.h stdtrit -- cdft2_wrap: dd->d -- cdf_wrappers.h stdtridf -- cdft3_wrap: dd->d -- cdf_wrappers.h nctdtr -- cdftnc1_wrap: ddd->d -- cdf_wrappers.h nctdtrit -- cdftnc2_wrap: ddd->d -- cdf_wrappers.h nctdtridf -- cdftnc3_wrap: ddd->d -- cdf_wrappers.h nctdtrinc -- cdftnc4_wrap: ddd->d -- cdf_wrappers.h tklmbda -- tukeylambdacdf: dd->d -- cdf_wrappers.h mathieu_a -- cem_cva_wrap: dd->d -- specfun_wrappers.h mathieu_b -- sem_cva_wrap: dd->d -- specfun_wrappers.h mathieu_cem -- cem_wrap: ddd*dd->*i -- specfun_wrappers.h mathieu_sem -- sem_wrap: ddd*dd->*i -- specfun_wrappers.h mathieu_modcem1 -- mcm1_wrap: ddd*dd->*i -- specfun_wrappers.h mathieu_modcem2 -- mcm2_wrap: ddd*dd->*i -- specfun_wrappers.h mathieu_modsem1 -- msm1_wrap: ddd*dd->*i -- specfun_wrappers.h mathieu_modsem2 -- msm2_wrap: ddd*dd->*i -- specfun_wrappers.h lpmv -- pmv_wrap: ddd->d -- specfun_wrappers.h pbwa -- pbwa_wrap: dd*dd->*i -- specfun_wrappers.h pbdv -- pbdv_wrap: dd*dd->*i -- specfun_wrappers.h pbvv -- pbvv_wrap: dd*dd->*i -- specfun_wrappers.h pro_cv -- prolate_segv_wrap: ddd->d -- specfun_wrappers.h obl_cv -- oblate_segv_wrap: ddd->d -- specfun_wrappers.h pro_ang1_cv -- prolate_aswfa_wrap: ddddd*dd->*i -- specfun_wrappers.h pro_rad1_cv -- prolate_radial1_wrap: ddddd*dd->*i -- specfun_wrappers.h pro_rad2_cv -- prolate_radial2_wrap: ddddd*dd->*i -- specfun_wrappers.h obl_ang1_cv -- oblate_aswfa_wrap: ddddd*dd->*i -- specfun_wrappers.h obl_rad1_cv -- oblate_radial1_wrap: ddddd*dd->*i -- specfun_wrappers.h obl_rad2_cv -- oblate_radial2_wrap: ddddd*dd->*i -- specfun_wrappers.h pro_ang1 -- prolate_aswfa_nocv_wrap: dddd*d->d -- specfun_wrappers.h pro_rad1 -- prolate_radial1_nocv_wrap: dddd*d->d -- specfun_wrappers.h pro_rad2 -- prolate_radial2_nocv_wrap: dddd*d->d -- specfun_wrappers.h obl_ang1 -- oblate_aswfa_nocv_wrap: dddd*d->d -- specfun_wrappers.h obl_rad1 -- oblate_radial1_nocv_wrap: dddd*d->d -- specfun_wrappers.h obl_rad2 -- oblate_radial2_nocv_wrap: dddd*d->d -- specfun_wrappers.h modfresnelp -- modified_fresnel_plus_wrap: d*DD->*i -- specfun_wrappers.h modfresnelm -- modified_fresnel_minus_wrap: d*DD->*i -- specfun_wrappers.h wofz -- faddeeva_w: D->D -- _faddeeva.h++ erfc -- erfc: d->d, faddeeva_erfc: D->D -- cephes.h, _faddeeva.h++ erf -- erf: d->d, faddeeva_erf: D->D -- cephes.h, _faddeeva.h++ dawsn -- faddeeva_dawsn: d->d, faddeeva_dawsn_complex: D->D -- _faddeeva.h++ erfcx -- faddeeva_erfcx: d->d, faddeeva_erfcx_complex: D->D -- _faddeeva.h++ erfi -- faddeeva_erfi: d->d, faddeeva_erfi_complex: D->D -- _faddeeva.h++ xlogy -- xlogy[double]: dd->d, xlogy[double_complex]: DD->D -- _xlogy.pxd xlog1py -- xlog1py[double]: dd->d, xlog1py[double_complex]: DD->D -- _xlogy.pxd poch -- poch: dd->d -- c_misc/misc.h boxcox -- boxcox: dd->d -- _boxcox.pxd boxcox1p -- boxcox1p: dd->d -- _boxcox.pxd inv_boxcox -- inv_boxcox: dd->d -- _boxcox.pxd inv_boxcox1p -- inv_boxcox1p: dd->d -- _boxcox.pxd entr -- entr: d->d -- _convex_analysis.pxd kl_div -- kl_div: dd->d -- _convex_analysis.pxd rel_entr -- rel_entr: dd->d -- _convex_analysis.pxd huber -- huber: dd->d -- _convex_analysis.pxd pseudo_huber -- pseudo_huber: dd->d -- _convex_analysis.pxd exprel -- exprel: d->d -- _exprel.pxd _spherical_yn -- spherical_yn_real: ld->d, spherical_yn_complex: lD->D -- _spherical_bessel.pxd _spherical_jn -- spherical_jn_real: ld->d, spherical_jn_complex: lD->D -- _spherical_bessel.pxd _spherical_in -- spherical_in_real: ld->d, spherical_in_complex: lD->D -- _spherical_bessel.pxd _spherical_kn -- spherical_kn_real: ld->d, spherical_kn_complex: lD->D -- _spherical_bessel.pxd _spherical_yn_d -- spherical_yn_d_real: ld->d, spherical_yn_d_complex: lD->D -- _spherical_bessel.pxd _spherical_jn_d -- spherical_jn_d_real: ld->d, spherical_jn_d_complex: lD->D -- _spherical_bessel.pxd _spherical_in_d -- spherical_in_d_real: ld->d, spherical_in_d_complex: lD->D -- _spherical_bessel.pxd _spherical_kn_d -- spherical_kn_d_real: ld->d, spherical_kn_d_complex: lD->D -- _spherical_bessel.pxd loggamma -- loggamma: D->D -- _loggamma.pxd _sinpi -- dsinpi: d->d, csinpi: D->D -- _trig.pxd _cospi -- dcospi: d->d, ccospi: D->D -- _trig.pxd _lgam1p -- lgam1p: d->d -- cephes.h _lanczos_sum_expg_scaled -- lanczos_sum_expg_scaled: d->d -- cephes.h _log1pmx -- log1pmx: d->d -- cephes.h _igam_fac -- igam_fac: dd->d -- cephes.h s~# This file is automatically generated by generate_ufuncs.py. # Do not edit manually! include "_ufuncs_extra_code_common.pxi" s!include "_ufuncs_extra_code.pxi" s# # Aliases # jn = jv sV# This file is automatically generated by generate_ufuncs.py. # Do not edit manually! s6# This file is automatically generated by generate_ufuncs.py. # Do not edit manually! """ .. highlight:: cython ================================ Cython API for Special Functions ================================ Scalar, typed versions of many of the functions in ``scipy.special`` can be accessed directly from Cython; the complete list is given below. Functions are overloaded using Cython fused types so their names match their ufunc counterpart. The module follows the following conventions: - If a function's ufunc counterpart returns multiple values, then the function returns its outputs via pointers in the final arguments - If a function's ufunc counterpart returns a single value, then the function's output is returned directly. The module is usable from Cython via:: cimport scipy.special.cython_special Error Handling ============== Functions can indicate an error by returning ``nan``; however they cannot emit warnings like their counterparts in ``scipy.special``. Available Functions =================== FUNCLIST """ include "_cython_special.pxi" Nt add_newdocstfloattftdoubletds long doubletgs float complextFsdouble complextDslong double complextGtinttitlongtltvoidtvt npy_floatt npy_doubletnpy_longdoublet npy_cfloatt npy_cdoubletnpy_clongdoubletnpy_inttnpy_longt NPY_FLOATt NPY_DOUBLEtNPY_LONGDOUBLEt NPY_CFLOATt NPY_CDOUBLEtNPY_CLONGDOUBLEtNPY_INTtNPY_LONGsd*ddddsD*DDDDtairytddtbetaterftexpreltgammatdDtjvtloggammatlogittpsicC`s|jddS(Nt t_(treplace(targ((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt underscorescC`s g|D]}dj|^qS(NtilfdgFDG(tindex(tctx((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt cast_orderstNPY_NANt 0xbad0bad0c C`st|t|kr'tdnt|t|krw|dkodt|dt|k rwtdnd|||||f}d|}|d7}|d7}|d 7}x.tt|D]}|d ||f7}qWx8tt|D]$}|d ||t|f7}qWg}g} g} xYtt|D]E}|jt||| jd t||t|||fqIWt|dt|krd} | j||d t|f7}nd} xmt|D]_\}} |dt| || f7}|jdt| | jd|| | j| qW|d7}t|dt|kr~d} nd} d| t|dj|dj| f}g}xutt|D]a}||||ftkrdt||t|||t|||f}|j|qqW|r|ddj|7}|d|7}|d7}|d7}xFt| D]+\}} |d|t| t| f7}qzWn ||7}xtt || D]\}\} }|| ftkrj|d|t| |f7}|dt| |t| |f7}|d7}|d 7}|d!t| |t| t| f7}q|d"t| |t| |f7}qWx.tt|D]}|d#||f7}qWx8tt|D]$}|d$||t|f7}qW|d%7}||fS(&s Generate a UFunc loop function that calls a function given as its data parameter with the specified input and output arguments and return value. This function can be passed to PyUFunc_FromFuncAndData. Parameters ---------- func_inputs, func_outputs, func_retval : str Signature of the function to call, given as type codes of the input, output and return value arguments. These 1-character codes are given according to the CY_TYPES and TYPE_NAMES lists above. The corresponding C function signature to be called is: retval func(intype1 iv1, intype2 iv2, ..., outtype1 *ov1, ...); If len(ufunc_outputs) == len(func_outputs)+1, the return value is treated as the first output argument. Otherwise, the return value is ignored. ufunc_inputs, ufunc_outputs : str Ufunc input and output signature. This does not have to exactly match the function signature, as long as the type casts work out on the C level. Returns ------- loop_name Name of the generated loop function. loop_body Generated C code for the loop. s2Function and ufunc have different number of inputsRis-Function retval and ufunc outputs don't matchsloop_%s_%s_%s_As_%s_%ssTcdef void %s(char **args, np.npy_intp *dims, np.npy_intp *steps, void *data) nogil: s$ cdef np.npy_intp i, n = dims[0] s( cdef void *func = (data)[0] s4 cdef char *func_name = (data)[1] s cdef char *ip%d = args[%d] s cdef char *op%d = args[%d] s<%s>(<%s*>ip%d)[0]s cdef %s ov0 is cdef %s ov%d s%s *s&ov%ds for i in range(n): sov0 = ts& %s(<%s(*)(%s) nogil>func)(%s) s, s$<%s>(<%s*>ip%d)[0] == (<%s*>ip%d)[0]s if %s: s and s s else: sQ sf_error.error(func_name, sf_error.DOMAIN, "invalid input argument") s ov%d = <%s>%s s if ov%d == <%s>ov%d: s' (<%s *>op%d)[0] = <%s>ov%d sI sf_error.error(func_name, sf_error.DOMAIN, "invalid output") s% (<%s *>op%d)[0] = <%s>%s s# (<%s *>op%d)[0] = <%s>ov%d s ip%d += steps[%d] s op%d += steps[%d] s" sf_error.check_fpe(func_name) ( tlent ValueErrortrangetappendtCY_TYPESt enumeratetjointDANGEROUS_DOWNCASTt NAN_VALUEtzip(t func_inputst func_outputst func_retvalt ufunc_inputst ufunc_outputstnametbodytjtftypestfvarst outtypecodest func_jofftouttypetrvtfuncallt input_checkstchktfouttype((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt generate_loops')    "     )    (  # " cC`shtd|}|d}d|dg}x|D]}|jd|q4Wdj|}||fS(s Generate name of and cython code for a fused type. Parameters ---------- typecodes : str Valid inputs to CY_TYPES (i.e. f, d, g, ...). cS`st|S(N(R>(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytrst _number_tsctypedef fused t:s s (tmapR=R@(tcodestcytypesRIt declarationtcytype((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytgenerate_fused_typehs   cC`s2d}ggg}}}|jd}t|dkrR|d}|d}n|d}d}gg} } xSt|D]E\} } dj| } | j| | jdjt| | q|Wd jd j| }|j||j|d g}xft|D]X\} } d j| } |jd j| djt| | }|j||qW|j|d|jd|ddj||d}dj||d}dj||}dj|||}|j|d|d j| }|j||d j| |}||fS(NR-it*iiR9sx{}s{} {}sdef {{}}(int N, {}):s, s cdef int nsy{}s&{}s cdef {} {}sfor n in range(N):is{}({})s_bench_{}_{}_{}tpytcys s_ufuncs.s (tsplitR:R?tformatR=R>R@(RIR[ttabttoptmiddletendttmptincodestoutcodestinargstinargs_and_typestntcodeR0tlinetoutargstpyfunctcyfunct pytemplatet cytemplatetpybenchtcybench((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytgenerate_bench{s>     ! ""c C`sd}dj|g}x|D]}|jd\}}|jd}ttd|d}t|dkrtd |d}|j|nt|} d j| |d j|} |jd || qWd j|d|dR`cS`st|S(N(R>(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWsiicS`sdjt|S(Ns{} *(RdR>(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWss {} {}({})s, is{} is s ( RdRctlistRZR:textendR>R@R=( RItspecsRetdoctspecRjRktintypesttypesRPRp((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt generate_docs  cC`sdj|}|S(s0Cast a cython double complex to a numpy cdouble.s1_complexstuff.npy_cdouble_from_double_complex({})(Rd(tvartres((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytnpy_cdouble_from_double_complexscC`sdj|}|S(s0Cast a numpy cdouble to a cython double complex.s1_complexstuff.double_complex_from_npy_cdouble({})(Rd(RR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytdouble_complex_from_npy_cdoublesc c`sdg}d|kpd|ksU|g|D] \}}|d|df^q+}nxo|D]g\}}|}|}xAt||D]0\}}|j||}|j||}qW||fVq\WdS(s. Generate variants of UFunc signatures, by changing variable types, within the limitation that the corresponding C types casts still work out. This does not generate all possibilities, just the ones required for the ufunc to work properly with the most common data types. Parameters ---------- inputs, outputs : str UFunc input and output signature strings Yields ------ new_input, new_output : str Modified input and output strings. Also the original input/output pair is yielded. R RR(tfFN(R R(RCR/( tinputstoutputstmapstatbtsrctdstt new_inputst new_outputs((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt iter_variantss 4tFunccB`sVeZdZdZdZdZedZededZ e dZ RS(s. Base class for Ufunc and FusedFunc. cC`s+||_|j|||_i|_dS(N(RIt_parse_signaturest signaturestfunction_name_overrides(tselfRIRtheaders((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt__init__s cC`sg|jdD]}|jr|j^q}g|jdD]}|jrD|j^qD}t|dkr|t|}nt|t|krtd|j||fngt||D]"\}}|j||f^qS(Nt,is<%s: Number of headers and signatures doesn't match: %r -- %r(RctstripR:R;RIRCt_parse_signature(Rtsigs_strt headers_strR5tsigsRth((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs44cC`stjd|}|rg|jD]}|j^q%\}}}}|jddkrztd|j|fn||||fStjd|}|rg|jD]}|j^q\}}}||d|fStd|j|fdS(NsG\s*(.*):\s*([fdgFDGil]*)\s*\*\s*([fdgFDGil]*)\s*->\s*([*fdgFDGil]*)\s*$R`is%s: Invalid signature: %rs1\s*(.*):\s*([fdgFDGil]*)\s*->\s*([fdgFDGil]?)\s*$R9(tretmatchtgroupsRtcountR;RI(RtsigtmR5tfunctinargtoutargtret((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs1.c C`s,g}x|jD]\}}}}}|jdd}g|D]}t|^q>g|D]}t|d^qX} g|D]}t|^qzg|D]}t|d^q} dt|dj| f} |jdr|r| d} ndt|dj| f} |j|| | |fqW|S( NR`R9s *s %s (*)(%s)s, Rtnogils%s (*)(%s) nogil(RR/tC_TYPESR>R@tendswithR=( Rt nptypes_for_ht prototypest func_nameRRRtheaderR5tc_argstcy_argstc_prototcy_proto((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytget_prototypess"" t_func_cC`s|r+||jkr+|j|}d}ntjd|}|rX|j\}}n |d}}|rd|||jddfSd||fSdS(NR9s^(.*?)(\[.*\])$s%s%s%sR-R.s%s%s(RRRRR/(Rtc_namet specializedtprefixtoverrideRt c_base_namet fused_part((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytcython_func_names   cC`sg}|j}|jx|D]}}|j}|sAq#ntjd|}|sltd|n|j||jd|jd|jdq#W|S(Ns$^([a-z0-9_]+)\s*--\s*(.*?)\s*--(.*)$sUnparseable line %riii(t splitlinestsortRRRR;R=tgroup(tclst ufunc_strtufuncstlinesRpR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt parse_all1s    8( t__name__t __module__t__doc__RRRtFalseRtTrueRt classmethodR(((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs   tUfunccB`s)eZdZdZdZdZRS(s Ufunc signature, restricted format suitable for special functions. Parameters ---------- name Name of the ufunc to create signature String of form 'func: fff*ff->f, func2: ddd->*i' describing the C-level functions and types of their input arguments and return values. The syntax is 'function_name: inputparams*outputparams->output_retval*ignored_retval' Attributes ---------- name : str Python name for the Ufunc signatures : list of (func_name, inarg_spec, outarg_spec, ret_spec, header_name) List of parsed signatures doc : str Docstring, obtained from add_newdocs function_name_overrides : dict of str->str Overrides for the function names in signatures cC`sstt|j|||tjd||_|jdkrTtd|ntj |jj |_dS(Nsscipy.special.sNo docstring for ufunc %r( tsuperRRRtgetR|tNoneR;ttextwraptdedentR(RRIRR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyR^s c `sddtgfd}xjD]\}}}}}tjdd||}|jdd}dkrt|t|ntt||d\} }|||||| |qCWxjD]y\}}}}}tjdd||}|jdd}x6t||D]%\} }|||||| |q@WqWj ddfS(Nc `s|krdSj||||f}d|krTtdj|fnt|ksxt|krtdj|fnt|||||\}}||<j||||fdS(NRs%s: void signature %rs6%s: signature %r does not have %d/%d input/output args(taddR;RIR:RVR=( RRRRtinptoutpRt loop_nametloop(t all_loopst inarg_numt outarg_numtseenRtvariants(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt add_variantls   $  s\*.*R9R`itkeycS`st|dS(Ni(R6(R((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWs( RtsetRRtsubR/R:RyRR( RRRRRRRRRR((RRRRRRs</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt_get_signatures_and_loopses(   !c C`sd}|j|\}}}g}g}g}xx|D]p\} } } } x| D]} |jt| qPWx| D]} |jt| qrW|j| |j| q7W|d|jt|f7}|d|jdt|f7}|d|jt|f7}|d|jt|f7}|d|j|jjdd jd d jd d f7}x4t|D]&\}}|d|j||f7}qnWx4t|D]&\}}|d|j||f7}qWx`t|D]R\}}|d|j||j|dt f7}|d|j||jf7}qWx:t|D],\}}|d|j||j|f7}q?W|dt|||||fjd|j7}|S(NR9s2cdef np.PyUFuncGenericFunction ufunc_%s_loops[%d] scdef void *ufunc_%s_ptr[%d] iscdef void *ufunc_%s_data[%d] scdef char ufunc_%s_types[%d] s&cdef char *ufunc_%s_doc = ( "%s") s\s\\t"s\"s s \n" "s3ufunc_%s_loops[%d] = %s sufunc_%s_types[%d] = %s sufunc_%s_ptr[2*%d] = %s Rs,ufunc_%s_ptr[2*%d+1] = ("%s") s(ufunc_%s_data[%d] = &ufunc_%s_ptr[2*%d] so@ = np.PyUFunc_FromFuncAndData(ufunc_@_loops, ufunc_@_data, ufunc_@_types, %d, %d, %d, 0, "@", ufunc_@_doc, 0) t@( RR=t TYPE_NAMESRIR:R|R/R?RR(RRttoplevelRRRtloopstfuncsRRRRRR5RKtfunctionttypeR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytgeneratesH   !2 (RRRRRR(((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRCs  1t FusedFunccB`seZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zd ZdZdZRS(s_ Generate code for a fused-type special function that can be cimported in cython. cC`stt|j|||d|j|_|j\|_|_t|_ |j |j\|_ }|j j ||j |j\|_ }|j j ||j\|_|_dS(Ns(See the documentation for scipy.special.(RRRRIR|t _get_codesRjRkRt fused_typest _get_typesR~tupdatetouttypest _get_varstinvarstoutvars(RRIRRt infused_typestoutfused_types((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs c`sd\}}gg}}x|jD]\}}}}}tjdd||} |dkrxt|}t| }ntt|| d\} } |j| |j| q#Wg} xRt|D]Dt t fd|} | j | jdj | qWg} xRt|D]Dt t fd|} | j | jdj | q#Wt | t | fS(Ns\*.*R9ic`s|S(N((R5(Rn(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWsc`s|S(N((R5(Rn(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWs(NN(RRRRR:RyRR=R<tuniqueRZRR@ttuple(RRRtall_inptall_outpR.RRRRRRjR[Rk((Rns</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs*       cC`sg}t}xl|D]d}t|dkrH|jt||fqt|\}}|j||j||fqW||fS(Ni(RR:R=R>R_R(RR[t all_typesRRot fused_typetdec((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs   cC`sg}x3tt|jD]}|jdj|qWg}x3tt|jD]}|jdj|qXW||fS(Nsx{}sy{}(R<R:R~R=RdR(RRRnR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRsc C`sg}t}xzt||D]i\\}}}t|dkrIqn||kr|jdj|tt||j|qqW|r|dkrdj|dj|} n!|r|dkrd} nd} | S(s_Generate an if/elif/else clause that selects a specialization of fused types. is{} is {}telses{} {}:s and selse:N( RRCR:R=RdR1R>RR@R( RRR[tadverbtclausesRttypttypcodeRoRp((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt_get_conditionals " # cC`sbg}xUt|D]G\}}|j|}|rM|dkrMt|}n|j|qW|S(swGenerate pure input variables to a specialization, i.e. variables that aren't used to return a value. sdouble complex(R?RRR=(RR~R4t incallvarsRntintypeR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt_get_incallvarss c C`sggg}}}t|jt|}|j|}xtt||D]\}\} } |r| dkrdj|} |j| |jdj| t| } |jdj| | qP|jdj| qPW|||fS(snGenerate output variables to a specialization, i.e. pointers that are used to return values. sdouble complexstmp{}s&{}s {}[0] = {}s{}(R:RR?RCRdR=R( RRR4t outcallvarsttmpvarstcastststartRRnRRPRittmpcast((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt_get_outcallvars*s (  c`sd gdg}}t}xt|j|j|jD]\}}}t|dkrdj|t|}|jd|q?|}|\} } | |kr?|j||j | q?q?W|s|Sg} x'|D]}|\} } | j| qWt | } t t d| } t j| d}x@ttj| D])\}} gg}}xt|D]\}}| |}|jtt|xVt|jD]E\}}|j||krdj|t|}|j|qqWqzW|dkrd }n||kr.d }nd }|j|| |}|jd||jt fd |qTW|S(snSet all variables to nan for specializations of fused types for which don't have signatures. R-iselse:is {}[0] = {}icS`s t|S(N(R:(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRW\sitifRtelifc`s d|S(Ni((R5(Re(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWoss (RRCRRRkR:RdRBR=RRRyRZtnumpytproductR?t itertoolsR1R>RRz(RRRRtoutvarRPRoRpRRIR.t all_codesR[tcodelenstlastRt fused_codestdecsRnRtcond((Res</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt _get_nan_decs?sN +      "      #cC`sWd}t|}|jg}x.|D]&}dj|}|j||q)W|S(sQGenerate the declarations of any necessary temporary variables. R-iscdef npy_cdouble {}s (RyRRdR=(Rt all_tmpvarsReRttmpdecsttmpvarRp((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt _get_tmp_decsrs   c C`sdd }gg}}xBt|j|jD]+\\}}}|jdj||q)Wdj|jdj|}|j|xLt|j|jD]5\\}}} dj|| }|j||qWt d|j} dj|jdj|jdj| }|j||d jdj|j}|j||d j|}|S( saGenerate a python wrapper for functions which pass their arguments as pointers. R-is{} {}sdef _{}_pywrap({}):s, s cdef {} {}cS`s dj|S(Ns&{}(Rd(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWss {}({}, {})s return {}s s ( RCR~RR=RdRIR@RRRZ( RReRJtcallvarsRR.tinvarRpRPRt addr_outvars((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt_get_python_wraps" ( (c C`sd}|\}}}}}|jdd}|jdd}|jdd}|jdrit} nt} |jdrt} nt} ttd|} ttd|} tjd d |}|sd }nt|} | r%d j |}d j | dj | | |}n|j |dt}|dkrOd}nd}|j |j ||}|r||g}d|}n g}|}||||| | | | ||f S(s4Generate code common to all the _generate_* methods.R-iR RRs++cS`st|S(N(R>(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWscS`st|S(N(R>(R5((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWss\*.*R9Rs$scipy.special._ufuncs_cxx._export_{}s(<{}(*)({}) nogil>{})s, RiR R is (R/RRRRyRZRRR>RdR@RRR~(RtsignumRReRRjRktretcodeRR4tcppR~RtrettypeRRRtsp((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt _get_commonsB         c C`sd}gg}}xt|jD]\}}|j||\ }}}} } } } } }}|j||j| | }dj|dj|}| r| dkrt|}n|dj|}|j|dj|| }|j|q#Wt |dkr%|j|d |j d \}}t |d}t |dkrdjt |}|jd ||q%xt|D]\}}|d krd }n||krd }nd}|j |j ||}|jd ||djt |}|jd||qWngg}}xCt|j D]2\}\}}|jdj||j|qBW|j d \}}dj||jdj|}|j|d|j|dj|jdj||}|||fS(NR-is{}({})s, sdouble complexs return {}s{}->{}iselse:iiR RR is{} {}scpdef {} {}({}) nogilRYs"""{}"""s s (R?RR#RzRRdR@RR=R:RRBRR~RRIR|(RReR{RJRRRRjRkRR~RR!R4RR"RtcallRpRPRRnRoRRtheadRR.RR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt$_generate_from_return_and_no_outargssN 0       "$!c `sd}t}gg}}x?t|jD].\}}|j||\ }}} } } } } }}|j||j| |}|j| |\}}}|j||j|dj|dj |}|j ||jt fd|t | dkr;dj|| }|j |q,dj|| }|j |q,Wt |dkr|j }|j|nt |jdkrd j|jd }|j ||ngg}}xBt|j|jD]+\}\}}|j d j||qWt |jdkrzxEt|j|jD]+\}\}}|j d j||qHWnt |jdkr|jd \}}d j||jdj |}ndj|jdj |}|j |d|j |dj|jt |jdkrm|jd }|jd \}}dj||}|j ||n|j|j|dj ||}|||fS(NR-is{}({})s, c`s|S(N((R5(R"(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRW sis{}->{}s{}*{}->vs return {}[0]is{} {}s{} *{}scpdef {} {}({}) nogilscdef void {}({}) nogilRYs"""{}"""s cdef {} {}s s (RR?RR#RzRRRRdR@R=RZR:RRRCRR~RRIR|R(RReRR{RJRRRRjRkRR~RR!R4RRRRRR$RpR%RRR.RRPRR((R"s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt$_generate_from_outargs_and_no_returnsX  0     (( $ c `sd}t}gg}}xIt|jD]8\}}|j||\ }}} } } } } }}|j||j| |}|j| |\}}}|j||j|dj|dj |}|r| dkrt |}ndj|j d|}|j ||jt fd|d j|| | }|j |q,Wt|d kr|j}|j|ngg}}xBt|j|jD]+\}\}}|j d j||qWxBt|j |jD]+\}\}}|j d j||qWd j|jdj |}|j |d|j |dj|j|j|j|dj ||}|||fS(NR-is{}({})s, sdouble complexs {}[0] = {}ic`s|S(N((R5(R"(s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRWGss{}*{}->vis{} {}s{} *{}scdef void {}({}) nogilRYs"""{}"""s s (RR?RR#RzRRRRdR@RRR=RZR:RRCRR~RRIR|R(RReRR{RJRRRRjRkRR~RR!R4RRRRRR$R%RRR.RRPRR((R"s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyt!_generate_from_outargs_and_return4s@  0     ((cC`s.|jd\}}}}}tjdd|}|s@d}nt|dkrv|dkrv|j\}}}nxt|dkr|dkr|j\}}}nBt|dkr|dkr|j\}}}n tdt|jdkr|j }nd}||||j |fS(Nis\*.*R9RsInvalid signaturei( RRRR:R&R'R(R;RRRR(RR.RkRRRR{twrap((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyR]s  (RRRRRRRRRRRRRR#R&R'R(R(((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyRs      3   1 3 9 )c C`sPg}g}|jddjddjdd}|jdr|jd|d |j|dd |j|fd |} |jd |jd | |jd | | |j|dtfnud|j||f} |jd||jd|jd | |jd||jd|jd ||||fS(s Construct a Cython declaration of a function coming either from a pxd or a header file. Do sufficient tricks to enable compile-time type checking against the signature expected by the ufunc. t[R.t]R-s.pxdsfrom %s cimport %s as %siRR9s _proto_%s_ts ctypedef %ss(*)scdef %s *%s_var = &%sRs%s "%s"scdef extern from "%s":s cdef %ss #include "%s"s%s;(R/RR=RR( tufuncRRRRtproto_h_filenametdefstdefs_htvar_namet proto_nametnew_name((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytget_declarationts"*   #c C`sc|d}|d}|d}|d}|d}d}g} g} i} g} ddg} g}|jddx=|D]5}|j}x|D]\}}}}|jd rR|d }t||||||\}}}| j||j|| jd ||j|d td tf| jd|fd||j |%sRRscdef void *_export_%ss"scipy.special._ufuncs_cxx._export_s tws0#ifndef UFUNCS_PROTO_H #define UFUNCS_PROTO_H 1 s #endif s # distutils: language = c++ (RRRR3RzR=RRRRRR@tsortedtvaluestopentwritetUFUNCS_EXTRA_CODE_COMMONtUFUNCS_EXTRA_CODEtUFUNCS_EXTRA_CODE_BOTTOMR(t fn_prefixt cxx_fn_prefixRtfilenameR-tcxx_proto_h_filenametcxx_pyx_filenametcxx_pxd_filenameRR.R/Rtcxx_defst cxx_pxd_defst cxx_defs_hR,tcfuncsRRRRt item_defst item_defs_hR0R.ttR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytgenerate_ufuncssr            # &          c C`stjjt}tjj||d}tjj||d}tjj||d}g}g}g} t} g} g} xa|D]Y} | jjdrqn| j\}}}}}|j ||j ||r|j |n| j || j dt }x`|D]X\}}}}|j drAqnt| |||||\}}}| j|qW| j t| j|| jtkrxCt| jD]1}t| j|\}}| j||gqWqqWt| } | jt|dZ}|jt|jd|jd j| |jd |jdj|WdQXt|d}t}|jd dj| }|j||jd|jdj| |jd |jd j||jd |jd j| WdQXdS( Ns.pxds.pyxs_defs.hR.Rs++R5s s tFUNCLIST(tostpathtdirnamet__file__R@RRIt startswithRR=RRRRR3RzRtCYTHON_SPECIAL_BENCHFUNCSRxRyRR8R9tCYTHON_SPECIAL_PXDtCYTHON_SPECIAL_PYXR/(tmodnametufunc_fn_prefixt fused_funcstpwdtpxdfiletpyxfileR-tsourcest declarationst bench_auxRR|R.RRRR{tfunc_fused_typesR)RFRRRRRGR.R[RvRwR((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytgenerate_fused_funcssd               cC`sPt}g}x:|D]2}||kr.qn|j||j|qW|S(sb Return a list without repeated entries (first occurrence is kept), preserving order. (RRR=(tlstRtnew_lsttitem((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pyR4s    cC`stjdtj}|j\}}t|dkrL|jdntjt }t dd|t jt }t dd|dS(Ntusageisinvalid number of argumentst_ufuncst _ufuncs_cxxtcython_special( toptparset OptionParserRRt parse_argsR:terrorRRtFUNCSRJRR^(tptoptionstargsRRV((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytmainCst__main__(-Rt __future__RRRRjR:R;R<RRRSRLRfRRR R t __import__RR>RRRQR1R6RRARBRVR_RxRRRRtobjectRRRR3RJR^RRnR(((s</tmp/pip-build-7oUkmx/scipy/scipy/special/generate_ufuncs.pytGs '                    ---   '    .U " W G