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