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