#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18727 = 0.5f; float r18728 = re; float r18729 = cos(r18728); float r18730 = r18727 * r18729; float r18731 = 0.0f; float r18732 = im; float r18733 = r18731 - r18732; float r18734 = exp(r18733); float r18735 = exp(r18732); float r18736 = r18734 - r18735; float r18737 = r18730 * r18736; return r18737; } double f_id(double re, double im) { double r18738 = 0.5; double r18739 = re; double r18740 = cos(r18739); double r18741 = r18738 * r18740; double r18742 = 0.0; double r18743 = im; double r18744 = r18742 - r18743; double r18745 = exp(r18744); double r18746 = exp(r18743); double r18747 = r18745 - r18746; double r18748 = r18741 * r18747; return r18748; } double f_of(float re, float im) { float r18749 = 0.01666666753590107f; float r18750 = im; float r18751 = 5.0f; float r18752 = pow(r18750, r18751); float r18753 = r18749 * r18752; float r18754 = 2.0f; float r18755 = r18754 * r18750; float r18756 = 0.3333333432674408f; float r18757 = r18750 * (r18750 * r18750); float r18758 = r18756 * r18757; float r18759 = r18755 + r18758; float r18760 = r18753 + r18759; float r18761 = -r18760; float r18762 = re; float r18763 = cos(r18762); float r18764 = 0.5f; float r18765 = r18763 * r18764; float r18766 = r18761 * r18765; return r18766; } double f_od(double re, double im) { double r18767 = 0.01666666753590107; double r18768 = im; double r18769 = 5.0; double r18770 = pow(r18768, r18769); double r18771 = r18767 * r18770; double r18772 = 2.0; double r18773 = r18772 * r18768; double r18774 = 0.3333333432674408; double r18775 = r18768 * (r18768 * r18768); double r18776 = r18774 * r18775; double r18777 = r18773 + r18776; double r18778 = r18771 + r18777; double r18779 = -r18778; double r18780 = re; double r18781 = cos(r18780); double r18782 = 0.5; double r18783 = r18781 * r18782; double r18784 = r18779 * r18783; return r18784; } 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 r18785, r18786, r18787, r18788, r18789, r18790, r18791, r18792, r18793, r18794, r18795; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18785, "0.5", 10, MPFR_RNDN); mpfr_init(r18786); mpfr_init(r18787); mpfr_init(r18788); mpfr_init_set_str(r18789, "0", 10, MPFR_RNDN); mpfr_init(r18790); mpfr_init(r18791); mpfr_init(r18792); mpfr_init(r18793); mpfr_init(r18794); mpfr_init(r18795); } double f_im(double re, double im) { ; mpfr_set_d(r18786, re, MPFR_RNDN); mpfr_cos(r18787, r18786, MPFR_RNDN); mpfr_mul(r18788, r18785, r18787, MPFR_RNDN); ; mpfr_set_d(r18790, im, MPFR_RNDN); mpfr_sub(r18791, r18789, r18790, MPFR_RNDN); mpfr_exp(r18792, r18791, MPFR_RNDN); mpfr_exp(r18793, r18790, MPFR_RNDN); mpfr_sub(r18794, r18792, r18793, MPFR_RNDN); mpfr_mul(r18795, r18788, r18794, MPFR_RNDN); return mpfr_get_d(r18795, MPFR_RNDN); } static mpfr_t r18796, r18797, r18798, r18799, r18800, r18801, r18802, r18803, r18804, r18805, r18806, r18807, r18808, r18809, r18810, r18811, r18812, r18813; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18796, "1/60", 10, MPFR_RNDN); mpfr_init(r18797); mpfr_init_set_str(r18798, "5", 10, MPFR_RNDN); mpfr_init(r18799); mpfr_init(r18800); mpfr_init_set_str(r18801, "2", 10, MPFR_RNDN); mpfr_init(r18802); mpfr_init_set_str(r18803, "1/3", 10, MPFR_RNDN); mpfr_init(r18804); mpfr_init(r18805); mpfr_init(r18806); mpfr_init(r18807); mpfr_init(r18808); mpfr_init(r18809); mpfr_init(r18810); mpfr_init_set_str(r18811, "0.5", 10, MPFR_RNDN); mpfr_init(r18812); mpfr_init(r18813); } double f_fm(double re, double im) { ; mpfr_set_d(r18797, im, MPFR_RNDN); ; mpfr_pow(r18799, r18797, r18798, MPFR_RNDN); mpfr_mul(r18800, r18796, r18799, MPFR_RNDN); ; mpfr_mul(r18802, r18801, r18797, MPFR_RNDN); ; mpfr_mul(r18804, r18797, r18797, MPFR_RNDN); mpfr_mul(r18804, r18804, r18797, MPFR_RNDN); mpfr_mul(r18805, r18803, r18804, MPFR_RNDN); mpfr_add(r18806, r18802, r18805, MPFR_RNDN); mpfr_add(r18807, r18800, r18806, MPFR_RNDN); mpfr_neg(r18808, r18807, MPFR_RNDN); mpfr_set_d(r18809, re, MPFR_RNDN); mpfr_cos(r18810, r18809, MPFR_RNDN); ; mpfr_mul(r18812, r18810, r18811, MPFR_RNDN); mpfr_mul(r18813, r18808, 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; void setup_mpfr_f_dm() { 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(r18822); mpfr_init(r18823); mpfr_init(r18824); mpfr_init(r18825); mpfr_init(r18826); mpfr_init(r18827); mpfr_init(r18828); mpfr_init_set_str(r18829, "0.5", 10, MPFR_RNDN); mpfr_init(r18830); mpfr_init(r18831); } double f_dm(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_mul(r18822, r18815, r18815, MPFR_RNDN); mpfr_mul(r18822, r18822, r18815, MPFR_RNDN); mpfr_mul(r18823, r18821, r18822, MPFR_RNDN); mpfr_add(r18824, r18820, r18823, MPFR_RNDN); mpfr_add(r18825, r18818, r18824, MPFR_RNDN); mpfr_neg(r18826, r18825, MPFR_RNDN); mpfr_set_d(r18827, re, MPFR_RNDN); mpfr_cos(r18828, r18827, MPFR_RNDN); ; mpfr_mul(r18830, r18828, r18829, MPFR_RNDN); mpfr_mul(r18831, r18826, r18830, MPFR_RNDN); return mpfr_get_d(r18831, MPFR_RNDN); }