#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r26125 = x; float r26126 = exp(r26125); float r26127 = 1; float r26128 = r26126 - r26127; float r26129 = r26128 / r26125; return r26129; } double f_id(double x) { double r26130 = x; double r26131 = exp(r26130); double r26132 = 1; double r26133 = r26131 - r26132; double r26134 = r26133 / r26130; return r26134; } double f_of(float x) { float r26135 = x; float r26136 = -0.00010281873558823207; bool r26137 = r26135 <= r26136; float r26138 = r26135 + r26135; float r26139 = exp(r26138); float r26140 = 1; float r26141 = r26139 - r26140; float r26142 = exp(r26135); float r26143 = r26142 + r26140; float r26144 = r26141 / r26143; float r26145 = r26144 / r26135; float r26146 = 1/2; float r26147 = 2; float r26148 = pow(r26135, r26147); float r26149 = r26146 * r26148; float r26150 = 1/6; float r26151 = 3; float r26152 = pow(r26135, r26151); float r26153 = r26150 * r26152; float r26154 = r26153 + r26135; float r26155 = r26149 + r26154; float r26156 = r26155 / r26135; float r26157 = r26137 ? r26145 : r26156; return r26157; } double f_od(double x) { double r26158 = x; double r26159 = -0.00010281873558823207; bool r26160 = r26158 <= r26159; double r26161 = r26158 + r26158; double r26162 = exp(r26161); double r26163 = 1; double r26164 = r26162 - r26163; double r26165 = exp(r26158); double r26166 = r26165 + r26163; double r26167 = r26164 / r26166; double r26168 = r26167 / r26158; double r26169 = 1/2; double r26170 = 2; double r26171 = pow(r26158, r26170); double r26172 = r26169 * r26171; double r26173 = 1/6; double r26174 = 3; double r26175 = pow(r26158, r26174); double r26176 = r26173 * r26175; double r26177 = r26176 + r26158; double r26178 = r26172 + r26177; double r26179 = r26178 / r26158; double r26180 = r26160 ? r26168 : r26179; return r26180; } 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 r26181, r26182, r26183, r26184, r26185; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26181); mpfr_init(r26182); mpfr_init_set_str(r26183, "1", 10, MPFR_RNDN); mpfr_init(r26184); mpfr_init(r26185); } double f_im(double x) { mpfr_set_d(r26181, x, MPFR_RNDN); mpfr_exp(r26182, r26181, MPFR_RNDN); ; mpfr_sub(r26184, r26182, r26183, MPFR_RNDN); mpfr_div(r26185, r26184, r26181, MPFR_RNDN); return mpfr_get_d(r26185, MPFR_RNDN); } static mpfr_t r26186, r26187, r26188, r26189, r26190, r26191, r26192, r26193, r26194, r26195, r26196, r26197, r26198, r26199, r26200, r26201, r26202, r26203, r26204, r26205, r26206, r26207, r26208; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26186); mpfr_init_set_str(r26187, "-0.00010281873558823207", 10, MPFR_RNDN); mpfr_init(r26188); mpfr_init(r26189); mpfr_init(r26190); mpfr_init_set_str(r26191, "1", 10, MPFR_RNDN); mpfr_init(r26192); mpfr_init(r26193); mpfr_init(r26194); mpfr_init(r26195); mpfr_init(r26196); mpfr_init_set_str(r26197, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26198, "2", 10, MPFR_RNDN); mpfr_init(r26199); mpfr_init(r26200); mpfr_init_set_str(r26201, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r26202, "3", 10, MPFR_RNDN); mpfr_init(r26203); mpfr_init(r26204); mpfr_init(r26205); mpfr_init(r26206); mpfr_init(r26207); mpfr_init(r26208); } double f_fm(double x) { mpfr_set_d(r26186, x, MPFR_RNDN); ; mpfr_set_si(r26188, mpfr_cmp(r26186, r26187) <= 0, MPFR_RNDN); mpfr_add(r26189, r26186, r26186, MPFR_RNDN); mpfr_exp(r26190, r26189, MPFR_RNDN); ; mpfr_sub(r26192, r26190, r26191, MPFR_RNDN); mpfr_exp(r26193, r26186, MPFR_RNDN); mpfr_add(r26194, r26193, r26191, MPFR_RNDN); mpfr_div(r26195, r26192, r26194, MPFR_RNDN); mpfr_div(r26196, r26195, r26186, MPFR_RNDN); ; ; mpfr_pow(r26199, r26186, r26198, MPFR_RNDN); mpfr_mul(r26200, r26197, r26199, MPFR_RNDN); ; ; mpfr_pow(r26203, r26186, r26202, MPFR_RNDN); mpfr_mul(r26204, r26201, r26203, MPFR_RNDN); mpfr_add(r26205, r26204, r26186, MPFR_RNDN); mpfr_add(r26206, r26200, r26205, MPFR_RNDN); mpfr_div(r26207, r26206, r26186, MPFR_RNDN); if (mpfr_get_si(r26188, MPFR_RNDN)) { mpfr_set(r26208, r26196, MPFR_RNDN); } else { mpfr_set(r26208, r26207, MPFR_RNDN); }; return mpfr_get_d(r26208, MPFR_RNDN); } static mpfr_t r26209, r26210, r26211, r26212, r26213, r26214, r26215, r26216, r26217, r26218, r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229, r26230, r26231; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26209); mpfr_init_set_str(r26210, "-0.00010281873558823207", 10, MPFR_RNDN); mpfr_init(r26211); mpfr_init(r26212); mpfr_init(r26213); mpfr_init_set_str(r26214, "1", 10, MPFR_RNDN); mpfr_init(r26215); mpfr_init(r26216); mpfr_init(r26217); mpfr_init(r26218); mpfr_init(r26219); mpfr_init_set_str(r26220, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26221, "2", 10, MPFR_RNDN); mpfr_init(r26222); mpfr_init(r26223); mpfr_init_set_str(r26224, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r26225, "3", 10, MPFR_RNDN); mpfr_init(r26226); mpfr_init(r26227); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init(r26231); } double f_dm(double x) { mpfr_set_d(r26209, x, MPFR_RNDN); ; mpfr_set_si(r26211, mpfr_cmp(r26209, r26210) <= 0, MPFR_RNDN); mpfr_add(r26212, r26209, r26209, MPFR_RNDN); mpfr_exp(r26213, r26212, MPFR_RNDN); ; mpfr_sub(r26215, r26213, r26214, MPFR_RNDN); mpfr_exp(r26216, r26209, MPFR_RNDN); mpfr_add(r26217, r26216, r26214, MPFR_RNDN); mpfr_div(r26218, r26215, r26217, MPFR_RNDN); mpfr_div(r26219, r26218, r26209, MPFR_RNDN); ; ; mpfr_pow(r26222, r26209, r26221, MPFR_RNDN); mpfr_mul(r26223, r26220, r26222, MPFR_RNDN); ; ; mpfr_pow(r26226, r26209, r26225, MPFR_RNDN); mpfr_mul(r26227, r26224, r26226, MPFR_RNDN); mpfr_add(r26228, r26227, r26209, MPFR_RNDN); mpfr_add(r26229, r26223, r26228, MPFR_RNDN); mpfr_div(r26230, r26229, r26209, MPFR_RNDN); if (mpfr_get_si(r26211, MPFR_RNDN)) { mpfr_set(r26231, r26219, MPFR_RNDN); } else { mpfr_set(r26231, r26230, MPFR_RNDN); }; return mpfr_get_d(r26231, MPFR_RNDN); }