#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (24)"; double f_if(float a, float b) { float r25442 = a; float r25443 = r25442 * r25442; float r25444 = b; float r25445 = r25444 * r25444; float r25446 = r25443 + r25445; float r25447 = 2; float r25448 = pow(r25446, r25447); float r25449 = 4; float r25450 = 1; float r25451 = r25450 - r25442; float r25452 = r25443 * r25451; float r25453 = 3; float r25454 = r25453 + r25442; float r25455 = r25445 * r25454; float r25456 = r25452 + r25455; float r25457 = r25449 * r25456; float r25458 = r25448 + r25457; float r25459 = r25458 - r25450; return r25459; } double f_id(double a, double b) { double r25460 = a; double r25461 = r25460 * r25460; double r25462 = b; double r25463 = r25462 * r25462; double r25464 = r25461 + r25463; double r25465 = 2; double r25466 = pow(r25464, r25465); double r25467 = 4; double r25468 = 1; double r25469 = r25468 - r25460; double r25470 = r25461 * r25469; double r25471 = 3; double r25472 = r25471 + r25460; double r25473 = r25463 * r25472; double r25474 = r25470 + r25473; double r25475 = r25467 * r25474; double r25476 = r25466 + r25475; double r25477 = r25476 - r25468; return r25477; } double f_of(float a, float b) { float r25478 = a; float r25479 = r25478 * r25478; float r25480 = b; float r25481 = r25480 * r25480; float r25482 = r25479 + r25481; float r25483 = 2; float r25484 = pow(r25482, r25483); float r25485 = 4; float r25486 = 1; float r25487 = r25486 - r25479; float r25488 = r25478 + r25486; float r25489 = r25488 / r25479; float r25490 = r25487 / r25489; float r25491 = cbrt(r25490); float r25492 = r25491 * r25491; float r25493 = r25492 * r25491; float r25494 = 3; float r25495 = r25494 + r25478; float r25496 = r25481 * r25495; float r25497 = r25493 + r25496; float r25498 = r25485 * r25497; float r25499 = r25484 + r25498; float r25500 = r25499 - r25486; return r25500; } double f_od(double a, double b) { double r25501 = a; double r25502 = r25501 * r25501; double r25503 = b; double r25504 = r25503 * r25503; double r25505 = r25502 + r25504; double r25506 = 2; double r25507 = pow(r25505, r25506); double r25508 = 4; double r25509 = 1; double r25510 = r25509 - r25502; double r25511 = r25501 + r25509; double r25512 = r25511 / r25502; double r25513 = r25510 / r25512; double r25514 = cbrt(r25513); double r25515 = r25514 * r25514; double r25516 = r25515 * r25514; double r25517 = 3; double r25518 = r25517 + r25501; double r25519 = r25504 * r25518; double r25520 = r25516 + r25519; double r25521 = r25508 * r25520; double r25522 = r25507 + r25521; double r25523 = r25522 - r25509; return r25523; } 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 r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init_set_str(r25529, "2", 10, MPFR_RNDN); mpfr_init(r25530); mpfr_init_set_str(r25531, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25532, "1", 10, MPFR_RNDN); mpfr_init(r25533); mpfr_init(r25534); mpfr_init_set_str(r25535, "3", 10, MPFR_RNDN); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); } double f_im(double a, double b) { mpfr_set_d(r25524, a, MPFR_RNDN); mpfr_mul(r25525, r25524, r25524, MPFR_RNDN); mpfr_set_d(r25526, b, MPFR_RNDN); mpfr_mul(r25527, r25526, r25526, MPFR_RNDN); mpfr_add(r25528, r25525, r25527, MPFR_RNDN); ; mpfr_pow(r25530, r25528, r25529, MPFR_RNDN); ; ; mpfr_sub(r25533, r25532, r25524, MPFR_RNDN); mpfr_mul(r25534, r25525, r25533, MPFR_RNDN); ; mpfr_add(r25536, r25535, r25524, MPFR_RNDN); mpfr_mul(r25537, r25527, r25536, MPFR_RNDN); mpfr_add(r25538, r25534, r25537, MPFR_RNDN); mpfr_mul(r25539, r25531, r25538, MPFR_RNDN); mpfr_add(r25540, r25530, r25539, MPFR_RNDN); mpfr_sub(r25541, r25540, r25532, MPFR_RNDN); return mpfr_get_d(r25541, MPFR_RNDN); } static mpfr_t r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init_set_str(r25547, "2", 10, MPFR_RNDN); mpfr_init(r25548); mpfr_init_set_str(r25549, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25550, "1", 10, MPFR_RNDN); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init_set_str(r25558, "3", 10, MPFR_RNDN); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); } double f_fm(double a, double b) { mpfr_set_d(r25542, a, MPFR_RNDN); mpfr_mul(r25543, r25542, r25542, MPFR_RNDN); mpfr_set_d(r25544, b, MPFR_RNDN); mpfr_mul(r25545, r25544, r25544, MPFR_RNDN); mpfr_add(r25546, r25543, r25545, MPFR_RNDN); ; mpfr_pow(r25548, r25546, r25547, MPFR_RNDN); ; ; mpfr_sub(r25551, r25550, r25543, MPFR_RNDN); mpfr_add(r25552, r25542, r25550, MPFR_RNDN); mpfr_div(r25553, r25552, r25543, MPFR_RNDN); mpfr_div(r25554, r25551, r25553, MPFR_RNDN); mpfr_cbrt(r25555, r25554, MPFR_RNDN); mpfr_mul(r25556, r25555, r25555, MPFR_RNDN); mpfr_mul(r25557, r25556, r25555, MPFR_RNDN); ; mpfr_add(r25559, r25558, r25542, MPFR_RNDN); mpfr_mul(r25560, r25545, r25559, MPFR_RNDN); mpfr_add(r25561, r25557, r25560, MPFR_RNDN); mpfr_mul(r25562, r25549, r25561, MPFR_RNDN); mpfr_add(r25563, r25548, r25562, MPFR_RNDN); mpfr_sub(r25564, r25563, r25550, MPFR_RNDN); return mpfr_get_d(r25564, MPFR_RNDN); } static mpfr_t r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579, r25580, r25581, r25582, r25583, r25584, r25585, r25586, r25587; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init_set_str(r25570, "2", 10, MPFR_RNDN); mpfr_init(r25571); mpfr_init_set_str(r25572, "4", 10, MPFR_RNDN); mpfr_init_set_str(r25573, "1", 10, MPFR_RNDN); mpfr_init(r25574); mpfr_init(r25575); mpfr_init(r25576); mpfr_init(r25577); mpfr_init(r25578); mpfr_init(r25579); mpfr_init(r25580); mpfr_init_set_str(r25581, "3", 10, MPFR_RNDN); mpfr_init(r25582); mpfr_init(r25583); mpfr_init(r25584); mpfr_init(r25585); mpfr_init(r25586); mpfr_init(r25587); } double f_dm(double a, double b) { mpfr_set_d(r25565, a, MPFR_RNDN); mpfr_mul(r25566, r25565, r25565, MPFR_RNDN); mpfr_set_d(r25567, b, MPFR_RNDN); mpfr_mul(r25568, r25567, r25567, MPFR_RNDN); mpfr_add(r25569, r25566, r25568, MPFR_RNDN); ; mpfr_pow(r25571, r25569, r25570, MPFR_RNDN); ; ; mpfr_sub(r25574, r25573, r25566, MPFR_RNDN); mpfr_add(r25575, r25565, r25573, MPFR_RNDN); mpfr_div(r25576, r25575, r25566, MPFR_RNDN); mpfr_div(r25577, r25574, r25576, MPFR_RNDN); mpfr_cbrt(r25578, r25577, MPFR_RNDN); mpfr_mul(r25579, r25578, r25578, MPFR_RNDN); mpfr_mul(r25580, r25579, r25578, MPFR_RNDN); ; mpfr_add(r25582, r25581, r25565, MPFR_RNDN); mpfr_mul(r25583, r25568, r25582, MPFR_RNDN); mpfr_add(r25584, r25580, r25583, MPFR_RNDN); mpfr_mul(r25585, r25572, r25584, MPFR_RNDN); mpfr_add(r25586, r25571, r25585, MPFR_RNDN); mpfr_sub(r25587, r25586, r25573, MPFR_RNDN); return mpfr_get_d(r25587, MPFR_RNDN); }