#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r26711 = 4; float r26712 = 3; float r26713 = atan2(1.0, 0.0); float r26714 = r26712 * r26713; float r26715 = 1; float r26716 = v; float r26717 = r26716 * r26716; float r26718 = r26715 - r26717; float r26719 = r26714 * r26718; float r26720 = 2; float r26721 = 6; float r26722 = r26721 * r26717; float r26723 = r26720 - r26722; float r26724 = sqrt(r26723); float r26725 = r26719 * r26724; float r26726 = r26711 / r26725; return r26726; } double f_id(double v) { double r26727 = 4; double r26728 = 3; double r26729 = atan2(1.0, 0.0); double r26730 = r26728 * r26729; double r26731 = 1; double r26732 = v; double r26733 = r26732 * r26732; double r26734 = r26731 - r26733; double r26735 = r26730 * r26734; double r26736 = 2; double r26737 = 6; double r26738 = r26737 * r26733; double r26739 = r26736 - r26738; double r26740 = sqrt(r26739); double r26741 = r26735 * r26740; double r26742 = r26727 / r26741; return r26742; } double f_of(float v) { float r26743 = 4; float r26744 = atan2(1.0, 0.0); float r26745 = 3; float r26746 = r26744 * r26745; float r26747 = r26743 / r26746; float r26748 = 1; float r26749 = v; float r26750 = r26749 * r26749; float r26751 = r26748 - r26750; float r26752 = r26747 / r26751; float r26753 = 6; float r26754 = -r26753; float r26755 = 2; float r26756 = fma(r26754, r26750, r26755); float r26757 = sqrt(r26756); float r26758 = r26752 / r26757; float r26759 = pow(r26758, r26745); float r26760 = cbrt(r26759); return r26760; } double f_od(double v) { double r26761 = 4; double r26762 = atan2(1.0, 0.0); double r26763 = 3; double r26764 = r26762 * r26763; double r26765 = r26761 / r26764; double r26766 = 1; double r26767 = v; double r26768 = r26767 * r26767; double r26769 = r26766 - r26768; double r26770 = r26765 / r26769; double r26771 = 6; double r26772 = -r26771; double r26773 = 2; double r26774 = fma(r26772, r26768, r26773); double r26775 = sqrt(r26774); double r26776 = r26770 / r26775; double r26777 = pow(r26776, r26763); double r26778 = cbrt(r26777); return r26778; } 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 r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792, r26793, r26794; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r26779, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26780, "3", 10, MPFR_RNDN); mpfr_init(r26781); mpfr_init(r26782); mpfr_init_set_str(r26783, "1", 10, MPFR_RNDN); mpfr_init(r26784); mpfr_init(r26785); mpfr_init(r26786); mpfr_init(r26787); mpfr_init_set_str(r26788, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26789, "6", 10, MPFR_RNDN); mpfr_init(r26790); mpfr_init(r26791); mpfr_init(r26792); mpfr_init(r26793); mpfr_init(r26794); } double f_im(double v) { ; ; mpfr_const_pi(r26781, MPFR_RNDN); mpfr_mul(r26782, r26780, r26781, MPFR_RNDN); ; mpfr_set_d(r26784, v, MPFR_RNDN); mpfr_mul(r26785, r26784, r26784, MPFR_RNDN); mpfr_sub(r26786, r26783, r26785, MPFR_RNDN); mpfr_mul(r26787, r26782, r26786, MPFR_RNDN); ; ; mpfr_mul(r26790, r26789, r26785, MPFR_RNDN); mpfr_sub(r26791, r26788, r26790, MPFR_RNDN); mpfr_sqrt(r26792, r26791, MPFR_RNDN); mpfr_mul(r26793, r26787, r26792, MPFR_RNDN); mpfr_div(r26794, r26779, r26793, MPFR_RNDN); return mpfr_get_d(r26794, MPFR_RNDN); } static mpfr_t r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26795, "4", 10, MPFR_RNDN); mpfr_init(r26796); mpfr_init_set_str(r26797, "3", 10, MPFR_RNDN); mpfr_init(r26798); mpfr_init(r26799); mpfr_init_set_str(r26800, "1", 10, MPFR_RNDN); mpfr_init(r26801); mpfr_init(r26802); mpfr_init(r26803); mpfr_init(r26804); mpfr_init_set_str(r26805, "6", 10, MPFR_RNDN); mpfr_init(r26806); mpfr_init_set_str(r26807, "2", 10, MPFR_RNDN); mpfr_init(r26808); mpfr_init(r26809); mpfr_init(r26810); mpfr_init(r26811); mpfr_init(r26812); } double f_fm(double v) { ; mpfr_const_pi(r26796, MPFR_RNDN); ; mpfr_mul(r26798, r26796, r26797, MPFR_RNDN); mpfr_div(r26799, r26795, r26798, MPFR_RNDN); ; mpfr_set_d(r26801, v, MPFR_RNDN); mpfr_mul(r26802, r26801, r26801, MPFR_RNDN); mpfr_sub(r26803, r26800, r26802, MPFR_RNDN); mpfr_div(r26804, r26799, r26803, MPFR_RNDN); ; mpfr_neg(r26806, r26805, MPFR_RNDN); ; mpfr_fma(r26808, r26806, r26802, r26807, MPFR_RNDN); mpfr_sqrt(r26809, r26808, MPFR_RNDN); mpfr_div(r26810, r26804, r26809, MPFR_RNDN); mpfr_pow(r26811, r26810, r26797, MPFR_RNDN); mpfr_cbrt(r26812, r26811, MPFR_RNDN); return mpfr_get_d(r26812, MPFR_RNDN); } static mpfr_t r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827, r26828, r26829, r26830; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26813, "4", 10, MPFR_RNDN); mpfr_init(r26814); mpfr_init_set_str(r26815, "3", 10, MPFR_RNDN); mpfr_init(r26816); mpfr_init(r26817); mpfr_init_set_str(r26818, "1", 10, MPFR_RNDN); mpfr_init(r26819); mpfr_init(r26820); mpfr_init(r26821); mpfr_init(r26822); mpfr_init_set_str(r26823, "6", 10, MPFR_RNDN); mpfr_init(r26824); mpfr_init_set_str(r26825, "2", 10, MPFR_RNDN); mpfr_init(r26826); mpfr_init(r26827); mpfr_init(r26828); mpfr_init(r26829); mpfr_init(r26830); } double f_dm(double v) { ; mpfr_const_pi(r26814, MPFR_RNDN); ; mpfr_mul(r26816, r26814, r26815, MPFR_RNDN); mpfr_div(r26817, r26813, r26816, MPFR_RNDN); ; mpfr_set_d(r26819, v, MPFR_RNDN); mpfr_mul(r26820, r26819, r26819, MPFR_RNDN); mpfr_sub(r26821, r26818, r26820, MPFR_RNDN); mpfr_div(r26822, r26817, r26821, MPFR_RNDN); ; mpfr_neg(r26824, r26823, MPFR_RNDN); ; mpfr_fma(r26826, r26824, r26820, r26825, MPFR_RNDN); mpfr_sqrt(r26827, r26826, MPFR_RNDN); mpfr_div(r26828, r26822, r26827, MPFR_RNDN); mpfr_pow(r26829, r26828, r26815, MPFR_RNDN); mpfr_cbrt(r26830, r26829, MPFR_RNDN); return mpfr_get_d(r26830, MPFR_RNDN); }