#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (26)"; double f_if(float a, float b) { float r25202 = a; float r25203 = r25202 * r25202; float r25204 = b; float r25205 = r25204 * r25204; float r25206 = r25203 + r25205; float r25207 = 2; float r25208 = pow(r25206, r25207); float r25209 = 4; float r25210 = r25209 * r25205; float r25211 = r25208 + r25210; float r25212 = 1; float r25213 = r25211 - r25212; return r25213; } double f_id(double a, double b) { double r25214 = a; double r25215 = r25214 * r25214; double r25216 = b; double r25217 = r25216 * r25216; double r25218 = r25215 + r25217; double r25219 = 2; double r25220 = pow(r25218, r25219); double r25221 = 4; double r25222 = r25221 * r25217; double r25223 = r25220 + r25222; double r25224 = 1; double r25225 = r25223 - r25224; return r25225; } double f_of(float a, float b) { float r25226 = b; float r25227 = a; float r25228 = r25227 * r25227; float r25229 = fma(r25226, r25226, r25228); float r25230 = cbrt(r25229); float r25231 = r25229 * r25230; float r25232 = r25226 * r25226; float r25233 = r25228 + r25232; float r25234 = cbrt(r25233); float r25235 = 2; float r25236 = pow(r25234, r25235); float r25237 = r25231 * r25236; float r25238 = 4; float r25239 = r25238 * r25232; float r25240 = r25237 + r25239; float r25241 = 1; float r25242 = r25240 - r25241; return r25242; } double f_od(double a, double b) { double r25243 = b; double r25244 = a; double r25245 = r25244 * r25244; double r25246 = fma(r25243, r25243, r25245); double r25247 = cbrt(r25246); double r25248 = r25246 * r25247; double r25249 = r25243 * r25243; double r25250 = r25245 + r25249; double r25251 = cbrt(r25250); double r25252 = 2; double r25253 = pow(r25251, r25252); double r25254 = r25248 * r25253; double r25255 = 4; double r25256 = r25255 * r25249; double r25257 = r25254 + r25256; double r25258 = 1; double r25259 = r25257 - r25258; return r25259; } 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 r25260, r25261, r25262, r25263, r25264, r25265, r25266, r25267, r25268, r25269, r25270, r25271; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init(r25264); mpfr_init_set_str(r25265, "2", 10, MPFR_RNDN); mpfr_init(r25266); mpfr_init_set_str(r25267, "4", 10, MPFR_RNDN); mpfr_init(r25268); mpfr_init(r25269); mpfr_init_set_str(r25270, "1", 10, MPFR_RNDN); mpfr_init(r25271); } double f_im(double a, double b) { mpfr_set_d(r25260, a, MPFR_RNDN); mpfr_mul(r25261, r25260, r25260, MPFR_RNDN); mpfr_set_d(r25262, b, MPFR_RNDN); mpfr_mul(r25263, r25262, r25262, MPFR_RNDN); mpfr_add(r25264, r25261, r25263, MPFR_RNDN); ; mpfr_pow(r25266, r25264, r25265, MPFR_RNDN); ; mpfr_mul(r25268, r25267, r25263, MPFR_RNDN); mpfr_add(r25269, r25266, r25268, MPFR_RNDN); ; mpfr_sub(r25271, r25269, r25270, MPFR_RNDN); return mpfr_get_d(r25271, MPFR_RNDN); } static mpfr_t r25272, r25273, r25274, r25275, r25276, r25277, r25278, r25279, r25280, r25281, r25282, r25283, r25284, r25285, r25286, r25287, r25288; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25272); mpfr_init(r25273); mpfr_init(r25274); mpfr_init(r25275); mpfr_init(r25276); mpfr_init(r25277); mpfr_init(r25278); mpfr_init(r25279); mpfr_init(r25280); mpfr_init_set_str(r25281, "2", 10, MPFR_RNDN); mpfr_init(r25282); mpfr_init(r25283); mpfr_init_set_str(r25284, "4", 10, MPFR_RNDN); mpfr_init(r25285); mpfr_init(r25286); mpfr_init_set_str(r25287, "1", 10, MPFR_RNDN); mpfr_init(r25288); } double f_fm(double a, double b) { mpfr_set_d(r25272, b, MPFR_RNDN); mpfr_set_d(r25273, a, MPFR_RNDN); mpfr_mul(r25274, r25273, r25273, MPFR_RNDN); mpfr_fma(r25275, r25272, r25272, r25274, MPFR_RNDN); mpfr_cbrt(r25276, r25275, MPFR_RNDN); mpfr_mul(r25277, r25275, r25276, MPFR_RNDN); mpfr_mul(r25278, r25272, r25272, MPFR_RNDN); mpfr_add(r25279, r25274, r25278, MPFR_RNDN); mpfr_cbrt(r25280, r25279, MPFR_RNDN); ; mpfr_pow(r25282, r25280, r25281, MPFR_RNDN); mpfr_mul(r25283, r25277, r25282, MPFR_RNDN); ; mpfr_mul(r25285, r25284, r25278, MPFR_RNDN); mpfr_add(r25286, r25283, r25285, MPFR_RNDN); ; mpfr_sub(r25288, r25286, r25287, MPFR_RNDN); return mpfr_get_d(r25288, MPFR_RNDN); } static mpfr_t r25289, r25290, r25291, r25292, r25293, r25294, r25295, r25296, r25297, r25298, r25299, r25300, r25301, r25302, r25303, r25304, r25305; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25289); mpfr_init(r25290); mpfr_init(r25291); mpfr_init(r25292); mpfr_init(r25293); mpfr_init(r25294); mpfr_init(r25295); mpfr_init(r25296); mpfr_init(r25297); mpfr_init_set_str(r25298, "2", 10, MPFR_RNDN); mpfr_init(r25299); mpfr_init(r25300); mpfr_init_set_str(r25301, "4", 10, MPFR_RNDN); mpfr_init(r25302); mpfr_init(r25303); mpfr_init_set_str(r25304, "1", 10, MPFR_RNDN); mpfr_init(r25305); } double f_dm(double a, double b) { mpfr_set_d(r25289, b, MPFR_RNDN); mpfr_set_d(r25290, a, MPFR_RNDN); mpfr_mul(r25291, r25290, r25290, MPFR_RNDN); mpfr_fma(r25292, r25289, r25289, r25291, MPFR_RNDN); mpfr_cbrt(r25293, r25292, MPFR_RNDN); mpfr_mul(r25294, r25292, r25293, MPFR_RNDN); mpfr_mul(r25295, r25289, r25289, MPFR_RNDN); mpfr_add(r25296, r25291, r25295, MPFR_RNDN); mpfr_cbrt(r25297, r25296, MPFR_RNDN); ; mpfr_pow(r25299, r25297, r25298, MPFR_RNDN); mpfr_mul(r25300, r25294, r25299, MPFR_RNDN); ; mpfr_mul(r25302, r25301, r25295, MPFR_RNDN); mpfr_add(r25303, r25300, r25302, MPFR_RNDN); ; mpfr_sub(r25305, r25303, r25304, MPFR_RNDN); return mpfr_get_d(r25305, MPFR_RNDN); }