npy_cdouble cexp1_wrap(npy_cdouble z) { npy_cdouble outz; F_FUNC(e1z,E1Z)(&z, &outz); ZCONVINF("cexp1", outz); return outz; }
npy_cdouble cexpi_wrap(npy_cdouble z) { npy_cdouble outz; F_FUNC(eixz,EIXZ)(&z, &outz); ZCONVINF("cexpi", outz); return outz; }
Py_complex cexpi_wrap(Py_complex z) { Py_complex outz; F_FUNC(eixz,EIXZ)(&z, &outz); ZCONVINF(outz); return outz; }
Py_complex cexp1_wrap(Py_complex z) { Py_complex outz; F_FUNC(e1z,E1Z)(&z, &outz); ZCONVINF(outz); return outz; }
double keip_wrap(double x) { npy_cdouble Be, Ke, Bep, Kep; if (x<0) return NPY_NAN; F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF("keip", Kep); return IMAG(Kep); }
double bei_wrap(double x) { npy_cdouble Be, Ke, Bep, Kep; if (x<0) x=-x; F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF("bei", Be); return IMAG(Be); }
double kerp_wrap(double x) { Py_complex Be, Ke, Bep, Kep; if (x<0) return NPY_NAN; F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF(Kep); return REAL(Kep); }
double ber_wrap(double x) { Py_complex Be, Ke, Bep, Kep; if (x<0) x=-x; F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF(Be); return REAL(Be); }
int kelvin_wrap(double x, npy_cdouble *Be, npy_cdouble *Ke, npy_cdouble *Bep, npy_cdouble *Kep) { int flag = 0; if (x<0) {x=-x; flag=1;} F_FUNC(klvna,KLVNA)(&x, F2C_CST(Be), F2C_CST(Ke), F2C_CST(Bep), F2C_CST(Kep)); ZCONVINF("klvna", *Be); ZCONVINF("klvna", *Ke); ZCONVINF("klvna", *Bep); ZCONVINF("klvna", *Kep); if (flag) { REAL(*Bep) = -REAL(*Bep); IMAG(*Bep) = -IMAG(*Bep); REAL(*Ke) = NPY_NAN; IMAG(*Ke) = NPY_NAN; REAL(*Kep) = NPY_NAN; IMAG(*Kep) = NPY_NAN; } return 0; }
double beip_wrap(double x) { npy_cdouble Be, Ke, Bep, Kep; int flag = 0; if (x<0) {x=-x; flag=1;} F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF("beip", Bep); if (flag) return -IMAG(Bep); return IMAG(Bep); }
double berp_wrap(double x) { Py_complex Be, Ke, Bep, Kep; int flag = 0; if (x<0) {x=-x; flag=1;} F_FUNC(klvna,KLVNA)(&x, CADDR(Be), CADDR(Ke), CADDR(Bep), CADDR(Kep)); ZCONVINF(Bep); if (flag) return -REAL(Bep); return REAL(Bep); }