#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r25269 = a; float r25270 = r25269 * r25269; float r25271 = b; float r25272 = r25271 * r25271; float r25273 = r25270 + r25272; float r25274 = 2; float r25275 = pow(r25273, r25274); float r25276 = 4; float r25277 = 1; float r25278 = r25277 + r25269; float r25279 = r25270 * r25278; float r25280 = 3; float r25281 = r25280 * r25269; float r25282 = r25277 - r25281; float r25283 = r25272 * r25282; float r25284 = r25279 + r25283; float r25285 = r25276 * r25284; float r25286 = r25275 + r25285; float r25287 = r25286 - r25277; return r25287; } double f_id(double a, double b) { double r25288 = a; double r25289 = r25288 * r25288; double r25290 = b; double r25291 = r25290 * r25290; double r25292 = r25289 + r25291; double r25293 = 2; double r25294 = pow(r25292, r25293); double r25295 = 4; double r25296 = 1; double r25297 = r25296 + r25288; double r25298 = r25289 * r25297; double r25299 = 3; double r25300 = r25299 * r25288; double r25301 = r25296 - r25300; double r25302 = r25291 * r25301; double r25303 = r25298 + r25302; double r25304 = r25295 * r25303; double r25305 = r25294 + r25304; double r25306 = r25305 - r25296; return r25306; } double f_of(float a, float b) { float r25307 = a; float r25308 = r25307 * r25307; float r25309 = b; float r25310 = r25309 * r25309; float r25311 = r25308 + r25310; float r25312 = 2; float r25313 = pow(r25311, r25312); float r25314 = 4; float r25315 = 1; float r25316 = r25315 + r25307; float r25317 = cbrt(r25316); float r25318 = r25317 * r25317; float r25319 = r25308 * r25318; float r25320 = r25319 * r25317; float r25321 = 3; float r25322 = r25321 * r25307; float r25323 = r25315 - r25322; float r25324 = r25310 * r25323; float r25325 = r25320 + r25324; float r25326 = r25314 * r25325; float r25327 = r25313 + r25326; float r25328 = r25327 - r25315; return r25328; } double f_od(double a, double b) { double r25329 = a; double r25330 = r25329 * r25329; double r25331 = b; double r25332 = r25331 * r25331; double r25333 = r25330 + r25332; double r25334 = 2; double r25335 = pow(r25333, r25334); double r25336 = 4; double r25337 = 1; double r25338 = r25337 + r25329; double r25339 = cbrt(r25338); double r25340 = r25339 * r25339; double r25341 = r25330 * r25340; double r25342 = r25341 * r25339; double r25343 = 3; double r25344 = r25343 * r25329; double r25345 = r25337 - r25344; double r25346 = r25332 * r25345; double r25347 = r25342 + r25346; double r25348 = r25336 * r25347; double r25349 = r25335 + r25348; double r25350 = r25349 - r25337; return r25350; } 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 r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364, r25365, r25366, r25367, r25368, r25369; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init_set_str(r25356, "2", 10, MPFR_RNDN); mpfr_init(r25357); mpfr_init_set_str(r25358, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25359, "1", 10, MPFR_RNDN); mpfr_init(r25360); mpfr_init(r25361); mpfr_init_set_str(r25362, "3", 10, MPFR_RNDN); mpfr_init(r25363); mpfr_init(r25364); mpfr_init(r25365); mpfr_init(r25366); mpfr_init(r25367); mpfr_init(r25368); mpfr_init(r25369); } double f_im(double a, double b) { mpfr_set_d(r25351, a, MPFR_RNDN); mpfr_mul(r25352, r25351, r25351, MPFR_RNDN); mpfr_set_d(r25353, b, MPFR_RNDN); mpfr_mul(r25354, r25353, r25353, MPFR_RNDN); mpfr_add(r25355, r25352, r25354, MPFR_RNDN); ; mpfr_pow(r25357, r25355, r25356, MPFR_RNDN); ; ; mpfr_add(r25360, r25359, r25351, MPFR_RNDN); mpfr_mul(r25361, r25352, r25360, MPFR_RNDN); ; mpfr_mul(r25363, r25362, r25351, MPFR_RNDN); mpfr_sub(r25364, r25359, r25363, MPFR_RNDN); mpfr_mul(r25365, r25354, r25364, MPFR_RNDN); mpfr_add(r25366, r25361, r25365, MPFR_RNDN); mpfr_mul(r25367, r25358, r25366, MPFR_RNDN); mpfr_add(r25368, r25357, r25367, MPFR_RNDN); mpfr_sub(r25369, r25368, r25359, MPFR_RNDN); return mpfr_get_d(r25369, MPFR_RNDN); } static mpfr_t r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r25370); mpfr_init(r25371); mpfr_init(r25372); mpfr_init(r25373); mpfr_init(r25374); mpfr_init_set_str(r25375, "2", 10, MPFR_RNDN); mpfr_init(r25376); mpfr_init_set_str(r25377, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25378, "1", 10, MPFR_RNDN); mpfr_init(r25379); mpfr_init(r25380); mpfr_init(r25381); mpfr_init(r25382); mpfr_init(r25383); mpfr_init_set_str(r25384, "3", 10, MPFR_RNDN); mpfr_init(r25385); mpfr_init(r25386); mpfr_init(r25387); mpfr_init(r25388); mpfr_init(r25389); mpfr_init(r25390); mpfr_init(r25391); } double f_fm(double a, double b) { mpfr_set_d(r25370, a, MPFR_RNDN); mpfr_mul(r25371, r25370, r25370, MPFR_RNDN); mpfr_set_d(r25372, b, MPFR_RNDN); mpfr_mul(r25373, r25372, r25372, MPFR_RNDN); mpfr_add(r25374, r25371, r25373, MPFR_RNDN); ; mpfr_pow(r25376, r25374, r25375, MPFR_RNDN); ; ; mpfr_add(r25379, r25378, r25370, MPFR_RNDN); mpfr_cbrt(r25380, r25379, MPFR_RNDN); mpfr_mul(r25381, r25380, r25380, MPFR_RNDN); mpfr_mul(r25382, r25371, r25381, MPFR_RNDN); mpfr_mul(r25383, r25382, r25380, MPFR_RNDN); ; mpfr_mul(r25385, r25384, r25370, MPFR_RNDN); mpfr_sub(r25386, r25378, r25385, MPFR_RNDN); mpfr_mul(r25387, r25373, r25386, MPFR_RNDN); mpfr_add(r25388, r25383, r25387, MPFR_RNDN); mpfr_mul(r25389, r25377, r25388, MPFR_RNDN); mpfr_add(r25390, r25376, r25389, MPFR_RNDN); mpfr_sub(r25391, r25390, r25378, MPFR_RNDN); return mpfr_get_d(r25391, MPFR_RNDN); } static mpfr_t r25392, r25393, r25394, r25395, r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r25392); mpfr_init(r25393); mpfr_init(r25394); mpfr_init(r25395); mpfr_init(r25396); mpfr_init_set_str(r25397, "2", 10, MPFR_RNDN); mpfr_init(r25398); mpfr_init_set_str(r25399, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25400, "1", 10, MPFR_RNDN); mpfr_init(r25401); mpfr_init(r25402); mpfr_init(r25403); mpfr_init(r25404); mpfr_init(r25405); mpfr_init_set_str(r25406, "3", 10, MPFR_RNDN); mpfr_init(r25407); mpfr_init(r25408); mpfr_init(r25409); mpfr_init(r25410); mpfr_init(r25411); mpfr_init(r25412); mpfr_init(r25413); } double f_dm(double a, double b) { mpfr_set_d(r25392, a, MPFR_RNDN); mpfr_mul(r25393, r25392, r25392, MPFR_RNDN); mpfr_set_d(r25394, b, MPFR_RNDN); mpfr_mul(r25395, r25394, r25394, MPFR_RNDN); mpfr_add(r25396, r25393, r25395, MPFR_RNDN); ; mpfr_pow(r25398, r25396, r25397, MPFR_RNDN); ; ; mpfr_add(r25401, r25400, r25392, MPFR_RNDN); mpfr_cbrt(r25402, r25401, MPFR_RNDN); mpfr_mul(r25403, r25402, r25402, MPFR_RNDN); mpfr_mul(r25404, r25393, r25403, MPFR_RNDN); mpfr_mul(r25405, r25404, r25402, MPFR_RNDN); ; mpfr_mul(r25407, r25406, r25392, MPFR_RNDN); mpfr_sub(r25408, r25400, r25407, MPFR_RNDN); mpfr_mul(r25409, r25395, r25408, MPFR_RNDN); mpfr_add(r25410, r25405, r25409, MPFR_RNDN); mpfr_mul(r25411, r25399, r25410, MPFR_RNDN); mpfr_add(r25412, r25398, r25411, MPFR_RNDN); mpfr_sub(r25413, r25412, r25400, MPFR_RNDN); return mpfr_get_d(r25413, MPFR_RNDN); }