#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r19767 = x; float r19768 = 1; float r19769 = B; float r19770 = tan(r19769); float r19771 = r19768 / r19770; float r19772 = r19767 * r19771; float r19773 = -r19772; float r19774 = F; float r19775 = sin(r19769); float r19776 = r19774 / r19775; float r19777 = r19774 * r19774; float r19778 = 2; float r19779 = r19777 + r19778; float r19780 = r19778 * r19767; float r19781 = r19779 + r19780; float r19782 = r19768 / r19778; float r19783 = -r19782; float r19784 = pow(r19781, r19783); float r19785 = r19776 * r19784; float r19786 = r19773 + r19785; return r19786; } double f_id(double F, double B, double x) { double r19787 = x; double r19788 = 1; double r19789 = B; double r19790 = tan(r19789); double r19791 = r19788 / r19790; double r19792 = r19787 * r19791; double r19793 = -r19792; double r19794 = F; double r19795 = sin(r19789); double r19796 = r19794 / r19795; double r19797 = r19794 * r19794; double r19798 = 2; double r19799 = r19797 + r19798; double r19800 = r19798 * r19787; double r19801 = r19799 + r19800; double r19802 = r19788 / r19798; double r19803 = -r19802; double r19804 = pow(r19801, r19803); double r19805 = r19796 * r19804; double r19806 = r19793 + r19805; return r19806; } double f_of(float F, float B, float x) { float r19807 = x; float r19808 = -r19807; float r19809 = B; float r19810 = tan(r19809); float r19811 = r19808 / r19810; float r19812 = r19807 + r19807; float r19813 = 2; float r19814 = F; float r19815 = r19814 * r19814; float r19816 = r19813 + r19815; float r19817 = r19812 + r19816; float r19818 = 1; float r19819 = r19818 / r19813; float r19820 = -r19819; float r19821 = pow(r19817, r19820); float r19822 = r19821 * r19814; float r19823 = sin(r19809); float r19824 = r19818 / r19823; float r19825 = r19822 * r19824; float r19826 = r19811 + r19825; return r19826; } double f_od(double F, double B, double x) { double r19827 = x; double r19828 = -r19827; double r19829 = B; double r19830 = tan(r19829); double r19831 = r19828 / r19830; double r19832 = r19827 + r19827; double r19833 = 2; double r19834 = F; double r19835 = r19834 * r19834; double r19836 = r19833 + r19835; double r19837 = r19832 + r19836; double r19838 = 1; double r19839 = r19838 / r19833; double r19840 = -r19839; double r19841 = pow(r19837, r19840); double r19842 = r19841 * r19834; double r19843 = sin(r19829); double r19844 = r19838 / r19843; double r19845 = r19842 * r19844; double r19846 = r19831 + r19845; return r19846; } 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 r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r19847); mpfr_init_set_str(r19848, "1", 10, MPFR_RNDN); mpfr_init(r19849); mpfr_init(r19850); mpfr_init(r19851); mpfr_init(r19852); mpfr_init(r19853); mpfr_init(r19854); mpfr_init(r19855); mpfr_init(r19856); mpfr_init(r19857); mpfr_init_set_str(r19858, "2", 10, MPFR_RNDN); mpfr_init(r19859); mpfr_init(r19860); mpfr_init(r19861); mpfr_init(r19862); mpfr_init(r19863); mpfr_init(r19864); mpfr_init(r19865); mpfr_init(r19866); } double f_im(double F, double B, double x) { mpfr_set_d(r19847, x, MPFR_RNDN); ; mpfr_set_d(r19849, B, MPFR_RNDN); mpfr_tan(r19850, r19849, MPFR_RNDN); mpfr_div(r19851, r19848, r19850, MPFR_RNDN); mpfr_mul(r19852, r19847, r19851, MPFR_RNDN); mpfr_neg(r19853, r19852, MPFR_RNDN); mpfr_set_d(r19854, F, MPFR_RNDN); mpfr_sin(r19855, r19849, MPFR_RNDN); mpfr_div(r19856, r19854, r19855, MPFR_RNDN); mpfr_mul(r19857, r19854, r19854, MPFR_RNDN); ; mpfr_add(r19859, r19857, r19858, MPFR_RNDN); mpfr_mul(r19860, r19858, r19847, MPFR_RNDN); mpfr_add(r19861, r19859, r19860, MPFR_RNDN); mpfr_div(r19862, r19848, r19858, MPFR_RNDN); mpfr_neg(r19863, r19862, MPFR_RNDN); mpfr_pow(r19864, r19861, r19863, MPFR_RNDN); mpfr_mul(r19865, r19856, r19864, MPFR_RNDN); mpfr_add(r19866, r19853, r19865, MPFR_RNDN); return mpfr_get_d(r19866, MPFR_RNDN); } static mpfr_t r19867, r19868, r19869, r19870, r19871, r19872, r19873, r19874, r19875, r19876, r19877, r19878, r19879, r19880, r19881, r19882, r19883, r19884, r19885, r19886; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r19867); mpfr_init(r19868); mpfr_init(r19869); mpfr_init(r19870); mpfr_init(r19871); mpfr_init(r19872); mpfr_init_set_str(r19873, "2", 10, MPFR_RNDN); mpfr_init(r19874); mpfr_init(r19875); mpfr_init(r19876); mpfr_init(r19877); mpfr_init_set_str(r19878, "1", 10, MPFR_RNDN); mpfr_init(r19879); mpfr_init(r19880); mpfr_init(r19881); mpfr_init(r19882); mpfr_init(r19883); mpfr_init(r19884); mpfr_init(r19885); mpfr_init(r19886); } double f_fm(double F, double B, double x) { mpfr_set_d(r19867, x, MPFR_RNDN); mpfr_neg(r19868, r19867, MPFR_RNDN); mpfr_set_d(r19869, B, MPFR_RNDN); mpfr_tan(r19870, r19869, MPFR_RNDN); mpfr_div(r19871, r19868, r19870, MPFR_RNDN); mpfr_add(r19872, r19867, r19867, MPFR_RNDN); ; mpfr_set_d(r19874, F, MPFR_RNDN); mpfr_mul(r19875, r19874, r19874, MPFR_RNDN); mpfr_add(r19876, r19873, r19875, MPFR_RNDN); mpfr_add(r19877, r19872, r19876, MPFR_RNDN); ; mpfr_div(r19879, r19878, r19873, MPFR_RNDN); mpfr_neg(r19880, r19879, MPFR_RNDN); mpfr_pow(r19881, r19877, r19880, MPFR_RNDN); mpfr_mul(r19882, r19881, r19874, MPFR_RNDN); mpfr_sin(r19883, r19869, MPFR_RNDN); mpfr_div(r19884, r19878, r19883, MPFR_RNDN); mpfr_mul(r19885, r19882, r19884, MPFR_RNDN); mpfr_add(r19886, r19871, r19885, MPFR_RNDN); return mpfr_get_d(r19886, MPFR_RNDN); } static mpfr_t r19887, r19888, r19889, r19890, r19891, r19892, r19893, r19894, r19895, r19896, r19897, r19898, r19899, r19900, r19901, r19902, r19903, r19904, r19905, r19906; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r19887); mpfr_init(r19888); mpfr_init(r19889); mpfr_init(r19890); mpfr_init(r19891); mpfr_init(r19892); mpfr_init_set_str(r19893, "2", 10, MPFR_RNDN); mpfr_init(r19894); mpfr_init(r19895); mpfr_init(r19896); mpfr_init(r19897); mpfr_init_set_str(r19898, "1", 10, MPFR_RNDN); mpfr_init(r19899); mpfr_init(r19900); mpfr_init(r19901); mpfr_init(r19902); mpfr_init(r19903); mpfr_init(r19904); mpfr_init(r19905); mpfr_init(r19906); } double f_dm(double F, double B, double x) { mpfr_set_d(r19887, x, MPFR_RNDN); mpfr_neg(r19888, r19887, MPFR_RNDN); mpfr_set_d(r19889, B, MPFR_RNDN); mpfr_tan(r19890, r19889, MPFR_RNDN); mpfr_div(r19891, r19888, r19890, MPFR_RNDN); mpfr_add(r19892, r19887, r19887, MPFR_RNDN); ; mpfr_set_d(r19894, F, MPFR_RNDN); mpfr_mul(r19895, r19894, r19894, MPFR_RNDN); mpfr_add(r19896, r19893, r19895, MPFR_RNDN); mpfr_add(r19897, r19892, r19896, MPFR_RNDN); ; mpfr_div(r19899, r19898, r19893, MPFR_RNDN); mpfr_neg(r19900, r19899, MPFR_RNDN); mpfr_pow(r19901, r19897, r19900, MPFR_RNDN); mpfr_mul(r19902, r19901, r19894, MPFR_RNDN); mpfr_sin(r19903, r19889, MPFR_RNDN); mpfr_div(r19904, r19898, r19903, MPFR_RNDN); mpfr_mul(r19905, r19902, r19904, MPFR_RNDN); mpfr_add(r19906, r19891, r19905, MPFR_RNDN); return mpfr_get_d(r19906, MPFR_RNDN); }