#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r18529 = 0.5f; float r18530 = re; float r18531 = sin(r18530); float r18532 = r18529 * r18531; float r18533 = im; float r18534 = -r18533; float r18535 = exp(r18534); float r18536 = exp(r18533); float r18537 = r18535 - r18536; float r18538 = r18532 * r18537; return r18538; } double f_id(double re, double im) { double r18539 = 0.5; double r18540 = re; double r18541 = sin(r18540); double r18542 = r18539 * r18541; double r18543 = im; double r18544 = -r18543; double r18545 = exp(r18544); double r18546 = exp(r18543); double r18547 = r18545 - r18546; double r18548 = r18542 * r18547; return r18548; } double f_of(float re, float im) { float r18549 = 0.5f; float r18550 = re; float r18551 = sin(r18550); float r18552 = r18549 * r18551; float r18553 = 0.01666666753590107f; float r18554 = im; float r18555 = 5.0f; float r18556 = pow(r18554, r18555); float r18557 = r18553 * r18556; float r18558 = 2.0f; float r18559 = r18558 * r18554; float r18560 = 0.3333333432674408f; float r18561 = 3.0f; float r18562 = pow(r18554, r18561); float r18563 = r18560 * r18562; float r18564 = r18559 + r18563; float r18565 = r18557 + r18564; float r18566 = -r18565; float r18567 = r18552 * r18566; return r18567; } double f_od(double re, double im) { double r18568 = 0.5; double r18569 = re; double r18570 = sin(r18569); double r18571 = r18568 * r18570; double r18572 = 0.01666666753590107; double r18573 = im; double r18574 = 5.0; double r18575 = pow(r18573, r18574); double r18576 = r18572 * r18575; double r18577 = 2.0; double r18578 = r18577 * r18573; double r18579 = 0.3333333432674408; double r18580 = 3.0; double r18581 = pow(r18573, r18580); double r18582 = r18579 * r18581; double r18583 = r18578 + r18582; double r18584 = r18576 + r18583; double r18585 = -r18584; double r18586 = r18571 * r18585; return r18586; } 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 r18587, r18588, r18589, r18590, r18591, r18592, r18593, r18594, r18595, r18596; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18587, "0.5", 10, MPFR_RNDN); mpfr_init(r18588); mpfr_init(r18589); mpfr_init(r18590); mpfr_init(r18591); mpfr_init(r18592); mpfr_init(r18593); mpfr_init(r18594); mpfr_init(r18595); mpfr_init(r18596); } double f_im(double re, double im) { ; mpfr_set_d(r18588, re, MPFR_RNDN); mpfr_sin(r18589, r18588, MPFR_RNDN); mpfr_mul(r18590, r18587, r18589, MPFR_RNDN); mpfr_set_d(r18591, im, MPFR_RNDN); mpfr_neg(r18592, r18591, MPFR_RNDN); mpfr_exp(r18593, r18592, MPFR_RNDN); mpfr_exp(r18594, r18591, MPFR_RNDN); mpfr_sub(r18595, r18593, r18594, MPFR_RNDN); mpfr_mul(r18596, r18590, r18595, MPFR_RNDN); return mpfr_get_d(r18596, MPFR_RNDN); } static mpfr_t r18597, r18598, r18599, r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609, r18610, r18611, r18612, r18613, r18614, r18615; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18597, "0.5", 10, MPFR_RNDN); mpfr_init(r18598); mpfr_init(r18599); mpfr_init(r18600); mpfr_init_set_str(r18601, "1/60", 10, MPFR_RNDN); mpfr_init(r18602); mpfr_init_set_str(r18603, "5", 10, MPFR_RNDN); mpfr_init(r18604); mpfr_init(r18605); mpfr_init_set_str(r18606, "2", 10, MPFR_RNDN); mpfr_init(r18607); mpfr_init_set_str(r18608, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18609, "3", 10, MPFR_RNDN); mpfr_init(r18610); mpfr_init(r18611); mpfr_init(r18612); mpfr_init(r18613); mpfr_init(r18614); mpfr_init(r18615); } double f_fm(double re, double im) { ; mpfr_set_d(r18598, re, MPFR_RNDN); mpfr_sin(r18599, r18598, MPFR_RNDN); mpfr_mul(r18600, r18597, r18599, MPFR_RNDN); ; mpfr_set_d(r18602, im, MPFR_RNDN); ; mpfr_pow(r18604, r18602, r18603, MPFR_RNDN); mpfr_mul(r18605, r18601, r18604, MPFR_RNDN); ; mpfr_mul(r18607, r18606, r18602, MPFR_RNDN); ; ; mpfr_pow(r18610, r18602, r18609, MPFR_RNDN); mpfr_mul(r18611, r18608, r18610, MPFR_RNDN); mpfr_add(r18612, r18607, r18611, MPFR_RNDN); mpfr_add(r18613, r18605, r18612, MPFR_RNDN); mpfr_neg(r18614, r18613, MPFR_RNDN); mpfr_mul(r18615, r18600, r18614, MPFR_RNDN); return mpfr_get_d(r18615, MPFR_RNDN); } static mpfr_t r18616, r18617, r18618, r18619, r18620, r18621, r18622, r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18616, "0.5", 10, MPFR_RNDN); mpfr_init(r18617); mpfr_init(r18618); mpfr_init(r18619); mpfr_init_set_str(r18620, "1/60", 10, MPFR_RNDN); mpfr_init(r18621); mpfr_init_set_str(r18622, "5", 10, MPFR_RNDN); mpfr_init(r18623); mpfr_init(r18624); mpfr_init_set_str(r18625, "2", 10, MPFR_RNDN); mpfr_init(r18626); mpfr_init_set_str(r18627, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18628, "3", 10, MPFR_RNDN); mpfr_init(r18629); mpfr_init(r18630); mpfr_init(r18631); mpfr_init(r18632); mpfr_init(r18633); mpfr_init(r18634); } double f_dm(double re, double im) { ; mpfr_set_d(r18617, re, MPFR_RNDN); mpfr_sin(r18618, r18617, MPFR_RNDN); mpfr_mul(r18619, r18616, r18618, MPFR_RNDN); ; mpfr_set_d(r18621, im, MPFR_RNDN); ; mpfr_pow(r18623, r18621, r18622, MPFR_RNDN); mpfr_mul(r18624, r18620, r18623, MPFR_RNDN); ; mpfr_mul(r18626, r18625, r18621, MPFR_RNDN); ; ; mpfr_pow(r18629, r18621, r18628, MPFR_RNDN); mpfr_mul(r18630, r18627, r18629, MPFR_RNDN); mpfr_add(r18631, r18626, r18630, MPFR_RNDN); mpfr_add(r18632, r18624, r18631, MPFR_RNDN); mpfr_neg(r18633, r18632, MPFR_RNDN); mpfr_mul(r18634, r18619, r18633, MPFR_RNDN); return mpfr_get_d(r18634, MPFR_RNDN); }