#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r10604 = a; float r10605 = r10604 * r10604; float r10606 = b; float r10607 = r10606 * r10606; float r10608 = r10605 + r10607; float r10609 = 2.0f; float r10610 = pow(r10608, r10609); float r10611 = 4.0f; float r10612 = 1.0f; float r10613 = r10612 + r10604; float r10614 = r10605 * r10613; float r10615 = 3.0f; float r10616 = r10615 * r10604; float r10617 = r10612 - r10616; float r10618 = r10607 * r10617; float r10619 = r10614 + r10618; float r10620 = r10611 * r10619; float r10621 = r10610 + r10620; float r10622 = r10621 - r10612; return r10622; } double f_id(double a, double b) { double r10623 = a; double r10624 = r10623 * r10623; double r10625 = b; double r10626 = r10625 * r10625; double r10627 = r10624 + r10626; double r10628 = 2.0; double r10629 = pow(r10627, r10628); double r10630 = 4.0; double r10631 = 1.0; double r10632 = r10631 + r10623; double r10633 = r10624 * r10632; double r10634 = 3.0; double r10635 = r10634 * r10623; double r10636 = r10631 - r10635; double r10637 = r10626 * r10636; double r10638 = r10633 + r10637; double r10639 = r10630 * r10638; double r10640 = r10629 + r10639; double r10641 = r10640 - r10631; return r10641; } double f_of(float a, float b) { float r10642 = a; float r10643 = r10642 * r10642; float r10644 = b; float r10645 = r10644 * r10644; float r10646 = r10643 + r10645; float r10647 = 2.0f; float r10648 = pow(r10646, r10647); float r10649 = 1.0f; float r10650 = r10642 + r10649; float r10651 = r10643 * r10650; float r10652 = 3.0f; float r10653 = r10652 * r10642; float r10654 = r10649 - r10653; float r10655 = r10645 * r10654; float r10656 = r10651 + r10655; float r10657 = 4.0f; float r10658 = r10656 * r10657; float r10659 = r10648 + r10658; float r10660 = sqrt(r10659); float r10661 = sqrt(r10660); float r10662 = r10661 * r10661; float r10663 = r10662 * r10660; float r10664 = r10663 - r10649; return r10664; } double f_od(double a, double b) { double r10665 = a; double r10666 = r10665 * r10665; double r10667 = b; double r10668 = r10667 * r10667; double r10669 = r10666 + r10668; double r10670 = 2.0; double r10671 = pow(r10669, r10670); double r10672 = 1.0; double r10673 = r10665 + r10672; double r10674 = r10666 * r10673; double r10675 = 3.0; double r10676 = r10675 * r10665; double r10677 = r10672 - r10676; double r10678 = r10668 * r10677; double r10679 = r10674 + r10678; double r10680 = 4.0; double r10681 = r10679 * r10680; double r10682 = r10671 + r10681; double r10683 = sqrt(r10682); double r10684 = sqrt(r10683); double r10685 = r10684 * r10684; double r10686 = r10685 * r10683; double r10687 = r10686 - r10672; return r10687; } 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 r10688, r10689, r10690, r10691, r10692, r10693, r10694, r10695, r10696, r10697, r10698, r10699, r10700, r10701, r10702, r10703, r10704, r10705, r10706; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10688); mpfr_init(r10689); mpfr_init(r10690); mpfr_init(r10691); mpfr_init(r10692); mpfr_init_set_str(r10693, "2", 10, MPFR_RNDN); mpfr_init(r10694); mpfr_init_set_str(r10695, "4", 10, MPFR_RNDN); mpfr_init_set_str(r10696, "1", 10, MPFR_RNDN); mpfr_init(r10697); mpfr_init(r10698); mpfr_init_set_str(r10699, "3", 10, MPFR_RNDN); mpfr_init(r10700); mpfr_init(r10701); mpfr_init(r10702); mpfr_init(r10703); mpfr_init(r10704); mpfr_init(r10705); mpfr_init(r10706); } double f_im(double a, double b) { mpfr_set_d(r10688, a, MPFR_RNDN); mpfr_mul(r10689, r10688, r10688, MPFR_RNDN); mpfr_set_d(r10690, b, MPFR_RNDN); mpfr_mul(r10691, r10690, r10690, MPFR_RNDN); mpfr_add(r10692, r10689, r10691, MPFR_RNDN); ; mpfr_pow(r10694, r10692, r10693, MPFR_RNDN); ; ; mpfr_add(r10697, r10696, r10688, MPFR_RNDN); mpfr_mul(r10698, r10689, r10697, MPFR_RNDN); ; mpfr_mul(r10700, r10699, r10688, MPFR_RNDN); mpfr_sub(r10701, r10696, r10700, MPFR_RNDN); mpfr_mul(r10702, r10691, r10701, MPFR_RNDN); mpfr_add(r10703, r10698, r10702, MPFR_RNDN); mpfr_mul(r10704, r10695, r10703, MPFR_RNDN); mpfr_add(r10705, r10694, r10704, MPFR_RNDN); mpfr_sub(r10706, r10705, r10696, MPFR_RNDN); return mpfr_get_d(r10706, MPFR_RNDN); } static mpfr_t r10707, r10708, r10709, r10710, r10711, r10712, r10713, r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721, r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10707); mpfr_init(r10708); mpfr_init(r10709); mpfr_init(r10710); mpfr_init(r10711); mpfr_init_set_str(r10712, "2", 10, MPFR_RNDN); mpfr_init(r10713); mpfr_init_set_str(r10714, "1", 10, MPFR_RNDN); mpfr_init(r10715); mpfr_init(r10716); mpfr_init_set_str(r10717, "3", 10, MPFR_RNDN); mpfr_init(r10718); mpfr_init(r10719); mpfr_init(r10720); mpfr_init(r10721); mpfr_init_set_str(r10722, "4", 10, MPFR_RNDN); mpfr_init(r10723); mpfr_init(r10724); mpfr_init(r10725); mpfr_init(r10726); mpfr_init(r10727); mpfr_init(r10728); mpfr_init(r10729); } double f_fm(double a, double b) { mpfr_set_d(r10707, a, MPFR_RNDN); mpfr_mul(r10708, r10707, r10707, MPFR_RNDN); mpfr_set_d(r10709, b, MPFR_RNDN); mpfr_mul(r10710, r10709, r10709, MPFR_RNDN); mpfr_add(r10711, r10708, r10710, MPFR_RNDN); ; mpfr_pow(r10713, r10711, r10712, MPFR_RNDN); ; mpfr_add(r10715, r10707, r10714, MPFR_RNDN); mpfr_mul(r10716, r10708, r10715, MPFR_RNDN); ; mpfr_mul(r10718, r10717, r10707, MPFR_RNDN); mpfr_sub(r10719, r10714, r10718, MPFR_RNDN); mpfr_mul(r10720, r10710, r10719, MPFR_RNDN); mpfr_add(r10721, r10716, r10720, MPFR_RNDN); ; mpfr_mul(r10723, r10721, r10722, MPFR_RNDN); mpfr_add(r10724, r10713, r10723, MPFR_RNDN); mpfr_sqrt(r10725, r10724, MPFR_RNDN); mpfr_sqrt(r10726, r10725, MPFR_RNDN); mpfr_mul(r10727, r10726, r10726, MPFR_RNDN); mpfr_mul(r10728, r10727, r10725, MPFR_RNDN); mpfr_sub(r10729, r10728, r10714, MPFR_RNDN); return mpfr_get_d(r10729, MPFR_RNDN); } static mpfr_t r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739, r10740, r10741, r10742, r10743, r10744, r10745, r10746, r10747, r10748, r10749, r10750, r10751, r10752; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10730); mpfr_init(r10731); mpfr_init(r10732); mpfr_init(r10733); mpfr_init(r10734); mpfr_init_set_str(r10735, "2", 10, MPFR_RNDN); mpfr_init(r10736); mpfr_init_set_str(r10737, "1", 10, MPFR_RNDN); mpfr_init(r10738); mpfr_init(r10739); mpfr_init_set_str(r10740, "3", 10, MPFR_RNDN); mpfr_init(r10741); mpfr_init(r10742); mpfr_init(r10743); mpfr_init(r10744); mpfr_init_set_str(r10745, "4", 10, MPFR_RNDN); mpfr_init(r10746); mpfr_init(r10747); mpfr_init(r10748); mpfr_init(r10749); mpfr_init(r10750); mpfr_init(r10751); mpfr_init(r10752); } double f_dm(double a, double b) { mpfr_set_d(r10730, a, MPFR_RNDN); mpfr_mul(r10731, r10730, r10730, MPFR_RNDN); mpfr_set_d(r10732, b, MPFR_RNDN); mpfr_mul(r10733, r10732, r10732, MPFR_RNDN); mpfr_add(r10734, r10731, r10733, MPFR_RNDN); ; mpfr_pow(r10736, r10734, r10735, MPFR_RNDN); ; mpfr_add(r10738, r10730, r10737, MPFR_RNDN); mpfr_mul(r10739, r10731, r10738, MPFR_RNDN); ; mpfr_mul(r10741, r10740, r10730, MPFR_RNDN); mpfr_sub(r10742, r10737, r10741, MPFR_RNDN); mpfr_mul(r10743, r10733, r10742, MPFR_RNDN); mpfr_add(r10744, r10739, r10743, MPFR_RNDN); ; mpfr_mul(r10746, r10744, r10745, MPFR_RNDN); mpfr_add(r10747, r10736, r10746, MPFR_RNDN); mpfr_sqrt(r10748, r10747, MPFR_RNDN); mpfr_sqrt(r10749, r10748, MPFR_RNDN); mpfr_mul(r10750, r10749, r10749, MPFR_RNDN); mpfr_mul(r10751, r10750, r10748, MPFR_RNDN); mpfr_sub(r10752, r10751, r10737, MPFR_RNDN); return mpfr_get_d(r10752, MPFR_RNDN); }