#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r26669 = a; float r26670 = r26669 * r26669; float r26671 = b; float r26672 = r26671 * r26671; float r26673 = r26670 + r26672; float r26674 = 2; float r26675 = pow(r26673, r26674); float r26676 = 4; float r26677 = 1; float r26678 = r26677 + r26669; float r26679 = r26670 * r26678; float r26680 = 3; float r26681 = r26680 * r26669; float r26682 = r26677 - r26681; float r26683 = r26672 * r26682; float r26684 = r26679 + r26683; float r26685 = r26676 * r26684; float r26686 = r26675 + r26685; float r26687 = r26686 - r26677; return r26687; } double f_id(double a, double b) { double r26688 = a; double r26689 = r26688 * r26688; double r26690 = b; double r26691 = r26690 * r26690; double r26692 = r26689 + r26691; double r26693 = 2; double r26694 = pow(r26692, r26693); double r26695 = 4; double r26696 = 1; double r26697 = r26696 + r26688; double r26698 = r26689 * r26697; double r26699 = 3; double r26700 = r26699 * r26688; double r26701 = r26696 - r26700; double r26702 = r26691 * r26701; double r26703 = r26698 + r26702; double r26704 = r26695 * r26703; double r26705 = r26694 + r26704; double r26706 = r26705 - r26696; return r26706; } double f_of(float a, float b) { float r26707 = b; float r26708 = a; float r26709 = r26708 * r26708; float r26710 = fma(r26707, r26707, r26709); float r26711 = fma(r26708, r26708, r26708); float r26712 = r26707 * r26707; float r26713 = fma(r26708, r26711, r26712); float r26714 = 4; float r26715 = 1; float r26716 = -r26715; float r26717 = fma(r26713, r26714, r26716); float r26718 = fma(r26710, r26710, r26717); return r26718; } double f_od(double a, double b) { double r26719 = b; double r26720 = a; double r26721 = r26720 * r26720; double r26722 = fma(r26719, r26719, r26721); double r26723 = fma(r26720, r26720, r26720); double r26724 = r26719 * r26719; double r26725 = fma(r26720, r26723, r26724); double r26726 = 4; double r26727 = 1; double r26728 = -r26727; double r26729 = fma(r26725, r26726, r26728); double r26730 = fma(r26722, r26722, r26729); return r26730; } 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 r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26731); mpfr_init(r26732); mpfr_init(r26733); mpfr_init(r26734); mpfr_init(r26735); mpfr_init_set_str(r26736, "2", 10, MPFR_RNDN); mpfr_init(r26737); mpfr_init_set_str(r26738, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26739, "1", 10, MPFR_RNDN); mpfr_init(r26740); mpfr_init(r26741); mpfr_init_set_str(r26742, "3", 10, MPFR_RNDN); mpfr_init(r26743); mpfr_init(r26744); mpfr_init(r26745); mpfr_init(r26746); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); } double f_im(double a, double b) { mpfr_set_d(r26731, a, MPFR_RNDN); mpfr_mul(r26732, r26731, r26731, MPFR_RNDN); mpfr_set_d(r26733, b, MPFR_RNDN); mpfr_mul(r26734, r26733, r26733, MPFR_RNDN); mpfr_add(r26735, r26732, r26734, MPFR_RNDN); ; mpfr_pow(r26737, r26735, r26736, MPFR_RNDN); ; ; mpfr_add(r26740, r26739, r26731, MPFR_RNDN); mpfr_mul(r26741, r26732, r26740, MPFR_RNDN); ; mpfr_mul(r26743, r26742, r26731, MPFR_RNDN); mpfr_sub(r26744, r26739, r26743, MPFR_RNDN); mpfr_mul(r26745, r26734, r26744, MPFR_RNDN); mpfr_add(r26746, r26741, r26745, MPFR_RNDN); mpfr_mul(r26747, r26738, r26746, MPFR_RNDN); mpfr_add(r26748, r26737, r26747, MPFR_RNDN); mpfr_sub(r26749, r26748, r26739, MPFR_RNDN); return mpfr_get_d(r26749, MPFR_RNDN); } static mpfr_t r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757, r26758, r26759, r26760, r26761; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); mpfr_init(r26754); mpfr_init(r26755); mpfr_init(r26756); mpfr_init_set_str(r26757, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26758, "1", 10, MPFR_RNDN); mpfr_init(r26759); mpfr_init(r26760); mpfr_init(r26761); } double f_fm(double a, double b) { mpfr_set_d(r26750, b, MPFR_RNDN); mpfr_set_d(r26751, a, MPFR_RNDN); mpfr_mul(r26752, r26751, r26751, MPFR_RNDN); mpfr_fma(r26753, r26750, r26750, r26752, MPFR_RNDN); mpfr_fma(r26754, r26751, r26751, r26751, MPFR_RNDN); mpfr_mul(r26755, r26750, r26750, MPFR_RNDN); mpfr_fma(r26756, r26751, r26754, r26755, MPFR_RNDN); ; ; mpfr_neg(r26759, r26758, MPFR_RNDN); mpfr_fma(r26760, r26756, r26757, r26759, MPFR_RNDN); mpfr_fma(r26761, r26753, r26753, r26760, MPFR_RNDN); return mpfr_get_d(r26761, MPFR_RNDN); } static mpfr_t r26762, r26763, r26764, r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r26762); mpfr_init(r26763); mpfr_init(r26764); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init_set_str(r26769, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26770, "1", 10, MPFR_RNDN); mpfr_init(r26771); mpfr_init(r26772); mpfr_init(r26773); } double f_dm(double a, double b) { mpfr_set_d(r26762, b, MPFR_RNDN); mpfr_set_d(r26763, a, MPFR_RNDN); mpfr_mul(r26764, r26763, r26763, MPFR_RNDN); mpfr_fma(r26765, r26762, r26762, r26764, MPFR_RNDN); mpfr_fma(r26766, r26763, r26763, r26763, MPFR_RNDN); mpfr_mul(r26767, r26762, r26762, MPFR_RNDN); mpfr_fma(r26768, r26763, r26766, r26767, MPFR_RNDN); ; ; mpfr_neg(r26771, r26770, MPFR_RNDN); mpfr_fma(r26772, r26768, r26769, r26771, MPFR_RNDN); mpfr_fma(r26773, r26765, r26765, r26772, MPFR_RNDN); return mpfr_get_d(r26773, MPFR_RNDN); }