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