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