#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r25164 = a; float r25165 = r25164 * r25164; float r25166 = b; float r25167 = r25166 * r25166; float r25168 = r25165 + r25167; float r25169 = 2; float r25170 = pow(r25168, r25169); float r25171 = 4; float r25172 = 1; float r25173 = r25172 + r25164; float r25174 = r25165 * r25173; float r25175 = 3; float r25176 = r25175 * r25164; float r25177 = r25172 - r25176; float r25178 = r25167 * r25177; float r25179 = r25174 + r25178; float r25180 = r25171 * r25179; float r25181 = r25170 + r25180; float r25182 = r25181 - r25172; return r25182; } double f_id(double a, double b) { double r25183 = a; double r25184 = r25183 * r25183; double r25185 = b; double r25186 = r25185 * r25185; double r25187 = r25184 + r25186; double r25188 = 2; double r25189 = pow(r25187, r25188); double r25190 = 4; double r25191 = 1; double r25192 = r25191 + r25183; double r25193 = r25184 * r25192; double r25194 = 3; double r25195 = r25194 * r25183; double r25196 = r25191 - r25195; double r25197 = r25186 * r25196; double r25198 = r25193 + r25197; double r25199 = r25190 * r25198; double r25200 = r25189 + r25199; double r25201 = r25200 - r25191; return r25201; } double f_of(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 = 1; float r25211 = r25210 + r25202; float r25212 = cbrt(r25211); float r25213 = r25212 * r25212; float r25214 = r25203 * r25213; float r25215 = r25214 * r25212; float r25216 = 3; float r25217 = r25216 * r25202; float r25218 = r25210 - r25217; float r25219 = r25205 * r25218; float r25220 = r25215 + r25219; float r25221 = r25209 * r25220; float r25222 = r25208 + r25221; float r25223 = r25222 - r25210; return r25223; } double f_od(double a, double b) { double r25224 = a; double r25225 = r25224 * r25224; double r25226 = b; double r25227 = r25226 * r25226; double r25228 = r25225 + r25227; double r25229 = 2; double r25230 = pow(r25228, r25229); double r25231 = 4; double r25232 = 1; double r25233 = r25232 + r25224; double r25234 = cbrt(r25233); double r25235 = r25234 * r25234; double r25236 = r25225 * r25235; double r25237 = r25236 * r25234; double r25238 = 3; double r25239 = r25238 * r25224; double r25240 = r25232 - r25239; double r25241 = r25227 * r25240; double r25242 = r25237 + r25241; double r25243 = r25231 * r25242; double r25244 = r25230 + r25243; double r25245 = r25244 - r25232; return r25245; } 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 r25246, r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r25246); mpfr_init(r25247); mpfr_init(r25248); mpfr_init(r25249); mpfr_init(r25250); mpfr_init_set_str(r25251, "2", 10, MPFR_RNDN); mpfr_init(r25252); mpfr_init_set_str(r25253, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25254, "1", 10, MPFR_RNDN); mpfr_init(r25255); mpfr_init(r25256); mpfr_init_set_str(r25257, "3", 10, MPFR_RNDN); mpfr_init(r25258); mpfr_init(r25259); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init(r25264); } double f_im(double a, double b) { mpfr_set_d(r25246, a, MPFR_RNDN); mpfr_mul(r25247, r25246, r25246, MPFR_RNDN); mpfr_set_d(r25248, b, MPFR_RNDN); mpfr_mul(r25249, r25248, r25248, MPFR_RNDN); mpfr_add(r25250, r25247, r25249, MPFR_RNDN); ; mpfr_pow(r25252, r25250, r25251, MPFR_RNDN); ; ; mpfr_add(r25255, r25254, r25246, MPFR_RNDN); mpfr_mul(r25256, r25247, r25255, MPFR_RNDN); ; mpfr_mul(r25258, r25257, r25246, MPFR_RNDN); mpfr_sub(r25259, r25254, r25258, MPFR_RNDN); mpfr_mul(r25260, r25249, r25259, MPFR_RNDN); mpfr_add(r25261, r25256, r25260, MPFR_RNDN); mpfr_mul(r25262, r25253, r25261, MPFR_RNDN); mpfr_add(r25263, r25252, r25262, MPFR_RNDN); mpfr_sub(r25264, r25263, r25254, MPFR_RNDN); return mpfr_get_d(r25264, MPFR_RNDN); } static mpfr_t r25265, r25266, r25267, r25268, r25269, r25270, r25271, r25272, r25273, r25274, r25275, r25276, r25277, r25278, r25279, r25280, r25281, r25282, r25283, r25284, r25285, r25286; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r25265); mpfr_init(r25266); mpfr_init(r25267); mpfr_init(r25268); mpfr_init(r25269); mpfr_init_set_str(r25270, "2", 10, MPFR_RNDN); mpfr_init(r25271); mpfr_init_set_str(r25272, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25273, "1", 10, MPFR_RNDN); mpfr_init(r25274); mpfr_init(r25275); mpfr_init(r25276); mpfr_init(r25277); mpfr_init(r25278); mpfr_init_set_str(r25279, "3", 10, MPFR_RNDN); mpfr_init(r25280); mpfr_init(r25281); mpfr_init(r25282); mpfr_init(r25283); mpfr_init(r25284); mpfr_init(r25285); mpfr_init(r25286); } double f_fm(double a, double b) { mpfr_set_d(r25265, a, MPFR_RNDN); mpfr_mul(r25266, r25265, r25265, MPFR_RNDN); mpfr_set_d(r25267, b, MPFR_RNDN); mpfr_mul(r25268, r25267, r25267, MPFR_RNDN); mpfr_add(r25269, r25266, r25268, MPFR_RNDN); ; mpfr_pow(r25271, r25269, r25270, MPFR_RNDN); ; ; mpfr_add(r25274, r25273, r25265, MPFR_RNDN); mpfr_cbrt(r25275, r25274, MPFR_RNDN); mpfr_mul(r25276, r25275, r25275, MPFR_RNDN); mpfr_mul(r25277, r25266, r25276, MPFR_RNDN); mpfr_mul(r25278, r25277, r25275, MPFR_RNDN); ; mpfr_mul(r25280, r25279, r25265, MPFR_RNDN); mpfr_sub(r25281, r25273, r25280, MPFR_RNDN); mpfr_mul(r25282, r25268, r25281, MPFR_RNDN); mpfr_add(r25283, r25278, r25282, MPFR_RNDN); mpfr_mul(r25284, r25272, r25283, MPFR_RNDN); mpfr_add(r25285, r25271, r25284, MPFR_RNDN); mpfr_sub(r25286, r25285, r25273, MPFR_RNDN); return mpfr_get_d(r25286, MPFR_RNDN); } static mpfr_t r25287, r25288, r25289, r25290, r25291, r25292, r25293, r25294, r25295, r25296, r25297, r25298, r25299, r25300, r25301, r25302, r25303, r25304, r25305, r25306, r25307, r25308; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r25287); mpfr_init(r25288); mpfr_init(r25289); mpfr_init(r25290); mpfr_init(r25291); mpfr_init_set_str(r25292, "2", 10, MPFR_RNDN); mpfr_init(r25293); mpfr_init_set_str(r25294, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25295, "1", 10, MPFR_RNDN); mpfr_init(r25296); mpfr_init(r25297); mpfr_init(r25298); mpfr_init(r25299); mpfr_init(r25300); mpfr_init_set_str(r25301, "3", 10, MPFR_RNDN); mpfr_init(r25302); mpfr_init(r25303); mpfr_init(r25304); mpfr_init(r25305); mpfr_init(r25306); mpfr_init(r25307); mpfr_init(r25308); } double f_dm(double a, double b) { mpfr_set_d(r25287, a, MPFR_RNDN); mpfr_mul(r25288, r25287, r25287, MPFR_RNDN); mpfr_set_d(r25289, b, MPFR_RNDN); mpfr_mul(r25290, r25289, r25289, MPFR_RNDN); mpfr_add(r25291, r25288, r25290, MPFR_RNDN); ; mpfr_pow(r25293, r25291, r25292, MPFR_RNDN); ; ; mpfr_add(r25296, r25295, r25287, MPFR_RNDN); mpfr_cbrt(r25297, r25296, MPFR_RNDN); mpfr_mul(r25298, r25297, r25297, MPFR_RNDN); mpfr_mul(r25299, r25288, r25298, MPFR_RNDN); mpfr_mul(r25300, r25299, r25297, MPFR_RNDN); ; mpfr_mul(r25302, r25301, r25287, MPFR_RNDN); mpfr_sub(r25303, r25295, r25302, MPFR_RNDN); mpfr_mul(r25304, r25290, r25303, MPFR_RNDN); mpfr_add(r25305, r25300, r25304, MPFR_RNDN); mpfr_mul(r25306, r25294, r25305, MPFR_RNDN); mpfr_add(r25307, r25293, r25306, MPFR_RNDN); mpfr_sub(r25308, r25307, r25295, MPFR_RNDN); return mpfr_get_d(r25308, MPFR_RNDN); }