#include #include #include #include #include char *name = "NMSE problem 3.3.1"; double f_if(float x) { float r15397 = 1.0f; float r15398 = x; float r15399 = r15398 + r15397; float r15400 = r15397 / r15399; float r15401 = r15397 / r15398; float r15402 = r15400 - r15401; return r15402; } double f_id(double x) { double r15403 = 1.0; double r15404 = x; double r15405 = r15404 + r15403; double r15406 = r15403 / r15405; double r15407 = r15403 / r15404; double r15408 = r15406 - r15407; return r15408; } double f_of(float x) { float r15409 = x; float r15410 = -840391.875f; bool r15411 = r15409 <= r15410; float r15412 = 1.0f; float r15413 = r15412 / r15409; float r15414 = r15409 * r15409; float r15415 = r15413 / r15414; float r15416 = r15413 / r15409; float r15417 = r15415 - r15416; float r15418 = 15032.4169921875f; bool r15419 = r15409 <= r15418; float r15420 = r15412 + r15409; float r15421 = r15409 - r15420; float r15422 = r15409 + r15412; float r15423 = r15422 * r15409; float r15424 = r15421 / r15423; float r15425 = r15419 ? r15424 : r15417; float r15426 = r15411 ? r15417 : r15425; return r15426; } double f_od(double x) { double r15427 = x; double r15428 = -840391.875; bool r15429 = r15427 <= r15428; double r15430 = 1.0; double r15431 = r15430 / r15427; double r15432 = r15427 * r15427; double r15433 = r15431 / r15432; double r15434 = r15431 / r15427; double r15435 = r15433 - r15434; double r15436 = 15032.4169921875; bool r15437 = r15427 <= r15436; double r15438 = r15430 + r15427; double r15439 = r15427 - r15438; double r15440 = r15427 + r15430; double r15441 = r15440 * r15427; double r15442 = r15439 / r15441; double r15443 = r15437 ? r15442 : r15435; double r15444 = r15429 ? r15435 : r15443; return r15444; } 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 r15445, r15446, r15447, r15448, r15449, r15450; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r15445, "1", 10, MPFR_RNDN); mpfr_init(r15446); mpfr_init(r15447); mpfr_init(r15448); mpfr_init(r15449); mpfr_init(r15450); } double f_im(double x) { ; mpfr_set_d(r15446, x, MPFR_RNDN); mpfr_add(r15447, r15446, r15445, MPFR_RNDN); mpfr_div(r15448, r15445, r15447, MPFR_RNDN); mpfr_div(r15449, r15445, r15446, MPFR_RNDN); mpfr_sub(r15450, r15448, r15449, MPFR_RNDN); return mpfr_get_d(r15450, MPFR_RNDN); } static mpfr_t r15451, r15452, r15453, r15454, r15455, r15456, r15457, r15458, r15459, r15460, r15461, r15462, r15463, r15464, r15465, r15466, r15467, r15468; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15451); mpfr_init_set_str(r15452, "-840391.9f0", 10, MPFR_RNDN); mpfr_init(r15453); mpfr_init_set_str(r15454, "1", 10, MPFR_RNDN); mpfr_init(r15455); mpfr_init(r15456); mpfr_init(r15457); mpfr_init(r15458); mpfr_init(r15459); mpfr_init_set_str(r15460, "15032.417f0", 10, MPFR_RNDN); mpfr_init(r15461); mpfr_init(r15462); mpfr_init(r15463); mpfr_init(r15464); mpfr_init(r15465); mpfr_init(r15466); mpfr_init(r15467); mpfr_init(r15468); } double f_fm(double x) { mpfr_set_d(r15451, x, MPFR_RNDN); ; mpfr_set_si(r15453, mpfr_cmp(r15451, r15452) <= 0, MPFR_RNDN); ; mpfr_div(r15455, r15454, r15451, MPFR_RNDN); mpfr_mul(r15456, r15451, r15451, MPFR_RNDN); mpfr_div(r15457, r15455, r15456, MPFR_RNDN); mpfr_div(r15458, r15455, r15451, MPFR_RNDN); mpfr_sub(r15459, r15457, r15458, MPFR_RNDN); ; mpfr_set_si(r15461, mpfr_cmp(r15451, r15460) <= 0, MPFR_RNDN); mpfr_add(r15462, r15454, r15451, MPFR_RNDN); mpfr_sub(r15463, r15451, r15462, MPFR_RNDN); mpfr_add(r15464, r15451, r15454, MPFR_RNDN); mpfr_mul(r15465, r15464, r15451, MPFR_RNDN); mpfr_div(r15466, r15463, r15465, MPFR_RNDN); if (mpfr_get_si(r15461, MPFR_RNDN)) { mpfr_set(r15467, r15466, MPFR_RNDN); } else { mpfr_set(r15467, r15459, MPFR_RNDN); }; if (mpfr_get_si(r15453, MPFR_RNDN)) { mpfr_set(r15468, r15459, MPFR_RNDN); } else { mpfr_set(r15468, r15467, MPFR_RNDN); }; return mpfr_get_d(r15468, MPFR_RNDN); } static mpfr_t r15469, r15470, r15471, r15472, r15473, r15474, r15475, r15476, r15477, r15478, r15479, r15480, r15481, r15482, r15483, r15484, r15485, r15486; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15469); mpfr_init_set_str(r15470, "-840391.9f0", 10, MPFR_RNDN); mpfr_init(r15471); mpfr_init_set_str(r15472, "1", 10, MPFR_RNDN); mpfr_init(r15473); mpfr_init(r15474); mpfr_init(r15475); mpfr_init(r15476); mpfr_init(r15477); mpfr_init_set_str(r15478, "15032.417f0", 10, MPFR_RNDN); mpfr_init(r15479); mpfr_init(r15480); mpfr_init(r15481); mpfr_init(r15482); mpfr_init(r15483); mpfr_init(r15484); mpfr_init(r15485); mpfr_init(r15486); } double f_dm(double x) { mpfr_set_d(r15469, x, MPFR_RNDN); ; mpfr_set_si(r15471, mpfr_cmp(r15469, r15470) <= 0, MPFR_RNDN); ; mpfr_div(r15473, r15472, r15469, MPFR_RNDN); mpfr_mul(r15474, r15469, r15469, MPFR_RNDN); mpfr_div(r15475, r15473, r15474, MPFR_RNDN); mpfr_div(r15476, r15473, r15469, MPFR_RNDN); mpfr_sub(r15477, r15475, r15476, MPFR_RNDN); ; mpfr_set_si(r15479, mpfr_cmp(r15469, r15478) <= 0, MPFR_RNDN); mpfr_add(r15480, r15472, r15469, MPFR_RNDN); mpfr_sub(r15481, r15469, r15480, MPFR_RNDN); mpfr_add(r15482, r15469, r15472, MPFR_RNDN); mpfr_mul(r15483, r15482, r15469, MPFR_RNDN); mpfr_div(r15484, r15481, r15483, MPFR_RNDN); if (mpfr_get_si(r15479, MPFR_RNDN)) { mpfr_set(r15485, r15484, MPFR_RNDN); } else { mpfr_set(r15485, r15477, MPFR_RNDN); }; if (mpfr_get_si(r15471, MPFR_RNDN)) { mpfr_set(r15486, r15477, MPFR_RNDN); } else { mpfr_set(r15486, r15485, MPFR_RNDN); }; return mpfr_get_d(r15486, MPFR_RNDN); }