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