#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r26619 = 0.5; float r26620 = re; float r26621 = cos(r26620); float r26622 = r26619 * r26621; float r26623 = 0; float r26624 = im; float r26625 = r26623 - r26624; float r26626 = exp(r26625); float r26627 = exp(r26624); float r26628 = r26626 - r26627; float r26629 = r26622 * r26628; return r26629; } double f_id(double re, double im) { double r26630 = 0.5; double r26631 = re; double r26632 = cos(r26631); double r26633 = r26630 * r26632; double r26634 = 0; double r26635 = im; double r26636 = r26634 - r26635; double r26637 = exp(r26636); double r26638 = exp(r26635); double r26639 = r26637 - r26638; double r26640 = r26633 * r26639; return r26640; } double f_of(float re, float im) { float r26641 = 0.5; float r26642 = re; float r26643 = cos(r26642); float r26644 = r26641 * r26643; float r26645 = 1/60; float r26646 = im; float r26647 = 5; float r26648 = pow(r26646, r26647); float r26649 = r26645 * r26648; float r26650 = 2; float r26651 = r26650 * r26646; float r26652 = 1/3; float r26653 = 3; float r26654 = pow(r26646, r26653); float r26655 = r26652 * r26654; float r26656 = r26651 + r26655; float r26657 = r26649 + r26656; float r26658 = -r26657; float r26659 = r26644 * r26658; return r26659; } double f_od(double re, double im) { double r26660 = 0.5; double r26661 = re; double r26662 = cos(r26661); double r26663 = r26660 * r26662; double r26664 = 1/60; double r26665 = im; double r26666 = 5; double r26667 = pow(r26665, r26666); double r26668 = r26664 * r26667; double r26669 = 2; double r26670 = r26669 * r26665; double r26671 = 1/3; double r26672 = 3; double r26673 = pow(r26665, r26672); double r26674 = r26671 * r26673; double r26675 = r26670 + r26674; double r26676 = r26668 + r26675; double r26677 = -r26676; double r26678 = r26663 * r26677; return r26678; } 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 r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26679, "0.5", 10, MPFR_RNDN); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); mpfr_init_set_str(r26683, "0", 10, MPFR_RNDN); mpfr_init(r26684); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); } double f_im(double re, double im) { ; mpfr_set_d(r26680, re, MPFR_RNDN); mpfr_cos(r26681, r26680, MPFR_RNDN); mpfr_mul(r26682, r26679, r26681, MPFR_RNDN); ; mpfr_set_d(r26684, im, MPFR_RNDN); mpfr_sub(r26685, r26683, r26684, MPFR_RNDN); mpfr_exp(r26686, r26685, MPFR_RNDN); mpfr_exp(r26687, r26684, MPFR_RNDN); mpfr_sub(r26688, r26686, r26687, MPFR_RNDN); mpfr_mul(r26689, r26682, r26688, MPFR_RNDN); return mpfr_get_d(r26689, MPFR_RNDN); } static mpfr_t r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26690, "0.5", 10, MPFR_RNDN); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init_set_str(r26694, "1/60", 10, MPFR_RNDN); mpfr_init(r26695); mpfr_init_set_str(r26696, "5", 10, MPFR_RNDN); mpfr_init(r26697); mpfr_init(r26698); mpfr_init_set_str(r26699, "2", 10, MPFR_RNDN); mpfr_init(r26700); mpfr_init_set_str(r26701, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r26702, "3", 10, MPFR_RNDN); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); } double f_fm(double re, double im) { ; mpfr_set_d(r26691, re, MPFR_RNDN); mpfr_cos(r26692, r26691, MPFR_RNDN); mpfr_mul(r26693, r26690, r26692, MPFR_RNDN); ; mpfr_set_d(r26695, im, MPFR_RNDN); ; mpfr_pow(r26697, r26695, r26696, MPFR_RNDN); mpfr_mul(r26698, r26694, r26697, MPFR_RNDN); ; mpfr_mul(r26700, r26699, r26695, MPFR_RNDN); ; ; mpfr_pow(r26703, r26695, r26702, MPFR_RNDN); mpfr_mul(r26704, r26701, r26703, MPFR_RNDN); mpfr_add(r26705, r26700, r26704, MPFR_RNDN); mpfr_add(r26706, r26698, r26705, MPFR_RNDN); mpfr_neg(r26707, r26706, MPFR_RNDN); mpfr_mul(r26708, r26693, r26707, MPFR_RNDN); return mpfr_get_d(r26708, MPFR_RNDN); } static mpfr_t r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725, r26726, r26727; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r26709, "0.5", 10, MPFR_RNDN); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); mpfr_init_set_str(r26713, "1/60", 10, MPFR_RNDN); mpfr_init(r26714); mpfr_init_set_str(r26715, "5", 10, MPFR_RNDN); mpfr_init(r26716); mpfr_init(r26717); mpfr_init_set_str(r26718, "2", 10, MPFR_RNDN); mpfr_init(r26719); mpfr_init_set_str(r26720, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r26721, "3", 10, MPFR_RNDN); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); mpfr_init(r26725); mpfr_init(r26726); mpfr_init(r26727); } double f_dm(double re, double im) { ; mpfr_set_d(r26710, re, MPFR_RNDN); mpfr_cos(r26711, r26710, MPFR_RNDN); mpfr_mul(r26712, r26709, r26711, MPFR_RNDN); ; mpfr_set_d(r26714, im, MPFR_RNDN); ; mpfr_pow(r26716, r26714, r26715, MPFR_RNDN); mpfr_mul(r26717, r26713, r26716, MPFR_RNDN); ; mpfr_mul(r26719, r26718, r26714, MPFR_RNDN); ; ; mpfr_pow(r26722, r26714, r26721, MPFR_RNDN); mpfr_mul(r26723, r26720, r26722, MPFR_RNDN); mpfr_add(r26724, r26719, r26723, MPFR_RNDN); mpfr_add(r26725, r26717, r26724, MPFR_RNDN); mpfr_neg(r26726, r26725, MPFR_RNDN); mpfr_mul(r26727, r26712, r26726, MPFR_RNDN); return mpfr_get_d(r26727, MPFR_RNDN); }