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