#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19635 = 0.5f; float r19636 = re; float r19637 = cos(r19636); float r19638 = r19635 * r19637; float r19639 = 0.0f; float r19640 = im; float r19641 = r19639 - r19640; float r19642 = exp(r19641); float r19643 = exp(r19640); float r19644 = r19642 - r19643; float r19645 = r19638 * r19644; return r19645; } double f_id(double re, double im) { double r19646 = 0.5; double r19647 = re; double r19648 = cos(r19647); double r19649 = r19646 * r19648; double r19650 = 0.0; double r19651 = im; double r19652 = r19650 - r19651; double r19653 = exp(r19652); double r19654 = exp(r19651); double r19655 = r19653 - r19654; double r19656 = r19649 * r19655; return r19656; } double f_of(float re, float im) { float r19657 = 0.016666666666666666f; float r19658 = im; float r19659 = 5.0f; float r19660 = pow(r19658, r19659); float r19661 = r19657 * r19660; float r19662 = 2.0f; float r19663 = r19662 * r19658; float r19664 = 0.3333333333333333f; float r19665 = 3.0f; float r19666 = pow(r19658, r19665); float r19667 = r19664 * r19666; float r19668 = r19663 + r19667; float r19669 = r19661 + r19668; float r19670 = -r19669; float r19671 = re; float r19672 = cos(r19671); float r19673 = 0.5f; float r19674 = r19672 * r19673; float r19675 = r19670 * r19674; return r19675; } double f_od(double re, double im) { double r19676 = 0.016666666666666666; double r19677 = im; double r19678 = 5.0; double r19679 = pow(r19677, r19678); double r19680 = r19676 * r19679; double r19681 = 2.0; double r19682 = r19681 * r19677; double r19683 = 0.3333333333333333; double r19684 = 3.0; double r19685 = pow(r19677, r19684); double r19686 = r19683 * r19685; double r19687 = r19682 + r19686; double r19688 = r19680 + r19687; double r19689 = -r19688; double r19690 = re; double r19691 = cos(r19690); double r19692 = 0.5; double r19693 = r19691 * r19692; double r19694 = r19689 * r19693; return r19694; } 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 r19695, r19696, r19697, r19698, r19699, r19700, r19701, r19702, r19703, r19704, r19705; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19695, "0.5", 10, MPFR_RNDN); mpfr_init(r19696); mpfr_init(r19697); mpfr_init(r19698); mpfr_init_set_str(r19699, "0", 10, MPFR_RNDN); mpfr_init(r19700); mpfr_init(r19701); mpfr_init(r19702); mpfr_init(r19703); mpfr_init(r19704); mpfr_init(r19705); } double f_im(double re, double im) { ; mpfr_set_d(r19696, re, MPFR_RNDN); mpfr_cos(r19697, r19696, MPFR_RNDN); mpfr_mul(r19698, r19695, r19697, MPFR_RNDN); ; mpfr_set_d(r19700, im, MPFR_RNDN); mpfr_sub(r19701, r19699, r19700, MPFR_RNDN); mpfr_exp(r19702, r19701, MPFR_RNDN); mpfr_exp(r19703, r19700, MPFR_RNDN); mpfr_sub(r19704, r19702, r19703, MPFR_RNDN); mpfr_mul(r19705, r19698, r19704, MPFR_RNDN); return mpfr_get_d(r19705, MPFR_RNDN); } static mpfr_t r19706, r19707, r19708, r19709, r19710, r19711, r19712, r19713, r19714, r19715, r19716, r19717, r19718, r19719, r19720, r19721, r19722, r19723, r19724; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19706, "1/60", 10, MPFR_RNDN); mpfr_init(r19707); mpfr_init_set_str(r19708, "5", 10, MPFR_RNDN); mpfr_init(r19709); mpfr_init(r19710); mpfr_init_set_str(r19711, "2", 10, MPFR_RNDN); mpfr_init(r19712); mpfr_init_set_str(r19713, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19714, "3", 10, MPFR_RNDN); mpfr_init(r19715); mpfr_init(r19716); mpfr_init(r19717); mpfr_init(r19718); mpfr_init(r19719); mpfr_init(r19720); mpfr_init(r19721); mpfr_init_set_str(r19722, "0.5", 10, MPFR_RNDN); mpfr_init(r19723); mpfr_init(r19724); } double f_fm(double re, double im) { ; mpfr_set_d(r19707, im, MPFR_RNDN); ; mpfr_pow(r19709, r19707, r19708, MPFR_RNDN); mpfr_mul(r19710, r19706, r19709, MPFR_RNDN); ; mpfr_mul(r19712, r19711, r19707, MPFR_RNDN); ; ; mpfr_pow(r19715, r19707, r19714, MPFR_RNDN); mpfr_mul(r19716, r19713, r19715, MPFR_RNDN); mpfr_add(r19717, r19712, r19716, MPFR_RNDN); mpfr_add(r19718, r19710, r19717, MPFR_RNDN); mpfr_neg(r19719, r19718, MPFR_RNDN); mpfr_set_d(r19720, re, MPFR_RNDN); mpfr_cos(r19721, r19720, MPFR_RNDN); ; mpfr_mul(r19723, r19721, r19722, MPFR_RNDN); mpfr_mul(r19724, r19719, r19723, MPFR_RNDN); return mpfr_get_d(r19724, MPFR_RNDN); } static mpfr_t r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19725, "1/60", 10, MPFR_RNDN); mpfr_init(r19726); mpfr_init_set_str(r19727, "5", 10, MPFR_RNDN); mpfr_init(r19728); mpfr_init(r19729); mpfr_init_set_str(r19730, "2", 10, MPFR_RNDN); mpfr_init(r19731); mpfr_init_set_str(r19732, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19733, "3", 10, MPFR_RNDN); mpfr_init(r19734); mpfr_init(r19735); mpfr_init(r19736); mpfr_init(r19737); mpfr_init(r19738); mpfr_init(r19739); mpfr_init(r19740); mpfr_init_set_str(r19741, "0.5", 10, MPFR_RNDN); mpfr_init(r19742); mpfr_init(r19743); } double f_dm(double re, double im) { ; mpfr_set_d(r19726, im, MPFR_RNDN); ; mpfr_pow(r19728, r19726, r19727, MPFR_RNDN); mpfr_mul(r19729, r19725, r19728, MPFR_RNDN); ; mpfr_mul(r19731, r19730, r19726, MPFR_RNDN); ; ; mpfr_pow(r19734, r19726, r19733, MPFR_RNDN); mpfr_mul(r19735, r19732, r19734, MPFR_RNDN); mpfr_add(r19736, r19731, r19735, MPFR_RNDN); mpfr_add(r19737, r19729, r19736, MPFR_RNDN); mpfr_neg(r19738, r19737, MPFR_RNDN); mpfr_set_d(r19739, re, MPFR_RNDN); mpfr_cos(r19740, r19739, MPFR_RNDN); ; mpfr_mul(r19742, r19740, r19741, MPFR_RNDN); mpfr_mul(r19743, r19738, r19742, MPFR_RNDN); return mpfr_get_d(r19743, MPFR_RNDN); }