#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19785 = 0.5f; float r19786 = re; float r19787 = sin(r19786); float r19788 = r19785 * r19787; float r19789 = im; float r19790 = -r19789; float r19791 = exp(r19790); float r19792 = exp(r19789); float r19793 = r19791 - r19792; float r19794 = r19788 * r19793; return r19794; } double f_id(double re, double im) { double r19795 = 0.5; double r19796 = re; double r19797 = sin(r19796); double r19798 = r19795 * r19797; double r19799 = im; double r19800 = -r19799; double r19801 = exp(r19800); double r19802 = exp(r19799); double r19803 = r19801 - r19802; double r19804 = r19798 * r19803; return r19804; } double f_of(float re, float im) { float r19805 = 0.5f; float r19806 = re; float r19807 = sin(r19806); float r19808 = r19805 * r19807; float r19809 = 0.016666666666666666f; float r19810 = im; float r19811 = 5.0f; float r19812 = pow(r19810, r19811); float r19813 = r19809 * r19812; float r19814 = 2.0f; float r19815 = r19814 * r19810; float r19816 = 0.3333333333333333f; float r19817 = 3.0f; float r19818 = pow(r19810, r19817); float r19819 = r19816 * r19818; float r19820 = r19815 + r19819; float r19821 = r19813 + r19820; float r19822 = -r19821; float r19823 = r19808 * r19822; return r19823; } double f_od(double re, double im) { double r19824 = 0.5; double r19825 = re; double r19826 = sin(r19825); double r19827 = r19824 * r19826; double r19828 = 0.016666666666666666; double r19829 = im; double r19830 = 5.0; double r19831 = pow(r19829, r19830); double r19832 = r19828 * r19831; double r19833 = 2.0; double r19834 = r19833 * r19829; double r19835 = 0.3333333333333333; double r19836 = 3.0; double r19837 = pow(r19829, r19836); double r19838 = r19835 * r19837; double r19839 = r19834 + r19838; double r19840 = r19832 + r19839; double r19841 = -r19840; double r19842 = r19827 * r19841; return r19842; } 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 r19843, r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851, r19852; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19843, "0.5", 10, MPFR_RNDN); mpfr_init(r19844); mpfr_init(r19845); mpfr_init(r19846); mpfr_init(r19847); mpfr_init(r19848); mpfr_init(r19849); mpfr_init(r19850); mpfr_init(r19851); mpfr_init(r19852); } double f_im(double re, double im) { ; mpfr_set_d(r19844, re, MPFR_RNDN); mpfr_sin(r19845, r19844, MPFR_RNDN); mpfr_mul(r19846, r19843, r19845, MPFR_RNDN); mpfr_set_d(r19847, im, MPFR_RNDN); mpfr_neg(r19848, r19847, MPFR_RNDN); mpfr_exp(r19849, r19848, MPFR_RNDN); mpfr_exp(r19850, r19847, MPFR_RNDN); mpfr_sub(r19851, r19849, r19850, MPFR_RNDN); mpfr_mul(r19852, r19846, r19851, MPFR_RNDN); return mpfr_get_d(r19852, MPFR_RNDN); } static mpfr_t r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866, r19867, r19868, r19869, r19870, r19871; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19853, "0.5", 10, MPFR_RNDN); mpfr_init(r19854); mpfr_init(r19855); mpfr_init(r19856); mpfr_init_set_str(r19857, "1/60", 10, MPFR_RNDN); mpfr_init(r19858); mpfr_init_set_str(r19859, "5", 10, MPFR_RNDN); mpfr_init(r19860); mpfr_init(r19861); mpfr_init_set_str(r19862, "2", 10, MPFR_RNDN); mpfr_init(r19863); mpfr_init_set_str(r19864, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19865, "3", 10, MPFR_RNDN); mpfr_init(r19866); mpfr_init(r19867); mpfr_init(r19868); mpfr_init(r19869); mpfr_init(r19870); mpfr_init(r19871); } double f_fm(double re, double im) { ; mpfr_set_d(r19854, re, MPFR_RNDN); mpfr_sin(r19855, r19854, MPFR_RNDN); mpfr_mul(r19856, r19853, r19855, MPFR_RNDN); ; mpfr_set_d(r19858, im, MPFR_RNDN); ; mpfr_pow(r19860, r19858, r19859, MPFR_RNDN); mpfr_mul(r19861, r19857, r19860, MPFR_RNDN); ; mpfr_mul(r19863, r19862, r19858, MPFR_RNDN); ; ; mpfr_pow(r19866, r19858, r19865, MPFR_RNDN); mpfr_mul(r19867, r19864, r19866, MPFR_RNDN); mpfr_add(r19868, r19863, r19867, MPFR_RNDN); mpfr_add(r19869, r19861, r19868, MPFR_RNDN); mpfr_neg(r19870, r19869, MPFR_RNDN); mpfr_mul(r19871, r19856, r19870, MPFR_RNDN); return mpfr_get_d(r19871, MPFR_RNDN); } static mpfr_t r19872, r19873, r19874, r19875, r19876, r19877, r19878, r19879, r19880, r19881, r19882, r19883, r19884, r19885, r19886, r19887, r19888, r19889, r19890; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19872, "0.5", 10, MPFR_RNDN); mpfr_init(r19873); mpfr_init(r19874); mpfr_init(r19875); mpfr_init_set_str(r19876, "1/60", 10, MPFR_RNDN); mpfr_init(r19877); mpfr_init_set_str(r19878, "5", 10, MPFR_RNDN); mpfr_init(r19879); mpfr_init(r19880); mpfr_init_set_str(r19881, "2", 10, MPFR_RNDN); mpfr_init(r19882); mpfr_init_set_str(r19883, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19884, "3", 10, MPFR_RNDN); mpfr_init(r19885); mpfr_init(r19886); mpfr_init(r19887); mpfr_init(r19888); mpfr_init(r19889); mpfr_init(r19890); } double f_dm(double re, double im) { ; mpfr_set_d(r19873, re, MPFR_RNDN); mpfr_sin(r19874, r19873, MPFR_RNDN); mpfr_mul(r19875, r19872, r19874, MPFR_RNDN); ; mpfr_set_d(r19877, im, MPFR_RNDN); ; mpfr_pow(r19879, r19877, r19878, MPFR_RNDN); mpfr_mul(r19880, r19876, r19879, MPFR_RNDN); ; mpfr_mul(r19882, r19881, r19877, MPFR_RNDN); ; ; mpfr_pow(r19885, r19877, r19884, MPFR_RNDN); mpfr_mul(r19886, r19883, r19885, MPFR_RNDN); mpfr_add(r19887, r19882, r19886, MPFR_RNDN); mpfr_add(r19888, r19880, r19887, MPFR_RNDN); mpfr_neg(r19889, r19888, MPFR_RNDN); mpfr_mul(r19890, r19875, r19889, MPFR_RNDN); return mpfr_get_d(r19890, MPFR_RNDN); }