#include #include #include #include #include char *name = "expax (section 3.5)"; double f_if(float a, float x) { float r24153 = a; float r24154 = x; float r24155 = r24153 * r24154; float r24156 = exp(r24155); float r24157 = 1; float r24158 = r24156 - r24157; return r24158; } double f_id(double a, double x) { double r24159 = a; double r24160 = x; double r24161 = r24159 * r24160; double r24162 = exp(r24161); double r24163 = 1; double r24164 = r24162 - r24163; return r24164; } double f_of(float a, float x) { float r24165 = a; float r24166 = x; float r24167 = r24165 * r24166; float r24168 = exp(r24167); float r24169 = 1; float r24170 = r24168 - r24169; float r24171 = sqrt(r24170); float r24172 = r24171 * r24171; float r24173 = 0.0054314674489082245; bool r24174 = r24172 <= r24173; float r24175 = r24167 * r24167; float r24176 = 1/6; float r24177 = r24167 * r24176; float r24178 = 1/2; float r24179 = r24177 + r24178; float r24180 = r24175 * r24179; float r24181 = r24180 + r24167; float r24182 = exp(r24170); float r24183 = log(r24182); float r24184 = r24174 ? r24181 : r24183; return r24184; } double f_od(double a, double x) { double r24185 = a; double r24186 = x; double r24187 = r24185 * r24186; double r24188 = exp(r24187); double r24189 = 1; double r24190 = r24188 - r24189; double r24191 = sqrt(r24190); double r24192 = r24191 * r24191; double r24193 = 0.0054314674489082245; bool r24194 = r24192 <= r24193; double r24195 = r24187 * r24187; double r24196 = 1/6; double r24197 = r24187 * r24196; double r24198 = 1/2; double r24199 = r24197 + r24198; double r24200 = r24195 * r24199; double r24201 = r24200 + r24187; double r24202 = exp(r24190); double r24203 = log(r24202); double r24204 = r24194 ? r24201 : r24203; return r24204; } 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 r24205, r24206, r24207, r24208, r24209, r24210; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init_set_str(r24209, "1", 10, MPFR_RNDN); mpfr_init(r24210); } double f_im(double a, double x) { mpfr_set_d(r24205, a, MPFR_RNDN); mpfr_set_d(r24206, x, MPFR_RNDN); mpfr_mul(r24207, r24205, r24206, MPFR_RNDN); mpfr_exp(r24208, r24207, MPFR_RNDN); ; mpfr_sub(r24210, r24208, r24209, MPFR_RNDN); return mpfr_get_d(r24210, MPFR_RNDN); } static mpfr_t r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218, r24219, r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init_set_str(r24215, "1", 10, MPFR_RNDN); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); mpfr_init_set_str(r24219, "0.0054314674489082245", 10, MPFR_RNDN); mpfr_init(r24220); mpfr_init(r24221); mpfr_init_set_str(r24222, "1/6", 10, MPFR_RNDN); mpfr_init(r24223); mpfr_init_set_str(r24224, "1/2", 10, MPFR_RNDN); mpfr_init(r24225); mpfr_init(r24226); mpfr_init(r24227); mpfr_init(r24228); mpfr_init(r24229); mpfr_init(r24230); } double f_fm(double a, double x) { mpfr_set_d(r24211, a, MPFR_RNDN); mpfr_set_d(r24212, x, MPFR_RNDN); mpfr_mul(r24213, r24211, r24212, MPFR_RNDN); mpfr_exp(r24214, r24213, MPFR_RNDN); ; mpfr_sub(r24216, r24214, r24215, MPFR_RNDN); mpfr_sqrt(r24217, r24216, MPFR_RNDN); mpfr_mul(r24218, r24217, r24217, MPFR_RNDN); ; mpfr_set_si(r24220, mpfr_cmp(r24218, r24219) <= 0, MPFR_RNDN); mpfr_mul(r24221, r24213, r24213, MPFR_RNDN); ; mpfr_mul(r24223, r24213, r24222, MPFR_RNDN); ; mpfr_add(r24225, r24223, r24224, MPFR_RNDN); mpfr_mul(r24226, r24221, r24225, MPFR_RNDN); mpfr_add(r24227, r24226, r24213, MPFR_RNDN); mpfr_exp(r24228, r24216, MPFR_RNDN); mpfr_log(r24229, r24228, MPFR_RNDN); if (mpfr_get_si(r24220, MPFR_RNDN)) { mpfr_set(r24230, r24227, MPFR_RNDN); } else { mpfr_set(r24230, r24229, MPFR_RNDN); }; return mpfr_get_d(r24230, MPFR_RNDN); } static mpfr_t r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242, r24243, r24244, r24245, r24246, r24247, r24248, r24249, r24250; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); 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.0054314674489082245", 10, MPFR_RNDN); mpfr_init(r24240); mpfr_init(r24241); mpfr_init_set_str(r24242, "1/6", 10, MPFR_RNDN); mpfr_init(r24243); mpfr_init_set_str(r24244, "1/2", 10, MPFR_RNDN); mpfr_init(r24245); mpfr_init(r24246); mpfr_init(r24247); mpfr_init(r24248); mpfr_init(r24249); mpfr_init(r24250); } double f_dm(double a, double x) { mpfr_set_d(r24231, a, MPFR_RNDN); mpfr_set_d(r24232, x, MPFR_RNDN); mpfr_mul(r24233, r24231, r24232, MPFR_RNDN); mpfr_exp(r24234, r24233, MPFR_RNDN); ; mpfr_sub(r24236, r24234, r24235, MPFR_RNDN); mpfr_sqrt(r24237, r24236, MPFR_RNDN); mpfr_mul(r24238, r24237, r24237, MPFR_RNDN); ; mpfr_set_si(r24240, mpfr_cmp(r24238, r24239) <= 0, MPFR_RNDN); mpfr_mul(r24241, r24233, r24233, MPFR_RNDN); ; mpfr_mul(r24243, r24233, r24242, MPFR_RNDN); ; mpfr_add(r24245, r24243, r24244, MPFR_RNDN); mpfr_mul(r24246, r24241, r24245, MPFR_RNDN); mpfr_add(r24247, r24246, r24233, MPFR_RNDN); mpfr_exp(r24248, r24236, MPFR_RNDN); mpfr_log(r24249, r24248, MPFR_RNDN); if (mpfr_get_si(r24240, MPFR_RNDN)) { mpfr_set(r24250, r24247, MPFR_RNDN); } else { mpfr_set(r24250, r24249, MPFR_RNDN); }; return mpfr_get_d(r24250, MPFR_RNDN); }