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