#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18623 = 0.5f; float r18624 = re; float r18625 = cos(r18624); float r18626 = r18623 * r18625; float r18627 = 0.0f; float r18628 = im; float r18629 = r18627 - r18628; float r18630 = exp(r18629); float r18631 = exp(r18628); float r18632 = r18630 - r18631; float r18633 = r18626 * r18632; return r18633; } double f_id(double re, double im) { double r18634 = 0.5; double r18635 = re; double r18636 = cos(r18635); double r18637 = r18634 * r18636; double r18638 = 0.0; double r18639 = im; double r18640 = r18638 - r18639; double r18641 = exp(r18640); double r18642 = exp(r18639); double r18643 = r18641 - r18642; double r18644 = r18637 * r18643; return r18644; } double f_of(float re, float im) { float r18645 = 0.01666666753590107f; float r18646 = im; float r18647 = 5.0f; float r18648 = pow(r18646, r18647); float r18649 = r18645 * r18648; float r18650 = 2.0f; float r18651 = r18650 * r18646; float r18652 = 0.3333333432674408f; float r18653 = 3.0f; float r18654 = pow(r18646, r18653); float r18655 = r18652 * r18654; float r18656 = r18651 + r18655; float r18657 = r18649 + r18656; float r18658 = -r18657; float r18659 = re; float r18660 = cos(r18659); float r18661 = 0.5f; float r18662 = r18660 * r18661; float r18663 = r18658 * r18662; return r18663; } double f_od(double re, double im) { double r18664 = 0.01666666753590107; double r18665 = im; double r18666 = 5.0; double r18667 = pow(r18665, r18666); double r18668 = r18664 * r18667; double r18669 = 2.0; double r18670 = r18669 * r18665; double r18671 = 0.3333333432674408; double r18672 = 3.0; double r18673 = pow(r18665, r18672); double r18674 = r18671 * r18673; double r18675 = r18670 + r18674; double r18676 = r18668 + r18675; double r18677 = -r18676; double r18678 = re; double r18679 = cos(r18678); double r18680 = 0.5; double r18681 = r18679 * r18680; double r18682 = r18677 * r18681; return r18682; } 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 r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18683, "0.5", 10, MPFR_RNDN); mpfr_init(r18684); mpfr_init(r18685); mpfr_init(r18686); mpfr_init_set_str(r18687, "0", 10, MPFR_RNDN); mpfr_init(r18688); mpfr_init(r18689); mpfr_init(r18690); mpfr_init(r18691); mpfr_init(r18692); mpfr_init(r18693); } double f_im(double re, double im) { ; mpfr_set_d(r18684, re, MPFR_RNDN); mpfr_cos(r18685, r18684, MPFR_RNDN); mpfr_mul(r18686, r18683, r18685, MPFR_RNDN); ; mpfr_set_d(r18688, im, MPFR_RNDN); mpfr_sub(r18689, r18687, r18688, MPFR_RNDN); mpfr_exp(r18690, r18689, MPFR_RNDN); mpfr_exp(r18691, r18688, MPFR_RNDN); mpfr_sub(r18692, r18690, r18691, MPFR_RNDN); mpfr_mul(r18693, r18686, r18692, MPFR_RNDN); return mpfr_get_d(r18693, MPFR_RNDN); } static mpfr_t r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703, r18704, r18705, r18706, r18707, r18708, r18709, r18710, r18711, r18712; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18694, "1/60", 10, MPFR_RNDN); mpfr_init(r18695); mpfr_init_set_str(r18696, "5", 10, MPFR_RNDN); mpfr_init(r18697); mpfr_init(r18698); mpfr_init_set_str(r18699, "2", 10, MPFR_RNDN); mpfr_init(r18700); mpfr_init_set_str(r18701, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18702, "3", 10, MPFR_RNDN); mpfr_init(r18703); mpfr_init(r18704); mpfr_init(r18705); mpfr_init(r18706); mpfr_init(r18707); mpfr_init(r18708); mpfr_init(r18709); mpfr_init_set_str(r18710, "0.5", 10, MPFR_RNDN); mpfr_init(r18711); mpfr_init(r18712); } double f_fm(double re, double im) { ; mpfr_set_d(r18695, im, MPFR_RNDN); ; mpfr_pow(r18697, r18695, r18696, MPFR_RNDN); mpfr_mul(r18698, r18694, r18697, MPFR_RNDN); ; mpfr_mul(r18700, r18699, r18695, MPFR_RNDN); ; ; mpfr_pow(r18703, r18695, r18702, MPFR_RNDN); mpfr_mul(r18704, r18701, r18703, MPFR_RNDN); mpfr_add(r18705, r18700, r18704, MPFR_RNDN); mpfr_add(r18706, r18698, r18705, MPFR_RNDN); mpfr_neg(r18707, r18706, MPFR_RNDN); mpfr_set_d(r18708, re, MPFR_RNDN); mpfr_cos(r18709, r18708, MPFR_RNDN); ; mpfr_mul(r18711, r18709, r18710, MPFR_RNDN); mpfr_mul(r18712, r18707, r18711, MPFR_RNDN); return mpfr_get_d(r18712, MPFR_RNDN); } static mpfr_t r18713, r18714, r18715, r18716, r18717, r18718, r18719, r18720, r18721, r18722, r18723, r18724, r18725, r18726, r18727, r18728, r18729, r18730, r18731; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18713, "1/60", 10, MPFR_RNDN); mpfr_init(r18714); mpfr_init_set_str(r18715, "5", 10, MPFR_RNDN); mpfr_init(r18716); mpfr_init(r18717); mpfr_init_set_str(r18718, "2", 10, MPFR_RNDN); mpfr_init(r18719); mpfr_init_set_str(r18720, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18721, "3", 10, MPFR_RNDN); mpfr_init(r18722); mpfr_init(r18723); mpfr_init(r18724); mpfr_init(r18725); mpfr_init(r18726); mpfr_init(r18727); mpfr_init(r18728); mpfr_init_set_str(r18729, "0.5", 10, MPFR_RNDN); mpfr_init(r18730); mpfr_init(r18731); } double f_dm(double re, double im) { ; mpfr_set_d(r18714, im, MPFR_RNDN); ; mpfr_pow(r18716, r18714, r18715, MPFR_RNDN); mpfr_mul(r18717, r18713, r18716, MPFR_RNDN); ; mpfr_mul(r18719, r18718, r18714, MPFR_RNDN); ; ; mpfr_pow(r18722, r18714, r18721, MPFR_RNDN); mpfr_mul(r18723, r18720, r18722, MPFR_RNDN); mpfr_add(r18724, r18719, r18723, MPFR_RNDN); mpfr_add(r18725, r18717, r18724, MPFR_RNDN); mpfr_neg(r18726, r18725, MPFR_RNDN); mpfr_set_d(r18727, re, MPFR_RNDN); mpfr_cos(r18728, r18727, MPFR_RNDN); ; mpfr_mul(r18730, r18728, r18729, MPFR_RNDN); mpfr_mul(r18731, r18726, r18730, MPFR_RNDN); return mpfr_get_d(r18731, MPFR_RNDN); }