/// Irregular modified cylindrical Bessel functions. double cyl_bessel_k(double nu, double x) { gsl_sf_result result; int stat = gsl_sf_bessel_Knu_e(nu, x, &result); if (stat != GSL_SUCCESS) { std::ostringstream msg("Error in cyl_bessel_k:"); msg << " nu=" << nu << " x=" << x; throw std::runtime_error(msg.str()); } else return result.val; }
double gsl_sf_bessel_Knu(const double nu, const double x) { EVAL_RESULT(gsl_sf_bessel_Knu_e(nu, x, &result)); }