#include #include #include #include #include char *name = "expq2 (section 3.11)"; double f_if(float x) { float r26715 = x; float r26716 = exp(r26715); float r26717 = 1; float r26718 = r26716 - r26717; float r26719 = r26716 / r26718; return r26719; } double f_id(double x) { double r26720 = x; double r26721 = exp(r26720); double r26722 = 1; double r26723 = r26721 - r26722; double r26724 = r26721 / r26723; return r26724; } double f_of(float x) { float r26725 = x; float r26726 = exp(r26725); float r26727 = 0.0012148373248771097; bool r26728 = r26726 <= r26727; float r26729 = 1.0000000000051388; bool r26730 = r26726 <= r26729; float r26731 = !r26730; bool r26732 = r26728 || r26731; float r26733 = 1; float r26734 = -r26725; float r26735 = exp(r26734); float r26736 = r26733 - r26735; float r26737 = r26733 / r26736; float r26738 = r26733 / r26725; float r26739 = 1/12; float r26740 = r26739 * r26725; float r26741 = r26738 + r26740; float r26742 = 1/2; float r26743 = r26741 + r26742; float r26744 = r26732 ? r26737 : r26743; return r26744; } double f_od(double x) { double r26745 = x; double r26746 = exp(r26745); double r26747 = 0.0012148373248771097; bool r26748 = r26746 <= r26747; double r26749 = 1.0000000000051388; bool r26750 = r26746 <= r26749; double r26751 = !r26750; bool r26752 = r26748 || r26751; double r26753 = 1; double r26754 = -r26745; double r26755 = exp(r26754); double r26756 = r26753 - r26755; double r26757 = r26753 / r26756; double r26758 = r26753 / r26745; double r26759 = 1/12; double r26760 = r26759 * r26745; double r26761 = r26758 + r26760; double r26762 = 1/2; double r26763 = r26761 + r26762; double r26764 = r26752 ? r26757 : r26763; return r26764; } 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 r26765, r26766, r26767, r26768, r26769; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26765); mpfr_init(r26766); mpfr_init_set_str(r26767, "1", 10, MPFR_RNDN); mpfr_init(r26768); mpfr_init(r26769); } double f_im(double x) { mpfr_set_d(r26765, x, MPFR_RNDN); mpfr_exp(r26766, r26765, MPFR_RNDN); ; mpfr_sub(r26768, r26766, r26767, MPFR_RNDN); mpfr_div(r26769, r26766, r26768, MPFR_RNDN); return mpfr_get_d(r26769, MPFR_RNDN); } static mpfr_t r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778, r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26770); mpfr_init(r26771); mpfr_init_set_str(r26772, "0.0012148373248771097", 10, MPFR_RNDN); mpfr_init(r26773); mpfr_init_set_str(r26774, "1.0000000000051388", 10, MPFR_RNDN); mpfr_init(r26775); mpfr_init(r26776); mpfr_init(r26777); mpfr_init_set_str(r26778, "1", 10, MPFR_RNDN); mpfr_init(r26779); mpfr_init(r26780); mpfr_init(r26781); mpfr_init(r26782); mpfr_init(r26783); mpfr_init_set_str(r26784, "1/12", 10, MPFR_RNDN); mpfr_init(r26785); mpfr_init(r26786); mpfr_init_set_str(r26787, "1/2", 10, MPFR_RNDN); mpfr_init(r26788); mpfr_init(r26789); } double f_fm(double x) { mpfr_set_d(r26770, x, MPFR_RNDN); mpfr_exp(r26771, r26770, MPFR_RNDN); ; mpfr_set_si(r26773, mpfr_cmp(r26771, r26772) <= 0, MPFR_RNDN); ; mpfr_set_si(r26775, mpfr_cmp(r26771, r26774) <= 0, MPFR_RNDN); mpfr_set_si(r26776, !mpfr_get_si(r26775, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26777, mpfr_get_si(r26773, MPFR_RNDN) || mpfr_get_si(r26776, MPFR_RNDN), MPFR_RNDN); ; mpfr_neg(r26779, r26770, MPFR_RNDN); mpfr_exp(r26780, r26779, MPFR_RNDN); mpfr_sub(r26781, r26778, r26780, MPFR_RNDN); mpfr_div(r26782, r26778, r26781, MPFR_RNDN); mpfr_div(r26783, r26778, r26770, MPFR_RNDN); ; mpfr_mul(r26785, r26784, r26770, MPFR_RNDN); mpfr_add(r26786, r26783, r26785, MPFR_RNDN); ; mpfr_add(r26788, r26786, r26787, MPFR_RNDN); if (mpfr_get_si(r26777, MPFR_RNDN)) { mpfr_set(r26789, r26782, MPFR_RNDN); } else { mpfr_set(r26789, r26788, MPFR_RNDN); }; return mpfr_get_d(r26789, MPFR_RNDN); } static mpfr_t r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26790); mpfr_init(r26791); mpfr_init_set_str(r26792, "0.0012148373248771097", 10, MPFR_RNDN); mpfr_init(r26793); mpfr_init_set_str(r26794, "1.0000000000051388", 10, MPFR_RNDN); mpfr_init(r26795); mpfr_init(r26796); mpfr_init(r26797); mpfr_init_set_str(r26798, "1", 10, MPFR_RNDN); mpfr_init(r26799); mpfr_init(r26800); mpfr_init(r26801); mpfr_init(r26802); mpfr_init(r26803); mpfr_init_set_str(r26804, "1/12", 10, MPFR_RNDN); mpfr_init(r26805); mpfr_init(r26806); mpfr_init_set_str(r26807, "1/2", 10, MPFR_RNDN); mpfr_init(r26808); mpfr_init(r26809); } double f_dm(double x) { mpfr_set_d(r26790, x, MPFR_RNDN); mpfr_exp(r26791, r26790, MPFR_RNDN); ; mpfr_set_si(r26793, mpfr_cmp(r26791, r26792) <= 0, MPFR_RNDN); ; mpfr_set_si(r26795, mpfr_cmp(r26791, r26794) <= 0, MPFR_RNDN); mpfr_set_si(r26796, !mpfr_get_si(r26795, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26797, mpfr_get_si(r26793, MPFR_RNDN) || mpfr_get_si(r26796, MPFR_RNDN), MPFR_RNDN); ; mpfr_neg(r26799, r26790, MPFR_RNDN); mpfr_exp(r26800, r26799, MPFR_RNDN); mpfr_sub(r26801, r26798, r26800, MPFR_RNDN); mpfr_div(r26802, r26798, r26801, MPFR_RNDN); mpfr_div(r26803, r26798, r26790, MPFR_RNDN); ; mpfr_mul(r26805, r26804, r26790, MPFR_RNDN); mpfr_add(r26806, r26803, r26805, MPFR_RNDN); ; mpfr_add(r26808, r26806, r26807, MPFR_RNDN); if (mpfr_get_si(r26797, MPFR_RNDN)) { mpfr_set(r26809, r26802, MPFR_RNDN); } else { mpfr_set(r26809, r26808, MPFR_RNDN); }; return mpfr_get_d(r26809, MPFR_RNDN); }