#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r22501 = x; float r22502 = exp(r22501); float r22503 = 1; float r22504 = r22502 - r22503; float r22505 = r22504 / r22501; return r22505; } double f_id(double x) { double r22506 = x; double r22507 = exp(r22506); double r22508 = 1; double r22509 = r22507 - r22508; double r22510 = r22509 / r22506; return r22510; } double f_of(float x) { float r22511 = 1/6; float r22512 = x; float r22513 = 2; float r22514 = pow(r22512, r22513); float r22515 = r22511 * r22514; float r22516 = 1; float r22517 = 1/2; float r22518 = r22517 * r22512; float r22519 = r22516 + r22518; float r22520 = r22515 + r22519; float r22521 = 1.0000018458005175; bool r22522 = r22520 <= r22521; float r22523 = r22512 + r22512; float r22524 = exp(r22523); float r22525 = r22524 - r22516; float r22526 = exp(r22512); float r22527 = r22526 + r22516; float r22528 = r22525 / r22527; float r22529 = r22528 / r22512; float r22530 = r22522 ? r22520 : r22529; return r22530; } double f_od(double x) { double r22531 = 1/6; double r22532 = x; double r22533 = 2; double r22534 = pow(r22532, r22533); double r22535 = r22531 * r22534; double r22536 = 1; double r22537 = 1/2; double r22538 = r22537 * r22532; double r22539 = r22536 + r22538; double r22540 = r22535 + r22539; double r22541 = 1.0000018458005175; bool r22542 = r22540 <= r22541; double r22543 = r22532 + r22532; double r22544 = exp(r22543); double r22545 = r22544 - r22536; double r22546 = exp(r22532); double r22547 = r22546 + r22536; double r22548 = r22545 / r22547; double r22549 = r22548 / r22532; double r22550 = r22542 ? r22540 : r22549; return r22550; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r22551, r22552, r22553, r22554, r22555; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22551); mpfr_init(r22552); mpfr_init_set_str(r22553, "1", 10, MPFR_RNDN); mpfr_init(r22554); mpfr_init(r22555); } double f_im(double x) { mpfr_set_d(r22551, x, MPFR_RNDN); mpfr_exp(r22552, r22551, MPFR_RNDN); ; mpfr_sub(r22554, r22552, r22553, MPFR_RNDN); mpfr_div(r22555, r22554, r22551, MPFR_RNDN); return mpfr_get_d(r22555, MPFR_RNDN); } static mpfr_t r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r22556, "1/6", 10, MPFR_RNDN); mpfr_init(r22557); mpfr_init_set_str(r22558, "2", 10, MPFR_RNDN); mpfr_init(r22559); mpfr_init(r22560); mpfr_init_set_str(r22561, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22562, "1/2", 10, MPFR_RNDN); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init_set_str(r22566, "1.0000018458005175", 10, MPFR_RNDN); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); } double f_fm(double x) { ; mpfr_set_d(r22557, x, MPFR_RNDN); ; mpfr_pow(r22559, r22557, r22558, MPFR_RNDN); mpfr_mul(r22560, r22556, r22559, MPFR_RNDN); ; ; mpfr_mul(r22563, r22562, r22557, MPFR_RNDN); mpfr_add(r22564, r22561, r22563, MPFR_RNDN); mpfr_add(r22565, r22560, r22564, MPFR_RNDN); ; mpfr_set_si(r22567, mpfr_cmp(r22565, r22566) <= 0, MPFR_RNDN); mpfr_add(r22568, r22557, r22557, MPFR_RNDN); mpfr_exp(r22569, r22568, MPFR_RNDN); mpfr_sub(r22570, r22569, r22561, MPFR_RNDN); mpfr_exp(r22571, r22557, MPFR_RNDN); mpfr_add(r22572, r22571, r22561, MPFR_RNDN); mpfr_div(r22573, r22570, r22572, MPFR_RNDN); mpfr_div(r22574, r22573, r22557, MPFR_RNDN); if (mpfr_get_si(r22567, MPFR_RNDN)) { mpfr_set(r22575, r22565, MPFR_RNDN); } else { mpfr_set(r22575, r22574, MPFR_RNDN); }; return mpfr_get_d(r22575, MPFR_RNDN); } static mpfr_t r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r22576, "1/6", 10, MPFR_RNDN); mpfr_init(r22577); mpfr_init_set_str(r22578, "2", 10, MPFR_RNDN); mpfr_init(r22579); mpfr_init(r22580); mpfr_init_set_str(r22581, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22582, "1/2", 10, MPFR_RNDN); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init_set_str(r22586, "1.0000018458005175", 10, MPFR_RNDN); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init(r22595); } double f_dm(double x) { ; mpfr_set_d(r22577, x, MPFR_RNDN); ; mpfr_pow(r22579, r22577, r22578, MPFR_RNDN); mpfr_mul(r22580, r22576, r22579, MPFR_RNDN); ; ; mpfr_mul(r22583, r22582, r22577, MPFR_RNDN); mpfr_add(r22584, r22581, r22583, MPFR_RNDN); mpfr_add(r22585, r22580, r22584, MPFR_RNDN); ; mpfr_set_si(r22587, mpfr_cmp(r22585, r22586) <= 0, MPFR_RNDN); mpfr_add(r22588, r22577, r22577, MPFR_RNDN); mpfr_exp(r22589, r22588, MPFR_RNDN); mpfr_sub(r22590, r22589, r22581, MPFR_RNDN); mpfr_exp(r22591, r22577, MPFR_RNDN); mpfr_add(r22592, r22591, r22581, MPFR_RNDN); mpfr_div(r22593, r22590, r22592, MPFR_RNDN); mpfr_div(r22594, r22593, r22577, MPFR_RNDN); if (mpfr_get_si(r22587, MPFR_RNDN)) { mpfr_set(r22595, r22585, MPFR_RNDN); } else { mpfr_set(r22595, r22594, MPFR_RNDN); }; return mpfr_get_d(r22595, MPFR_RNDN); }