npy_cdouble cpsi_wrap( npy_cdouble z) { npy_cdouble cy; if (IMAG(z)==0.0) { REAL(cy) = cephes_psi(REAL(z)); IMAG(cy) = 0.0; } else { F_FUNC(cpsi,CPSI)(CADDR(z), CADDR(cy)); } return cy; }
Py_complex cpsi_wrap( Py_complex z) { Py_complex cy; if (IMAG(z)==0.0) { REAL(cy) = cephes_psi(REAL(z)); IMAG(cy) = 0.0; } else { F_FUNC(cpsi,CPSI)(CADDR(z), CADDR(cy)); } return cy; }