#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r22320 = x; float r22321 = exp(r22320); float r22322 = 1; float r22323 = r22321 - r22322; float r22324 = r22323 / r22320; return r22324; } double f_id(double x) { double r22325 = x; double r22326 = exp(r22325); double r22327 = 1; double r22328 = r22326 - r22327; double r22329 = r22328 / r22325; return r22329; } double f_of(float x) { float r22330 = x; float r22331 = -0.0001429957783623544; bool r22332 = r22330 <= r22331; float r22333 = r22330 + r22330; float r22334 = exp(r22333); float r22335 = 1; float r22336 = r22334 - r22335; float r22337 = exp(r22330); float r22338 = r22337 + r22335; float r22339 = r22336 / r22338; float r22340 = r22339 / r22330; float r22341 = 1/6; float r22342 = 2; float r22343 = pow(r22330, r22342); float r22344 = r22341 * r22343; float r22345 = 1/2; float r22346 = r22345 * r22330; float r22347 = r22335 + r22346; float r22348 = r22344 + r22347; float r22349 = r22332 ? r22340 : r22348; return r22349; } double f_od(double x) { double r22350 = x; double r22351 = -0.0001429957783623544; bool r22352 = r22350 <= r22351; double r22353 = r22350 + r22350; double r22354 = exp(r22353); double r22355 = 1; double r22356 = r22354 - r22355; double r22357 = exp(r22350); double r22358 = r22357 + r22355; double r22359 = r22356 / r22358; double r22360 = r22359 / r22350; double r22361 = 1/6; double r22362 = 2; double r22363 = pow(r22350, r22362); double r22364 = r22361 * r22363; double r22365 = 1/2; double r22366 = r22365 * r22350; double r22367 = r22355 + r22366; double r22368 = r22364 + r22367; double r22369 = r22352 ? r22360 : r22368; return r22369; } 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 r22370, r22371, r22372, r22373, r22374; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22370); mpfr_init(r22371); mpfr_init_set_str(r22372, "1", 10, MPFR_RNDN); mpfr_init(r22373); mpfr_init(r22374); } double f_im(double x) { mpfr_set_d(r22370, x, MPFR_RNDN); mpfr_exp(r22371, r22370, MPFR_RNDN); ; mpfr_sub(r22373, r22371, r22372, MPFR_RNDN); mpfr_div(r22374, r22373, r22370, MPFR_RNDN); return mpfr_get_d(r22374, MPFR_RNDN); } static mpfr_t r22375, r22376, r22377, r22378, r22379, r22380, r22381, r22382, r22383, r22384, r22385, r22386, r22387, r22388, r22389, r22390, r22391, r22392, r22393, r22394; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22375); mpfr_init_set_str(r22376, "-0.0001429957783623544", 10, MPFR_RNDN); mpfr_init(r22377); mpfr_init(r22378); mpfr_init(r22379); mpfr_init_set_str(r22380, "1", 10, MPFR_RNDN); mpfr_init(r22381); mpfr_init(r22382); mpfr_init(r22383); mpfr_init(r22384); mpfr_init(r22385); mpfr_init_set_str(r22386, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r22387, "2", 10, MPFR_RNDN); mpfr_init(r22388); mpfr_init(r22389); mpfr_init_set_str(r22390, "1/2", 10, MPFR_RNDN); mpfr_init(r22391); mpfr_init(r22392); mpfr_init(r22393); mpfr_init(r22394); } double f_fm(double x) { mpfr_set_d(r22375, x, MPFR_RNDN); ; mpfr_set_si(r22377, mpfr_cmp(r22375, r22376) <= 0, MPFR_RNDN); mpfr_add(r22378, r22375, r22375, MPFR_RNDN); mpfr_exp(r22379, r22378, MPFR_RNDN); ; mpfr_sub(r22381, r22379, r22380, MPFR_RNDN); mpfr_exp(r22382, r22375, MPFR_RNDN); mpfr_add(r22383, r22382, r22380, MPFR_RNDN); mpfr_div(r22384, r22381, r22383, MPFR_RNDN); mpfr_div(r22385, r22384, r22375, MPFR_RNDN); ; ; mpfr_pow(r22388, r22375, r22387, MPFR_RNDN); mpfr_mul(r22389, r22386, r22388, MPFR_RNDN); ; mpfr_mul(r22391, r22390, r22375, MPFR_RNDN); mpfr_add(r22392, r22380, r22391, MPFR_RNDN); mpfr_add(r22393, r22389, r22392, MPFR_RNDN); if (mpfr_get_si(r22377, MPFR_RNDN)) { mpfr_set(r22394, r22385, MPFR_RNDN); } else { mpfr_set(r22394, r22393, MPFR_RNDN); }; return mpfr_get_d(r22394, MPFR_RNDN); } static mpfr_t r22395, r22396, r22397, r22398, r22399, r22400, r22401, r22402, r22403, r22404, r22405, r22406, r22407, r22408, r22409, r22410, r22411, r22412, r22413, r22414; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22395); mpfr_init_set_str(r22396, "-0.0001429957783623544", 10, MPFR_RNDN); mpfr_init(r22397); mpfr_init(r22398); mpfr_init(r22399); mpfr_init_set_str(r22400, "1", 10, MPFR_RNDN); mpfr_init(r22401); mpfr_init(r22402); mpfr_init(r22403); mpfr_init(r22404); mpfr_init(r22405); mpfr_init_set_str(r22406, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r22407, "2", 10, MPFR_RNDN); mpfr_init(r22408); mpfr_init(r22409); mpfr_init_set_str(r22410, "1/2", 10, MPFR_RNDN); mpfr_init(r22411); mpfr_init(r22412); mpfr_init(r22413); mpfr_init(r22414); } double f_dm(double x) { mpfr_set_d(r22395, x, MPFR_RNDN); ; mpfr_set_si(r22397, mpfr_cmp(r22395, r22396) <= 0, MPFR_RNDN); mpfr_add(r22398, r22395, r22395, MPFR_RNDN); mpfr_exp(r22399, r22398, MPFR_RNDN); ; mpfr_sub(r22401, r22399, r22400, MPFR_RNDN); mpfr_exp(r22402, r22395, MPFR_RNDN); mpfr_add(r22403, r22402, r22400, MPFR_RNDN); mpfr_div(r22404, r22401, r22403, MPFR_RNDN); mpfr_div(r22405, r22404, r22395, MPFR_RNDN); ; ; mpfr_pow(r22408, r22395, r22407, MPFR_RNDN); mpfr_mul(r22409, r22406, r22408, MPFR_RNDN); ; mpfr_mul(r22411, r22410, r22395, MPFR_RNDN); mpfr_add(r22412, r22400, r22411, MPFR_RNDN); mpfr_add(r22413, r22409, r22412, MPFR_RNDN); if (mpfr_get_si(r22397, MPFR_RNDN)) { mpfr_set(r22414, r22405, MPFR_RNDN); } else { mpfr_set(r22414, r22413, MPFR_RNDN); }; return mpfr_get_d(r22414, MPFR_RNDN); }