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