#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r21701 = x; float r21702 = 1.0f; float r21703 = B; float r21704 = tan(r21703); float r21705 = r21702 / r21704; float r21706 = r21701 * r21705; float r21707 = -r21706; float r21708 = F; float r21709 = sin(r21703); float r21710 = r21708 / r21709; float r21711 = r21708 * r21708; float r21712 = 2.0f; float r21713 = r21711 + r21712; float r21714 = r21712 * r21701; float r21715 = r21713 + r21714; float r21716 = r21702 / r21712; float r21717 = -r21716; float r21718 = pow(r21715, r21717); float r21719 = r21710 * r21718; float r21720 = r21707 + r21719; return r21720; } double f_id(double F, double B, double x) { double r21721 = x; double r21722 = 1.0; double r21723 = B; double r21724 = tan(r21723); double r21725 = r21722 / r21724; double r21726 = r21721 * r21725; double r21727 = -r21726; double r21728 = F; double r21729 = sin(r21723); double r21730 = r21728 / r21729; double r21731 = r21728 * r21728; double r21732 = 2.0; double r21733 = r21731 + r21732; double r21734 = r21732 * r21721; double r21735 = r21733 + r21734; double r21736 = r21722 / r21732; double r21737 = -r21736; double r21738 = pow(r21735, r21737); double r21739 = r21730 * r21738; double r21740 = r21727 + r21739; return r21740; } double f_of(float F, float B, float x) { float r21741 = x; float r21742 = -r21741; float r21743 = B; float r21744 = tan(r21743); float r21745 = r21742 / r21744; float r21746 = F; float r21747 = sin(r21743); float r21748 = r21741 + r21741; float r21749 = r21746 * r21746; float r21750 = 2.0f; float r21751 = r21749 + r21750; float r21752 = r21748 + r21751; float r21753 = 1.0f; float r21754 = r21753 / r21750; float r21755 = pow(r21752, r21754); float r21756 = r21747 * r21755; float r21757 = r21746 / r21756; float r21758 = r21745 + r21757; return r21758; } double f_od(double F, double B, double x) { double r21759 = x; double r21760 = -r21759; double r21761 = B; double r21762 = tan(r21761); double r21763 = r21760 / r21762; double r21764 = F; double r21765 = sin(r21761); double r21766 = r21759 + r21759; double r21767 = r21764 * r21764; double r21768 = 2.0; double r21769 = r21767 + r21768; double r21770 = r21766 + r21769; double r21771 = 1.0; double r21772 = r21771 / r21768; double r21773 = pow(r21770, r21772); double r21774 = r21765 * r21773; double r21775 = r21764 / r21774; double r21776 = r21763 + r21775; return r21776; } 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 r21777, r21778, r21779, r21780, r21781, r21782, r21783, r21784, r21785, r21786, r21787, r21788, r21789, r21790, r21791, r21792, r21793, r21794, r21795, r21796; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r21777); mpfr_init_set_str(r21778, "1", 10, MPFR_RNDN); mpfr_init(r21779); mpfr_init(r21780); mpfr_init(r21781); mpfr_init(r21782); mpfr_init(r21783); mpfr_init(r21784); mpfr_init(r21785); mpfr_init(r21786); mpfr_init(r21787); mpfr_init_set_str(r21788, "2", 10, MPFR_RNDN); mpfr_init(r21789); mpfr_init(r21790); mpfr_init(r21791); mpfr_init(r21792); mpfr_init(r21793); mpfr_init(r21794); mpfr_init(r21795); mpfr_init(r21796); } double f_im(double F, double B, double x) { mpfr_set_d(r21777, x, MPFR_RNDN); ; mpfr_set_d(r21779, B, MPFR_RNDN); mpfr_tan(r21780, r21779, MPFR_RNDN); mpfr_div(r21781, r21778, r21780, MPFR_RNDN); mpfr_mul(r21782, r21777, r21781, MPFR_RNDN); mpfr_neg(r21783, r21782, MPFR_RNDN); mpfr_set_d(r21784, F, MPFR_RNDN); mpfr_sin(r21785, r21779, MPFR_RNDN); mpfr_div(r21786, r21784, r21785, MPFR_RNDN); mpfr_sqr(r21787, r21784, MPFR_RNDN); ; mpfr_add(r21789, r21787, r21788, MPFR_RNDN); mpfr_mul(r21790, r21788, r21777, MPFR_RNDN); mpfr_add(r21791, r21789, r21790, MPFR_RNDN); mpfr_div(r21792, r21778, r21788, MPFR_RNDN); mpfr_neg(r21793, r21792, MPFR_RNDN); mpfr_pow(r21794, r21791, r21793, MPFR_RNDN); mpfr_mul(r21795, r21786, r21794, MPFR_RNDN); mpfr_add(r21796, r21783, r21795, MPFR_RNDN); return mpfr_get_d(r21796, MPFR_RNDN); } static mpfr_t r21797, r21798, r21799, r21800, r21801, r21802, r21803, r21804, r21805, r21806, r21807, r21808, r21809, r21810, r21811, r21812, r21813, r21814; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r21797); mpfr_init(r21798); mpfr_init(r21799); mpfr_init(r21800); mpfr_init(r21801); mpfr_init(r21802); mpfr_init(r21803); mpfr_init(r21804); mpfr_init(r21805); mpfr_init_set_str(r21806, "2", 10, MPFR_RNDN); mpfr_init(r21807); mpfr_init(r21808); mpfr_init_set_str(r21809, "1", 10, MPFR_RNDN); mpfr_init(r21810); mpfr_init(r21811); mpfr_init(r21812); mpfr_init(r21813); mpfr_init(r21814); } double f_fm(double F, double B, double x) { mpfr_set_d(r21797, x, MPFR_RNDN); mpfr_neg(r21798, r21797, MPFR_RNDN); mpfr_set_d(r21799, B, MPFR_RNDN); mpfr_tan(r21800, r21799, MPFR_RNDN); mpfr_div(r21801, r21798, r21800, MPFR_RNDN); mpfr_set_d(r21802, F, MPFR_RNDN); mpfr_sin(r21803, r21799, MPFR_RNDN); mpfr_add(r21804, r21797, r21797, MPFR_RNDN); mpfr_sqr(r21805, r21802, MPFR_RNDN); ; mpfr_add(r21807, r21805, r21806, MPFR_RNDN); mpfr_add(r21808, r21804, r21807, MPFR_RNDN); ; mpfr_div(r21810, r21809, r21806, MPFR_RNDN); mpfr_pow(r21811, r21808, r21810, MPFR_RNDN); mpfr_mul(r21812, r21803, r21811, MPFR_RNDN); mpfr_div(r21813, r21802, r21812, MPFR_RNDN); mpfr_add(r21814, r21801, r21813, MPFR_RNDN); return mpfr_get_d(r21814, MPFR_RNDN); } static mpfr_t r21815, r21816, r21817, r21818, r21819, r21820, r21821, r21822, r21823, r21824, r21825, r21826, r21827, r21828, r21829, r21830, r21831, r21832; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r21815); mpfr_init(r21816); mpfr_init(r21817); mpfr_init(r21818); mpfr_init(r21819); mpfr_init(r21820); mpfr_init(r21821); mpfr_init(r21822); mpfr_init(r21823); mpfr_init_set_str(r21824, "2", 10, MPFR_RNDN); mpfr_init(r21825); mpfr_init(r21826); mpfr_init_set_str(r21827, "1", 10, MPFR_RNDN); mpfr_init(r21828); mpfr_init(r21829); mpfr_init(r21830); mpfr_init(r21831); mpfr_init(r21832); } double f_dm(double F, double B, double x) { mpfr_set_d(r21815, x, MPFR_RNDN); mpfr_neg(r21816, r21815, MPFR_RNDN); mpfr_set_d(r21817, B, MPFR_RNDN); mpfr_tan(r21818, r21817, MPFR_RNDN); mpfr_div(r21819, r21816, r21818, MPFR_RNDN); mpfr_set_d(r21820, F, MPFR_RNDN); mpfr_sin(r21821, r21817, MPFR_RNDN); mpfr_add(r21822, r21815, r21815, MPFR_RNDN); mpfr_sqr(r21823, r21820, MPFR_RNDN); ; mpfr_add(r21825, r21823, r21824, MPFR_RNDN); mpfr_add(r21826, r21822, r21825, MPFR_RNDN); ; mpfr_div(r21828, r21827, r21824, MPFR_RNDN); mpfr_pow(r21829, r21826, r21828, MPFR_RNDN); mpfr_mul(r21830, r21821, r21829, MPFR_RNDN); mpfr_div(r21831, r21820, r21830, MPFR_RNDN); mpfr_add(r21832, r21819, r21831, MPFR_RNDN); return mpfr_get_d(r21832, MPFR_RNDN); }