#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r25002 = x; float r25003 = r25002 * r25002; float r25004 = 1; float r25005 = r25003 + r25004; float r25006 = r25002 / r25005; return r25006; } double f_id(double x) { double r25007 = x; double r25008 = r25007 * r25007; double r25009 = 1; double r25010 = r25008 + r25009; double r25011 = r25007 / r25010; return r25011; } double f_of(float x) { float r25012 = x; float r25013 = -1631567455.901017; bool r25014 = r25012 <= r25013; float r25015 = 164813693.3988662; bool r25016 = r25012 <= r25015; float r25017 = !r25016; bool r25018 = r25014 || r25017; float r25019 = 1; float r25020 = 5; float r25021 = pow(r25012, r25020); float r25022 = r25019 / r25021; float r25023 = r25019 / r25012; float r25024 = r25022 + r25023; float r25025 = 3; float r25026 = pow(r25012, r25025); float r25027 = r25019 / r25026; float r25028 = r25024 - r25027; float r25029 = fma(r25012, r25012, r25019); float r25030 = r25019 / r25029; float r25031 = r25030 * r25012; float r25032 = r25018 ? r25028 : r25031; return r25032; } double f_od(double x) { double r25033 = x; double r25034 = -1631567455.901017; bool r25035 = r25033 <= r25034; double r25036 = 164813693.3988662; bool r25037 = r25033 <= r25036; double r25038 = !r25037; bool r25039 = r25035 || r25038; double r25040 = 1; double r25041 = 5; double r25042 = pow(r25033, r25041); double r25043 = r25040 / r25042; double r25044 = r25040 / r25033; double r25045 = r25043 + r25044; double r25046 = 3; double r25047 = pow(r25033, r25046); double r25048 = r25040 / r25047; double r25049 = r25045 - r25048; double r25050 = fma(r25033, r25033, r25040); double r25051 = r25040 / r25050; double r25052 = r25051 * r25033; double r25053 = r25039 ? r25049 : r25052; return r25053; } 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 r25054, r25055, r25056, r25057, r25058; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25054); mpfr_init(r25055); mpfr_init_set_str(r25056, "1", 10, MPFR_RNDN); mpfr_init(r25057); mpfr_init(r25058); } double f_im(double x) { mpfr_set_d(r25054, x, MPFR_RNDN); mpfr_mul(r25055, r25054, r25054, MPFR_RNDN); ; mpfr_add(r25057, r25055, r25056, MPFR_RNDN); mpfr_div(r25058, r25054, r25057, MPFR_RNDN); return mpfr_get_d(r25058, MPFR_RNDN); } static mpfr_t r25059, r25060, r25061, r25062, r25063, r25064, r25065, r25066, r25067, r25068, r25069, r25070, r25071, r25072, r25073, r25074, r25075, r25076, r25077, r25078, r25079; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25059); mpfr_init_set_str(r25060, "-1631567455.901017", 10, MPFR_RNDN); mpfr_init(r25061); mpfr_init_set_str(r25062, "164813693.3988662", 10, MPFR_RNDN); mpfr_init(r25063); mpfr_init(r25064); mpfr_init(r25065); mpfr_init_set_str(r25066, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25067, "5", 10, MPFR_RNDN); mpfr_init(r25068); mpfr_init(r25069); mpfr_init(r25070); mpfr_init(r25071); mpfr_init_set_str(r25072, "3", 10, MPFR_RNDN); mpfr_init(r25073); mpfr_init(r25074); mpfr_init(r25075); mpfr_init(r25076); mpfr_init(r25077); mpfr_init(r25078); mpfr_init(r25079); } double f_fm(double x) { mpfr_set_d(r25059, x, MPFR_RNDN); ; mpfr_set_si(r25061, mpfr_cmp(r25059, r25060) <= 0, MPFR_RNDN); ; mpfr_set_si(r25063, mpfr_cmp(r25059, r25062) <= 0, MPFR_RNDN); mpfr_set_si(r25064, !mpfr_get_si(r25063, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25065, mpfr_get_si(r25061, MPFR_RNDN) || mpfr_get_si(r25064, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25068, r25059, r25067, MPFR_RNDN); mpfr_div(r25069, r25066, r25068, MPFR_RNDN); mpfr_div(r25070, r25066, r25059, MPFR_RNDN); mpfr_add(r25071, r25069, r25070, MPFR_RNDN); ; mpfr_pow(r25073, r25059, r25072, MPFR_RNDN); mpfr_div(r25074, r25066, r25073, MPFR_RNDN); mpfr_sub(r25075, r25071, r25074, MPFR_RNDN); mpfr_fma(r25076, r25059, r25059, r25066, MPFR_RNDN); mpfr_div(r25077, r25066, r25076, MPFR_RNDN); mpfr_mul(r25078, r25077, r25059, MPFR_RNDN); if (mpfr_get_si(r25065, MPFR_RNDN)) { mpfr_set(r25079, r25075, MPFR_RNDN); } else { mpfr_set(r25079, r25078, MPFR_RNDN); }; return mpfr_get_d(r25079, MPFR_RNDN); } static mpfr_t r25080, r25081, r25082, r25083, r25084, r25085, r25086, r25087, r25088, r25089, r25090, r25091, r25092, r25093, r25094, r25095, r25096, r25097, r25098, r25099, r25100; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25080); mpfr_init_set_str(r25081, "-1631567455.901017", 10, MPFR_RNDN); mpfr_init(r25082); mpfr_init_set_str(r25083, "164813693.3988662", 10, MPFR_RNDN); mpfr_init(r25084); mpfr_init(r25085); mpfr_init(r25086); mpfr_init_set_str(r25087, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25088, "5", 10, MPFR_RNDN); mpfr_init(r25089); mpfr_init(r25090); mpfr_init(r25091); mpfr_init(r25092); mpfr_init_set_str(r25093, "3", 10, MPFR_RNDN); mpfr_init(r25094); mpfr_init(r25095); mpfr_init(r25096); mpfr_init(r25097); mpfr_init(r25098); mpfr_init(r25099); mpfr_init(r25100); } double f_dm(double x) { mpfr_set_d(r25080, x, MPFR_RNDN); ; mpfr_set_si(r25082, mpfr_cmp(r25080, r25081) <= 0, MPFR_RNDN); ; mpfr_set_si(r25084, mpfr_cmp(r25080, r25083) <= 0, MPFR_RNDN); mpfr_set_si(r25085, !mpfr_get_si(r25084, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25086, mpfr_get_si(r25082, MPFR_RNDN) || mpfr_get_si(r25085, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25089, r25080, r25088, MPFR_RNDN); mpfr_div(r25090, r25087, r25089, MPFR_RNDN); mpfr_div(r25091, r25087, r25080, MPFR_RNDN); mpfr_add(r25092, r25090, r25091, MPFR_RNDN); ; mpfr_pow(r25094, r25080, r25093, MPFR_RNDN); mpfr_div(r25095, r25087, r25094, MPFR_RNDN); mpfr_sub(r25096, r25092, r25095, MPFR_RNDN); mpfr_fma(r25097, r25080, r25080, r25087, MPFR_RNDN); mpfr_div(r25098, r25087, r25097, MPFR_RNDN); mpfr_mul(r25099, r25098, r25080, MPFR_RNDN); if (mpfr_get_si(r25086, MPFR_RNDN)) { mpfr_set(r25100, r25096, MPFR_RNDN); } else { mpfr_set(r25100, r25099, MPFR_RNDN); }; return mpfr_get_d(r25100, MPFR_RNDN); }