#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17473 = a; float r17474 = r17473 * r17473; float r17475 = b; float r17476 = r17475 * r17475; float r17477 = r17474 + r17476; float r17478 = r17477 * r17477; float r17479 = 4.0f; float r17480 = 1.0f; float r17481 = r17480 + r17473; float r17482 = r17474 * r17481; float r17483 = 3.0f; float r17484 = r17483 * r17473; float r17485 = r17480 - r17484; float r17486 = r17476 * r17485; float r17487 = r17482 + r17486; float r17488 = r17479 * r17487; float r17489 = r17478 + r17488; float r17490 = r17489 - r17480; return r17490; } double f_id(double a, double b) { double r17491 = a; double r17492 = r17491 * r17491; double r17493 = b; double r17494 = r17493 * r17493; double r17495 = r17492 + r17494; double r17496 = r17495 * r17495; double r17497 = 4.0; double r17498 = 1.0; double r17499 = r17498 + r17491; double r17500 = r17492 * r17499; double r17501 = 3.0; double r17502 = r17501 * r17491; double r17503 = r17498 - r17502; double r17504 = r17494 * r17503; double r17505 = r17500 + r17504; double r17506 = r17497 * r17505; double r17507 = r17496 + r17506; double r17508 = r17507 - r17498; return r17508; } double f_of(float a, float b) { float r17509 = 1.0f; float r17510 = a; float r17511 = 3.0f; float r17512 = r17510 * r17511; float r17513 = r17509 - r17512; float r17514 = b; float r17515 = r17514 * r17514; float r17516 = fma(r17510, r17510, r17510); float r17517 = r17516 * r17510; float r17518 = fma(r17513, r17515, r17517); float r17519 = 4.0f; float r17520 = 2.0f; float r17521 = r17514 * r17520; float r17522 = r17521 * r17514; float r17523 = r17510 * r17510; float r17524 = pow(r17510, r17519); float r17525 = pow(r17514, r17519); float r17526 = r17524 + r17525; float r17527 = fma(r17522, r17523, r17526); float r17528 = fma(r17518, r17519, r17527); float r17529 = r17528 - r17509; return r17529; } double f_od(double a, double b) { double r17530 = 1.0; double r17531 = a; double r17532 = 3.0; double r17533 = r17531 * r17532; double r17534 = r17530 - r17533; double r17535 = b; double r17536 = r17535 * r17535; double r17537 = fma(r17531, r17531, r17531); double r17538 = r17537 * r17531; double r17539 = fma(r17534, r17536, r17538); double r17540 = 4.0; double r17541 = 2.0; double r17542 = r17535 * r17541; double r17543 = r17542 * r17535; double r17544 = r17531 * r17531; double r17545 = pow(r17531, r17540); double r17546 = pow(r17535, r17540); double r17547 = r17545 + r17546; double r17548 = fma(r17543, r17544, r17547); double r17549 = fma(r17539, r17540, r17548); double r17550 = r17549 - r17530; return r17550; } 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 r17551, r17552, r17553, r17554, r17555, r17556, r17557, r17558, r17559, r17560, r17561, r17562, r17563, r17564, r17565, r17566, r17567, r17568; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17551); mpfr_init(r17552); mpfr_init(r17553); mpfr_init(r17554); mpfr_init(r17555); mpfr_init(r17556); mpfr_init_set_str(r17557, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17558, "1", 10, MPFR_RNDN); mpfr_init(r17559); mpfr_init(r17560); mpfr_init_set_str(r17561, "3", 10, MPFR_RNDN); mpfr_init(r17562); mpfr_init(r17563); mpfr_init(r17564); mpfr_init(r17565); mpfr_init(r17566); mpfr_init(r17567); mpfr_init(r17568); } double f_im(double a, double b) { mpfr_set_d(r17551, a, MPFR_RNDN); mpfr_sqr(r17552, r17551, MPFR_RNDN); mpfr_set_d(r17553, b, MPFR_RNDN); mpfr_sqr(r17554, r17553, MPFR_RNDN); mpfr_add(r17555, r17552, r17554, MPFR_RNDN); mpfr_sqr(r17556, r17555, MPFR_RNDN); ; ; mpfr_add(r17559, r17558, r17551, MPFR_RNDN); mpfr_mul(r17560, r17552, r17559, MPFR_RNDN); ; mpfr_mul(r17562, r17561, r17551, MPFR_RNDN); mpfr_sub(r17563, r17558, r17562, MPFR_RNDN); mpfr_mul(r17564, r17554, r17563, MPFR_RNDN); mpfr_add(r17565, r17560, r17564, MPFR_RNDN); mpfr_mul(r17566, r17557, r17565, MPFR_RNDN); mpfr_add(r17567, r17556, r17566, MPFR_RNDN); mpfr_sub(r17568, r17567, r17558, MPFR_RNDN); return mpfr_get_d(r17568, MPFR_RNDN); } static mpfr_t r17569, r17570, r17571, r17572, r17573, r17574, r17575, r17576, r17577, r17578, r17579, r17580, r17581, r17582, r17583, r17584, r17585, r17586, r17587, r17588, r17589; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17569, "1", 10, MPFR_RNDN); mpfr_init(r17570); mpfr_init_set_str(r17571, "3", 10, MPFR_RNDN); mpfr_init(r17572); mpfr_init(r17573); mpfr_init(r17574); mpfr_init(r17575); mpfr_init(r17576); mpfr_init(r17577); mpfr_init(r17578); mpfr_init_set_str(r17579, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17580, "2", 10, MPFR_RNDN); mpfr_init(r17581); mpfr_init(r17582); mpfr_init(r17583); mpfr_init(r17584); mpfr_init(r17585); mpfr_init(r17586); mpfr_init(r17587); mpfr_init(r17588); mpfr_init(r17589); } double f_fm(double a, double b) { ; mpfr_set_d(r17570, a, MPFR_RNDN); ; mpfr_mul(r17572, r17570, r17571, MPFR_RNDN); mpfr_sub(r17573, r17569, r17572, MPFR_RNDN); mpfr_set_d(r17574, b, MPFR_RNDN); mpfr_sqr(r17575, r17574, MPFR_RNDN); mpfr_fma(r17576, r17570, r17570, r17570, MPFR_RNDN); mpfr_mul(r17577, r17576, r17570, MPFR_RNDN); mpfr_fma(r17578, r17573, r17575, r17577, MPFR_RNDN); ; ; mpfr_mul(r17581, r17574, r17580, MPFR_RNDN); mpfr_mul(r17582, r17581, r17574, MPFR_RNDN); mpfr_mul(r17583, r17570, r17570, MPFR_RNDN); mpfr_pow(r17584, r17570, r17579, MPFR_RNDN); mpfr_pow(r17585, r17574, r17579, MPFR_RNDN); mpfr_add(r17586, r17584, r17585, MPFR_RNDN); mpfr_fma(r17587, r17582, r17583, r17586, MPFR_RNDN); mpfr_fma(r17588, r17578, r17579, r17587, MPFR_RNDN); mpfr_sub(r17589, r17588, r17569, MPFR_RNDN); return mpfr_get_d(r17589, MPFR_RNDN); } static mpfr_t r17590, r17591, r17592, r17593, r17594, r17595, r17596, r17597, r17598, r17599, r17600, r17601, r17602, r17603, r17604, r17605, r17606, r17607, r17608, r17609, r17610; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17590, "1", 10, MPFR_RNDN); mpfr_init(r17591); mpfr_init_set_str(r17592, "3", 10, MPFR_RNDN); mpfr_init(r17593); mpfr_init(r17594); mpfr_init(r17595); mpfr_init(r17596); mpfr_init(r17597); mpfr_init(r17598); mpfr_init(r17599); mpfr_init_set_str(r17600, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17601, "2", 10, MPFR_RNDN); mpfr_init(r17602); mpfr_init(r17603); mpfr_init(r17604); mpfr_init(r17605); mpfr_init(r17606); mpfr_init(r17607); mpfr_init(r17608); mpfr_init(r17609); mpfr_init(r17610); } double f_dm(double a, double b) { ; mpfr_set_d(r17591, a, MPFR_RNDN); ; mpfr_mul(r17593, r17591, r17592, MPFR_RNDN); mpfr_sub(r17594, r17590, r17593, MPFR_RNDN); mpfr_set_d(r17595, b, MPFR_RNDN); mpfr_sqr(r17596, r17595, MPFR_RNDN); mpfr_fma(r17597, r17591, r17591, r17591, MPFR_RNDN); mpfr_mul(r17598, r17597, r17591, MPFR_RNDN); mpfr_fma(r17599, r17594, r17596, r17598, MPFR_RNDN); ; ; mpfr_mul(r17602, r17595, r17601, MPFR_RNDN); mpfr_mul(r17603, r17602, r17595, MPFR_RNDN); mpfr_mul(r17604, r17591, r17591, MPFR_RNDN); mpfr_pow(r17605, r17591, r17600, MPFR_RNDN); mpfr_pow(r17606, r17595, r17600, MPFR_RNDN); mpfr_add(r17607, r17605, r17606, MPFR_RNDN); mpfr_fma(r17608, r17603, r17604, r17607, MPFR_RNDN); mpfr_fma(r17609, r17599, r17600, r17608, MPFR_RNDN); mpfr_sub(r17610, r17609, r17590, MPFR_RNDN); return mpfr_get_d(r17610, MPFR_RNDN); }