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