#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r22561 = x; float r22562 = exp(r22561); float r22563 = 1; float r22564 = r22562 - r22563; float r22565 = r22564 / r22561; return r22565; } double f_id(double x) { double r22566 = x; double r22567 = exp(r22566); double r22568 = 1; double r22569 = r22567 - r22568; double r22570 = r22569 / r22566; return r22570; } double f_of(float x) { float r22571 = x; float r22572 = -0.00012975755436313898; bool r22573 = r22571 <= r22572; float r22574 = r22571 + r22571; float r22575 = exp(r22574); float r22576 = 1; float r22577 = r22575 - r22576; float r22578 = exp(r22571); float r22579 = r22578 + r22576; float r22580 = r22577 / r22579; float r22581 = r22580 / r22571; float r22582 = 1/6; float r22583 = r22582 * r22571; float r22584 = 1/2; float r22585 = r22583 + r22584; float r22586 = r22571 * r22571; float r22587 = r22585 * r22586; float r22588 = r22587 + r22571; float r22589 = r22588 / r22571; float r22590 = r22573 ? r22581 : r22589; return r22590; } double f_od(double x) { double r22591 = x; double r22592 = -0.00012975755436313898; bool r22593 = r22591 <= r22592; double r22594 = r22591 + r22591; double r22595 = exp(r22594); double r22596 = 1; double r22597 = r22595 - r22596; double r22598 = exp(r22591); double r22599 = r22598 + r22596; double r22600 = r22597 / r22599; double r22601 = r22600 / r22591; double r22602 = 1/6; double r22603 = r22602 * r22591; double r22604 = 1/2; double r22605 = r22603 + r22604; double r22606 = r22591 * r22591; double r22607 = r22605 * r22606; double r22608 = r22607 + r22591; double r22609 = r22608 / r22591; double r22610 = r22593 ? r22601 : r22609; return r22610; } 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 r22611, r22612, r22613, r22614, r22615; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22611); mpfr_init(r22612); mpfr_init_set_str(r22613, "1", 10, MPFR_RNDN); mpfr_init(r22614); mpfr_init(r22615); } double f_im(double x) { mpfr_set_d(r22611, x, MPFR_RNDN); mpfr_exp(r22612, r22611, MPFR_RNDN); ; mpfr_sub(r22614, r22612, r22613, MPFR_RNDN); mpfr_div(r22615, r22614, r22611, MPFR_RNDN); return mpfr_get_d(r22615, MPFR_RNDN); } static mpfr_t r22616, r22617, r22618, r22619, r22620, r22621, r22622, r22623, r22624, r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22616); mpfr_init_set_str(r22617, "-0.00012975755436313898", 10, MPFR_RNDN); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); mpfr_init_set_str(r22621, "1", 10, MPFR_RNDN); mpfr_init(r22622); mpfr_init(r22623); mpfr_init(r22624); mpfr_init(r22625); mpfr_init(r22626); mpfr_init_set_str(r22627, "1/6", 10, MPFR_RNDN); mpfr_init(r22628); mpfr_init_set_str(r22629, "1/2", 10, MPFR_RNDN); mpfr_init(r22630); mpfr_init(r22631); mpfr_init(r22632); mpfr_init(r22633); mpfr_init(r22634); mpfr_init(r22635); } double f_fm(double x) { mpfr_set_d(r22616, x, MPFR_RNDN); ; mpfr_set_si(r22618, mpfr_cmp(r22616, r22617) <= 0, MPFR_RNDN); mpfr_add(r22619, r22616, r22616, MPFR_RNDN); mpfr_exp(r22620, r22619, MPFR_RNDN); ; mpfr_sub(r22622, r22620, r22621, MPFR_RNDN); mpfr_exp(r22623, r22616, MPFR_RNDN); mpfr_add(r22624, r22623, r22621, MPFR_RNDN); mpfr_div(r22625, r22622, r22624, MPFR_RNDN); mpfr_div(r22626, r22625, r22616, MPFR_RNDN); ; mpfr_mul(r22628, r22627, r22616, MPFR_RNDN); ; mpfr_add(r22630, r22628, r22629, MPFR_RNDN); mpfr_mul(r22631, r22616, r22616, MPFR_RNDN); mpfr_mul(r22632, r22630, r22631, MPFR_RNDN); mpfr_add(r22633, r22632, r22616, MPFR_RNDN); mpfr_div(r22634, r22633, r22616, MPFR_RNDN); if (mpfr_get_si(r22618, MPFR_RNDN)) { mpfr_set(r22635, r22626, MPFR_RNDN); } else { mpfr_set(r22635, r22634, MPFR_RNDN); }; return mpfr_get_d(r22635, MPFR_RNDN); } static mpfr_t r22636, r22637, r22638, r22639, r22640, r22641, r22642, r22643, r22644, r22645, r22646, r22647, r22648, r22649, r22650, r22651, r22652, r22653, r22654, r22655; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22636); mpfr_init_set_str(r22637, "-0.00012975755436313898", 10, MPFR_RNDN); mpfr_init(r22638); mpfr_init(r22639); mpfr_init(r22640); mpfr_init_set_str(r22641, "1", 10, MPFR_RNDN); mpfr_init(r22642); mpfr_init(r22643); mpfr_init(r22644); mpfr_init(r22645); mpfr_init(r22646); mpfr_init_set_str(r22647, "1/6", 10, MPFR_RNDN); mpfr_init(r22648); mpfr_init_set_str(r22649, "1/2", 10, MPFR_RNDN); mpfr_init(r22650); mpfr_init(r22651); mpfr_init(r22652); mpfr_init(r22653); mpfr_init(r22654); mpfr_init(r22655); } double f_dm(double x) { mpfr_set_d(r22636, x, MPFR_RNDN); ; mpfr_set_si(r22638, mpfr_cmp(r22636, r22637) <= 0, MPFR_RNDN); mpfr_add(r22639, r22636, r22636, MPFR_RNDN); mpfr_exp(r22640, r22639, MPFR_RNDN); ; mpfr_sub(r22642, r22640, r22641, MPFR_RNDN); mpfr_exp(r22643, r22636, MPFR_RNDN); mpfr_add(r22644, r22643, r22641, MPFR_RNDN); mpfr_div(r22645, r22642, r22644, MPFR_RNDN); mpfr_div(r22646, r22645, r22636, MPFR_RNDN); ; mpfr_mul(r22648, r22647, r22636, MPFR_RNDN); ; mpfr_add(r22650, r22648, r22649, MPFR_RNDN); mpfr_mul(r22651, r22636, r22636, MPFR_RNDN); mpfr_mul(r22652, r22650, r22651, MPFR_RNDN); mpfr_add(r22653, r22652, r22636, MPFR_RNDN); mpfr_div(r22654, r22653, r22636, MPFR_RNDN); if (mpfr_get_si(r22638, MPFR_RNDN)) { mpfr_set(r22655, r22646, MPFR_RNDN); } else { mpfr_set(r22655, r22654, MPFR_RNDN); }; return mpfr_get_d(r22655, MPFR_RNDN); }