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