#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r9109 = x; float r9110 = r9109 * r9109; float r9111 = 1; float r9112 = r9110 + r9111; float r9113 = r9109 / r9112; return r9113; } double f_id(double x) { double r9114 = x; double r9115 = r9114 * r9114; double r9116 = 1; double r9117 = r9115 + r9116; double r9118 = r9114 / r9117; return r9118; } double f_of(float x) { float r9119 = x; float r9120 = r9119 * r9119; float r9121 = 1; float r9122 = r9120 + r9121; float r9123 = r9119 / r9122; float r9124 = -4.3207196525733e-310; bool r9125 = r9123 <= r9124; float r9126 = sqrt(r9122); float r9127 = r9119 / r9126; float r9128 = r9127 / r9126; float r9129 = 4.21129123229683e-310; bool r9130 = r9123 <= r9129; float r9131 = 5; float r9132 = pow(r9119, r9131); float r9133 = r9121 / r9132; float r9134 = r9121 / r9119; float r9135 = r9133 + r9134; float r9136 = 3; float r9137 = pow(r9119, r9136); float r9138 = r9121 / r9137; float r9139 = r9135 - r9138; float r9140 = r9130 ? r9139 : r9128; float r9141 = r9125 ? r9128 : r9140; return r9141; } double f_od(double x) { double r9142 = x; double r9143 = r9142 * r9142; double r9144 = 1; double r9145 = r9143 + r9144; double r9146 = r9142 / r9145; double r9147 = -4.3207196525733e-310; bool r9148 = r9146 <= r9147; double r9149 = sqrt(r9145); double r9150 = r9142 / r9149; double r9151 = r9150 / r9149; double r9152 = 4.21129123229683e-310; bool r9153 = r9146 <= r9152; double r9154 = 5; double r9155 = pow(r9142, r9154); double r9156 = r9144 / r9155; double r9157 = r9144 / r9142; double r9158 = r9156 + r9157; double r9159 = 3; double r9160 = pow(r9142, r9159); double r9161 = r9144 / r9160; double r9162 = r9158 - r9161; double r9163 = r9153 ? r9162 : r9151; double r9164 = r9148 ? r9151 : r9163; return r9164; } 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 r9165, r9166, r9167, r9168, r9169; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r9165); mpfr_init(r9166); mpfr_init_set_str(r9167, "1", 10, MPFR_RNDN); mpfr_init(r9168); mpfr_init(r9169); } double f_im(double x) { mpfr_set_d(r9165, x, MPFR_RNDN); mpfr_mul(r9166, r9165, r9165, MPFR_RNDN); ; mpfr_add(r9168, r9166, r9167, MPFR_RNDN); mpfr_div(r9169, r9165, r9168, MPFR_RNDN); return mpfr_get_d(r9169, MPFR_RNDN); } static mpfr_t r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185, r9186, r9187, r9188, r9189, r9190, r9191, r9192; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r9170); mpfr_init(r9171); mpfr_init_set_str(r9172, "1", 10, MPFR_RNDN); mpfr_init(r9173); mpfr_init(r9174); mpfr_init_set_str(r9175, "-4.3207196525733e-310", 10, MPFR_RNDN); mpfr_init(r9176); mpfr_init(r9177); mpfr_init(r9178); mpfr_init(r9179); mpfr_init_set_str(r9180, "4.21129123229683e-310", 10, MPFR_RNDN); mpfr_init(r9181); mpfr_init_set_str(r9182, "5", 10, MPFR_RNDN); mpfr_init(r9183); mpfr_init(r9184); mpfr_init(r9185); mpfr_init(r9186); mpfr_init_set_str(r9187, "3", 10, MPFR_RNDN); mpfr_init(r9188); mpfr_init(r9189); mpfr_init(r9190); mpfr_init(r9191); mpfr_init(r9192); } double f_fm(double x) { mpfr_set_d(r9170, x, MPFR_RNDN); mpfr_mul(r9171, r9170, r9170, MPFR_RNDN); ; mpfr_add(r9173, r9171, r9172, MPFR_RNDN); mpfr_div(r9174, r9170, r9173, MPFR_RNDN); ; mpfr_set_si(r9176, mpfr_cmp(r9174, r9175) <= 0, MPFR_RNDN); mpfr_sqrt(r9177, r9173, MPFR_RNDN); mpfr_div(r9178, r9170, r9177, MPFR_RNDN); mpfr_div(r9179, r9178, r9177, MPFR_RNDN); ; mpfr_set_si(r9181, mpfr_cmp(r9174, r9180) <= 0, MPFR_RNDN); ; mpfr_pow(r9183, r9170, r9182, MPFR_RNDN); mpfr_div(r9184, r9172, r9183, MPFR_RNDN); mpfr_div(r9185, r9172, r9170, MPFR_RNDN); mpfr_add(r9186, r9184, r9185, MPFR_RNDN); ; mpfr_pow(r9188, r9170, r9187, MPFR_RNDN); mpfr_div(r9189, r9172, r9188, MPFR_RNDN); mpfr_sub(r9190, r9186, r9189, MPFR_RNDN); if (mpfr_get_si(r9181, MPFR_RNDN)) { mpfr_set(r9191, r9190, MPFR_RNDN); } else { mpfr_set(r9191, r9179, MPFR_RNDN); }; if (mpfr_get_si(r9176, MPFR_RNDN)) { mpfr_set(r9192, r9179, MPFR_RNDN); } else { mpfr_set(r9192, r9191, MPFR_RNDN); }; return mpfr_get_d(r9192, MPFR_RNDN); } static mpfr_t r9193, r9194, r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204, r9205, r9206, r9207, r9208, r9209, r9210, r9211, r9212, r9213, r9214, r9215; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r9193); mpfr_init(r9194); mpfr_init_set_str(r9195, "1", 10, MPFR_RNDN); mpfr_init(r9196); mpfr_init(r9197); mpfr_init_set_str(r9198, "-4.3207196525733e-310", 10, MPFR_RNDN); mpfr_init(r9199); mpfr_init(r9200); mpfr_init(r9201); mpfr_init(r9202); mpfr_init_set_str(r9203, "4.21129123229683e-310", 10, MPFR_RNDN); mpfr_init(r9204); mpfr_init_set_str(r9205, "5", 10, MPFR_RNDN); mpfr_init(r9206); mpfr_init(r9207); mpfr_init(r9208); mpfr_init(r9209); mpfr_init_set_str(r9210, "3", 10, MPFR_RNDN); mpfr_init(r9211); mpfr_init(r9212); mpfr_init(r9213); mpfr_init(r9214); mpfr_init(r9215); } double f_dm(double x) { mpfr_set_d(r9193, x, MPFR_RNDN); mpfr_mul(r9194, r9193, r9193, MPFR_RNDN); ; mpfr_add(r9196, r9194, r9195, MPFR_RNDN); mpfr_div(r9197, r9193, r9196, MPFR_RNDN); ; mpfr_set_si(r9199, mpfr_cmp(r9197, r9198) <= 0, MPFR_RNDN); mpfr_sqrt(r9200, r9196, MPFR_RNDN); mpfr_div(r9201, r9193, r9200, MPFR_RNDN); mpfr_div(r9202, r9201, r9200, MPFR_RNDN); ; mpfr_set_si(r9204, mpfr_cmp(r9197, r9203) <= 0, MPFR_RNDN); ; mpfr_pow(r9206, r9193, r9205, MPFR_RNDN); mpfr_div(r9207, r9195, r9206, MPFR_RNDN); mpfr_div(r9208, r9195, r9193, MPFR_RNDN); mpfr_add(r9209, r9207, r9208, MPFR_RNDN); ; mpfr_pow(r9211, r9193, r9210, MPFR_RNDN); mpfr_div(r9212, r9195, r9211, MPFR_RNDN); mpfr_sub(r9213, r9209, r9212, MPFR_RNDN); if (mpfr_get_si(r9204, MPFR_RNDN)) { mpfr_set(r9214, r9213, MPFR_RNDN); } else { mpfr_set(r9214, r9202, MPFR_RNDN); }; if (mpfr_get_si(r9199, MPFR_RNDN)) { mpfr_set(r9215, r9202, MPFR_RNDN); } else { mpfr_set(r9215, r9214, MPFR_RNDN); }; return mpfr_get_d(r9215, MPFR_RNDN); }