#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (26)"; double f_if(float a, float b) { float r25369 = a; float r25370 = r25369 * r25369; float r25371 = b; float r25372 = r25371 * r25371; float r25373 = r25370 + r25372; float r25374 = 2; float r25375 = pow(r25373, r25374); float r25376 = 4; float r25377 = r25376 * r25372; float r25378 = r25375 + r25377; float r25379 = 1; float r25380 = r25378 - r25379; return r25380; } double f_id(double a, double b) { double r25381 = a; double r25382 = r25381 * r25381; double r25383 = b; double r25384 = r25383 * r25383; double r25385 = r25382 + r25384; double r25386 = 2; double r25387 = pow(r25385, r25386); double r25388 = 4; double r25389 = r25388 * r25384; double r25390 = r25387 + r25389; double r25391 = 1; double r25392 = r25390 - r25391; return r25392; } double f_of(float a, float b) { float r25393 = b; float r25394 = r25393 * r25393; float r25395 = a; float r25396 = r25395 * r25395; float r25397 = r25394 + r25396; float r25398 = cbrt(r25397); float r25399 = r25397 * r25398; float r25400 = r25396 + r25394; float r25401 = cbrt(r25400); float r25402 = 2; float r25403 = pow(r25401, r25402); float r25404 = r25399 * r25403; float r25405 = 4; float r25406 = r25405 * r25394; float r25407 = r25404 + r25406; float r25408 = 1; float r25409 = r25407 - r25408; return r25409; } double f_od(double a, double b) { double r25410 = b; double r25411 = r25410 * r25410; double r25412 = a; double r25413 = r25412 * r25412; double r25414 = r25411 + r25413; double r25415 = cbrt(r25414); double r25416 = r25414 * r25415; double r25417 = r25413 + r25411; double r25418 = cbrt(r25417); double r25419 = 2; double r25420 = pow(r25418, r25419); double r25421 = r25416 * r25420; double r25422 = 4; double r25423 = r25422 * r25411; double r25424 = r25421 + r25423; double r25425 = 1; double r25426 = r25424 - r25425; return r25426; } 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 r25427, r25428, r25429, r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25427); mpfr_init(r25428); mpfr_init(r25429); mpfr_init(r25430); mpfr_init(r25431); mpfr_init_set_str(r25432, "2", 10, MPFR_RNDN); mpfr_init(r25433); mpfr_init_set_str(r25434, "4", 10, MPFR_RNDN); mpfr_init(r25435); mpfr_init(r25436); mpfr_init_set_str(r25437, "1", 10, MPFR_RNDN); mpfr_init(r25438); } double f_im(double a, double b) { mpfr_set_d(r25427, a, MPFR_RNDN); mpfr_mul(r25428, r25427, r25427, MPFR_RNDN); mpfr_set_d(r25429, b, MPFR_RNDN); mpfr_mul(r25430, r25429, r25429, MPFR_RNDN); mpfr_add(r25431, r25428, r25430, MPFR_RNDN); ; mpfr_pow(r25433, r25431, r25432, MPFR_RNDN); ; mpfr_mul(r25435, r25434, r25430, MPFR_RNDN); mpfr_add(r25436, r25433, r25435, MPFR_RNDN); ; mpfr_sub(r25438, r25436, r25437, MPFR_RNDN); return mpfr_get_d(r25438, MPFR_RNDN); } static mpfr_t r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25439); mpfr_init(r25440); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init_set_str(r25448, "2", 10, MPFR_RNDN); mpfr_init(r25449); mpfr_init(r25450); mpfr_init_set_str(r25451, "4", 10, MPFR_RNDN); mpfr_init(r25452); mpfr_init(r25453); mpfr_init_set_str(r25454, "1", 10, MPFR_RNDN); mpfr_init(r25455); } double f_fm(double a, double b) { mpfr_set_d(r25439, b, MPFR_RNDN); mpfr_mul(r25440, r25439, r25439, MPFR_RNDN); mpfr_set_d(r25441, a, MPFR_RNDN); mpfr_mul(r25442, r25441, r25441, MPFR_RNDN); mpfr_add(r25443, r25440, r25442, MPFR_RNDN); mpfr_cbrt(r25444, r25443, MPFR_RNDN); mpfr_mul(r25445, r25443, r25444, MPFR_RNDN); mpfr_add(r25446, r25442, r25440, MPFR_RNDN); mpfr_cbrt(r25447, r25446, MPFR_RNDN); ; mpfr_pow(r25449, r25447, r25448, MPFR_RNDN); mpfr_mul(r25450, r25445, r25449, MPFR_RNDN); ; mpfr_mul(r25452, r25451, r25440, MPFR_RNDN); mpfr_add(r25453, r25450, r25452, MPFR_RNDN); ; mpfr_sub(r25455, r25453, r25454, MPFR_RNDN); return mpfr_get_d(r25455, MPFR_RNDN); } static mpfr_t r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25456); mpfr_init(r25457); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); mpfr_init(r25464); mpfr_init_set_str(r25465, "2", 10, MPFR_RNDN); mpfr_init(r25466); mpfr_init(r25467); mpfr_init_set_str(r25468, "4", 10, MPFR_RNDN); mpfr_init(r25469); mpfr_init(r25470); mpfr_init_set_str(r25471, "1", 10, MPFR_RNDN); mpfr_init(r25472); } double f_dm(double a, double b) { mpfr_set_d(r25456, b, MPFR_RNDN); mpfr_mul(r25457, r25456, r25456, MPFR_RNDN); mpfr_set_d(r25458, a, MPFR_RNDN); mpfr_mul(r25459, r25458, r25458, MPFR_RNDN); mpfr_add(r25460, r25457, r25459, MPFR_RNDN); mpfr_cbrt(r25461, r25460, MPFR_RNDN); mpfr_mul(r25462, r25460, r25461, MPFR_RNDN); mpfr_add(r25463, r25459, r25457, MPFR_RNDN); mpfr_cbrt(r25464, r25463, MPFR_RNDN); ; mpfr_pow(r25466, r25464, r25465, MPFR_RNDN); mpfr_mul(r25467, r25462, r25466, MPFR_RNDN); ; mpfr_mul(r25469, r25468, r25457, MPFR_RNDN); mpfr_add(r25470, r25467, r25469, MPFR_RNDN); ; mpfr_sub(r25472, r25470, r25471, MPFR_RNDN); return mpfr_get_d(r25472, MPFR_RNDN); }