#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r21547 = x; float r21548 = exp(r21547); float r21549 = 1; float r21550 = r21548 - r21549; float r21551 = r21550 / r21547; return r21551; } double f_id(double x) { double r21552 = x; double r21553 = exp(r21552); double r21554 = 1; double r21555 = r21553 - r21554; double r21556 = r21555 / r21552; return r21556; } double f_of(float x) { float r21557 = x; float r21558 = -9.72316615530548e-09; bool r21559 = r21557 <= r21558; float r21560 = r21557 + r21557; float r21561 = exp(r21560); float r21562 = 1; float r21563 = r21561 - r21562; float r21564 = exp(r21563); float r21565 = log(r21564); float r21566 = exp(r21557); float r21567 = r21566 + r21562; float r21568 = r21565 / r21567; float r21569 = r21568 / r21557; float r21570 = 1/6; float r21571 = 2; float r21572 = pow(r21557, r21571); float r21573 = r21570 * r21572; float r21574 = 1/2; float r21575 = r21574 * r21557; float r21576 = r21562 + r21575; float r21577 = r21573 + r21576; float r21578 = r21559 ? r21569 : r21577; return r21578; } double f_od(double x) { double r21579 = x; double r21580 = -9.72316615530548e-09; bool r21581 = r21579 <= r21580; double r21582 = r21579 + r21579; double r21583 = exp(r21582); double r21584 = 1; double r21585 = r21583 - r21584; double r21586 = exp(r21585); double r21587 = log(r21586); double r21588 = exp(r21579); double r21589 = r21588 + r21584; double r21590 = r21587 / r21589; double r21591 = r21590 / r21579; double r21592 = 1/6; double r21593 = 2; double r21594 = pow(r21579, r21593); double r21595 = r21592 * r21594; double r21596 = 1/2; double r21597 = r21596 * r21579; double r21598 = r21584 + r21597; double r21599 = r21595 + r21598; double r21600 = r21581 ? r21591 : r21599; return r21600; } 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 r21601, r21602, r21603, r21604, r21605; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r21601); mpfr_init(r21602); mpfr_init_set_str(r21603, "1", 10, MPFR_RNDN); mpfr_init(r21604); mpfr_init(r21605); } double f_im(double x) { mpfr_set_d(r21601, x, MPFR_RNDN); mpfr_exp(r21602, r21601, MPFR_RNDN); ; mpfr_sub(r21604, r21602, r21603, MPFR_RNDN); mpfr_div(r21605, r21604, r21601, MPFR_RNDN); return mpfr_get_d(r21605, MPFR_RNDN); } static mpfr_t r21606, r21607, r21608, r21609, r21610, r21611, r21612, r21613, r21614, r21615, r21616, r21617, r21618, r21619, r21620, r21621, r21622, r21623, r21624, r21625, r21626, r21627; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r21606); mpfr_init_set_str(r21607, "-9.72316615530548e-09", 10, MPFR_RNDN); mpfr_init(r21608); mpfr_init(r21609); mpfr_init(r21610); mpfr_init_set_str(r21611, "1", 10, MPFR_RNDN); mpfr_init(r21612); mpfr_init(r21613); mpfr_init(r21614); mpfr_init(r21615); mpfr_init(r21616); mpfr_init(r21617); mpfr_init(r21618); mpfr_init_set_str(r21619, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r21620, "2", 10, MPFR_RNDN); mpfr_init(r21621); mpfr_init(r21622); mpfr_init_set_str(r21623, "1/2", 10, MPFR_RNDN); mpfr_init(r21624); mpfr_init(r21625); mpfr_init(r21626); mpfr_init(r21627); } double f_fm(double x) { mpfr_set_d(r21606, x, MPFR_RNDN); ; mpfr_set_si(r21608, mpfr_cmp(r21606, r21607) <= 0, MPFR_RNDN); mpfr_add(r21609, r21606, r21606, MPFR_RNDN); mpfr_exp(r21610, r21609, MPFR_RNDN); ; mpfr_sub(r21612, r21610, r21611, MPFR_RNDN); mpfr_exp(r21613, r21612, MPFR_RNDN); mpfr_log(r21614, r21613, MPFR_RNDN); mpfr_exp(r21615, r21606, MPFR_RNDN); mpfr_add(r21616, r21615, r21611, MPFR_RNDN); mpfr_div(r21617, r21614, r21616, MPFR_RNDN); mpfr_div(r21618, r21617, r21606, MPFR_RNDN); ; ; mpfr_pow(r21621, r21606, r21620, MPFR_RNDN); mpfr_mul(r21622, r21619, r21621, MPFR_RNDN); ; mpfr_mul(r21624, r21623, r21606, MPFR_RNDN); mpfr_add(r21625, r21611, r21624, MPFR_RNDN); mpfr_add(r21626, r21622, r21625, MPFR_RNDN); if (mpfr_get_si(r21608, MPFR_RNDN)) { mpfr_set(r21627, r21618, MPFR_RNDN); } else { mpfr_set(r21627, r21626, MPFR_RNDN); }; return mpfr_get_d(r21627, MPFR_RNDN); } static mpfr_t r21628, r21629, r21630, r21631, r21632, r21633, r21634, r21635, r21636, r21637, r21638, r21639, r21640, r21641, r21642, r21643, r21644, r21645, r21646, r21647, r21648, r21649; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r21628); mpfr_init_set_str(r21629, "-9.72316615530548e-09", 10, MPFR_RNDN); mpfr_init(r21630); mpfr_init(r21631); mpfr_init(r21632); mpfr_init_set_str(r21633, "1", 10, MPFR_RNDN); mpfr_init(r21634); mpfr_init(r21635); mpfr_init(r21636); mpfr_init(r21637); mpfr_init(r21638); mpfr_init(r21639); mpfr_init(r21640); mpfr_init_set_str(r21641, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r21642, "2", 10, MPFR_RNDN); mpfr_init(r21643); mpfr_init(r21644); mpfr_init_set_str(r21645, "1/2", 10, MPFR_RNDN); mpfr_init(r21646); mpfr_init(r21647); mpfr_init(r21648); mpfr_init(r21649); } double f_dm(double x) { mpfr_set_d(r21628, x, MPFR_RNDN); ; mpfr_set_si(r21630, mpfr_cmp(r21628, r21629) <= 0, MPFR_RNDN); mpfr_add(r21631, r21628, r21628, MPFR_RNDN); mpfr_exp(r21632, r21631, MPFR_RNDN); ; mpfr_sub(r21634, r21632, r21633, MPFR_RNDN); mpfr_exp(r21635, r21634, MPFR_RNDN); mpfr_log(r21636, r21635, MPFR_RNDN); mpfr_exp(r21637, r21628, MPFR_RNDN); mpfr_add(r21638, r21637, r21633, MPFR_RNDN); mpfr_div(r21639, r21636, r21638, MPFR_RNDN); mpfr_div(r21640, r21639, r21628, MPFR_RNDN); ; ; mpfr_pow(r21643, r21628, r21642, MPFR_RNDN); mpfr_mul(r21644, r21641, r21643, MPFR_RNDN); ; mpfr_mul(r21646, r21645, r21628, MPFR_RNDN); mpfr_add(r21647, r21633, r21646, MPFR_RNDN); mpfr_add(r21648, r21644, r21647, MPFR_RNDN); if (mpfr_get_si(r21630, MPFR_RNDN)) { mpfr_set(r21649, r21640, MPFR_RNDN); } else { mpfr_set(r21649, r21648, MPFR_RNDN); }; return mpfr_get_d(r21649, MPFR_RNDN); }