npy_cdouble crgamma_wrap( npy_cdouble z) { int kf = 1; npy_cdouble cy; npy_cdouble cy2; double magsq; F_FUNC(cgama,CGAMA)(CADDR(z), &kf, CADDR(cy)); magsq = ABSQ(cy); REAL(cy2) = REAL(cy) / magsq; IMAG(cy2) = -IMAG(cy) / magsq; return cy2; }
Py_complex crgamma_wrap( Py_complex z) { int kf = 1; Py_complex cy; Py_complex cy2; double magsq; F_FUNC(cgama,CGAMA)(CADDR(z), &kf, CADDR(cy)); magsq = ABSQ(cy); REAL(cy2) = REAL(cy) / magsq; IMAG(cy2) = -IMAG(cy) / magsq; return cy2; }