#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r18565 = 0.5f; float r18566 = re; float r18567 = sin(r18566); float r18568 = r18565 * r18567; float r18569 = im; float r18570 = -r18569; float r18571 = exp(r18570); float r18572 = exp(r18569); float r18573 = r18571 - r18572; float r18574 = r18568 * r18573; return r18574; } double f_id(double re, double im) { double r18575 = 0.5; double r18576 = re; double r18577 = sin(r18576); double r18578 = r18575 * r18577; double r18579 = im; double r18580 = -r18579; double r18581 = exp(r18580); double r18582 = exp(r18579); double r18583 = r18581 - r18582; double r18584 = r18578 * r18583; return r18584; } double f_of(float re, float im) { float r18585 = 0.5f; float r18586 = re; float r18587 = sin(r18586); float r18588 = r18585 * r18587; float r18589 = 0.01666666753590107f; float r18590 = im; float r18591 = 5.0f; float r18592 = pow(r18590, r18591); float r18593 = r18589 * r18592; float r18594 = 2.0f; float r18595 = r18594 * r18590; float r18596 = 0.3333333432674408f; float r18597 = 3.0f; float r18598 = pow(r18590, r18597); float r18599 = r18596 * r18598; float r18600 = r18595 + r18599; float r18601 = r18593 + r18600; float r18602 = -r18601; float r18603 = r18588 * r18602; return r18603; } double f_od(double re, double im) { double r18604 = 0.5; double r18605 = re; double r18606 = sin(r18605); double r18607 = r18604 * r18606; double r18608 = 0.01666666753590107; double r18609 = im; double r18610 = 5.0; double r18611 = pow(r18609, r18610); double r18612 = r18608 * r18611; double r18613 = 2.0; double r18614 = r18613 * r18609; double r18615 = 0.3333333432674408; double r18616 = 3.0; double r18617 = pow(r18609, r18616); double r18618 = r18615 * r18617; double r18619 = r18614 + r18618; double r18620 = r18612 + r18619; double r18621 = -r18620; double r18622 = r18607 * r18621; return r18622; } 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 r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18623, "0.5", 10, MPFR_RNDN); mpfr_init(r18624); mpfr_init(r18625); mpfr_init(r18626); mpfr_init(r18627); mpfr_init(r18628); mpfr_init(r18629); mpfr_init(r18630); mpfr_init(r18631); mpfr_init(r18632); } double f_im(double re, double im) { ; mpfr_set_d(r18624, re, MPFR_RNDN); mpfr_sin(r18625, r18624, MPFR_RNDN); mpfr_mul(r18626, r18623, r18625, MPFR_RNDN); mpfr_set_d(r18627, im, MPFR_RNDN); mpfr_neg(r18628, r18627, MPFR_RNDN); mpfr_exp(r18629, r18628, MPFR_RNDN); mpfr_exp(r18630, r18627, MPFR_RNDN); mpfr_sub(r18631, r18629, r18630, MPFR_RNDN); mpfr_mul(r18632, r18626, r18631, MPFR_RNDN); return mpfr_get_d(r18632, MPFR_RNDN); } static mpfr_t r18633, r18634, r18635, r18636, r18637, r18638, r18639, r18640, r18641, r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18633, "0.5", 10, MPFR_RNDN); mpfr_init(r18634); mpfr_init(r18635); mpfr_init(r18636); mpfr_init_set_str(r18637, "1/60", 10, MPFR_RNDN); mpfr_init(r18638); mpfr_init_set_str(r18639, "5", 10, MPFR_RNDN); mpfr_init(r18640); mpfr_init(r18641); mpfr_init_set_str(r18642, "2", 10, MPFR_RNDN); mpfr_init(r18643); mpfr_init_set_str(r18644, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18645, "3", 10, MPFR_RNDN); mpfr_init(r18646); mpfr_init(r18647); mpfr_init(r18648); mpfr_init(r18649); mpfr_init(r18650); mpfr_init(r18651); } double f_fm(double re, double im) { ; mpfr_set_d(r18634, re, MPFR_RNDN); mpfr_sin(r18635, r18634, MPFR_RNDN); mpfr_mul(r18636, r18633, r18635, MPFR_RNDN); ; mpfr_set_d(r18638, im, MPFR_RNDN); ; mpfr_pow(r18640, r18638, r18639, MPFR_RNDN); mpfr_mul(r18641, r18637, r18640, MPFR_RNDN); ; mpfr_mul(r18643, r18642, r18638, MPFR_RNDN); ; ; mpfr_pow(r18646, r18638, r18645, MPFR_RNDN); mpfr_mul(r18647, r18644, r18646, MPFR_RNDN); mpfr_add(r18648, r18643, r18647, MPFR_RNDN); mpfr_add(r18649, r18641, r18648, MPFR_RNDN); mpfr_neg(r18650, r18649, MPFR_RNDN); mpfr_mul(r18651, r18636, r18650, MPFR_RNDN); return mpfr_get_d(r18651, MPFR_RNDN); } static mpfr_t r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667, r18668, r18669, r18670; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18652, "0.5", 10, MPFR_RNDN); mpfr_init(r18653); mpfr_init(r18654); mpfr_init(r18655); mpfr_init_set_str(r18656, "1/60", 10, MPFR_RNDN); mpfr_init(r18657); mpfr_init_set_str(r18658, "5", 10, MPFR_RNDN); mpfr_init(r18659); mpfr_init(r18660); mpfr_init_set_str(r18661, "2", 10, MPFR_RNDN); mpfr_init(r18662); mpfr_init_set_str(r18663, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18664, "3", 10, MPFR_RNDN); mpfr_init(r18665); mpfr_init(r18666); mpfr_init(r18667); mpfr_init(r18668); mpfr_init(r18669); mpfr_init(r18670); } double f_dm(double re, double im) { ; mpfr_set_d(r18653, re, MPFR_RNDN); mpfr_sin(r18654, r18653, MPFR_RNDN); mpfr_mul(r18655, r18652, r18654, MPFR_RNDN); ; mpfr_set_d(r18657, im, MPFR_RNDN); ; mpfr_pow(r18659, r18657, r18658, MPFR_RNDN); mpfr_mul(r18660, r18656, r18659, MPFR_RNDN); ; mpfr_mul(r18662, r18661, r18657, MPFR_RNDN); ; ; mpfr_pow(r18665, r18657, r18664, MPFR_RNDN); mpfr_mul(r18666, r18663, r18665, MPFR_RNDN); mpfr_add(r18667, r18662, r18666, MPFR_RNDN); mpfr_add(r18668, r18660, r18667, MPFR_RNDN); mpfr_neg(r18669, r18668, MPFR_RNDN); mpfr_mul(r18670, r18655, r18669, MPFR_RNDN); return mpfr_get_d(r18670, MPFR_RNDN); }