#include #include #include #include #include char *name = "expax (section 3.5)"; double f_if(float a, float x) { float r24681 = a; float r24682 = x; float r24683 = r24681 * r24682; float r24684 = exp(r24683); float r24685 = 1; float r24686 = r24684 - r24685; return r24686; } double f_id(double a, double x) { double r24687 = a; double r24688 = x; double r24689 = r24687 * r24688; double r24690 = exp(r24689); double r24691 = 1; double r24692 = r24690 - r24691; return r24692; } double f_of(float a, float x) { float r24693 = a; float r24694 = x; float r24695 = r24693 * r24694; float r24696 = exp(r24695); float r24697 = 1; float r24698 = r24696 - r24697; float r24699 = sqrt(r24698); float r24700 = r24699 * r24699; float r24701 = 0.0005425738380245641; bool r24702 = r24700 <= r24701; float r24703 = r24695 * r24695; float r24704 = 1/6; float r24705 = r24695 * r24704; float r24706 = 1/2; float r24707 = r24705 + r24706; float r24708 = r24703 * r24707; float r24709 = r24708 + r24695; float r24710 = cbrt(r24698); float r24711 = r24710 * r24710; float r24712 = r24711 * r24710; float r24713 = r24702 ? r24709 : r24712; return r24713; } double f_od(double a, double x) { double r24714 = a; double r24715 = x; double r24716 = r24714 * r24715; double r24717 = exp(r24716); double r24718 = 1; double r24719 = r24717 - r24718; double r24720 = sqrt(r24719); double r24721 = r24720 * r24720; double r24722 = 0.0005425738380245641; bool r24723 = r24721 <= r24722; double r24724 = r24716 * r24716; double r24725 = 1/6; double r24726 = r24716 * r24725; double r24727 = 1/2; double r24728 = r24726 + r24727; double r24729 = r24724 * r24728; double r24730 = r24729 + r24716; double r24731 = cbrt(r24719); double r24732 = r24731 * r24731; double r24733 = r24732 * r24731; double r24734 = r24723 ? r24730 : r24733; return r24734; } 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 r24735, r24736, r24737, r24738, r24739, r24740; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24735); mpfr_init(r24736); mpfr_init(r24737); mpfr_init(r24738); mpfr_init_set_str(r24739, "1", 10, MPFR_RNDN); mpfr_init(r24740); } double f_im(double a, double x) { mpfr_set_d(r24735, a, MPFR_RNDN); mpfr_set_d(r24736, x, MPFR_RNDN); mpfr_mul(r24737, r24735, r24736, MPFR_RNDN); mpfr_exp(r24738, r24737, MPFR_RNDN); ; mpfr_sub(r24740, r24738, r24739, MPFR_RNDN); return mpfr_get_d(r24740, MPFR_RNDN); } static mpfr_t r24741, r24742, r24743, r24744, r24745, r24746, r24747, r24748, r24749, r24750, r24751, r24752, r24753, r24754, r24755, r24756, r24757, r24758, r24759, r24760, r24761; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24741); mpfr_init(r24742); mpfr_init(r24743); mpfr_init(r24744); mpfr_init_set_str(r24745, "1", 10, MPFR_RNDN); mpfr_init(r24746); mpfr_init(r24747); mpfr_init(r24748); mpfr_init_set_str(r24749, "0.0005425738380245641", 10, MPFR_RNDN); mpfr_init(r24750); mpfr_init(r24751); mpfr_init_set_str(r24752, "1/6", 10, MPFR_RNDN); mpfr_init(r24753); mpfr_init_set_str(r24754, "1/2", 10, MPFR_RNDN); mpfr_init(r24755); mpfr_init(r24756); mpfr_init(r24757); mpfr_init(r24758); mpfr_init(r24759); mpfr_init(r24760); mpfr_init(r24761); } double f_fm(double a, double x) { mpfr_set_d(r24741, a, MPFR_RNDN); mpfr_set_d(r24742, x, MPFR_RNDN); mpfr_mul(r24743, r24741, r24742, MPFR_RNDN); mpfr_exp(r24744, r24743, MPFR_RNDN); ; mpfr_sub(r24746, r24744, r24745, MPFR_RNDN); mpfr_sqrt(r24747, r24746, MPFR_RNDN); mpfr_mul(r24748, r24747, r24747, MPFR_RNDN); ; mpfr_set_si(r24750, mpfr_cmp(r24748, r24749) <= 0, MPFR_RNDN); mpfr_mul(r24751, r24743, r24743, MPFR_RNDN); ; mpfr_mul(r24753, r24743, r24752, MPFR_RNDN); ; mpfr_add(r24755, r24753, r24754, MPFR_RNDN); mpfr_mul(r24756, r24751, r24755, MPFR_RNDN); mpfr_add(r24757, r24756, r24743, MPFR_RNDN); mpfr_cbrt(r24758, r24746, MPFR_RNDN); mpfr_mul(r24759, r24758, r24758, MPFR_RNDN); mpfr_mul(r24760, r24759, r24758, MPFR_RNDN); if (mpfr_get_si(r24750, MPFR_RNDN)) { mpfr_set(r24761, r24757, MPFR_RNDN); } else { mpfr_set(r24761, r24760, MPFR_RNDN); }; return mpfr_get_d(r24761, MPFR_RNDN); } static mpfr_t r24762, r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773, r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24762); mpfr_init(r24763); mpfr_init(r24764); mpfr_init(r24765); mpfr_init_set_str(r24766, "1", 10, MPFR_RNDN); mpfr_init(r24767); mpfr_init(r24768); mpfr_init(r24769); mpfr_init_set_str(r24770, "0.0005425738380245641", 10, MPFR_RNDN); mpfr_init(r24771); mpfr_init(r24772); mpfr_init_set_str(r24773, "1/6", 10, MPFR_RNDN); mpfr_init(r24774); mpfr_init_set_str(r24775, "1/2", 10, MPFR_RNDN); mpfr_init(r24776); mpfr_init(r24777); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init(r24781); mpfr_init(r24782); } double f_dm(double a, double x) { mpfr_set_d(r24762, a, MPFR_RNDN); mpfr_set_d(r24763, x, MPFR_RNDN); mpfr_mul(r24764, r24762, r24763, MPFR_RNDN); mpfr_exp(r24765, r24764, MPFR_RNDN); ; mpfr_sub(r24767, r24765, r24766, MPFR_RNDN); mpfr_sqrt(r24768, r24767, MPFR_RNDN); mpfr_mul(r24769, r24768, r24768, MPFR_RNDN); ; mpfr_set_si(r24771, mpfr_cmp(r24769, r24770) <= 0, MPFR_RNDN); mpfr_mul(r24772, r24764, r24764, MPFR_RNDN); ; mpfr_mul(r24774, r24764, r24773, MPFR_RNDN); ; mpfr_add(r24776, r24774, r24775, MPFR_RNDN); mpfr_mul(r24777, r24772, r24776, MPFR_RNDN); mpfr_add(r24778, r24777, r24764, MPFR_RNDN); mpfr_cbrt(r24779, r24767, MPFR_RNDN); mpfr_mul(r24780, r24779, r24779, MPFR_RNDN); mpfr_mul(r24781, r24780, r24779, MPFR_RNDN); if (mpfr_get_si(r24771, MPFR_RNDN)) { mpfr_set(r24782, r24778, MPFR_RNDN); } else { mpfr_set(r24782, r24781, MPFR_RNDN); }; return mpfr_get_d(r24782, MPFR_RNDN); }