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