#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r22540 = x; float r22541 = exp(r22540); float r22542 = 1; float r22543 = r22541 - r22542; float r22544 = r22543 / r22540; return r22544; } double f_id(double x) { double r22545 = x; double r22546 = exp(r22545); double r22547 = 1; double r22548 = r22546 - r22547; double r22549 = r22548 / r22545; return r22549; } double f_of(float x) { float r22550 = x; float r22551 = exp(r22550); float r22552 = 1; float r22553 = r22551 - r22552; float r22554 = r22553 / r22550; float r22555 = -0.0; bool r22556 = r22554 <= r22555; float r22557 = 1/6; float r22558 = 2; float r22559 = pow(r22550, r22558); float r22560 = r22557 * r22559; float r22561 = 1/2; float r22562 = r22561 * r22550; float r22563 = r22552 + r22562; float r22564 = r22560 + r22563; float r22565 = 0.985336851063437; bool r22566 = r22554 <= r22565; float r22567 = exp(r22553); float r22568 = log(r22567); float r22569 = r22568 / r22550; float r22570 = r22566 ? r22569 : r22564; float r22571 = r22556 ? r22564 : r22570; return r22571; } double f_od(double x) { double r22572 = x; double r22573 = exp(r22572); double r22574 = 1; double r22575 = r22573 - r22574; double r22576 = r22575 / r22572; double r22577 = -0.0; bool r22578 = r22576 <= r22577; double r22579 = 1/6; double r22580 = 2; double r22581 = pow(r22572, r22580); double r22582 = r22579 * r22581; double r22583 = 1/2; double r22584 = r22583 * r22572; double r22585 = r22574 + r22584; double r22586 = r22582 + r22585; double r22587 = 0.985336851063437; bool r22588 = r22576 <= r22587; double r22589 = exp(r22575); double r22590 = log(r22589); double r22591 = r22590 / r22572; double r22592 = r22588 ? r22591 : r22586; double r22593 = r22578 ? r22586 : r22592; return r22593; } 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 r22594, r22595, r22596, r22597, r22598; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22594); mpfr_init(r22595); mpfr_init_set_str(r22596, "1", 10, MPFR_RNDN); mpfr_init(r22597); mpfr_init(r22598); } double f_im(double x) { mpfr_set_d(r22594, x, MPFR_RNDN); mpfr_exp(r22595, r22594, MPFR_RNDN); ; mpfr_sub(r22597, r22595, r22596, MPFR_RNDN); mpfr_div(r22598, r22597, r22594, MPFR_RNDN); return mpfr_get_d(r22598, MPFR_RNDN); } static mpfr_t r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608, r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619, r22620; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22599); mpfr_init(r22600); mpfr_init_set_str(r22601, "1", 10, MPFR_RNDN); mpfr_init(r22602); mpfr_init(r22603); mpfr_init_set_str(r22604, "-0.0", 10, MPFR_RNDN); mpfr_init(r22605); mpfr_init_set_str(r22606, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r22607, "2", 10, MPFR_RNDN); mpfr_init(r22608); mpfr_init(r22609); mpfr_init_set_str(r22610, "1/2", 10, MPFR_RNDN); mpfr_init(r22611); mpfr_init(r22612); mpfr_init(r22613); mpfr_init_set_str(r22614, "0.985336851063437", 10, MPFR_RNDN); mpfr_init(r22615); mpfr_init(r22616); mpfr_init(r22617); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); } double f_fm(double x) { mpfr_set_d(r22599, x, MPFR_RNDN); mpfr_exp(r22600, r22599, MPFR_RNDN); ; mpfr_sub(r22602, r22600, r22601, MPFR_RNDN); mpfr_div(r22603, r22602, r22599, MPFR_RNDN); ; mpfr_set_si(r22605, mpfr_cmp(r22603, r22604) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22608, r22599, r22607, MPFR_RNDN); mpfr_mul(r22609, r22606, r22608, MPFR_RNDN); ; mpfr_mul(r22611, r22610, r22599, MPFR_RNDN); mpfr_add(r22612, r22601, r22611, MPFR_RNDN); mpfr_add(r22613, r22609, r22612, MPFR_RNDN); ; mpfr_set_si(r22615, mpfr_cmp(r22603, r22614) <= 0, MPFR_RNDN); mpfr_exp(r22616, r22602, MPFR_RNDN); mpfr_log(r22617, r22616, MPFR_RNDN); mpfr_div(r22618, r22617, r22599, MPFR_RNDN); if (mpfr_get_si(r22615, MPFR_RNDN)) { mpfr_set(r22619, r22618, MPFR_RNDN); } else { mpfr_set(r22619, r22613, MPFR_RNDN); }; if (mpfr_get_si(r22605, MPFR_RNDN)) { mpfr_set(r22620, r22613, MPFR_RNDN); } else { mpfr_set(r22620, r22619, MPFR_RNDN); }; return mpfr_get_d(r22620, MPFR_RNDN); } static mpfr_t r22621, r22622, r22623, r22624, r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635, r22636, r22637, r22638, r22639, r22640, r22641, r22642; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22621); mpfr_init(r22622); mpfr_init_set_str(r22623, "1", 10, MPFR_RNDN); mpfr_init(r22624); mpfr_init(r22625); mpfr_init_set_str(r22626, "-0.0", 10, MPFR_RNDN); mpfr_init(r22627); mpfr_init_set_str(r22628, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r22629, "2", 10, MPFR_RNDN); mpfr_init(r22630); mpfr_init(r22631); mpfr_init_set_str(r22632, "1/2", 10, MPFR_RNDN); mpfr_init(r22633); mpfr_init(r22634); mpfr_init(r22635); mpfr_init_set_str(r22636, "0.985336851063437", 10, MPFR_RNDN); mpfr_init(r22637); mpfr_init(r22638); mpfr_init(r22639); mpfr_init(r22640); mpfr_init(r22641); mpfr_init(r22642); } double f_dm(double x) { mpfr_set_d(r22621, x, MPFR_RNDN); mpfr_exp(r22622, r22621, MPFR_RNDN); ; mpfr_sub(r22624, r22622, r22623, MPFR_RNDN); mpfr_div(r22625, r22624, r22621, MPFR_RNDN); ; mpfr_set_si(r22627, mpfr_cmp(r22625, r22626) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22630, r22621, r22629, MPFR_RNDN); mpfr_mul(r22631, r22628, r22630, MPFR_RNDN); ; mpfr_mul(r22633, r22632, r22621, MPFR_RNDN); mpfr_add(r22634, r22623, r22633, MPFR_RNDN); mpfr_add(r22635, r22631, r22634, MPFR_RNDN); ; mpfr_set_si(r22637, mpfr_cmp(r22625, r22636) <= 0, MPFR_RNDN); mpfr_exp(r22638, r22624, MPFR_RNDN); mpfr_log(r22639, r22638, MPFR_RNDN); mpfr_div(r22640, r22639, r22621, MPFR_RNDN); if (mpfr_get_si(r22637, MPFR_RNDN)) { mpfr_set(r22641, r22640, MPFR_RNDN); } else { mpfr_set(r22641, r22635, MPFR_RNDN); }; if (mpfr_get_si(r22627, MPFR_RNDN)) { mpfr_set(r22642, r22635, MPFR_RNDN); } else { mpfr_set(r22642, r22641, MPFR_RNDN); }; return mpfr_get_d(r22642, MPFR_RNDN); }