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