#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18743 = 0.5f; float r18744 = re; float r18745 = cos(r18744); float r18746 = r18743 * r18745; float r18747 = 0.0f; float r18748 = im; float r18749 = r18747 - r18748; float r18750 = exp(r18749); float r18751 = exp(r18748); float r18752 = r18750 - r18751; float r18753 = r18746 * r18752; return r18753; } double f_id(double re, double im) { double r18754 = 0.5; double r18755 = re; double r18756 = cos(r18755); double r18757 = r18754 * r18756; double r18758 = 0.0; double r18759 = im; double r18760 = r18758 - r18759; double r18761 = exp(r18760); double r18762 = exp(r18759); double r18763 = r18761 - r18762; double r18764 = r18757 * r18763; return r18764; } double f_of(float re, float im) { float r18765 = 0.01666666753590107f; float r18766 = im; float r18767 = 5.0f; float r18768 = pow(r18766, r18767); float r18769 = r18765 * r18768; float r18770 = 2.0f; float r18771 = r18770 * r18766; float r18772 = 0.3333333432674408f; float r18773 = 3.0f; float r18774 = pow(r18766, r18773); float r18775 = r18772 * r18774; float r18776 = r18771 + r18775; float r18777 = r18769 + r18776; float r18778 = -r18777; float r18779 = re; float r18780 = cos(r18779); float r18781 = 0.5f; float r18782 = r18780 * r18781; float r18783 = r18778 * r18782; return r18783; } double f_od(double re, double im) { double r18784 = 0.01666666753590107; double r18785 = im; double r18786 = 5.0; double r18787 = pow(r18785, r18786); double r18788 = r18784 * r18787; double r18789 = 2.0; double r18790 = r18789 * r18785; double r18791 = 0.3333333432674408; double r18792 = 3.0; double r18793 = pow(r18785, r18792); double r18794 = r18791 * r18793; double r18795 = r18790 + r18794; double r18796 = r18788 + r18795; double r18797 = -r18796; double r18798 = re; double r18799 = cos(r18798); double r18800 = 0.5; double r18801 = r18799 * r18800; double r18802 = r18797 * r18801; return r18802; } 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 r18803, r18804, r18805, r18806, r18807, r18808, r18809, r18810, r18811, r18812, r18813; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18803, "0.5", 10, MPFR_RNDN); mpfr_init(r18804); mpfr_init(r18805); mpfr_init(r18806); mpfr_init_set_str(r18807, "0", 10, MPFR_RNDN); mpfr_init(r18808); mpfr_init(r18809); mpfr_init(r18810); mpfr_init(r18811); mpfr_init(r18812); mpfr_init(r18813); } double f_im(double re, double im) { ; mpfr_set_d(r18804, re, MPFR_RNDN); mpfr_cos(r18805, r18804, MPFR_RNDN); mpfr_mul(r18806, r18803, r18805, MPFR_RNDN); ; mpfr_set_d(r18808, im, MPFR_RNDN); mpfr_sub(r18809, r18807, r18808, MPFR_RNDN); mpfr_exp(r18810, r18809, MPFR_RNDN); mpfr_exp(r18811, r18808, MPFR_RNDN); mpfr_sub(r18812, r18810, r18811, MPFR_RNDN); mpfr_mul(r18813, r18806, r18812, MPFR_RNDN); return mpfr_get_d(r18813, MPFR_RNDN); } static mpfr_t r18814, r18815, r18816, r18817, r18818, r18819, r18820, r18821, r18822, r18823, r18824, r18825, r18826, r18827, r18828, r18829, r18830, r18831, r18832; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18814, "1/60", 10, MPFR_RNDN); mpfr_init(r18815); mpfr_init_set_str(r18816, "5", 10, MPFR_RNDN); mpfr_init(r18817); mpfr_init(r18818); mpfr_init_set_str(r18819, "2", 10, MPFR_RNDN); mpfr_init(r18820); mpfr_init_set_str(r18821, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18822, "3", 10, MPFR_RNDN); mpfr_init(r18823); mpfr_init(r18824); mpfr_init(r18825); mpfr_init(r18826); mpfr_init(r18827); mpfr_init(r18828); mpfr_init(r18829); mpfr_init_set_str(r18830, "0.5", 10, MPFR_RNDN); mpfr_init(r18831); mpfr_init(r18832); } double f_fm(double re, double im) { ; mpfr_set_d(r18815, im, MPFR_RNDN); ; mpfr_pow(r18817, r18815, r18816, MPFR_RNDN); mpfr_mul(r18818, r18814, r18817, MPFR_RNDN); ; mpfr_mul(r18820, r18819, r18815, MPFR_RNDN); ; ; mpfr_pow(r18823, r18815, r18822, MPFR_RNDN); mpfr_mul(r18824, r18821, r18823, MPFR_RNDN); mpfr_add(r18825, r18820, r18824, MPFR_RNDN); mpfr_add(r18826, r18818, r18825, MPFR_RNDN); mpfr_neg(r18827, r18826, MPFR_RNDN); mpfr_set_d(r18828, re, MPFR_RNDN); mpfr_cos(r18829, r18828, MPFR_RNDN); ; mpfr_mul(r18831, r18829, r18830, MPFR_RNDN); mpfr_mul(r18832, r18827, r18831, MPFR_RNDN); return mpfr_get_d(r18832, MPFR_RNDN); } static mpfr_t r18833, r18834, r18835, r18836, r18837, r18838, r18839, r18840, r18841, r18842, r18843, r18844, r18845, r18846, r18847, r18848, r18849, r18850, r18851; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18833, "1/60", 10, MPFR_RNDN); mpfr_init(r18834); mpfr_init_set_str(r18835, "5", 10, MPFR_RNDN); mpfr_init(r18836); mpfr_init(r18837); mpfr_init_set_str(r18838, "2", 10, MPFR_RNDN); mpfr_init(r18839); mpfr_init_set_str(r18840, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18841, "3", 10, MPFR_RNDN); mpfr_init(r18842); mpfr_init(r18843); mpfr_init(r18844); mpfr_init(r18845); mpfr_init(r18846); mpfr_init(r18847); mpfr_init(r18848); mpfr_init_set_str(r18849, "0.5", 10, MPFR_RNDN); mpfr_init(r18850); mpfr_init(r18851); } double f_dm(double re, double im) { ; mpfr_set_d(r18834, im, MPFR_RNDN); ; mpfr_pow(r18836, r18834, r18835, MPFR_RNDN); mpfr_mul(r18837, r18833, r18836, MPFR_RNDN); ; mpfr_mul(r18839, r18838, r18834, MPFR_RNDN); ; ; mpfr_pow(r18842, r18834, r18841, MPFR_RNDN); mpfr_mul(r18843, r18840, r18842, MPFR_RNDN); mpfr_add(r18844, r18839, r18843, MPFR_RNDN); mpfr_add(r18845, r18837, r18844, MPFR_RNDN); mpfr_neg(r18846, r18845, MPFR_RNDN); mpfr_set_d(r18847, re, MPFR_RNDN); mpfr_cos(r18848, r18847, MPFR_RNDN); ; mpfr_mul(r18850, r18848, r18849, MPFR_RNDN); mpfr_mul(r18851, r18846, r18850, MPFR_RNDN); return mpfr_get_d(r18851, MPFR_RNDN); }