#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r16183 = x; float r16184 = exp(r16183); float r16185 = 1.0f; float r16186 = r16184 - r16185; float r16187 = r16186 / r16183; return r16187; } double f_id(double x) { double r16188 = x; double r16189 = exp(r16188); double r16190 = 1.0; double r16191 = r16189 - r16190; double r16192 = r16191 / r16188; return r16192; } double f_of(float x) { float r16193 = x; float r16194 = -3.2784791869173944e-06f; bool r16195 = r16193 <= r16194; float r16196 = exp(r16193); float r16197 = r16196 * (r16196 * r16196); float r16198 = 1.0f; float r16199 = r16197 - r16198; float r16200 = r16196 + r16198; float r16201 = r16193 + r16193; float r16202 = exp(r16201); float r16203 = r16200 + r16202; float r16204 = r16199 / r16203; float r16205 = r16204 / r16193; float r16206 = 0.5f; float r16207 = 0.16666666666666666f; float r16208 = r16207 * r16193; float r16209 = r16206 + r16208; float r16210 = r16193 * r16193; float r16211 = r16209 * r16210; float r16212 = r16211 + r16193; float r16213 = r16212 / r16193; float r16214 = r16195 ? r16205 : r16213; return r16214; } double f_od(double x) { double r16215 = x; double r16216 = -3.2784791869173944e-06; bool r16217 = r16215 <= r16216; double r16218 = exp(r16215); double r16219 = r16218 * (r16218 * r16218); double r16220 = 1.0; double r16221 = r16219 - r16220; double r16222 = r16218 + r16220; double r16223 = r16215 + r16215; double r16224 = exp(r16223); double r16225 = r16222 + r16224; double r16226 = r16221 / r16225; double r16227 = r16226 / r16215; double r16228 = 0.5; double r16229 = 0.16666666666666666; double r16230 = r16229 * r16215; double r16231 = r16228 + r16230; double r16232 = r16215 * r16215; double r16233 = r16231 * r16232; double r16234 = r16233 + r16215; double r16235 = r16234 / r16215; double r16236 = r16217 ? r16227 : r16235; return r16236; } 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 r16237, r16238, r16239, r16240, r16241; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16237); mpfr_init(r16238); mpfr_init_set_str(r16239, "1", 10, MPFR_RNDN); mpfr_init(r16240); mpfr_init(r16241); } double f_im(double x) { mpfr_set_d(r16237, x, MPFR_RNDN); mpfr_exp(r16238, r16237, MPFR_RNDN); ; mpfr_sub(r16240, r16238, r16239, MPFR_RNDN); mpfr_div(r16241, r16240, r16237, MPFR_RNDN); return mpfr_get_d(r16241, MPFR_RNDN); } static mpfr_t r16242, r16243, r16244, r16245, r16246, r16247, r16248, r16249, r16250, r16251, r16252, r16253, r16254, r16255, r16256, r16257, r16258, r16259, r16260, r16261, r16262, r16263; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16242); mpfr_init_set_str(r16243, "-3.2784791869173944e-06", 10, MPFR_RNDN); mpfr_init(r16244); mpfr_init(r16245); mpfr_init(r16246); mpfr_init_set_str(r16247, "1", 10, MPFR_RNDN); mpfr_init(r16248); mpfr_init(r16249); mpfr_init(r16250); mpfr_init(r16251); mpfr_init(r16252); mpfr_init(r16253); mpfr_init(r16254); mpfr_init_set_str(r16255, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r16256, "1/6", 10, MPFR_RNDN); mpfr_init(r16257); mpfr_init(r16258); mpfr_init(r16259); mpfr_init(r16260); mpfr_init(r16261); mpfr_init(r16262); mpfr_init(r16263); } double f_fm(double x) { mpfr_set_d(r16242, x, MPFR_RNDN); ; mpfr_set_si(r16244, mpfr_cmp(r16242, r16243) <= 0, MPFR_RNDN); mpfr_exp(r16245, r16242, MPFR_RNDN); mpfr_mul(r16246, r16245, r16245, MPFR_RNDN); mpfr_mul(r16246, r16246, r16245, MPFR_RNDN); ; mpfr_sub(r16248, r16246, r16247, MPFR_RNDN); mpfr_add(r16249, r16245, r16247, MPFR_RNDN); mpfr_add(r16250, r16242, r16242, MPFR_RNDN); mpfr_exp(r16251, r16250, MPFR_RNDN); mpfr_add(r16252, r16249, r16251, MPFR_RNDN); mpfr_div(r16253, r16248, r16252, MPFR_RNDN); mpfr_div(r16254, r16253, r16242, MPFR_RNDN); ; ; mpfr_mul(r16257, r16256, r16242, MPFR_RNDN); mpfr_add(r16258, r16255, r16257, MPFR_RNDN); mpfr_sqr(r16259, r16242, MPFR_RNDN); mpfr_mul(r16260, r16258, r16259, MPFR_RNDN); mpfr_add(r16261, r16260, r16242, MPFR_RNDN); mpfr_div(r16262, r16261, r16242, MPFR_RNDN); if (mpfr_get_si(r16244, MPFR_RNDN)) { mpfr_set(r16263, r16254, MPFR_RNDN); } else { mpfr_set(r16263, r16262, MPFR_RNDN); }; return mpfr_get_d(r16263, MPFR_RNDN); } static mpfr_t r16264, r16265, r16266, r16267, r16268, r16269, r16270, r16271, r16272, r16273, r16274, r16275, r16276, r16277, r16278, r16279, r16280, r16281, r16282, r16283, r16284, r16285; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16264); mpfr_init_set_str(r16265, "-3.2784791869173944e-06", 10, MPFR_RNDN); mpfr_init(r16266); mpfr_init(r16267); mpfr_init(r16268); mpfr_init_set_str(r16269, "1", 10, MPFR_RNDN); mpfr_init(r16270); mpfr_init(r16271); mpfr_init(r16272); mpfr_init(r16273); mpfr_init(r16274); mpfr_init(r16275); mpfr_init(r16276); mpfr_init_set_str(r16277, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r16278, "1/6", 10, MPFR_RNDN); mpfr_init(r16279); mpfr_init(r16280); mpfr_init(r16281); mpfr_init(r16282); mpfr_init(r16283); mpfr_init(r16284); mpfr_init(r16285); } double f_dm(double x) { mpfr_set_d(r16264, x, MPFR_RNDN); ; mpfr_set_si(r16266, mpfr_cmp(r16264, r16265) <= 0, MPFR_RNDN); mpfr_exp(r16267, r16264, MPFR_RNDN); mpfr_mul(r16268, r16267, r16267, MPFR_RNDN); mpfr_mul(r16268, r16268, r16267, MPFR_RNDN); ; mpfr_sub(r16270, r16268, r16269, MPFR_RNDN); mpfr_add(r16271, r16267, r16269, MPFR_RNDN); mpfr_add(r16272, r16264, r16264, MPFR_RNDN); mpfr_exp(r16273, r16272, MPFR_RNDN); mpfr_add(r16274, r16271, r16273, MPFR_RNDN); mpfr_div(r16275, r16270, r16274, MPFR_RNDN); mpfr_div(r16276, r16275, r16264, MPFR_RNDN); ; ; mpfr_mul(r16279, r16278, r16264, MPFR_RNDN); mpfr_add(r16280, r16277, r16279, MPFR_RNDN); mpfr_sqr(r16281, r16264, MPFR_RNDN); mpfr_mul(r16282, r16280, r16281, MPFR_RNDN); mpfr_add(r16283, r16282, r16264, MPFR_RNDN); mpfr_div(r16284, r16283, r16264, MPFR_RNDN); if (mpfr_get_si(r16266, MPFR_RNDN)) { mpfr_set(r16285, r16276, MPFR_RNDN); } else { mpfr_set(r16285, r16284, MPFR_RNDN); }; return mpfr_get_d(r16285, MPFR_RNDN); }