#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22363 = x; float r22364 = r22363 * r22363; float r22365 = 1; float r22366 = r22364 + r22365; float r22367 = r22363 / r22366; return r22367; } double f_id(double x) { double r22368 = x; double r22369 = r22368 * r22368; double r22370 = 1; double r22371 = r22369 + r22370; double r22372 = r22368 / r22371; return r22372; } double f_of(float x) { float r22373 = x; float r22374 = -19409101.03762743; bool r22375 = r22373 <= r22374; float r22376 = 1; float r22377 = 5; float r22378 = pow(r22373, r22377); float r22379 = r22376 / r22378; float r22380 = r22376 / r22373; float r22381 = r22379 + r22380; float r22382 = 3; float r22383 = pow(r22373, r22382); float r22384 = r22376 / r22383; float r22385 = r22381 - r22384; float r22386 = 144528077.66840905; bool r22387 = r22373 <= r22386; float r22388 = fma(r22373, r22373, r22376); float r22389 = r22373 / r22388; float r22390 = log1p(r22389); float r22391 = expm1(r22390); float r22392 = r22387 ? r22391 : r22385; float r22393 = r22375 ? r22385 : r22392; return r22393; } double f_od(double x) { double r22394 = x; double r22395 = -19409101.03762743; bool r22396 = r22394 <= r22395; double r22397 = 1; double r22398 = 5; double r22399 = pow(r22394, r22398); double r22400 = r22397 / r22399; double r22401 = r22397 / r22394; double r22402 = r22400 + r22401; double r22403 = 3; double r22404 = pow(r22394, r22403); double r22405 = r22397 / r22404; double r22406 = r22402 - r22405; double r22407 = 144528077.66840905; bool r22408 = r22394 <= r22407; double r22409 = fma(r22394, r22394, r22397); double r22410 = r22394 / r22409; double r22411 = log1p(r22410); double r22412 = expm1(r22411); double r22413 = r22408 ? r22412 : r22406; double r22414 = r22396 ? r22406 : r22413; return r22414; } 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 r22415, r22416, r22417, r22418, r22419; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22415); mpfr_init(r22416); mpfr_init_set_str(r22417, "1", 10, MPFR_RNDN); mpfr_init(r22418); mpfr_init(r22419); } double f_im(double x) { mpfr_set_d(r22415, x, MPFR_RNDN); mpfr_mul(r22416, r22415, r22415, MPFR_RNDN); ; mpfr_add(r22418, r22416, r22417, MPFR_RNDN); mpfr_div(r22419, r22415, r22418, MPFR_RNDN); return mpfr_get_d(r22419, MPFR_RNDN); } static mpfr_t r22420, r22421, r22422, r22423, r22424, r22425, r22426, r22427, r22428, r22429, r22430, r22431, r22432, r22433, r22434, r22435, r22436, r22437, r22438, r22439, r22440; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22420); mpfr_init_set_str(r22421, "-19409101.03762743", 10, MPFR_RNDN); mpfr_init(r22422); mpfr_init_set_str(r22423, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22424, "5", 10, MPFR_RNDN); mpfr_init(r22425); mpfr_init(r22426); mpfr_init(r22427); mpfr_init(r22428); mpfr_init_set_str(r22429, "3", 10, MPFR_RNDN); mpfr_init(r22430); mpfr_init(r22431); mpfr_init(r22432); mpfr_init_set_str(r22433, "144528077.66840905", 10, MPFR_RNDN); mpfr_init(r22434); mpfr_init(r22435); mpfr_init(r22436); mpfr_init(r22437); mpfr_init(r22438); mpfr_init(r22439); mpfr_init(r22440); } double f_fm(double x) { mpfr_set_d(r22420, x, MPFR_RNDN); ; mpfr_set_si(r22422, mpfr_cmp(r22420, r22421) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22425, r22420, r22424, MPFR_RNDN); mpfr_div(r22426, r22423, r22425, MPFR_RNDN); mpfr_div(r22427, r22423, r22420, MPFR_RNDN); mpfr_add(r22428, r22426, r22427, MPFR_RNDN); ; mpfr_pow(r22430, r22420, r22429, MPFR_RNDN); mpfr_div(r22431, r22423, r22430, MPFR_RNDN); mpfr_sub(r22432, r22428, r22431, MPFR_RNDN); ; mpfr_set_si(r22434, mpfr_cmp(r22420, r22433) <= 0, MPFR_RNDN); mpfr_fma(r22435, r22420, r22420, r22423, MPFR_RNDN); mpfr_div(r22436, r22420, r22435, MPFR_RNDN); mpfr_log1p(r22437, r22436, MPFR_RNDN); mpfr_expm1(r22438, r22437, MPFR_RNDN); if (mpfr_get_si(r22434, MPFR_RNDN)) { mpfr_set(r22439, r22438, MPFR_RNDN); } else { mpfr_set(r22439, r22432, MPFR_RNDN); }; if (mpfr_get_si(r22422, MPFR_RNDN)) { mpfr_set(r22440, r22432, MPFR_RNDN); } else { mpfr_set(r22440, r22439, MPFR_RNDN); }; return mpfr_get_d(r22440, MPFR_RNDN); } static mpfr_t r22441, r22442, r22443, r22444, r22445, r22446, r22447, r22448, r22449, r22450, r22451, r22452, r22453, r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22441); mpfr_init_set_str(r22442, "-19409101.03762743", 10, MPFR_RNDN); mpfr_init(r22443); mpfr_init_set_str(r22444, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22445, "5", 10, MPFR_RNDN); mpfr_init(r22446); mpfr_init(r22447); mpfr_init(r22448); mpfr_init(r22449); mpfr_init_set_str(r22450, "3", 10, MPFR_RNDN); mpfr_init(r22451); mpfr_init(r22452); mpfr_init(r22453); mpfr_init_set_str(r22454, "144528077.66840905", 10, MPFR_RNDN); mpfr_init(r22455); mpfr_init(r22456); mpfr_init(r22457); mpfr_init(r22458); mpfr_init(r22459); mpfr_init(r22460); mpfr_init(r22461); } double f_dm(double x) { mpfr_set_d(r22441, x, MPFR_RNDN); ; mpfr_set_si(r22443, mpfr_cmp(r22441, r22442) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22446, r22441, r22445, MPFR_RNDN); mpfr_div(r22447, r22444, r22446, MPFR_RNDN); mpfr_div(r22448, r22444, r22441, MPFR_RNDN); mpfr_add(r22449, r22447, r22448, MPFR_RNDN); ; mpfr_pow(r22451, r22441, r22450, MPFR_RNDN); mpfr_div(r22452, r22444, r22451, MPFR_RNDN); mpfr_sub(r22453, r22449, r22452, MPFR_RNDN); ; mpfr_set_si(r22455, mpfr_cmp(r22441, r22454) <= 0, MPFR_RNDN); mpfr_fma(r22456, r22441, r22441, r22444, MPFR_RNDN); mpfr_div(r22457, r22441, r22456, MPFR_RNDN); mpfr_log1p(r22458, r22457, MPFR_RNDN); mpfr_expm1(r22459, r22458, MPFR_RNDN); if (mpfr_get_si(r22455, MPFR_RNDN)) { mpfr_set(r22460, r22459, MPFR_RNDN); } else { mpfr_set(r22460, r22453, MPFR_RNDN); }; if (mpfr_get_si(r22443, MPFR_RNDN)) { mpfr_set(r22461, r22453, MPFR_RNDN); } else { mpfr_set(r22461, r22460, MPFR_RNDN); }; return mpfr_get_d(r22461, MPFR_RNDN); }