#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17691 = a; float r17692 = r17691 * r17691; float r17693 = b; float r17694 = r17693 * r17693; float r17695 = r17692 + r17694; float r17696 = r17695 * r17695; float r17697 = 4.0f; float r17698 = 1.0f; float r17699 = r17698 + r17691; float r17700 = r17692 * r17699; float r17701 = 3.0f; float r17702 = r17701 * r17691; float r17703 = r17698 - r17702; float r17704 = r17694 * r17703; float r17705 = r17700 + r17704; float r17706 = r17697 * r17705; float r17707 = r17696 + r17706; float r17708 = r17707 - r17698; return r17708; } double f_id(double a, double b) { double r17709 = a; double r17710 = r17709 * r17709; double r17711 = b; double r17712 = r17711 * r17711; double r17713 = r17710 + r17712; double r17714 = r17713 * r17713; double r17715 = 4.0; double r17716 = 1.0; double r17717 = r17716 + r17709; double r17718 = r17710 * r17717; double r17719 = 3.0; double r17720 = r17719 * r17709; double r17721 = r17716 - r17720; double r17722 = r17712 * r17721; double r17723 = r17718 + r17722; double r17724 = r17715 * r17723; double r17725 = r17714 + r17724; double r17726 = r17725 - r17716; return r17726; } double f_of(float a, float b) { float r17727 = 1.0f; float r17728 = a; float r17729 = 3.0f; float r17730 = r17728 * r17729; float r17731 = r17727 - r17730; float r17732 = b; float r17733 = r17732 * r17732; float r17734 = fma(r17728, r17728, r17728); float r17735 = r17734 * r17728; float r17736 = fma(r17731, r17733, r17735); float r17737 = 4.0f; float r17738 = 2.0f; float r17739 = r17732 * r17738; float r17740 = r17739 * r17732; float r17741 = r17728 * r17728; float r17742 = pow(r17728, r17737); float r17743 = pow(r17732, r17737); float r17744 = r17742 + r17743; float r17745 = fma(r17740, r17741, r17744); float r17746 = fma(r17736, r17737, r17745); float r17747 = r17746 - r17727; return r17747; } double f_od(double a, double b) { double r17748 = 1.0; double r17749 = a; double r17750 = 3.0; double r17751 = r17749 * r17750; double r17752 = r17748 - r17751; double r17753 = b; double r17754 = r17753 * r17753; double r17755 = fma(r17749, r17749, r17749); double r17756 = r17755 * r17749; double r17757 = fma(r17752, r17754, r17756); double r17758 = 4.0; double r17759 = 2.0; double r17760 = r17753 * r17759; double r17761 = r17760 * r17753; double r17762 = r17749 * r17749; double r17763 = pow(r17749, r17758); double r17764 = pow(r17753, r17758); double r17765 = r17763 + r17764; double r17766 = fma(r17761, r17762, r17765); double r17767 = fma(r17757, r17758, r17766); double r17768 = r17767 - r17748; return r17768; } 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 r17769, r17770, r17771, r17772, r17773, r17774, r17775, r17776, r17777, r17778, r17779, r17780, r17781, r17782, r17783, r17784, r17785, r17786; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17769); mpfr_init(r17770); mpfr_init(r17771); mpfr_init(r17772); mpfr_init(r17773); mpfr_init(r17774); mpfr_init_set_str(r17775, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17776, "1", 10, MPFR_RNDN); mpfr_init(r17777); mpfr_init(r17778); mpfr_init_set_str(r17779, "3", 10, MPFR_RNDN); mpfr_init(r17780); mpfr_init(r17781); mpfr_init(r17782); mpfr_init(r17783); mpfr_init(r17784); mpfr_init(r17785); mpfr_init(r17786); } double f_im(double a, double b) { mpfr_set_d(r17769, a, MPFR_RNDN); mpfr_sqr(r17770, r17769, MPFR_RNDN); mpfr_set_d(r17771, b, MPFR_RNDN); mpfr_sqr(r17772, r17771, MPFR_RNDN); mpfr_add(r17773, r17770, r17772, MPFR_RNDN); mpfr_sqr(r17774, r17773, MPFR_RNDN); ; ; mpfr_add(r17777, r17776, r17769, MPFR_RNDN); mpfr_mul(r17778, r17770, r17777, MPFR_RNDN); ; mpfr_mul(r17780, r17779, r17769, MPFR_RNDN); mpfr_sub(r17781, r17776, r17780, MPFR_RNDN); mpfr_mul(r17782, r17772, r17781, MPFR_RNDN); mpfr_add(r17783, r17778, r17782, MPFR_RNDN); mpfr_mul(r17784, r17775, r17783, MPFR_RNDN); mpfr_add(r17785, r17774, r17784, MPFR_RNDN); mpfr_sub(r17786, r17785, r17776, MPFR_RNDN); return mpfr_get_d(r17786, MPFR_RNDN); } static mpfr_t r17787, r17788, r17789, r17790, r17791, r17792, r17793, r17794, r17795, r17796, r17797, r17798, r17799, r17800, r17801, r17802, r17803, r17804, r17805, r17806, r17807; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17787, "1", 10, MPFR_RNDN); mpfr_init(r17788); mpfr_init_set_str(r17789, "3", 10, MPFR_RNDN); mpfr_init(r17790); mpfr_init(r17791); mpfr_init(r17792); mpfr_init(r17793); mpfr_init(r17794); mpfr_init(r17795); mpfr_init(r17796); mpfr_init_set_str(r17797, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17798, "2", 10, MPFR_RNDN); mpfr_init(r17799); mpfr_init(r17800); mpfr_init(r17801); mpfr_init(r17802); mpfr_init(r17803); mpfr_init(r17804); mpfr_init(r17805); mpfr_init(r17806); mpfr_init(r17807); } double f_fm(double a, double b) { ; mpfr_set_d(r17788, a, MPFR_RNDN); ; mpfr_mul(r17790, r17788, r17789, MPFR_RNDN); mpfr_sub(r17791, r17787, r17790, MPFR_RNDN); mpfr_set_d(r17792, b, MPFR_RNDN); mpfr_sqr(r17793, r17792, MPFR_RNDN); mpfr_fma(r17794, r17788, r17788, r17788, MPFR_RNDN); mpfr_mul(r17795, r17794, r17788, MPFR_RNDN); mpfr_fma(r17796, r17791, r17793, r17795, MPFR_RNDN); ; ; mpfr_mul(r17799, r17792, r17798, MPFR_RNDN); mpfr_mul(r17800, r17799, r17792, MPFR_RNDN); mpfr_mul(r17801, r17788, r17788, MPFR_RNDN); mpfr_pow(r17802, r17788, r17797, MPFR_RNDN); mpfr_pow(r17803, r17792, r17797, MPFR_RNDN); mpfr_add(r17804, r17802, r17803, MPFR_RNDN); mpfr_fma(r17805, r17800, r17801, r17804, MPFR_RNDN); mpfr_fma(r17806, r17796, r17797, r17805, MPFR_RNDN); mpfr_sub(r17807, r17806, r17787, MPFR_RNDN); return mpfr_get_d(r17807, MPFR_RNDN); } static mpfr_t r17808, r17809, r17810, r17811, r17812, r17813, r17814, r17815, r17816, r17817, r17818, r17819, r17820, r17821, r17822, r17823, r17824, r17825, r17826, r17827, r17828; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17808, "1", 10, MPFR_RNDN); mpfr_init(r17809); mpfr_init_set_str(r17810, "3", 10, MPFR_RNDN); mpfr_init(r17811); mpfr_init(r17812); mpfr_init(r17813); mpfr_init(r17814); mpfr_init(r17815); mpfr_init(r17816); mpfr_init(r17817); mpfr_init_set_str(r17818, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17819, "2", 10, MPFR_RNDN); mpfr_init(r17820); mpfr_init(r17821); mpfr_init(r17822); mpfr_init(r17823); mpfr_init(r17824); mpfr_init(r17825); mpfr_init(r17826); mpfr_init(r17827); mpfr_init(r17828); } double f_dm(double a, double b) { ; mpfr_set_d(r17809, a, MPFR_RNDN); ; mpfr_mul(r17811, r17809, r17810, MPFR_RNDN); mpfr_sub(r17812, r17808, r17811, MPFR_RNDN); mpfr_set_d(r17813, b, MPFR_RNDN); mpfr_sqr(r17814, r17813, MPFR_RNDN); mpfr_fma(r17815, r17809, r17809, r17809, MPFR_RNDN); mpfr_mul(r17816, r17815, r17809, MPFR_RNDN); mpfr_fma(r17817, r17812, r17814, r17816, MPFR_RNDN); ; ; mpfr_mul(r17820, r17813, r17819, MPFR_RNDN); mpfr_mul(r17821, r17820, r17813, MPFR_RNDN); mpfr_mul(r17822, r17809, r17809, MPFR_RNDN); mpfr_pow(r17823, r17809, r17818, MPFR_RNDN); mpfr_pow(r17824, r17813, r17818, MPFR_RNDN); mpfr_add(r17825, r17823, r17824, MPFR_RNDN); mpfr_fma(r17826, r17821, r17822, r17825, MPFR_RNDN); mpfr_fma(r17827, r17817, r17818, r17826, MPFR_RNDN); mpfr_sub(r17828, r17827, r17808, MPFR_RNDN); return mpfr_get_d(r17828, MPFR_RNDN); }