#include #include #include #include #include char *name = "expq3 (problem 3.4.2)"; double f_if(float a, float b, float eps) { float r24123 = eps; float r24124 = a; float r24125 = b; float r24126 = r24124 + r24125; float r24127 = r24126 * r24123; float r24128 = exp(r24127); float r24129 = 1; float r24130 = r24128 - r24129; float r24131 = r24123 * r24130; float r24132 = r24124 * r24123; float r24133 = exp(r24132); float r24134 = r24133 - r24129; float r24135 = r24125 * r24123; float r24136 = exp(r24135); float r24137 = r24136 - r24129; float r24138 = r24134 * r24137; float r24139 = r24131 / r24138; return r24139; } double f_id(double a, double b, double eps) { double r24140 = eps; double r24141 = a; double r24142 = b; double r24143 = r24141 + r24142; double r24144 = r24143 * r24140; double r24145 = exp(r24144); double r24146 = 1; double r24147 = r24145 - r24146; double r24148 = r24140 * r24147; double r24149 = r24141 * r24140; double r24150 = exp(r24149); double r24151 = r24150 - r24146; double r24152 = r24142 * r24140; double r24153 = exp(r24152); double r24154 = r24153 - r24146; double r24155 = r24151 * r24154; double r24156 = r24148 / r24155; return r24156; } double f_of(float a, float b, float eps) { float r24157 = eps; float r24158 = a; float r24159 = r24158 * r24157; float r24160 = expm1(r24159); float r24161 = r24157 / r24160; float r24162 = b; float r24163 = r24162 + r24158; float r24164 = r24163 * r24157; float r24165 = expm1(r24164); float r24166 = r24157 * r24162; float r24167 = expm1(r24166); float r24168 = r24165 / r24167; float r24169 = r24161 * r24168; float r24170 = -9.407038744126801e-44; bool r24171 = r24169 <= r24170; float r24172 = 1; float r24173 = r24172 / r24162; float r24174 = r24172 / r24158; float r24175 = r24173 + r24174; float r24176 = 0.005120215871996153; bool r24177 = r24169 <= r24176; float r24178 = r24177 ? r24169 : r24175; float r24179 = r24171 ? r24175 : r24178; return r24179; } double f_od(double a, double b, double eps) { double r24180 = eps; double r24181 = a; double r24182 = r24181 * r24180; double r24183 = expm1(r24182); double r24184 = r24180 / r24183; double r24185 = b; double r24186 = r24185 + r24181; double r24187 = r24186 * r24180; double r24188 = expm1(r24187); double r24189 = r24180 * r24185; double r24190 = expm1(r24189); double r24191 = r24188 / r24190; double r24192 = r24184 * r24191; double r24193 = -9.407038744126801e-44; bool r24194 = r24192 <= r24193; double r24195 = 1; double r24196 = r24195 / r24185; double r24197 = r24195 / r24181; double r24198 = r24196 + r24197; double r24199 = 0.005120215871996153; bool r24200 = r24192 <= r24199; double r24201 = r24200 ? r24192 : r24198; double r24202 = r24194 ? r24198 : r24201; return r24202; } 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 r24203, r24204, r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218, r24219; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r24203); mpfr_init(r24204); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init_set_str(r24209, "1", 10, MPFR_RNDN); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init(r24215); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); mpfr_init(r24219); } double f_im(double a, double b, double eps) { mpfr_set_d(r24203, eps, MPFR_RNDN); mpfr_set_d(r24204, a, MPFR_RNDN); mpfr_set_d(r24205, b, MPFR_RNDN); mpfr_add(r24206, r24204, r24205, MPFR_RNDN); mpfr_mul(r24207, r24206, r24203, MPFR_RNDN); mpfr_exp(r24208, r24207, MPFR_RNDN); ; mpfr_sub(r24210, r24208, r24209, MPFR_RNDN); mpfr_mul(r24211, r24203, r24210, MPFR_RNDN); mpfr_mul(r24212, r24204, r24203, MPFR_RNDN); mpfr_exp(r24213, r24212, MPFR_RNDN); mpfr_sub(r24214, r24213, r24209, MPFR_RNDN); mpfr_mul(r24215, r24205, r24203, MPFR_RNDN); mpfr_exp(r24216, r24215, MPFR_RNDN); mpfr_sub(r24217, r24216, r24209, MPFR_RNDN); mpfr_mul(r24218, r24214, r24217, MPFR_RNDN); mpfr_div(r24219, r24211, r24218, MPFR_RNDN); return mpfr_get_d(r24219, MPFR_RNDN); } static mpfr_t r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r24220); mpfr_init(r24221); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init(r24225); mpfr_init(r24226); mpfr_init(r24227); mpfr_init(r24228); mpfr_init(r24229); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init_set_str(r24233, "-9.407038744126801e-44", 10, MPFR_RNDN); mpfr_init(r24234); mpfr_init_set_str(r24235, "1", 10, MPFR_RNDN); mpfr_init(r24236); mpfr_init(r24237); mpfr_init(r24238); mpfr_init_set_str(r24239, "0.005120215871996153", 10, MPFR_RNDN); mpfr_init(r24240); mpfr_init(r24241); mpfr_init(r24242); } double f_fm(double a, double b, double eps) { mpfr_set_d(r24220, eps, MPFR_RNDN); mpfr_set_d(r24221, a, MPFR_RNDN); mpfr_mul(r24222, r24221, r24220, MPFR_RNDN); mpfr_expm1(r24223, r24222, MPFR_RNDN); mpfr_div(r24224, r24220, r24223, MPFR_RNDN); mpfr_set_d(r24225, b, MPFR_RNDN); mpfr_add(r24226, r24225, r24221, MPFR_RNDN); mpfr_mul(r24227, r24226, r24220, MPFR_RNDN); mpfr_expm1(r24228, r24227, MPFR_RNDN); mpfr_mul(r24229, r24220, r24225, MPFR_RNDN); mpfr_expm1(r24230, r24229, MPFR_RNDN); mpfr_div(r24231, r24228, r24230, MPFR_RNDN); mpfr_mul(r24232, r24224, r24231, MPFR_RNDN); ; mpfr_set_si(r24234, mpfr_cmp(r24232, r24233) <= 0, MPFR_RNDN); ; mpfr_div(r24236, r24235, r24225, MPFR_RNDN); mpfr_div(r24237, r24235, r24221, MPFR_RNDN); mpfr_add(r24238, r24236, r24237, MPFR_RNDN); ; mpfr_set_si(r24240, mpfr_cmp(r24232, r24239) <= 0, MPFR_RNDN); if (mpfr_get_si(r24240, MPFR_RNDN)) { mpfr_set(r24241, r24232, MPFR_RNDN); } else { mpfr_set(r24241, r24238, MPFR_RNDN); }; if (mpfr_get_si(r24234, MPFR_RNDN)) { mpfr_set(r24242, r24238, MPFR_RNDN); } else { mpfr_set(r24242, r24241, MPFR_RNDN); }; return mpfr_get_d(r24242, MPFR_RNDN); } static mpfr_t r24243, r24244, r24245, r24246, r24247, r24248, r24249, r24250, r24251, r24252, r24253, r24254, r24255, r24256, r24257, r24258, r24259, r24260, r24261, r24262, r24263, r24264, r24265; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r24243); mpfr_init(r24244); mpfr_init(r24245); mpfr_init(r24246); mpfr_init(r24247); mpfr_init(r24248); mpfr_init(r24249); mpfr_init(r24250); mpfr_init(r24251); mpfr_init(r24252); mpfr_init(r24253); mpfr_init(r24254); mpfr_init(r24255); mpfr_init_set_str(r24256, "-9.407038744126801e-44", 10, MPFR_RNDN); mpfr_init(r24257); mpfr_init_set_str(r24258, "1", 10, MPFR_RNDN); mpfr_init(r24259); mpfr_init(r24260); mpfr_init(r24261); mpfr_init_set_str(r24262, "0.005120215871996153", 10, MPFR_RNDN); mpfr_init(r24263); mpfr_init(r24264); mpfr_init(r24265); } double f_dm(double a, double b, double eps) { mpfr_set_d(r24243, eps, MPFR_RNDN); mpfr_set_d(r24244, a, MPFR_RNDN); mpfr_mul(r24245, r24244, r24243, MPFR_RNDN); mpfr_expm1(r24246, r24245, MPFR_RNDN); mpfr_div(r24247, r24243, r24246, MPFR_RNDN); mpfr_set_d(r24248, b, MPFR_RNDN); mpfr_add(r24249, r24248, r24244, MPFR_RNDN); mpfr_mul(r24250, r24249, r24243, MPFR_RNDN); mpfr_expm1(r24251, r24250, MPFR_RNDN); mpfr_mul(r24252, r24243, r24248, MPFR_RNDN); mpfr_expm1(r24253, r24252, MPFR_RNDN); mpfr_div(r24254, r24251, r24253, MPFR_RNDN); mpfr_mul(r24255, r24247, r24254, MPFR_RNDN); ; mpfr_set_si(r24257, mpfr_cmp(r24255, r24256) <= 0, MPFR_RNDN); ; mpfr_div(r24259, r24258, r24248, MPFR_RNDN); mpfr_div(r24260, r24258, r24244, MPFR_RNDN); mpfr_add(r24261, r24259, r24260, MPFR_RNDN); ; mpfr_set_si(r24263, mpfr_cmp(r24255, r24262) <= 0, MPFR_RNDN); if (mpfr_get_si(r24263, MPFR_RNDN)) { mpfr_set(r24264, r24255, MPFR_RNDN); } else { mpfr_set(r24264, r24261, MPFR_RNDN); }; if (mpfr_get_si(r24257, MPFR_RNDN)) { mpfr_set(r24265, r24261, MPFR_RNDN); } else { mpfr_set(r24265, r24264, MPFR_RNDN); }; return mpfr_get_d(r24265, MPFR_RNDN); }