#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17254 = a; float r17255 = r17254 * r17254; float r17256 = b; float r17257 = r17256 * r17256; float r17258 = r17255 + r17257; float r17259 = r17258 * r17258; float r17260 = 4.0f; float r17261 = 1.0f; float r17262 = r17261 + r17254; float r17263 = r17255 * r17262; float r17264 = 3.0f; float r17265 = r17264 * r17254; float r17266 = r17261 - r17265; float r17267 = r17257 * r17266; float r17268 = r17263 + r17267; float r17269 = r17260 * r17268; float r17270 = r17259 + r17269; float r17271 = r17270 - r17261; return r17271; } double f_id(double a, double b) { double r17272 = a; double r17273 = r17272 * r17272; double r17274 = b; double r17275 = r17274 * r17274; double r17276 = r17273 + r17275; double r17277 = r17276 * r17276; double r17278 = 4.0; double r17279 = 1.0; double r17280 = r17279 + r17272; double r17281 = r17273 * r17280; double r17282 = 3.0; double r17283 = r17282 * r17272; double r17284 = r17279 - r17283; double r17285 = r17275 * r17284; double r17286 = r17281 + r17285; double r17287 = r17278 * r17286; double r17288 = r17277 + r17287; double r17289 = r17288 - r17279; return r17289; } double f_of(float a, float b) { float r17290 = 1.0f; float r17291 = a; float r17292 = 3.0f; float r17293 = r17291 * r17292; float r17294 = r17290 - r17293; float r17295 = b; float r17296 = r17295 * r17295; float r17297 = fma(r17291, r17291, r17291); float r17298 = r17297 * r17291; float r17299 = fma(r17294, r17296, r17298); float r17300 = 4.0f; float r17301 = 2.0f; float r17302 = r17295 * r17301; float r17303 = r17302 * r17295; float r17304 = r17291 * r17291; float r17305 = pow(r17291, r17300); float r17306 = pow(r17295, r17300); float r17307 = r17305 + r17306; float r17308 = fma(r17303, r17304, r17307); float r17309 = fma(r17299, r17300, r17308); float r17310 = r17309 - r17290; return r17310; } double f_od(double a, double b) { double r17311 = 1.0; double r17312 = a; double r17313 = 3.0; double r17314 = r17312 * r17313; double r17315 = r17311 - r17314; double r17316 = b; double r17317 = r17316 * r17316; double r17318 = fma(r17312, r17312, r17312); double r17319 = r17318 * r17312; double r17320 = fma(r17315, r17317, r17319); double r17321 = 4.0; double r17322 = 2.0; double r17323 = r17316 * r17322; double r17324 = r17323 * r17316; double r17325 = r17312 * r17312; double r17326 = pow(r17312, r17321); double r17327 = pow(r17316, r17321); double r17328 = r17326 + r17327; double r17329 = fma(r17324, r17325, r17328); double r17330 = fma(r17320, r17321, r17329); double r17331 = r17330 - r17311; return r17331; } 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 r17332, r17333, r17334, r17335, r17336, r17337, r17338, r17339, r17340, r17341, r17342, r17343, r17344, r17345, r17346, r17347, r17348, r17349; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17332); mpfr_init(r17333); mpfr_init(r17334); mpfr_init(r17335); mpfr_init(r17336); mpfr_init(r17337); mpfr_init_set_str(r17338, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17339, "1", 10, MPFR_RNDN); mpfr_init(r17340); mpfr_init(r17341); mpfr_init_set_str(r17342, "3", 10, MPFR_RNDN); mpfr_init(r17343); mpfr_init(r17344); mpfr_init(r17345); mpfr_init(r17346); mpfr_init(r17347); mpfr_init(r17348); mpfr_init(r17349); } double f_im(double a, double b) { mpfr_set_d(r17332, a, MPFR_RNDN); mpfr_sqr(r17333, r17332, MPFR_RNDN); mpfr_set_d(r17334, b, MPFR_RNDN); mpfr_sqr(r17335, r17334, MPFR_RNDN); mpfr_add(r17336, r17333, r17335, MPFR_RNDN); mpfr_sqr(r17337, r17336, MPFR_RNDN); ; ; mpfr_add(r17340, r17339, r17332, MPFR_RNDN); mpfr_mul(r17341, r17333, r17340, MPFR_RNDN); ; mpfr_mul(r17343, r17342, r17332, MPFR_RNDN); mpfr_sub(r17344, r17339, r17343, MPFR_RNDN); mpfr_mul(r17345, r17335, r17344, MPFR_RNDN); mpfr_add(r17346, r17341, r17345, MPFR_RNDN); mpfr_mul(r17347, r17338, r17346, MPFR_RNDN); mpfr_add(r17348, r17337, r17347, MPFR_RNDN); mpfr_sub(r17349, r17348, r17339, MPFR_RNDN); return mpfr_get_d(r17349, MPFR_RNDN); } static mpfr_t r17350, r17351, r17352, r17353, r17354, r17355, r17356, r17357, r17358, r17359, r17360, r17361, r17362, r17363, r17364, r17365, r17366, r17367, r17368, r17369, r17370; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17350, "1", 10, MPFR_RNDN); mpfr_init(r17351); mpfr_init_set_str(r17352, "3", 10, MPFR_RNDN); mpfr_init(r17353); mpfr_init(r17354); mpfr_init(r17355); mpfr_init(r17356); mpfr_init(r17357); mpfr_init(r17358); mpfr_init(r17359); mpfr_init_set_str(r17360, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17361, "2", 10, MPFR_RNDN); mpfr_init(r17362); mpfr_init(r17363); mpfr_init(r17364); mpfr_init(r17365); mpfr_init(r17366); mpfr_init(r17367); mpfr_init(r17368); mpfr_init(r17369); mpfr_init(r17370); } double f_fm(double a, double b) { ; mpfr_set_d(r17351, a, MPFR_RNDN); ; mpfr_mul(r17353, r17351, r17352, MPFR_RNDN); mpfr_sub(r17354, r17350, r17353, MPFR_RNDN); mpfr_set_d(r17355, b, MPFR_RNDN); mpfr_sqr(r17356, r17355, MPFR_RNDN); mpfr_fma(r17357, r17351, r17351, r17351, MPFR_RNDN); mpfr_mul(r17358, r17357, r17351, MPFR_RNDN); mpfr_fma(r17359, r17354, r17356, r17358, MPFR_RNDN); ; ; mpfr_mul(r17362, r17355, r17361, MPFR_RNDN); mpfr_mul(r17363, r17362, r17355, MPFR_RNDN); mpfr_mul(r17364, r17351, r17351, MPFR_RNDN); mpfr_pow(r17365, r17351, r17360, MPFR_RNDN); mpfr_pow(r17366, r17355, r17360, MPFR_RNDN); mpfr_add(r17367, r17365, r17366, MPFR_RNDN); mpfr_fma(r17368, r17363, r17364, r17367, MPFR_RNDN); mpfr_fma(r17369, r17359, r17360, r17368, MPFR_RNDN); mpfr_sub(r17370, r17369, r17350, MPFR_RNDN); return mpfr_get_d(r17370, MPFR_RNDN); } static mpfr_t r17371, r17372, r17373, r17374, r17375, r17376, r17377, r17378, r17379, r17380, r17381, r17382, r17383, r17384, r17385, r17386, r17387, r17388, r17389, r17390, r17391; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17371, "1", 10, MPFR_RNDN); mpfr_init(r17372); mpfr_init_set_str(r17373, "3", 10, MPFR_RNDN); mpfr_init(r17374); mpfr_init(r17375); mpfr_init(r17376); mpfr_init(r17377); mpfr_init(r17378); mpfr_init(r17379); mpfr_init(r17380); mpfr_init_set_str(r17381, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17382, "2", 10, MPFR_RNDN); mpfr_init(r17383); mpfr_init(r17384); mpfr_init(r17385); mpfr_init(r17386); mpfr_init(r17387); mpfr_init(r17388); mpfr_init(r17389); mpfr_init(r17390); mpfr_init(r17391); } double f_dm(double a, double b) { ; mpfr_set_d(r17372, a, MPFR_RNDN); ; mpfr_mul(r17374, r17372, r17373, MPFR_RNDN); mpfr_sub(r17375, r17371, r17374, MPFR_RNDN); mpfr_set_d(r17376, b, MPFR_RNDN); mpfr_sqr(r17377, r17376, MPFR_RNDN); mpfr_fma(r17378, r17372, r17372, r17372, MPFR_RNDN); mpfr_mul(r17379, r17378, r17372, MPFR_RNDN); mpfr_fma(r17380, r17375, r17377, r17379, MPFR_RNDN); ; ; mpfr_mul(r17383, r17376, r17382, MPFR_RNDN); mpfr_mul(r17384, r17383, r17376, MPFR_RNDN); mpfr_mul(r17385, r17372, r17372, MPFR_RNDN); mpfr_pow(r17386, r17372, r17381, MPFR_RNDN); mpfr_pow(r17387, r17376, r17381, MPFR_RNDN); mpfr_add(r17388, r17386, r17387, MPFR_RNDN); mpfr_fma(r17389, r17384, r17385, r17388, MPFR_RNDN); mpfr_fma(r17390, r17380, r17381, r17389, MPFR_RNDN); mpfr_sub(r17391, r17390, r17371, MPFR_RNDN); return mpfr_get_d(r17391, MPFR_RNDN); }