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