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