#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r19201 = x; float r19202 = r19201 * r19201; float r19203 = 1.0f; float r19204 = r19202 + r19203; float r19205 = r19201 / r19204; return r19205; } double f_id(double x) { double r19206 = x; double r19207 = r19206 * r19206; double r19208 = 1.0; double r19209 = r19207 + r19208; double r19210 = r19206 / r19209; return r19210; } double f_of(float x) { float r19211 = x; float r19212 = -4.248855072932943e+22f; bool r19213 = r19211 <= r19212; float r19214 = 1.0f; float r19215 = 5.0f; float r19216 = pow(r19211, r19215); float r19217 = r19214 / r19216; float r19218 = r19214 / r19211; float r19219 = r19217 + r19218; float r19220 = r19211 * (r19211 * r19211); float r19221 = r19214 / r19220; float r19222 = r19219 - r19221; float r19223 = 95035.66200477978f; bool r19224 = r19211 <= r19223; float r19225 = r19211 * r19211; float r19226 = r19225 + r19214; float r19227 = r19211 / r19226; float r19228 = r19224 ? r19227 : r19222; float r19229 = r19213 ? r19222 : r19228; return r19229; } double f_od(double x) { double r19230 = x; double r19231 = -4.248855072932943e+22; bool r19232 = r19230 <= r19231; double r19233 = 1.0; double r19234 = 5.0; double r19235 = pow(r19230, r19234); double r19236 = r19233 / r19235; double r19237 = r19233 / r19230; double r19238 = r19236 + r19237; double r19239 = r19230 * (r19230 * r19230); double r19240 = r19233 / r19239; double r19241 = r19238 - r19240; double r19242 = 95035.66200477978; bool r19243 = r19230 <= r19242; double r19244 = r19230 * r19230; double r19245 = r19244 + r19233; double r19246 = r19230 / r19245; double r19247 = r19243 ? r19246 : r19241; double r19248 = r19232 ? r19241 : r19247; return r19248; } 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 r19249, r19250, r19251, r19252, r19253; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r19249); mpfr_init(r19250); mpfr_init_set_str(r19251, "1", 10, MPFR_RNDN); mpfr_init(r19252); mpfr_init(r19253); } double f_im(double x) { mpfr_set_d(r19249, x, MPFR_RNDN); mpfr_sqr(r19250, r19249, MPFR_RNDN); ; mpfr_add(r19252, r19250, r19251, MPFR_RNDN); mpfr_div(r19253, r19249, r19252, MPFR_RNDN); return mpfr_get_d(r19253, MPFR_RNDN); } static mpfr_t r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265, r19266, r19267, r19268, r19269, r19270, r19271, r19272; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r19254); mpfr_init_set_str(r19255, "-4.248855072932943e+22", 10, MPFR_RNDN); mpfr_init(r19256); mpfr_init_set_str(r19257, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19258, "5", 10, MPFR_RNDN); mpfr_init(r19259); mpfr_init(r19260); mpfr_init(r19261); mpfr_init(r19262); mpfr_init(r19263); mpfr_init(r19264); mpfr_init(r19265); mpfr_init_set_str(r19266, "95035.66200477978", 10, MPFR_RNDN); mpfr_init(r19267); mpfr_init(r19268); mpfr_init(r19269); mpfr_init(r19270); mpfr_init(r19271); mpfr_init(r19272); } double f_fm(double x) { mpfr_set_d(r19254, x, MPFR_RNDN); ; mpfr_set_si(r19256, mpfr_cmp(r19254, r19255) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19259, r19254, r19258, MPFR_RNDN); mpfr_div(r19260, r19257, r19259, MPFR_RNDN); mpfr_div(r19261, r19257, r19254, MPFR_RNDN); mpfr_add(r19262, r19260, r19261, MPFR_RNDN); mpfr_mul(r19263, r19254, r19254, MPFR_RNDN); mpfr_mul(r19263, r19263, r19254, MPFR_RNDN); mpfr_div(r19264, r19257, r19263, MPFR_RNDN); mpfr_sub(r19265, r19262, r19264, MPFR_RNDN); ; mpfr_set_si(r19267, mpfr_cmp(r19254, r19266) <= 0, MPFR_RNDN); mpfr_sqr(r19268, r19254, MPFR_RNDN); mpfr_add(r19269, r19268, r19257, MPFR_RNDN); mpfr_div(r19270, r19254, r19269, MPFR_RNDN); if (mpfr_get_si(r19267, MPFR_RNDN)) { mpfr_set(r19271, r19270, MPFR_RNDN); } else { mpfr_set(r19271, r19265, MPFR_RNDN); }; if (mpfr_get_si(r19256, MPFR_RNDN)) { mpfr_set(r19272, r19265, MPFR_RNDN); } else { mpfr_set(r19272, r19271, MPFR_RNDN); }; return mpfr_get_d(r19272, MPFR_RNDN); } static mpfr_t r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284, r19285, r19286, r19287, r19288, r19289, r19290, r19291; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r19273); mpfr_init_set_str(r19274, "-4.248855072932943e+22", 10, MPFR_RNDN); mpfr_init(r19275); mpfr_init_set_str(r19276, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19277, "5", 10, MPFR_RNDN); mpfr_init(r19278); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init(r19282); mpfr_init(r19283); mpfr_init(r19284); mpfr_init_set_str(r19285, "95035.66200477978", 10, MPFR_RNDN); mpfr_init(r19286); mpfr_init(r19287); mpfr_init(r19288); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); } double f_dm(double x) { mpfr_set_d(r19273, x, MPFR_RNDN); ; mpfr_set_si(r19275, mpfr_cmp(r19273, r19274) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19278, r19273, r19277, MPFR_RNDN); mpfr_div(r19279, r19276, r19278, MPFR_RNDN); mpfr_div(r19280, r19276, r19273, MPFR_RNDN); mpfr_add(r19281, r19279, r19280, MPFR_RNDN); mpfr_mul(r19282, r19273, r19273, MPFR_RNDN); mpfr_mul(r19282, r19282, r19273, MPFR_RNDN); mpfr_div(r19283, r19276, r19282, MPFR_RNDN); mpfr_sub(r19284, r19281, r19283, MPFR_RNDN); ; mpfr_set_si(r19286, mpfr_cmp(r19273, r19285) <= 0, MPFR_RNDN); mpfr_sqr(r19287, r19273, MPFR_RNDN); mpfr_add(r19288, r19287, r19276, MPFR_RNDN); mpfr_div(r19289, r19273, r19288, MPFR_RNDN); if (mpfr_get_si(r19286, MPFR_RNDN)) { mpfr_set(r19290, r19289, MPFR_RNDN); } else { mpfr_set(r19290, r19284, MPFR_RNDN); }; if (mpfr_get_si(r19275, MPFR_RNDN)) { mpfr_set(r19291, r19284, MPFR_RNDN); } else { mpfr_set(r19291, r19290, MPFR_RNDN); }; return mpfr_get_d(r19291, MPFR_RNDN); }