#include #include #include #include #include char *name = "expq2 (section 3.11)"; double f_if(float x) { float r24259 = x; float r24260 = exp(r24259); float r24261 = 1; float r24262 = r24260 - r24261; float r24263 = r24260 / r24262; return r24263; } double f_id(double x) { double r24264 = x; double r24265 = exp(r24264); double r24266 = 1; double r24267 = r24265 - r24266; double r24268 = r24265 / r24267; return r24268; } double f_of(float x) { float r24269 = x; float r24270 = exp(r24269); float r24271 = 0.9991732492398515; bool r24272 = r24270 <= r24271; float r24273 = r24270 * r24270; float r24274 = 1; float r24275 = r24273 - r24274; float r24276 = 3; float r24277 = pow(r24275, r24276); float r24278 = cbrt(r24277); float r24279 = r24270 + r24274; float r24280 = r24278 / r24279; float r24281 = r24270 / r24280; float r24282 = 1/2; float r24283 = r24274 / r24269; float r24284 = 1/12; float r24285 = r24284 * r24269; float r24286 = r24283 + r24285; float r24287 = r24282 + r24286; float r24288 = r24272 ? r24281 : r24287; return r24288; } double f_od(double x) { double r24289 = x; double r24290 = exp(r24289); double r24291 = 0.9991732492398515; bool r24292 = r24290 <= r24291; double r24293 = r24290 * r24290; double r24294 = 1; double r24295 = r24293 - r24294; double r24296 = 3; double r24297 = pow(r24295, r24296); double r24298 = cbrt(r24297); double r24299 = r24290 + r24294; double r24300 = r24298 / r24299; double r24301 = r24290 / r24300; double r24302 = 1/2; double r24303 = r24294 / r24289; double r24304 = 1/12; double r24305 = r24304 * r24289; double r24306 = r24303 + r24305; double r24307 = r24302 + r24306; double r24308 = r24292 ? r24301 : r24307; return r24308; } 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 r24309, r24310, r24311, r24312, r24313; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24309); mpfr_init(r24310); mpfr_init_set_str(r24311, "1", 10, MPFR_RNDN); mpfr_init(r24312); mpfr_init(r24313); } double f_im(double x) { mpfr_set_d(r24309, x, MPFR_RNDN); mpfr_exp(r24310, r24309, MPFR_RNDN); ; mpfr_sub(r24312, r24310, r24311, MPFR_RNDN); mpfr_div(r24313, r24310, r24312, MPFR_RNDN); return mpfr_get_d(r24313, MPFR_RNDN); } static mpfr_t r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332, r24333; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24314); mpfr_init(r24315); mpfr_init_set_str(r24316, "0.9991732492398515", 10, MPFR_RNDN); mpfr_init(r24317); mpfr_init(r24318); mpfr_init_set_str(r24319, "1", 10, MPFR_RNDN); mpfr_init(r24320); mpfr_init_set_str(r24321, "3", 10, MPFR_RNDN); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init_set_str(r24327, "1/2", 10, MPFR_RNDN); mpfr_init(r24328); mpfr_init_set_str(r24329, "1/12", 10, MPFR_RNDN); mpfr_init(r24330); mpfr_init(r24331); mpfr_init(r24332); mpfr_init(r24333); } double f_fm(double x) { mpfr_set_d(r24314, x, MPFR_RNDN); mpfr_exp(r24315, r24314, MPFR_RNDN); ; mpfr_set_si(r24317, mpfr_cmp(r24315, r24316) <= 0, MPFR_RNDN); mpfr_mul(r24318, r24315, r24315, MPFR_RNDN); ; mpfr_sub(r24320, r24318, r24319, MPFR_RNDN); ; mpfr_pow(r24322, r24320, r24321, MPFR_RNDN); mpfr_cbrt(r24323, r24322, MPFR_RNDN); mpfr_add(r24324, r24315, r24319, MPFR_RNDN); mpfr_div(r24325, r24323, r24324, MPFR_RNDN); mpfr_div(r24326, r24315, r24325, MPFR_RNDN); ; mpfr_div(r24328, r24319, r24314, MPFR_RNDN); ; mpfr_mul(r24330, r24329, r24314, MPFR_RNDN); mpfr_add(r24331, r24328, r24330, MPFR_RNDN); mpfr_add(r24332, r24327, r24331, MPFR_RNDN); if (mpfr_get_si(r24317, MPFR_RNDN)) { mpfr_set(r24333, r24326, MPFR_RNDN); } else { mpfr_set(r24333, r24332, MPFR_RNDN); }; return mpfr_get_d(r24333, MPFR_RNDN); } static mpfr_t r24334, r24335, r24336, r24337, r24338, r24339, r24340, r24341, r24342, r24343, r24344, r24345, r24346, r24347, r24348, r24349, r24350, r24351, r24352, r24353; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24334); mpfr_init(r24335); mpfr_init_set_str(r24336, "0.9991732492398515", 10, MPFR_RNDN); mpfr_init(r24337); mpfr_init(r24338); mpfr_init_set_str(r24339, "1", 10, MPFR_RNDN); mpfr_init(r24340); mpfr_init_set_str(r24341, "3", 10, MPFR_RNDN); mpfr_init(r24342); mpfr_init(r24343); mpfr_init(r24344); mpfr_init(r24345); mpfr_init(r24346); mpfr_init_set_str(r24347, "1/2", 10, MPFR_RNDN); mpfr_init(r24348); mpfr_init_set_str(r24349, "1/12", 10, MPFR_RNDN); mpfr_init(r24350); mpfr_init(r24351); mpfr_init(r24352); mpfr_init(r24353); } double f_dm(double x) { mpfr_set_d(r24334, x, MPFR_RNDN); mpfr_exp(r24335, r24334, MPFR_RNDN); ; mpfr_set_si(r24337, mpfr_cmp(r24335, r24336) <= 0, MPFR_RNDN); mpfr_mul(r24338, r24335, r24335, MPFR_RNDN); ; mpfr_sub(r24340, r24338, r24339, MPFR_RNDN); ; mpfr_pow(r24342, r24340, r24341, MPFR_RNDN); mpfr_cbrt(r24343, r24342, MPFR_RNDN); mpfr_add(r24344, r24335, r24339, MPFR_RNDN); mpfr_div(r24345, r24343, r24344, MPFR_RNDN); mpfr_div(r24346, r24335, r24345, MPFR_RNDN); ; mpfr_div(r24348, r24339, r24334, MPFR_RNDN); ; mpfr_mul(r24350, r24349, r24334, MPFR_RNDN); mpfr_add(r24351, r24348, r24350, MPFR_RNDN); mpfr_add(r24352, r24347, r24351, MPFR_RNDN); if (mpfr_get_si(r24337, MPFR_RNDN)) { mpfr_set(r24353, r24346, MPFR_RNDN); } else { mpfr_set(r24353, r24352, MPFR_RNDN); }; return mpfr_get_d(r24353, MPFR_RNDN); }