#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19809 = 0.5f; float r19810 = re; float r19811 = sin(r19810); float r19812 = r19809 * r19811; float r19813 = im; float r19814 = -r19813; float r19815 = exp(r19814); float r19816 = exp(r19813); float r19817 = r19815 - r19816; float r19818 = r19812 * r19817; return r19818; } double f_id(double re, double im) { double r19819 = 0.5; double r19820 = re; double r19821 = sin(r19820); double r19822 = r19819 * r19821; double r19823 = im; double r19824 = -r19823; double r19825 = exp(r19824); double r19826 = exp(r19823); double r19827 = r19825 - r19826; double r19828 = r19822 * r19827; return r19828; } double f_of(float re, float im) { float r19829 = 0.5f; float r19830 = re; float r19831 = sin(r19830); float r19832 = r19829 * r19831; float r19833 = 0.01666666753590107f; float r19834 = im; float r19835 = 5.0f; float r19836 = pow(r19834, r19835); float r19837 = r19833 * r19836; float r19838 = 2.0f; float r19839 = r19838 * r19834; float r19840 = 0.3333333432674408f; float r19841 = 3.0f; float r19842 = pow(r19834, r19841); float r19843 = r19840 * r19842; float r19844 = r19839 + r19843; float r19845 = r19837 + r19844; float r19846 = -r19845; float r19847 = r19832 * r19846; return r19847; } double f_od(double re, double im) { double r19848 = 0.5; double r19849 = re; double r19850 = sin(r19849); double r19851 = r19848 * r19850; double r19852 = 0.01666666753590107; double r19853 = im; double r19854 = 5.0; double r19855 = pow(r19853, r19854); double r19856 = r19852 * r19855; double r19857 = 2.0; double r19858 = r19857 * r19853; double r19859 = 0.3333333432674408; double r19860 = 3.0; double r19861 = pow(r19853, r19860); double r19862 = r19859 * r19861; double r19863 = r19858 + r19862; double r19864 = r19856 + r19863; double r19865 = -r19864; double r19866 = r19851 * r19865; return r19866; } 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 r19867, r19868, r19869, r19870, r19871, r19872, r19873, r19874, r19875, r19876; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19867, "0.5", 10, MPFR_RNDN); mpfr_init(r19868); mpfr_init(r19869); mpfr_init(r19870); mpfr_init(r19871); mpfr_init(r19872); mpfr_init(r19873); mpfr_init(r19874); mpfr_init(r19875); mpfr_init(r19876); } double f_im(double re, double im) { ; mpfr_set_d(r19868, re, MPFR_RNDN); mpfr_sin(r19869, r19868, MPFR_RNDN); mpfr_mul(r19870, r19867, r19869, MPFR_RNDN); mpfr_set_d(r19871, im, MPFR_RNDN); mpfr_neg(r19872, r19871, MPFR_RNDN); mpfr_exp(r19873, r19872, MPFR_RNDN); mpfr_exp(r19874, r19871, MPFR_RNDN); mpfr_sub(r19875, r19873, r19874, MPFR_RNDN); mpfr_mul(r19876, r19870, r19875, MPFR_RNDN); return mpfr_get_d(r19876, MPFR_RNDN); } static mpfr_t r19877, r19878, r19879, r19880, r19881, r19882, r19883, r19884, r19885, r19886, r19887, r19888, r19889, r19890, r19891, r19892, r19893, r19894, r19895; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19877, "0.5", 10, MPFR_RNDN); mpfr_init(r19878); mpfr_init(r19879); mpfr_init(r19880); mpfr_init_set_str(r19881, "1/60", 10, MPFR_RNDN); mpfr_init(r19882); mpfr_init_set_str(r19883, "5", 10, MPFR_RNDN); mpfr_init(r19884); mpfr_init(r19885); mpfr_init_set_str(r19886, "2", 10, MPFR_RNDN); mpfr_init(r19887); mpfr_init_set_str(r19888, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19889, "3", 10, MPFR_RNDN); mpfr_init(r19890); mpfr_init(r19891); mpfr_init(r19892); mpfr_init(r19893); mpfr_init(r19894); mpfr_init(r19895); } double f_fm(double re, double im) { ; mpfr_set_d(r19878, re, MPFR_RNDN); mpfr_sin(r19879, r19878, MPFR_RNDN); mpfr_mul(r19880, r19877, r19879, MPFR_RNDN); ; mpfr_set_d(r19882, im, MPFR_RNDN); ; mpfr_pow(r19884, r19882, r19883, MPFR_RNDN); mpfr_mul(r19885, r19881, r19884, MPFR_RNDN); ; mpfr_mul(r19887, r19886, r19882, MPFR_RNDN); ; ; mpfr_pow(r19890, r19882, r19889, MPFR_RNDN); mpfr_mul(r19891, r19888, r19890, MPFR_RNDN); mpfr_add(r19892, r19887, r19891, MPFR_RNDN); mpfr_add(r19893, r19885, r19892, MPFR_RNDN); mpfr_neg(r19894, r19893, MPFR_RNDN); mpfr_mul(r19895, r19880, r19894, MPFR_RNDN); return mpfr_get_d(r19895, MPFR_RNDN); } static mpfr_t r19896, r19897, r19898, r19899, r19900, r19901, r19902, r19903, r19904, r19905, r19906, r19907, r19908, r19909, r19910, r19911, r19912, r19913, r19914; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19896, "0.5", 10, MPFR_RNDN); mpfr_init(r19897); mpfr_init(r19898); mpfr_init(r19899); mpfr_init_set_str(r19900, "1/60", 10, MPFR_RNDN); mpfr_init(r19901); mpfr_init_set_str(r19902, "5", 10, MPFR_RNDN); mpfr_init(r19903); mpfr_init(r19904); mpfr_init_set_str(r19905, "2", 10, MPFR_RNDN); mpfr_init(r19906); mpfr_init_set_str(r19907, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19908, "3", 10, MPFR_RNDN); mpfr_init(r19909); mpfr_init(r19910); mpfr_init(r19911); mpfr_init(r19912); mpfr_init(r19913); mpfr_init(r19914); } double f_dm(double re, double im) { ; mpfr_set_d(r19897, re, MPFR_RNDN); mpfr_sin(r19898, r19897, MPFR_RNDN); mpfr_mul(r19899, r19896, r19898, MPFR_RNDN); ; mpfr_set_d(r19901, im, MPFR_RNDN); ; mpfr_pow(r19903, r19901, r19902, MPFR_RNDN); mpfr_mul(r19904, r19900, r19903, MPFR_RNDN); ; mpfr_mul(r19906, r19905, r19901, MPFR_RNDN); ; ; mpfr_pow(r19909, r19901, r19908, MPFR_RNDN); mpfr_mul(r19910, r19907, r19909, MPFR_RNDN); mpfr_add(r19911, r19906, r19910, MPFR_RNDN); mpfr_add(r19912, r19904, r19911, MPFR_RNDN); mpfr_neg(r19913, r19912, MPFR_RNDN); mpfr_mul(r19914, r19899, r19913, MPFR_RNDN); return mpfr_get_d(r19914, MPFR_RNDN); }