#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17508 = a; float r17509 = r17508 * r17508; float r17510 = b; float r17511 = r17510 * r17510; float r17512 = r17509 + r17511; float r17513 = r17512 * r17512; float r17514 = 4.0f; float r17515 = 1.0f; float r17516 = r17515 + r17508; float r17517 = r17509 * r17516; float r17518 = 3.0f; float r17519 = r17518 * r17508; float r17520 = r17515 - r17519; float r17521 = r17511 * r17520; float r17522 = r17517 + r17521; float r17523 = r17514 * r17522; float r17524 = r17513 + r17523; float r17525 = r17524 - r17515; return r17525; } double f_id(double a, double b) { double r17526 = a; double r17527 = r17526 * r17526; double r17528 = b; double r17529 = r17528 * r17528; double r17530 = r17527 + r17529; double r17531 = r17530 * r17530; double r17532 = 4.0; double r17533 = 1.0; double r17534 = r17533 + r17526; double r17535 = r17527 * r17534; double r17536 = 3.0; double r17537 = r17536 * r17526; double r17538 = r17533 - r17537; double r17539 = r17529 * r17538; double r17540 = r17535 + r17539; double r17541 = r17532 * r17540; double r17542 = r17531 + r17541; double r17543 = r17542 - r17533; return r17543; } double f_of(float a, float b) { float r17544 = 1.0f; float r17545 = a; float r17546 = 3.0f; float r17547 = r17545 * r17546; float r17548 = r17544 - r17547; float r17549 = b; float r17550 = r17549 * r17549; float r17551 = fma(r17545, r17545, r17545); float r17552 = r17551 * r17545; float r17553 = fma(r17548, r17550, r17552); float r17554 = 4.0f; float r17555 = 2.0f; float r17556 = r17549 * r17555; float r17557 = r17556 * r17549; float r17558 = r17545 * r17545; float r17559 = pow(r17545, r17554); float r17560 = pow(r17549, r17554); float r17561 = r17559 + r17560; float r17562 = fma(r17557, r17558, r17561); float r17563 = fma(r17553, r17554, r17562); float r17564 = r17563 - r17544; return r17564; } double f_od(double a, double b) { double r17565 = 1.0; double r17566 = a; double r17567 = 3.0; double r17568 = r17566 * r17567; double r17569 = r17565 - r17568; double r17570 = b; double r17571 = r17570 * r17570; double r17572 = fma(r17566, r17566, r17566); double r17573 = r17572 * r17566; double r17574 = fma(r17569, r17571, r17573); double r17575 = 4.0; double r17576 = 2.0; double r17577 = r17570 * r17576; double r17578 = r17577 * r17570; double r17579 = r17566 * r17566; double r17580 = pow(r17566, r17575); double r17581 = pow(r17570, r17575); double r17582 = r17580 + r17581; double r17583 = fma(r17578, r17579, r17582); double r17584 = fma(r17574, r17575, r17583); double r17585 = r17584 - r17565; return r17585; } 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 r17586, r17587, r17588, r17589, r17590, r17591, r17592, r17593, r17594, r17595, r17596, r17597, r17598, r17599, r17600, r17601, r17602, r17603; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17586); mpfr_init(r17587); mpfr_init(r17588); mpfr_init(r17589); mpfr_init(r17590); mpfr_init(r17591); mpfr_init_set_str(r17592, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17593, "1", 10, MPFR_RNDN); mpfr_init(r17594); mpfr_init(r17595); mpfr_init_set_str(r17596, "3", 10, MPFR_RNDN); mpfr_init(r17597); mpfr_init(r17598); mpfr_init(r17599); mpfr_init(r17600); mpfr_init(r17601); mpfr_init(r17602); mpfr_init(r17603); } double f_im(double a, double b) { mpfr_set_d(r17586, a, MPFR_RNDN); mpfr_sqr(r17587, r17586, MPFR_RNDN); mpfr_set_d(r17588, b, MPFR_RNDN); mpfr_sqr(r17589, r17588, MPFR_RNDN); mpfr_add(r17590, r17587, r17589, MPFR_RNDN); mpfr_sqr(r17591, r17590, MPFR_RNDN); ; ; mpfr_add(r17594, r17593, r17586, MPFR_RNDN); mpfr_mul(r17595, r17587, r17594, MPFR_RNDN); ; mpfr_mul(r17597, r17596, r17586, MPFR_RNDN); mpfr_sub(r17598, r17593, r17597, MPFR_RNDN); mpfr_mul(r17599, r17589, r17598, MPFR_RNDN); mpfr_add(r17600, r17595, r17599, MPFR_RNDN); mpfr_mul(r17601, r17592, r17600, MPFR_RNDN); mpfr_add(r17602, r17591, r17601, MPFR_RNDN); mpfr_sub(r17603, r17602, r17593, MPFR_RNDN); return mpfr_get_d(r17603, MPFR_RNDN); } static mpfr_t r17604, r17605, r17606, r17607, r17608, r17609, r17610, r17611, r17612, r17613, r17614, r17615, r17616, r17617, r17618, r17619, r17620, r17621, r17622, r17623, r17624; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17604, "1", 10, MPFR_RNDN); mpfr_init(r17605); mpfr_init_set_str(r17606, "3", 10, MPFR_RNDN); mpfr_init(r17607); mpfr_init(r17608); mpfr_init(r17609); mpfr_init(r17610); mpfr_init(r17611); mpfr_init(r17612); mpfr_init(r17613); mpfr_init_set_str(r17614, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17615, "2", 10, MPFR_RNDN); mpfr_init(r17616); mpfr_init(r17617); mpfr_init(r17618); mpfr_init(r17619); mpfr_init(r17620); mpfr_init(r17621); mpfr_init(r17622); mpfr_init(r17623); mpfr_init(r17624); } double f_fm(double a, double b) { ; mpfr_set_d(r17605, a, MPFR_RNDN); ; mpfr_mul(r17607, r17605, r17606, MPFR_RNDN); mpfr_sub(r17608, r17604, r17607, MPFR_RNDN); mpfr_set_d(r17609, b, MPFR_RNDN); mpfr_sqr(r17610, r17609, MPFR_RNDN); mpfr_fma(r17611, r17605, r17605, r17605, MPFR_RNDN); mpfr_mul(r17612, r17611, r17605, MPFR_RNDN); mpfr_fma(r17613, r17608, r17610, r17612, MPFR_RNDN); ; ; mpfr_mul(r17616, r17609, r17615, MPFR_RNDN); mpfr_mul(r17617, r17616, r17609, MPFR_RNDN); mpfr_mul(r17618, r17605, r17605, MPFR_RNDN); mpfr_pow(r17619, r17605, r17614, MPFR_RNDN); mpfr_pow(r17620, r17609, r17614, MPFR_RNDN); mpfr_add(r17621, r17619, r17620, MPFR_RNDN); mpfr_fma(r17622, r17617, r17618, r17621, MPFR_RNDN); mpfr_fma(r17623, r17613, r17614, r17622, MPFR_RNDN); mpfr_sub(r17624, r17623, r17604, MPFR_RNDN); return mpfr_get_d(r17624, MPFR_RNDN); } static mpfr_t r17625, r17626, r17627, r17628, r17629, r17630, r17631, r17632, r17633, r17634, r17635, r17636, r17637, r17638, r17639, r17640, r17641, r17642, r17643, r17644, r17645; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17625, "1", 10, MPFR_RNDN); mpfr_init(r17626); mpfr_init_set_str(r17627, "3", 10, MPFR_RNDN); mpfr_init(r17628); mpfr_init(r17629); mpfr_init(r17630); mpfr_init(r17631); mpfr_init(r17632); mpfr_init(r17633); mpfr_init(r17634); mpfr_init_set_str(r17635, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17636, "2", 10, MPFR_RNDN); mpfr_init(r17637); mpfr_init(r17638); mpfr_init(r17639); mpfr_init(r17640); mpfr_init(r17641); mpfr_init(r17642); mpfr_init(r17643); mpfr_init(r17644); mpfr_init(r17645); } double f_dm(double a, double b) { ; mpfr_set_d(r17626, a, MPFR_RNDN); ; mpfr_mul(r17628, r17626, r17627, MPFR_RNDN); mpfr_sub(r17629, r17625, r17628, MPFR_RNDN); mpfr_set_d(r17630, b, MPFR_RNDN); mpfr_sqr(r17631, r17630, MPFR_RNDN); mpfr_fma(r17632, r17626, r17626, r17626, MPFR_RNDN); mpfr_mul(r17633, r17632, r17626, MPFR_RNDN); mpfr_fma(r17634, r17629, r17631, r17633, MPFR_RNDN); ; ; mpfr_mul(r17637, r17630, r17636, MPFR_RNDN); mpfr_mul(r17638, r17637, r17630, MPFR_RNDN); mpfr_mul(r17639, r17626, r17626, MPFR_RNDN); mpfr_pow(r17640, r17626, r17635, MPFR_RNDN); mpfr_pow(r17641, r17630, r17635, MPFR_RNDN); mpfr_add(r17642, r17640, r17641, MPFR_RNDN); mpfr_fma(r17643, r17638, r17639, r17642, MPFR_RNDN); mpfr_fma(r17644, r17634, r17635, r17643, MPFR_RNDN); mpfr_sub(r17645, r17644, r17625, MPFR_RNDN); return mpfr_get_d(r17645, MPFR_RNDN); }