#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r26745 = 4; float r26746 = 3; float r26747 = atan2(1.0, 0.0); float r26748 = r26746 * r26747; float r26749 = 1; float r26750 = v; float r26751 = r26750 * r26750; float r26752 = r26749 - r26751; float r26753 = r26748 * r26752; float r26754 = 2; float r26755 = 6; float r26756 = r26755 * r26751; float r26757 = r26754 - r26756; float r26758 = sqrt(r26757); float r26759 = r26753 * r26758; float r26760 = r26745 / r26759; return r26760; } double f_id(double v) { double r26761 = 4; double r26762 = 3; double r26763 = atan2(1.0, 0.0); double r26764 = r26762 * r26763; double r26765 = 1; double r26766 = v; double r26767 = r26766 * r26766; double r26768 = r26765 - r26767; double r26769 = r26764 * r26768; double r26770 = 2; double r26771 = 6; double r26772 = r26771 * r26767; double r26773 = r26770 - r26772; double r26774 = sqrt(r26773); double r26775 = r26769 * r26774; double r26776 = r26761 / r26775; return r26776; } double f_of(float v) { float r26777 = 4; float r26778 = 3; float r26779 = atan2(1.0, 0.0); float r26780 = r26778 * r26779; float r26781 = 1; float r26782 = v; float r26783 = r26782 * r26782; float r26784 = r26781 - r26783; float r26785 = r26780 * r26784; float r26786 = r26777 / r26785; float r26787 = 2; float r26788 = 6; float r26789 = r26788 * r26783; float r26790 = r26787 - r26789; float r26791 = sqrt(r26790); float r26792 = r26786 / r26791; return r26792; } double f_od(double v) { double r26793 = 4; double r26794 = 3; double r26795 = atan2(1.0, 0.0); double r26796 = r26794 * r26795; double r26797 = 1; double r26798 = v; double r26799 = r26798 * r26798; double r26800 = r26797 - r26799; double r26801 = r26796 * r26800; double r26802 = r26793 / r26801; double r26803 = 2; double r26804 = 6; double r26805 = r26804 * r26799; double r26806 = r26803 - r26805; double r26807 = sqrt(r26806); double r26808 = r26802 / r26807; return r26808; } 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 r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r26809, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26810, "3", 10, MPFR_RNDN); mpfr_init(r26811); mpfr_init(r26812); mpfr_init_set_str(r26813, "1", 10, MPFR_RNDN); mpfr_init(r26814); mpfr_init(r26815); mpfr_init(r26816); mpfr_init(r26817); mpfr_init_set_str(r26818, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26819, "6", 10, MPFR_RNDN); mpfr_init(r26820); mpfr_init(r26821); mpfr_init(r26822); mpfr_init(r26823); mpfr_init(r26824); } double f_im(double v) { ; ; mpfr_const_pi(r26811, MPFR_RNDN); mpfr_mul(r26812, r26810, r26811, MPFR_RNDN); ; mpfr_set_d(r26814, v, MPFR_RNDN); mpfr_mul(r26815, r26814, r26814, MPFR_RNDN); mpfr_sub(r26816, r26813, r26815, MPFR_RNDN); mpfr_mul(r26817, r26812, r26816, MPFR_RNDN); ; ; mpfr_mul(r26820, r26819, r26815, MPFR_RNDN); mpfr_sub(r26821, r26818, r26820, MPFR_RNDN); mpfr_sqrt(r26822, r26821, MPFR_RNDN); mpfr_mul(r26823, r26817, r26822, MPFR_RNDN); mpfr_div(r26824, r26809, r26823, MPFR_RNDN); return mpfr_get_d(r26824, MPFR_RNDN); } static mpfr_t r26825, r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834, r26835, r26836, r26837, r26838, r26839, r26840; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26825, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26826, "3", 10, MPFR_RNDN); mpfr_init(r26827); mpfr_init(r26828); mpfr_init_set_str(r26829, "1", 10, MPFR_RNDN); mpfr_init(r26830); mpfr_init(r26831); mpfr_init(r26832); mpfr_init(r26833); mpfr_init(r26834); mpfr_init_set_str(r26835, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26836, "6", 10, MPFR_RNDN); mpfr_init(r26837); mpfr_init(r26838); mpfr_init(r26839); mpfr_init(r26840); } double f_fm(double v) { ; ; mpfr_const_pi(r26827, MPFR_RNDN); mpfr_mul(r26828, r26826, r26827, MPFR_RNDN); ; mpfr_set_d(r26830, v, MPFR_RNDN); mpfr_mul(r26831, r26830, r26830, MPFR_RNDN); mpfr_sub(r26832, r26829, r26831, MPFR_RNDN); mpfr_mul(r26833, r26828, r26832, MPFR_RNDN); mpfr_div(r26834, r26825, r26833, MPFR_RNDN); ; ; mpfr_mul(r26837, r26836, r26831, MPFR_RNDN); mpfr_sub(r26838, r26835, r26837, MPFR_RNDN); mpfr_sqrt(r26839, r26838, MPFR_RNDN); mpfr_div(r26840, r26834, r26839, MPFR_RNDN); return mpfr_get_d(r26840, MPFR_RNDN); } static mpfr_t r26841, r26842, r26843, r26844, r26845, r26846, r26847, r26848, r26849, r26850, r26851, r26852, r26853, r26854, r26855, r26856; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26841, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26842, "3", 10, MPFR_RNDN); mpfr_init(r26843); mpfr_init(r26844); mpfr_init_set_str(r26845, "1", 10, MPFR_RNDN); mpfr_init(r26846); mpfr_init(r26847); mpfr_init(r26848); mpfr_init(r26849); mpfr_init(r26850); mpfr_init_set_str(r26851, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26852, "6", 10, MPFR_RNDN); mpfr_init(r26853); mpfr_init(r26854); mpfr_init(r26855); mpfr_init(r26856); } double f_dm(double v) { ; ; mpfr_const_pi(r26843, MPFR_RNDN); mpfr_mul(r26844, r26842, r26843, MPFR_RNDN); ; mpfr_set_d(r26846, v, MPFR_RNDN); mpfr_mul(r26847, r26846, r26846, MPFR_RNDN); mpfr_sub(r26848, r26845, r26847, MPFR_RNDN); mpfr_mul(r26849, r26844, r26848, MPFR_RNDN); mpfr_div(r26850, r26841, r26849, MPFR_RNDN); ; ; mpfr_mul(r26853, r26852, r26847, MPFR_RNDN); mpfr_sub(r26854, r26851, r26853, MPFR_RNDN); mpfr_sqrt(r26855, r26854, MPFR_RNDN); mpfr_div(r26856, r26850, r26855, MPFR_RNDN); return mpfr_get_d(r26856, MPFR_RNDN); }