#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18699 = 0.5f; float r18700 = re; float r18701 = cos(r18700); float r18702 = r18699 * r18701; float r18703 = 0.0f; float r18704 = im; float r18705 = r18703 - r18704; float r18706 = exp(r18705); float r18707 = exp(r18704); float r18708 = r18706 - r18707; float r18709 = r18702 * r18708; return r18709; } double f_id(double re, double im) { double r18710 = 0.5; double r18711 = re; double r18712 = cos(r18711); double r18713 = r18710 * r18712; double r18714 = 0.0; double r18715 = im; double r18716 = r18714 - r18715; double r18717 = exp(r18716); double r18718 = exp(r18715); double r18719 = r18717 - r18718; double r18720 = r18713 * r18719; return r18720; } double f_of(float re, float im) { float r18721 = 0.01666666753590107f; float r18722 = im; float r18723 = 5.0f; float r18724 = pow(r18722, r18723); float r18725 = r18721 * r18724; float r18726 = 2.0f; float r18727 = r18726 * r18722; float r18728 = 0.3333333432674408f; float r18729 = 3.0f; float r18730 = pow(r18722, r18729); float r18731 = r18728 * r18730; float r18732 = r18727 + r18731; float r18733 = r18725 + r18732; float r18734 = -r18733; float r18735 = re; float r18736 = cos(r18735); float r18737 = 0.5f; float r18738 = r18736 * r18737; float r18739 = r18734 * r18738; return r18739; } double f_od(double re, double im) { double r18740 = 0.01666666753590107; double r18741 = im; double r18742 = 5.0; double r18743 = pow(r18741, r18742); double r18744 = r18740 * r18743; double r18745 = 2.0; double r18746 = r18745 * r18741; double r18747 = 0.3333333432674408; double r18748 = 3.0; double r18749 = pow(r18741, r18748); double r18750 = r18747 * r18749; double r18751 = r18746 + r18750; double r18752 = r18744 + r18751; double r18753 = -r18752; double r18754 = re; double r18755 = cos(r18754); double r18756 = 0.5; double r18757 = r18755 * r18756; double r18758 = r18753 * r18757; return r18758; } 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 r18759, r18760, r18761, r18762, r18763, r18764, r18765, r18766, r18767, r18768, r18769; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18759, "0.5", 10, MPFR_RNDN); mpfr_init(r18760); mpfr_init(r18761); mpfr_init(r18762); mpfr_init_set_str(r18763, "0", 10, MPFR_RNDN); mpfr_init(r18764); mpfr_init(r18765); mpfr_init(r18766); mpfr_init(r18767); mpfr_init(r18768); mpfr_init(r18769); } double f_im(double re, double im) { ; mpfr_set_d(r18760, re, MPFR_RNDN); mpfr_cos(r18761, r18760, MPFR_RNDN); mpfr_mul(r18762, r18759, r18761, MPFR_RNDN); ; mpfr_set_d(r18764, im, MPFR_RNDN); mpfr_sub(r18765, r18763, r18764, MPFR_RNDN); mpfr_exp(r18766, r18765, MPFR_RNDN); mpfr_exp(r18767, r18764, MPFR_RNDN); mpfr_sub(r18768, r18766, r18767, MPFR_RNDN); mpfr_mul(r18769, r18762, r18768, MPFR_RNDN); return mpfr_get_d(r18769, MPFR_RNDN); } static mpfr_t r18770, r18771, r18772, r18773, r18774, r18775, r18776, r18777, r18778, r18779, r18780, r18781, r18782, r18783, r18784, r18785, r18786, r18787, r18788; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18770, "1/60", 10, MPFR_RNDN); mpfr_init(r18771); mpfr_init_set_str(r18772, "5", 10, MPFR_RNDN); mpfr_init(r18773); mpfr_init(r18774); mpfr_init_set_str(r18775, "2", 10, MPFR_RNDN); mpfr_init(r18776); mpfr_init_set_str(r18777, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18778, "3", 10, MPFR_RNDN); mpfr_init(r18779); mpfr_init(r18780); mpfr_init(r18781); mpfr_init(r18782); mpfr_init(r18783); mpfr_init(r18784); mpfr_init(r18785); mpfr_init_set_str(r18786, "0.5", 10, MPFR_RNDN); mpfr_init(r18787); mpfr_init(r18788); } double f_fm(double re, double im) { ; mpfr_set_d(r18771, im, MPFR_RNDN); ; mpfr_pow(r18773, r18771, r18772, MPFR_RNDN); mpfr_mul(r18774, r18770, r18773, MPFR_RNDN); ; mpfr_mul(r18776, r18775, r18771, MPFR_RNDN); ; ; mpfr_pow(r18779, r18771, r18778, MPFR_RNDN); mpfr_mul(r18780, r18777, r18779, MPFR_RNDN); mpfr_add(r18781, r18776, r18780, MPFR_RNDN); mpfr_add(r18782, r18774, r18781, MPFR_RNDN); mpfr_neg(r18783, r18782, MPFR_RNDN); mpfr_set_d(r18784, re, MPFR_RNDN); mpfr_cos(r18785, r18784, MPFR_RNDN); ; mpfr_mul(r18787, r18785, r18786, MPFR_RNDN); mpfr_mul(r18788, r18783, r18787, MPFR_RNDN); return mpfr_get_d(r18788, MPFR_RNDN); } static mpfr_t r18789, r18790, r18791, r18792, r18793, r18794, r18795, r18796, r18797, r18798, r18799, r18800, r18801, r18802, r18803, r18804, r18805, r18806, r18807; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18789, "1/60", 10, MPFR_RNDN); mpfr_init(r18790); mpfr_init_set_str(r18791, "5", 10, MPFR_RNDN); mpfr_init(r18792); mpfr_init(r18793); mpfr_init_set_str(r18794, "2", 10, MPFR_RNDN); mpfr_init(r18795); mpfr_init_set_str(r18796, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18797, "3", 10, MPFR_RNDN); mpfr_init(r18798); mpfr_init(r18799); mpfr_init(r18800); mpfr_init(r18801); mpfr_init(r18802); mpfr_init(r18803); mpfr_init(r18804); mpfr_init_set_str(r18805, "0.5", 10, MPFR_RNDN); mpfr_init(r18806); mpfr_init(r18807); } double f_dm(double re, double im) { ; mpfr_set_d(r18790, im, MPFR_RNDN); ; mpfr_pow(r18792, r18790, r18791, MPFR_RNDN); mpfr_mul(r18793, r18789, r18792, MPFR_RNDN); ; mpfr_mul(r18795, r18794, r18790, MPFR_RNDN); ; ; mpfr_pow(r18798, r18790, r18797, MPFR_RNDN); mpfr_mul(r18799, r18796, r18798, MPFR_RNDN); mpfr_add(r18800, r18795, r18799, MPFR_RNDN); mpfr_add(r18801, r18793, r18800, MPFR_RNDN); mpfr_neg(r18802, r18801, MPFR_RNDN); mpfr_set_d(r18803, re, MPFR_RNDN); mpfr_cos(r18804, r18803, MPFR_RNDN); ; mpfr_mul(r18806, r18804, r18805, MPFR_RNDN); mpfr_mul(r18807, r18802, r18806, MPFR_RNDN); return mpfr_get_d(r18807, MPFR_RNDN); }