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