#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25248 = 0.5; float r25249 = re; float r25250 = cos(r25249); float r25251 = r25248 * r25250; float r25252 = 0; float r25253 = im; float r25254 = r25252 - r25253; float r25255 = exp(r25254); float r25256 = exp(r25253); float r25257 = r25255 - r25256; float r25258 = r25251 * r25257; return r25258; } double f_id(double re, double im) { double r25259 = 0.5; double r25260 = re; double r25261 = cos(r25260); double r25262 = r25259 * r25261; double r25263 = 0; double r25264 = im; double r25265 = r25263 - r25264; double r25266 = exp(r25265); double r25267 = exp(r25264); double r25268 = r25266 - r25267; double r25269 = r25262 * r25268; return r25269; } double f_of(float re, float im) { float r25270 = 0.5; float r25271 = re; float r25272 = cos(r25271); float r25273 = r25270 * r25272; float r25274 = 1/60; float r25275 = im; float r25276 = 5; float r25277 = pow(r25275, r25276); float r25278 = r25274 * r25277; float r25279 = 2; float r25280 = r25279 * r25275; float r25281 = 1/3; float r25282 = 3; float r25283 = pow(r25275, r25282); float r25284 = r25281 * r25283; float r25285 = r25280 + r25284; float r25286 = r25278 + r25285; float r25287 = -r25286; float r25288 = r25273 * r25287; return r25288; } double f_od(double re, double im) { double r25289 = 0.5; double r25290 = re; double r25291 = cos(r25290); double r25292 = r25289 * r25291; double r25293 = 1/60; double r25294 = im; double r25295 = 5; double r25296 = pow(r25294, r25295); double r25297 = r25293 * r25296; double r25298 = 2; double r25299 = r25298 * r25294; double r25300 = 1/3; double r25301 = 3; double r25302 = pow(r25294, r25301); double r25303 = r25300 * r25302; double r25304 = r25299 + r25303; double r25305 = r25297 + r25304; double r25306 = -r25305; double r25307 = r25292 * r25306; return r25307; } 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 r25308, r25309, r25310, r25311, r25312, r25313, r25314, r25315, r25316, r25317, r25318; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25308, "0.5", 10, MPFR_RNDN); mpfr_init(r25309); mpfr_init(r25310); mpfr_init(r25311); mpfr_init_set_str(r25312, "0", 10, MPFR_RNDN); mpfr_init(r25313); mpfr_init(r25314); mpfr_init(r25315); mpfr_init(r25316); mpfr_init(r25317); mpfr_init(r25318); } double f_im(double re, double im) { ; mpfr_set_d(r25309, re, MPFR_RNDN); mpfr_cos(r25310, r25309, MPFR_RNDN); mpfr_mul(r25311, r25308, r25310, MPFR_RNDN); ; mpfr_set_d(r25313, im, MPFR_RNDN); mpfr_sub(r25314, r25312, r25313, MPFR_RNDN); mpfr_exp(r25315, r25314, MPFR_RNDN); mpfr_exp(r25316, r25313, MPFR_RNDN); mpfr_sub(r25317, r25315, r25316, MPFR_RNDN); mpfr_mul(r25318, r25311, r25317, MPFR_RNDN); return mpfr_get_d(r25318, MPFR_RNDN); } static mpfr_t r25319, r25320, r25321, r25322, r25323, r25324, r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25319, "0.5", 10, MPFR_RNDN); mpfr_init(r25320); mpfr_init(r25321); mpfr_init(r25322); mpfr_init_set_str(r25323, "1/60", 10, MPFR_RNDN); mpfr_init(r25324); mpfr_init_set_str(r25325, "5", 10, MPFR_RNDN); mpfr_init(r25326); mpfr_init(r25327); mpfr_init_set_str(r25328, "2", 10, MPFR_RNDN); mpfr_init(r25329); mpfr_init_set_str(r25330, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25331, "3", 10, MPFR_RNDN); mpfr_init(r25332); mpfr_init(r25333); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init(r25337); } double f_fm(double re, double im) { ; mpfr_set_d(r25320, re, MPFR_RNDN); mpfr_cos(r25321, r25320, MPFR_RNDN); mpfr_mul(r25322, r25319, r25321, MPFR_RNDN); ; mpfr_set_d(r25324, im, MPFR_RNDN); ; mpfr_pow(r25326, r25324, r25325, MPFR_RNDN); mpfr_mul(r25327, r25323, r25326, MPFR_RNDN); ; mpfr_mul(r25329, r25328, r25324, MPFR_RNDN); ; ; mpfr_pow(r25332, r25324, r25331, MPFR_RNDN); mpfr_mul(r25333, r25330, r25332, MPFR_RNDN); mpfr_add(r25334, r25329, r25333, MPFR_RNDN); mpfr_add(r25335, r25327, r25334, MPFR_RNDN); mpfr_neg(r25336, r25335, MPFR_RNDN); mpfr_mul(r25337, r25322, r25336, MPFR_RNDN); return mpfr_get_d(r25337, MPFR_RNDN); } static mpfr_t r25338, r25339, r25340, r25341, r25342, r25343, r25344, r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25338, "0.5", 10, MPFR_RNDN); mpfr_init(r25339); mpfr_init(r25340); mpfr_init(r25341); mpfr_init_set_str(r25342, "1/60", 10, MPFR_RNDN); mpfr_init(r25343); mpfr_init_set_str(r25344, "5", 10, MPFR_RNDN); mpfr_init(r25345); mpfr_init(r25346); mpfr_init_set_str(r25347, "2", 10, MPFR_RNDN); mpfr_init(r25348); mpfr_init_set_str(r25349, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25350, "3", 10, MPFR_RNDN); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init(r25356); } double f_dm(double re, double im) { ; mpfr_set_d(r25339, re, MPFR_RNDN); mpfr_cos(r25340, r25339, MPFR_RNDN); mpfr_mul(r25341, r25338, r25340, MPFR_RNDN); ; mpfr_set_d(r25343, im, MPFR_RNDN); ; mpfr_pow(r25345, r25343, r25344, MPFR_RNDN); mpfr_mul(r25346, r25342, r25345, MPFR_RNDN); ; mpfr_mul(r25348, r25347, r25343, MPFR_RNDN); ; ; mpfr_pow(r25351, r25343, r25350, MPFR_RNDN); mpfr_mul(r25352, r25349, r25351, MPFR_RNDN); mpfr_add(r25353, r25348, r25352, MPFR_RNDN); mpfr_add(r25354, r25346, r25353, MPFR_RNDN); mpfr_neg(r25355, r25354, MPFR_RNDN); mpfr_mul(r25356, r25341, r25355, MPFR_RNDN); return mpfr_get_d(r25356, MPFR_RNDN); }