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