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