#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25658 = 0.5; float r25659 = re; float r25660 = cos(r25659); float r25661 = r25658 * r25660; float r25662 = 0; float r25663 = im; float r25664 = r25662 - r25663; float r25665 = exp(r25664); float r25666 = exp(r25663); float r25667 = r25665 - r25666; float r25668 = r25661 * r25667; return r25668; } double f_id(double re, double im) { double r25669 = 0.5; double r25670 = re; double r25671 = cos(r25670); double r25672 = r25669 * r25671; double r25673 = 0; double r25674 = im; double r25675 = r25673 - r25674; double r25676 = exp(r25675); double r25677 = exp(r25674); double r25678 = r25676 - r25677; double r25679 = r25672 * r25678; return r25679; } double f_of(float re, float im) { float r25680 = 0.5; float r25681 = re; float r25682 = cos(r25681); float r25683 = r25680 * r25682; float r25684 = 1/60; float r25685 = im; float r25686 = 5; float r25687 = pow(r25685, r25686); float r25688 = r25684 * r25687; float r25689 = 2; float r25690 = r25689 * r25685; float r25691 = 1/3; float r25692 = 3; float r25693 = pow(r25685, r25692); float r25694 = r25691 * r25693; float r25695 = r25690 + r25694; float r25696 = r25688 + r25695; float r25697 = -r25696; float r25698 = r25683 * r25697; return r25698; } double f_od(double re, double im) { double r25699 = 0.5; double r25700 = re; double r25701 = cos(r25700); double r25702 = r25699 * r25701; double r25703 = 1/60; double r25704 = im; double r25705 = 5; double r25706 = pow(r25704, r25705); double r25707 = r25703 * r25706; double r25708 = 2; double r25709 = r25708 * r25704; double r25710 = 1/3; double r25711 = 3; double r25712 = pow(r25704, r25711); double r25713 = r25710 * r25712; double r25714 = r25709 + r25713; double r25715 = r25707 + r25714; double r25716 = -r25715; double r25717 = r25702 * r25716; return r25717; } 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 r25718, r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25718, "0.5", 10, MPFR_RNDN); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init_set_str(r25722, "0", 10, MPFR_RNDN); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init(r25727); mpfr_init(r25728); } double f_im(double re, double im) { ; mpfr_set_d(r25719, re, MPFR_RNDN); mpfr_cos(r25720, r25719, MPFR_RNDN); mpfr_mul(r25721, r25718, r25720, MPFR_RNDN); ; mpfr_set_d(r25723, im, MPFR_RNDN); mpfr_sub(r25724, r25722, r25723, MPFR_RNDN); mpfr_exp(r25725, r25724, MPFR_RNDN); mpfr_exp(r25726, r25723, MPFR_RNDN); mpfr_sub(r25727, r25725, r25726, MPFR_RNDN); mpfr_mul(r25728, r25721, r25727, MPFR_RNDN); return mpfr_get_d(r25728, MPFR_RNDN); } static mpfr_t r25729, r25730, r25731, r25732, r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25729, "0.5", 10, MPFR_RNDN); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); mpfr_init_set_str(r25733, "1/60", 10, MPFR_RNDN); mpfr_init(r25734); mpfr_init_set_str(r25735, "5", 10, MPFR_RNDN); mpfr_init(r25736); mpfr_init(r25737); mpfr_init_set_str(r25738, "2", 10, MPFR_RNDN); mpfr_init(r25739); mpfr_init_set_str(r25740, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25741, "3", 10, MPFR_RNDN); mpfr_init(r25742); mpfr_init(r25743); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init(r25747); } double f_fm(double re, double im) { ; mpfr_set_d(r25730, re, MPFR_RNDN); mpfr_cos(r25731, r25730, MPFR_RNDN); mpfr_mul(r25732, r25729, r25731, MPFR_RNDN); ; mpfr_set_d(r25734, im, MPFR_RNDN); ; mpfr_pow(r25736, r25734, r25735, MPFR_RNDN); mpfr_mul(r25737, r25733, r25736, MPFR_RNDN); ; mpfr_mul(r25739, r25738, r25734, MPFR_RNDN); ; ; mpfr_pow(r25742, r25734, r25741, MPFR_RNDN); mpfr_mul(r25743, r25740, r25742, MPFR_RNDN); mpfr_add(r25744, r25739, r25743, MPFR_RNDN); mpfr_add(r25745, r25737, r25744, MPFR_RNDN); mpfr_neg(r25746, r25745, MPFR_RNDN); mpfr_mul(r25747, r25732, r25746, MPFR_RNDN); return mpfr_get_d(r25747, MPFR_RNDN); } static mpfr_t r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755, r25756, r25757, r25758, r25759, r25760, r25761, r25762, r25763, r25764, r25765, r25766; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25748, "0.5", 10, MPFR_RNDN); mpfr_init(r25749); mpfr_init(r25750); mpfr_init(r25751); mpfr_init_set_str(r25752, "1/60", 10, MPFR_RNDN); mpfr_init(r25753); mpfr_init_set_str(r25754, "5", 10, MPFR_RNDN); mpfr_init(r25755); mpfr_init(r25756); mpfr_init_set_str(r25757, "2", 10, MPFR_RNDN); mpfr_init(r25758); mpfr_init_set_str(r25759, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25760, "3", 10, MPFR_RNDN); mpfr_init(r25761); mpfr_init(r25762); mpfr_init(r25763); mpfr_init(r25764); mpfr_init(r25765); mpfr_init(r25766); } double f_dm(double re, double im) { ; mpfr_set_d(r25749, re, MPFR_RNDN); mpfr_cos(r25750, r25749, MPFR_RNDN); mpfr_mul(r25751, r25748, r25750, MPFR_RNDN); ; mpfr_set_d(r25753, im, MPFR_RNDN); ; mpfr_pow(r25755, r25753, r25754, MPFR_RNDN); mpfr_mul(r25756, r25752, r25755, MPFR_RNDN); ; mpfr_mul(r25758, r25757, r25753, MPFR_RNDN); ; ; mpfr_pow(r25761, r25753, r25760, MPFR_RNDN); mpfr_mul(r25762, r25759, r25761, MPFR_RNDN); mpfr_add(r25763, r25758, r25762, MPFR_RNDN); mpfr_add(r25764, r25756, r25763, MPFR_RNDN); mpfr_neg(r25765, r25764, MPFR_RNDN); mpfr_mul(r25766, r25751, r25765, MPFR_RNDN); return mpfr_get_d(r25766, MPFR_RNDN); }