#include #include #include #include #include char *name = "expq2 (section 3.11)"; double f_if(float x) { float r24358 = x; float r24359 = exp(r24358); float r24360 = 1; float r24361 = r24359 - r24360; float r24362 = r24359 / r24361; return r24362; } double f_id(double x) { double r24363 = x; double r24364 = exp(r24363); double r24365 = 1; double r24366 = r24364 - r24365; double r24367 = r24364 / r24366; return r24367; } double f_of(float x) { float r24368 = x; float r24369 = exp(r24368); float r24370 = 0.9997312861954971; bool r24371 = r24369 <= r24370; float r24372 = 1; float r24373 = expm1(r24368); float r24374 = cbrt(r24373); float r24375 = r24374 * r24374; float r24376 = r24372 / r24375; float r24377 = r24369 / r24374; float r24378 = r24376 * r24377; float r24379 = r24372 / r24368; float r24380 = 1/12; float r24381 = 1/2; float r24382 = fma(r24380, r24368, r24381); float r24383 = 1/3; float r24384 = r24383 + r24383; float r24385 = pow(r24382, r24384); float r24386 = cbrt(r24382); float r24387 = r24385 * r24386; float r24388 = r24379 + r24387; float r24389 = r24371 ? r24378 : r24388; return r24389; } double f_od(double x) { double r24390 = x; double r24391 = exp(r24390); double r24392 = 0.9997312861954971; bool r24393 = r24391 <= r24392; double r24394 = 1; double r24395 = expm1(r24390); double r24396 = cbrt(r24395); double r24397 = r24396 * r24396; double r24398 = r24394 / r24397; double r24399 = r24391 / r24396; double r24400 = r24398 * r24399; double r24401 = r24394 / r24390; double r24402 = 1/12; double r24403 = 1/2; double r24404 = fma(r24402, r24390, r24403); double r24405 = 1/3; double r24406 = r24405 + r24405; double r24407 = pow(r24404, r24406); double r24408 = cbrt(r24404); double r24409 = r24407 * r24408; double r24410 = r24401 + r24409; double r24411 = r24393 ? r24400 : r24410; return r24411; } 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 r24412, r24413, r24414, r24415, r24416; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24412); mpfr_init(r24413); mpfr_init_set_str(r24414, "1", 10, MPFR_RNDN); mpfr_init(r24415); mpfr_init(r24416); } double f_im(double x) { mpfr_set_d(r24412, x, MPFR_RNDN); mpfr_exp(r24413, r24412, MPFR_RNDN); ; mpfr_sub(r24415, r24413, r24414, MPFR_RNDN); mpfr_div(r24416, r24413, r24415, MPFR_RNDN); return mpfr_get_d(r24416, MPFR_RNDN); } static mpfr_t r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24417); mpfr_init(r24418); mpfr_init_set_str(r24419, "0.9997312861954971", 10, MPFR_RNDN); mpfr_init(r24420); mpfr_init_set_str(r24421, "1", 10, MPFR_RNDN); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init_set_str(r24429, "1/12", 10, MPFR_RNDN); mpfr_init_set_str(r24430, "1/2", 10, MPFR_RNDN); mpfr_init(r24431); mpfr_init_set_str(r24432, "1/3", 10, MPFR_RNDN); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); } double f_fm(double x) { mpfr_set_d(r24417, x, MPFR_RNDN); mpfr_exp(r24418, r24417, MPFR_RNDN); ; mpfr_set_si(r24420, mpfr_cmp(r24418, r24419) <= 0, MPFR_RNDN); ; mpfr_expm1(r24422, r24417, MPFR_RNDN); mpfr_cbrt(r24423, r24422, MPFR_RNDN); mpfr_mul(r24424, r24423, r24423, MPFR_RNDN); mpfr_div(r24425, r24421, r24424, MPFR_RNDN); mpfr_div(r24426, r24418, r24423, MPFR_RNDN); mpfr_mul(r24427, r24425, r24426, MPFR_RNDN); mpfr_div(r24428, r24421, r24417, MPFR_RNDN); ; ; mpfr_fma(r24431, r24429, r24417, r24430, MPFR_RNDN); ; mpfr_add(r24433, r24432, r24432, MPFR_RNDN); mpfr_pow(r24434, r24431, r24433, MPFR_RNDN); mpfr_cbrt(r24435, r24431, MPFR_RNDN); mpfr_mul(r24436, r24434, r24435, MPFR_RNDN); mpfr_add(r24437, r24428, r24436, MPFR_RNDN); if (mpfr_get_si(r24420, MPFR_RNDN)) { mpfr_set(r24438, r24427, MPFR_RNDN); } else { mpfr_set(r24438, r24437, MPFR_RNDN); }; return mpfr_get_d(r24438, MPFR_RNDN); } static mpfr_t r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24439); mpfr_init(r24440); mpfr_init_set_str(r24441, "0.9997312861954971", 10, MPFR_RNDN); mpfr_init(r24442); mpfr_init_set_str(r24443, "1", 10, MPFR_RNDN); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); mpfr_init_set_str(r24451, "1/12", 10, MPFR_RNDN); mpfr_init_set_str(r24452, "1/2", 10, MPFR_RNDN); mpfr_init(r24453); mpfr_init_set_str(r24454, "1/3", 10, MPFR_RNDN); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); } double f_dm(double x) { mpfr_set_d(r24439, x, MPFR_RNDN); mpfr_exp(r24440, r24439, MPFR_RNDN); ; mpfr_set_si(r24442, mpfr_cmp(r24440, r24441) <= 0, MPFR_RNDN); ; mpfr_expm1(r24444, r24439, MPFR_RNDN); mpfr_cbrt(r24445, r24444, MPFR_RNDN); mpfr_mul(r24446, r24445, r24445, MPFR_RNDN); mpfr_div(r24447, r24443, r24446, MPFR_RNDN); mpfr_div(r24448, r24440, r24445, MPFR_RNDN); mpfr_mul(r24449, r24447, r24448, MPFR_RNDN); mpfr_div(r24450, r24443, r24439, MPFR_RNDN); ; ; mpfr_fma(r24453, r24451, r24439, r24452, MPFR_RNDN); ; mpfr_add(r24455, r24454, r24454, MPFR_RNDN); mpfr_pow(r24456, r24453, r24455, MPFR_RNDN); mpfr_cbrt(r24457, r24453, MPFR_RNDN); mpfr_mul(r24458, r24456, r24457, MPFR_RNDN); mpfr_add(r24459, r24450, r24458, MPFR_RNDN); if (mpfr_get_si(r24442, MPFR_RNDN)) { mpfr_set(r24460, r24449, MPFR_RNDN); } else { mpfr_set(r24460, r24459, MPFR_RNDN); }; return mpfr_get_d(r24460, MPFR_RNDN); }