#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r23248 = 0.5f; float r23249 = re; float r23250 = sin(r23249); float r23251 = r23248 * r23250; float r23252 = im; float r23253 = -r23252; float r23254 = exp(r23253); float r23255 = exp(r23252); float r23256 = r23254 - r23255; float r23257 = r23251 * r23256; return r23257; } double f_id(double re, double im) { double r23258 = 0.5; double r23259 = re; double r23260 = sin(r23259); double r23261 = r23258 * r23260; double r23262 = im; double r23263 = -r23262; double r23264 = exp(r23263); double r23265 = exp(r23262); double r23266 = r23264 - r23265; double r23267 = r23261 * r23266; return r23267; } double f_of(float re, float im) { float r23268 = im; float r23269 = r23268 + r23268; float r23270 = 5.0f; float r23271 = pow(r23268, r23270); float r23272 = 0.016666666666666666f; float r23273 = r23271 * r23272; float r23274 = r23268 * (r23268 * r23268); float r23275 = 0.3333333333333333f; float r23276 = r23274 * r23275; float r23277 = r23273 + r23276; float r23278 = r23269 + r23277; float r23279 = re; float r23280 = sin(r23279); float r23281 = 0.5f; float r23282 = -r23281; float r23283 = r23280 * r23282; float r23284 = r23278 * r23283; return r23284; } double f_od(double re, double im) { double r23285 = im; double r23286 = r23285 + r23285; double r23287 = 5.0; double r23288 = pow(r23285, r23287); double r23289 = 0.016666666666666666; double r23290 = r23288 * r23289; double r23291 = r23285 * (r23285 * r23285); double r23292 = 0.3333333333333333; double r23293 = r23291 * r23292; double r23294 = r23290 + r23293; double r23295 = r23286 + r23294; double r23296 = re; double r23297 = sin(r23296); double r23298 = 0.5; double r23299 = -r23298; double r23300 = r23297 * r23299; double r23301 = r23295 * r23300; return r23301; } 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 r23302, r23303, r23304, r23305, r23306, r23307, r23308, r23309, r23310, r23311; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r23302, "0.5", 10, MPFR_RNDN); mpfr_init(r23303); mpfr_init(r23304); mpfr_init(r23305); mpfr_init(r23306); mpfr_init(r23307); mpfr_init(r23308); mpfr_init(r23309); mpfr_init(r23310); mpfr_init(r23311); } double f_im(double re, double im) { ; mpfr_set_d(r23303, re, MPFR_RNDN); mpfr_sin(r23304, r23303, MPFR_RNDN); mpfr_mul(r23305, r23302, r23304, MPFR_RNDN); mpfr_set_d(r23306, im, MPFR_RNDN); mpfr_neg(r23307, r23306, MPFR_RNDN); mpfr_exp(r23308, r23307, MPFR_RNDN); mpfr_exp(r23309, r23306, MPFR_RNDN); mpfr_sub(r23310, r23308, r23309, MPFR_RNDN); mpfr_mul(r23311, r23305, r23310, MPFR_RNDN); return mpfr_get_d(r23311, MPFR_RNDN); } static mpfr_t r23312, r23313, r23314, r23315, r23316, r23317, r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327, r23328; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23312); mpfr_init(r23313); mpfr_init_set_str(r23314, "5", 10, MPFR_RNDN); mpfr_init(r23315); mpfr_init_set_str(r23316, "1/60", 10, MPFR_RNDN); mpfr_init(r23317); mpfr_init(r23318); mpfr_init_set_str(r23319, "1/3", 10, MPFR_RNDN); mpfr_init(r23320); mpfr_init(r23321); mpfr_init(r23322); mpfr_init(r23323); mpfr_init(r23324); mpfr_init_set_str(r23325, "0.5", 10, MPFR_RNDN); mpfr_init(r23326); mpfr_init(r23327); mpfr_init(r23328); } double f_fm(double re, double im) { mpfr_set_d(r23312, im, MPFR_RNDN); mpfr_add(r23313, r23312, r23312, MPFR_RNDN); ; mpfr_pow(r23315, r23312, r23314, MPFR_RNDN); ; mpfr_mul(r23317, r23315, r23316, MPFR_RNDN); mpfr_mul(r23318, r23312, r23312, MPFR_RNDN); mpfr_mul(r23318, r23318, r23312, MPFR_RNDN); ; mpfr_mul(r23320, r23318, r23319, MPFR_RNDN); mpfr_add(r23321, r23317, r23320, MPFR_RNDN); mpfr_add(r23322, r23313, r23321, MPFR_RNDN); mpfr_set_d(r23323, re, MPFR_RNDN); mpfr_sin(r23324, r23323, MPFR_RNDN); ; mpfr_neg(r23326, r23325, MPFR_RNDN); mpfr_mul(r23327, r23324, r23326, MPFR_RNDN); mpfr_mul(r23328, r23322, r23327, MPFR_RNDN); return mpfr_get_d(r23328, MPFR_RNDN); } static mpfr_t r23329, r23330, r23331, r23332, r23333, r23334, r23335, r23336, r23337, r23338, r23339, r23340, r23341, r23342, r23343, r23344, r23345; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23329); mpfr_init(r23330); mpfr_init_set_str(r23331, "5", 10, MPFR_RNDN); mpfr_init(r23332); mpfr_init_set_str(r23333, "1/60", 10, MPFR_RNDN); mpfr_init(r23334); mpfr_init(r23335); mpfr_init_set_str(r23336, "1/3", 10, MPFR_RNDN); mpfr_init(r23337); mpfr_init(r23338); mpfr_init(r23339); mpfr_init(r23340); mpfr_init(r23341); mpfr_init_set_str(r23342, "0.5", 10, MPFR_RNDN); mpfr_init(r23343); mpfr_init(r23344); mpfr_init(r23345); } double f_dm(double re, double im) { mpfr_set_d(r23329, im, MPFR_RNDN); mpfr_add(r23330, r23329, r23329, MPFR_RNDN); ; mpfr_pow(r23332, r23329, r23331, MPFR_RNDN); ; mpfr_mul(r23334, r23332, r23333, MPFR_RNDN); mpfr_mul(r23335, r23329, r23329, MPFR_RNDN); mpfr_mul(r23335, r23335, r23329, MPFR_RNDN); ; mpfr_mul(r23337, r23335, r23336, MPFR_RNDN); mpfr_add(r23338, r23334, r23337, MPFR_RNDN); mpfr_add(r23339, r23330, r23338, MPFR_RNDN); mpfr_set_d(r23340, re, MPFR_RNDN); mpfr_sin(r23341, r23340, MPFR_RNDN); ; mpfr_neg(r23343, r23342, MPFR_RNDN); mpfr_mul(r23344, r23341, r23343, MPFR_RNDN); mpfr_mul(r23345, r23339, r23344, MPFR_RNDN); return mpfr_get_d(r23345, MPFR_RNDN); }