void hydrogenicR_1(double *Z, double *r, int *len, double *val, double *err, int *status) { int i; gsl_sf_result result; gsl_set_error_handler_off(); for(i = 0; i< *len ; i++){ status[i] = gsl_sf_hydrogenicR_1_e(Z[i], r[i], &result) ; val[i] = result.val; err[i] = result.err; } }
/** * C++ version of gsl_sf_hydrogenicR_1(). * Normalized hydrogenic bound states, radial dependence. * R_1 := 2Z sqrt(Z) exp(-Z r) * @param Z A real value * @param r A real value * @param result The result as a @c gsl::sf::result object * @return Error code on failure */ inline int hydrogenicR_1_e( double const Z, double const r, result& result ){ return gsl_sf_hydrogenicR_1_e( Z, r, &result ); }
double gsl_sf_hydrogenicR_1(const double Z, const double r) { EVAL_RESULT(gsl_sf_hydrogenicR_1_e(Z, r, &result)); }