#include #include #include #include #include char *name = "expax (section 3.5)"; double f_if(float a, float x) { float r27075 = a; float r27076 = x; float r27077 = r27075 * r27076; float r27078 = exp(r27077); float r27079 = 1; float r27080 = r27078 - r27079; return r27080; } double f_id(double a, double x) { double r27081 = a; double r27082 = x; double r27083 = r27081 * r27082; double r27084 = exp(r27083); double r27085 = 1; double r27086 = r27084 - r27085; return r27086; } double f_of(float a, float x) { float r27087 = a; float r27088 = x; float r27089 = r27087 * r27088; float r27090 = exp(r27089); float r27091 = 1; float r27092 = r27090 - r27091; float r27093 = -0.0008112531176381211; bool r27094 = r27092 <= r27093; float r27095 = r27088 * r27087; float r27096 = exp(r27095); float r27097 = r27096 - r27091; float r27098 = 3; float r27099 = pow(r27097, r27098); float r27100 = cbrt(r27099); float r27101 = r27089 * r27089; float r27102 = 1/6; float r27103 = r27089 * r27102; float r27104 = 1/2; float r27105 = r27103 + r27104; float r27106 = r27101 * r27105; float r27107 = r27106 + r27089; float r27108 = r27094 ? r27100 : r27107; return r27108; } double f_od(double a, double x) { double r27109 = a; double r27110 = x; double r27111 = r27109 * r27110; double r27112 = exp(r27111); double r27113 = 1; double r27114 = r27112 - r27113; double r27115 = -0.0008112531176381211; bool r27116 = r27114 <= r27115; double r27117 = r27110 * r27109; double r27118 = exp(r27117); double r27119 = r27118 - r27113; double r27120 = 3; double r27121 = pow(r27119, r27120); double r27122 = cbrt(r27121); double r27123 = r27111 * r27111; double r27124 = 1/6; double r27125 = r27111 * r27124; double r27126 = 1/2; double r27127 = r27125 + r27126; double r27128 = r27123 * r27127; double r27129 = r27128 + r27111; double r27130 = r27116 ? r27122 : r27129; return r27130; } 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 r27131, r27132, r27133, r27134, r27135, r27136; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init_set_str(r27135, "1", 10, MPFR_RNDN); mpfr_init(r27136); } double f_im(double a, double x) { mpfr_set_d(r27131, a, MPFR_RNDN); mpfr_set_d(r27132, x, MPFR_RNDN); mpfr_mul(r27133, r27131, r27132, MPFR_RNDN); mpfr_exp(r27134, r27133, MPFR_RNDN); ; mpfr_sub(r27136, r27134, r27135, MPFR_RNDN); return mpfr_get_d(r27136, MPFR_RNDN); } static mpfr_t r27137, r27138, r27139, r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27137); mpfr_init(r27138); mpfr_init(r27139); mpfr_init(r27140); mpfr_init_set_str(r27141, "1", 10, MPFR_RNDN); mpfr_init(r27142); mpfr_init_set_str(r27143, "-0.0008112531176381211", 10, MPFR_RNDN); mpfr_init(r27144); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init_set_str(r27148, "3", 10, MPFR_RNDN); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init_set_str(r27152, "1/6", 10, MPFR_RNDN); mpfr_init(r27153); mpfr_init_set_str(r27154, "1/2", 10, MPFR_RNDN); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); mpfr_init(r27158); } double f_fm(double a, double x) { mpfr_set_d(r27137, a, MPFR_RNDN); mpfr_set_d(r27138, x, MPFR_RNDN); mpfr_mul(r27139, r27137, r27138, MPFR_RNDN); mpfr_exp(r27140, r27139, MPFR_RNDN); ; mpfr_sub(r27142, r27140, r27141, MPFR_RNDN); ; mpfr_set_si(r27144, mpfr_cmp(r27142, r27143) <= 0, MPFR_RNDN); mpfr_mul(r27145, r27138, r27137, MPFR_RNDN); mpfr_exp(r27146, r27145, MPFR_RNDN); mpfr_sub(r27147, r27146, r27141, MPFR_RNDN); ; mpfr_pow(r27149, r27147, r27148, MPFR_RNDN); mpfr_cbrt(r27150, r27149, MPFR_RNDN); mpfr_mul(r27151, r27139, r27139, MPFR_RNDN); ; mpfr_mul(r27153, r27139, r27152, MPFR_RNDN); ; mpfr_add(r27155, r27153, r27154, MPFR_RNDN); mpfr_mul(r27156, r27151, r27155, MPFR_RNDN); mpfr_add(r27157, r27156, r27139, MPFR_RNDN); if (mpfr_get_si(r27144, MPFR_RNDN)) { mpfr_set(r27158, r27150, MPFR_RNDN); } else { mpfr_set(r27158, r27157, MPFR_RNDN); }; return mpfr_get_d(r27158, MPFR_RNDN); } static mpfr_t r27159, r27160, r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175, r27176, r27177, r27178, r27179, r27180; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); mpfr_init(r27162); mpfr_init_set_str(r27163, "1", 10, MPFR_RNDN); mpfr_init(r27164); mpfr_init_set_str(r27165, "-0.0008112531176381211", 10, MPFR_RNDN); mpfr_init(r27166); mpfr_init(r27167); mpfr_init(r27168); mpfr_init(r27169); mpfr_init_set_str(r27170, "3", 10, MPFR_RNDN); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init_set_str(r27174, "1/6", 10, MPFR_RNDN); mpfr_init(r27175); mpfr_init_set_str(r27176, "1/2", 10, MPFR_RNDN); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); mpfr_init(r27180); } double f_dm(double a, double x) { mpfr_set_d(r27159, a, MPFR_RNDN); mpfr_set_d(r27160, x, MPFR_RNDN); mpfr_mul(r27161, r27159, r27160, MPFR_RNDN); mpfr_exp(r27162, r27161, MPFR_RNDN); ; mpfr_sub(r27164, r27162, r27163, MPFR_RNDN); ; mpfr_set_si(r27166, mpfr_cmp(r27164, r27165) <= 0, MPFR_RNDN); mpfr_mul(r27167, r27160, r27159, MPFR_RNDN); mpfr_exp(r27168, r27167, MPFR_RNDN); mpfr_sub(r27169, r27168, r27163, MPFR_RNDN); ; mpfr_pow(r27171, r27169, r27170, MPFR_RNDN); mpfr_cbrt(r27172, r27171, MPFR_RNDN); mpfr_mul(r27173, r27161, r27161, MPFR_RNDN); ; mpfr_mul(r27175, r27161, r27174, MPFR_RNDN); ; mpfr_add(r27177, r27175, r27176, MPFR_RNDN); mpfr_mul(r27178, r27173, r27177, MPFR_RNDN); mpfr_add(r27179, r27178, r27161, MPFR_RNDN); if (mpfr_get_si(r27166, MPFR_RNDN)) { mpfr_set(r27180, r27172, MPFR_RNDN); } else { mpfr_set(r27180, r27179, MPFR_RNDN); }; return mpfr_get_d(r27180, MPFR_RNDN); }