#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r23543 = x; float r23544 = exp(r23543); float r23545 = 1; float r23546 = r23544 - r23545; float r23547 = r23546 / r23543; return r23547; } double f_id(double x) { double r23548 = x; double r23549 = exp(r23548); double r23550 = 1; double r23551 = r23549 - r23550; double r23552 = r23551 / r23548; return r23552; } double f_of(float x) { float r23553 = x; float r23554 = -0.00013140414695441723; bool r23555 = r23553 <= r23554; float r23556 = r23553 + r23553; float r23557 = exp(r23556); float r23558 = 1; float r23559 = r23557 - r23558; float r23560 = exp(r23553); float r23561 = r23560 + r23558; float r23562 = r23559 / r23561; float r23563 = r23562 / r23553; float r23564 = 1/2; float r23565 = 2; float r23566 = pow(r23553, r23565); float r23567 = r23564 * r23566; float r23568 = 1/6; float r23569 = 3; float r23570 = pow(r23553, r23569); float r23571 = r23568 * r23570; float r23572 = r23571 + r23553; float r23573 = r23567 + r23572; float r23574 = r23573 / r23553; float r23575 = r23555 ? r23563 : r23574; return r23575; } double f_od(double x) { double r23576 = x; double r23577 = -0.00013140414695441723; bool r23578 = r23576 <= r23577; double r23579 = r23576 + r23576; double r23580 = exp(r23579); double r23581 = 1; double r23582 = r23580 - r23581; double r23583 = exp(r23576); double r23584 = r23583 + r23581; double r23585 = r23582 / r23584; double r23586 = r23585 / r23576; double r23587 = 1/2; double r23588 = 2; double r23589 = pow(r23576, r23588); double r23590 = r23587 * r23589; double r23591 = 1/6; double r23592 = 3; double r23593 = pow(r23576, r23592); double r23594 = r23591 * r23593; double r23595 = r23594 + r23576; double r23596 = r23590 + r23595; double r23597 = r23596 / r23576; double r23598 = r23578 ? r23586 : r23597; return r23598; } 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 r23599, r23600, r23601, r23602, r23603; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23599); mpfr_init(r23600); mpfr_init_set_str(r23601, "1", 10, MPFR_RNDN); mpfr_init(r23602); mpfr_init(r23603); } double f_im(double x) { mpfr_set_d(r23599, x, MPFR_RNDN); mpfr_exp(r23600, r23599, MPFR_RNDN); ; mpfr_sub(r23602, r23600, r23601, MPFR_RNDN); mpfr_div(r23603, r23602, r23599, MPFR_RNDN); return mpfr_get_d(r23603, MPFR_RNDN); } static mpfr_t r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23604); mpfr_init_set_str(r23605, "-0.00013140414695441723", 10, MPFR_RNDN); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init_set_str(r23609, "1", 10, MPFR_RNDN); mpfr_init(r23610); mpfr_init(r23611); mpfr_init(r23612); mpfr_init(r23613); mpfr_init(r23614); mpfr_init_set_str(r23615, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r23616, "2", 10, MPFR_RNDN); mpfr_init(r23617); mpfr_init(r23618); mpfr_init_set_str(r23619, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r23620, "3", 10, MPFR_RNDN); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); mpfr_init(r23626); } double f_fm(double x) { mpfr_set_d(r23604, x, MPFR_RNDN); ; mpfr_set_si(r23606, mpfr_cmp(r23604, r23605) <= 0, MPFR_RNDN); mpfr_add(r23607, r23604, r23604, MPFR_RNDN); mpfr_exp(r23608, r23607, MPFR_RNDN); ; mpfr_sub(r23610, r23608, r23609, MPFR_RNDN); mpfr_exp(r23611, r23604, MPFR_RNDN); mpfr_add(r23612, r23611, r23609, MPFR_RNDN); mpfr_div(r23613, r23610, r23612, MPFR_RNDN); mpfr_div(r23614, r23613, r23604, MPFR_RNDN); ; ; mpfr_pow(r23617, r23604, r23616, MPFR_RNDN); mpfr_mul(r23618, r23615, r23617, MPFR_RNDN); ; ; mpfr_pow(r23621, r23604, r23620, MPFR_RNDN); mpfr_mul(r23622, r23619, r23621, MPFR_RNDN); mpfr_add(r23623, r23622, r23604, MPFR_RNDN); mpfr_add(r23624, r23618, r23623, MPFR_RNDN); mpfr_div(r23625, r23624, r23604, MPFR_RNDN); if (mpfr_get_si(r23606, MPFR_RNDN)) { mpfr_set(r23626, r23614, MPFR_RNDN); } else { mpfr_set(r23626, r23625, MPFR_RNDN); }; return mpfr_get_d(r23626, MPFR_RNDN); } static mpfr_t r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648, r23649; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23627); mpfr_init_set_str(r23628, "-0.00013140414695441723", 10, MPFR_RNDN); mpfr_init(r23629); mpfr_init(r23630); mpfr_init(r23631); mpfr_init_set_str(r23632, "1", 10, MPFR_RNDN); mpfr_init(r23633); mpfr_init(r23634); mpfr_init(r23635); mpfr_init(r23636); mpfr_init(r23637); mpfr_init_set_str(r23638, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r23639, "2", 10, MPFR_RNDN); mpfr_init(r23640); mpfr_init(r23641); mpfr_init_set_str(r23642, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r23643, "3", 10, MPFR_RNDN); mpfr_init(r23644); mpfr_init(r23645); mpfr_init(r23646); mpfr_init(r23647); mpfr_init(r23648); mpfr_init(r23649); } double f_dm(double x) { mpfr_set_d(r23627, x, MPFR_RNDN); ; mpfr_set_si(r23629, mpfr_cmp(r23627, r23628) <= 0, MPFR_RNDN); mpfr_add(r23630, r23627, r23627, MPFR_RNDN); mpfr_exp(r23631, r23630, MPFR_RNDN); ; mpfr_sub(r23633, r23631, r23632, MPFR_RNDN); mpfr_exp(r23634, r23627, MPFR_RNDN); mpfr_add(r23635, r23634, r23632, MPFR_RNDN); mpfr_div(r23636, r23633, r23635, MPFR_RNDN); mpfr_div(r23637, r23636, r23627, MPFR_RNDN); ; ; mpfr_pow(r23640, r23627, r23639, MPFR_RNDN); mpfr_mul(r23641, r23638, r23640, MPFR_RNDN); ; ; mpfr_pow(r23644, r23627, r23643, MPFR_RNDN); mpfr_mul(r23645, r23642, r23644, MPFR_RNDN); mpfr_add(r23646, r23645, r23627, MPFR_RNDN); mpfr_add(r23647, r23641, r23646, MPFR_RNDN); mpfr_div(r23648, r23647, r23627, MPFR_RNDN); if (mpfr_get_si(r23629, MPFR_RNDN)) { mpfr_set(r23649, r23637, MPFR_RNDN); } else { mpfr_set(r23649, r23648, MPFR_RNDN); }; return mpfr_get_d(r23649, MPFR_RNDN); }