#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22347 = x; float r22348 = r22347 * r22347; float r22349 = 1; float r22350 = r22348 + r22349; float r22351 = r22347 / r22350; return r22351; } double f_id(double x) { double r22352 = x; double r22353 = r22352 * r22352; double r22354 = 1; double r22355 = r22353 + r22354; double r22356 = r22352 / r22355; return r22356; } double f_of(float x) { float r22357 = x; float r22358 = -544269350.2683507; bool r22359 = r22357 <= r22358; float r22360 = 1; float r22361 = 5; float r22362 = pow(r22357, r22361); float r22363 = r22360 / r22362; float r22364 = r22360 / r22357; float r22365 = r22363 + r22364; float r22366 = 3; float r22367 = pow(r22357, r22366); float r22368 = r22360 / r22367; float r22369 = r22365 - r22368; float r22370 = 117645.7431062391; bool r22371 = r22357 <= r22370; float r22372 = r22357 * r22357; float r22373 = r22372 + r22360; float r22374 = r22357 / r22373; float r22375 = r22371 ? r22374 : r22369; float r22376 = r22359 ? r22369 : r22375; return r22376; } double f_od(double x) { double r22377 = x; double r22378 = -544269350.2683507; bool r22379 = r22377 <= r22378; double r22380 = 1; double r22381 = 5; double r22382 = pow(r22377, r22381); double r22383 = r22380 / r22382; double r22384 = r22380 / r22377; double r22385 = r22383 + r22384; double r22386 = 3; double r22387 = pow(r22377, r22386); double r22388 = r22380 / r22387; double r22389 = r22385 - r22388; double r22390 = 117645.7431062391; bool r22391 = r22377 <= r22390; double r22392 = r22377 * r22377; double r22393 = r22392 + r22380; double r22394 = r22377 / r22393; double r22395 = r22391 ? r22394 : r22389; double r22396 = r22379 ? r22389 : r22395; return r22396; } 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 r22397, r22398, r22399, r22400, r22401; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22397); mpfr_init(r22398); mpfr_init_set_str(r22399, "1", 10, MPFR_RNDN); mpfr_init(r22400); mpfr_init(r22401); } double f_im(double x) { mpfr_set_d(r22397, x, MPFR_RNDN); mpfr_mul(r22398, r22397, r22397, MPFR_RNDN); ; mpfr_add(r22400, r22398, r22399, MPFR_RNDN); mpfr_div(r22401, r22397, r22400, MPFR_RNDN); return mpfr_get_d(r22401, MPFR_RNDN); } static mpfr_t r22402, r22403, r22404, r22405, r22406, r22407, r22408, r22409, r22410, r22411, r22412, r22413, r22414, r22415, r22416, r22417, r22418, r22419, r22420, r22421; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22402); mpfr_init_set_str(r22403, "-544269350.2683507", 10, MPFR_RNDN); mpfr_init(r22404); mpfr_init_set_str(r22405, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22406, "5", 10, MPFR_RNDN); mpfr_init(r22407); mpfr_init(r22408); mpfr_init(r22409); mpfr_init(r22410); mpfr_init_set_str(r22411, "3", 10, MPFR_RNDN); mpfr_init(r22412); mpfr_init(r22413); mpfr_init(r22414); mpfr_init_set_str(r22415, "117645.7431062391", 10, MPFR_RNDN); mpfr_init(r22416); mpfr_init(r22417); mpfr_init(r22418); mpfr_init(r22419); mpfr_init(r22420); mpfr_init(r22421); } double f_fm(double x) { mpfr_set_d(r22402, x, MPFR_RNDN); ; mpfr_set_si(r22404, mpfr_cmp(r22402, r22403) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22407, r22402, r22406, MPFR_RNDN); mpfr_div(r22408, r22405, r22407, MPFR_RNDN); mpfr_div(r22409, r22405, r22402, MPFR_RNDN); mpfr_add(r22410, r22408, r22409, MPFR_RNDN); ; mpfr_pow(r22412, r22402, r22411, MPFR_RNDN); mpfr_div(r22413, r22405, r22412, MPFR_RNDN); mpfr_sub(r22414, r22410, r22413, MPFR_RNDN); ; mpfr_set_si(r22416, mpfr_cmp(r22402, r22415) <= 0, MPFR_RNDN); mpfr_mul(r22417, r22402, r22402, MPFR_RNDN); mpfr_add(r22418, r22417, r22405, MPFR_RNDN); mpfr_div(r22419, r22402, r22418, MPFR_RNDN); if (mpfr_get_si(r22416, MPFR_RNDN)) { mpfr_set(r22420, r22419, MPFR_RNDN); } else { mpfr_set(r22420, r22414, MPFR_RNDN); }; if (mpfr_get_si(r22404, MPFR_RNDN)) { mpfr_set(r22421, r22414, MPFR_RNDN); } else { mpfr_set(r22421, r22420, MPFR_RNDN); }; return mpfr_get_d(r22421, MPFR_RNDN); } static mpfr_t r22422, r22423, r22424, r22425, r22426, r22427, r22428, r22429, r22430, r22431, r22432, r22433, r22434, r22435, r22436, r22437, r22438, r22439, r22440, r22441; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22422); mpfr_init_set_str(r22423, "-544269350.2683507", 10, MPFR_RNDN); mpfr_init(r22424); mpfr_init_set_str(r22425, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22426, "5", 10, MPFR_RNDN); mpfr_init(r22427); mpfr_init(r22428); mpfr_init(r22429); mpfr_init(r22430); mpfr_init_set_str(r22431, "3", 10, MPFR_RNDN); mpfr_init(r22432); mpfr_init(r22433); mpfr_init(r22434); mpfr_init_set_str(r22435, "117645.7431062391", 10, MPFR_RNDN); mpfr_init(r22436); mpfr_init(r22437); mpfr_init(r22438); mpfr_init(r22439); mpfr_init(r22440); mpfr_init(r22441); } double f_dm(double x) { mpfr_set_d(r22422, x, MPFR_RNDN); ; mpfr_set_si(r22424, mpfr_cmp(r22422, r22423) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22427, r22422, r22426, MPFR_RNDN); mpfr_div(r22428, r22425, r22427, MPFR_RNDN); mpfr_div(r22429, r22425, r22422, MPFR_RNDN); mpfr_add(r22430, r22428, r22429, MPFR_RNDN); ; mpfr_pow(r22432, r22422, r22431, MPFR_RNDN); mpfr_div(r22433, r22425, r22432, MPFR_RNDN); mpfr_sub(r22434, r22430, r22433, MPFR_RNDN); ; mpfr_set_si(r22436, mpfr_cmp(r22422, r22435) <= 0, MPFR_RNDN); mpfr_mul(r22437, r22422, r22422, MPFR_RNDN); mpfr_add(r22438, r22437, r22425, MPFR_RNDN); mpfr_div(r22439, r22422, r22438, MPFR_RNDN); if (mpfr_get_si(r22436, MPFR_RNDN)) { mpfr_set(r22440, r22439, MPFR_RNDN); } else { mpfr_set(r22440, r22434, MPFR_RNDN); }; if (mpfr_get_si(r22424, MPFR_RNDN)) { mpfr_set(r22441, r22434, MPFR_RNDN); } else { mpfr_set(r22441, r22440, MPFR_RNDN); }; return mpfr_get_d(r22441, MPFR_RNDN); }