#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r28701 = 0.5; float r28702 = re; float r28703 = sin(r28702); float r28704 = r28701 * r28703; float r28705 = im; float r28706 = -r28705; float r28707 = exp(r28706); float r28708 = exp(r28705); float r28709 = r28707 - r28708; float r28710 = r28704 * r28709; return r28710; } double f_id(double re, double im) { double r28711 = 0.5; double r28712 = re; double r28713 = sin(r28712); double r28714 = r28711 * r28713; double r28715 = im; double r28716 = -r28715; double r28717 = exp(r28716); double r28718 = exp(r28715); double r28719 = r28717 - r28718; double r28720 = r28714 * r28719; return r28720; } double f_of(float re, float im) { float r28721 = 0.5; float r28722 = re; float r28723 = sin(r28722); float r28724 = r28721 * r28723; float r28725 = 1/60; float r28726 = im; float r28727 = 5; float r28728 = pow(r28726, r28727); float r28729 = r28725 * r28728; float r28730 = 2; float r28731 = r28730 * r28726; float r28732 = 1/3; float r28733 = 3; float r28734 = pow(r28726, r28733); float r28735 = r28732 * r28734; float r28736 = r28731 + r28735; float r28737 = r28729 + r28736; float r28738 = -r28737; float r28739 = r28724 * r28738; return r28739; } double f_od(double re, double im) { double r28740 = 0.5; double r28741 = re; double r28742 = sin(r28741); double r28743 = r28740 * r28742; double r28744 = 1/60; double r28745 = im; double r28746 = 5; double r28747 = pow(r28745, r28746); double r28748 = r28744 * r28747; double r28749 = 2; double r28750 = r28749 * r28745; double r28751 = 1/3; double r28752 = 3; double r28753 = pow(r28745, r28752); double r28754 = r28751 * r28753; double r28755 = r28750 + r28754; double r28756 = r28748 + r28755; double r28757 = -r28756; double r28758 = r28743 * r28757; return r28758; } 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 r28759, r28760, r28761, r28762, r28763, r28764, r28765, r28766, r28767, r28768; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r28759, "0.5", 10, MPFR_RNDN); mpfr_init(r28760); mpfr_init(r28761); mpfr_init(r28762); mpfr_init(r28763); mpfr_init(r28764); mpfr_init(r28765); mpfr_init(r28766); mpfr_init(r28767); mpfr_init(r28768); } double f_im(double re, double im) { ; mpfr_set_d(r28760, re, MPFR_RNDN); mpfr_sin(r28761, r28760, MPFR_RNDN); mpfr_mul(r28762, r28759, r28761, MPFR_RNDN); mpfr_set_d(r28763, im, MPFR_RNDN); mpfr_neg(r28764, r28763, MPFR_RNDN); mpfr_exp(r28765, r28764, MPFR_RNDN); mpfr_exp(r28766, r28763, MPFR_RNDN); mpfr_sub(r28767, r28765, r28766, MPFR_RNDN); mpfr_mul(r28768, r28762, r28767, MPFR_RNDN); return mpfr_get_d(r28768, MPFR_RNDN); } static mpfr_t r28769, r28770, r28771, r28772, r28773, r28774, r28775, r28776, r28777, r28778, r28779, r28780, r28781, r28782, r28783, r28784, r28785, r28786, r28787; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r28769, "0.5", 10, MPFR_RNDN); mpfr_init(r28770); mpfr_init(r28771); mpfr_init(r28772); mpfr_init_set_str(r28773, "1/60", 10, MPFR_RNDN); mpfr_init(r28774); mpfr_init_set_str(r28775, "5", 10, MPFR_RNDN); mpfr_init(r28776); mpfr_init(r28777); mpfr_init_set_str(r28778, "2", 10, MPFR_RNDN); mpfr_init(r28779); mpfr_init_set_str(r28780, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r28781, "3", 10, MPFR_RNDN); mpfr_init(r28782); mpfr_init(r28783); mpfr_init(r28784); mpfr_init(r28785); mpfr_init(r28786); mpfr_init(r28787); } double f_fm(double re, double im) { ; mpfr_set_d(r28770, re, MPFR_RNDN); mpfr_sin(r28771, r28770, MPFR_RNDN); mpfr_mul(r28772, r28769, r28771, MPFR_RNDN); ; mpfr_set_d(r28774, im, MPFR_RNDN); ; mpfr_pow(r28776, r28774, r28775, MPFR_RNDN); mpfr_mul(r28777, r28773, r28776, MPFR_RNDN); ; mpfr_mul(r28779, r28778, r28774, MPFR_RNDN); ; ; mpfr_pow(r28782, r28774, r28781, MPFR_RNDN); mpfr_mul(r28783, r28780, r28782, MPFR_RNDN); mpfr_add(r28784, r28779, r28783, MPFR_RNDN); mpfr_add(r28785, r28777, r28784, MPFR_RNDN); mpfr_neg(r28786, r28785, MPFR_RNDN); mpfr_mul(r28787, r28772, r28786, MPFR_RNDN); return mpfr_get_d(r28787, MPFR_RNDN); } static mpfr_t r28788, r28789, r28790, r28791, r28792, r28793, r28794, r28795, r28796, r28797, r28798, r28799, r28800, r28801, r28802, r28803, r28804, r28805, r28806; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r28788, "0.5", 10, MPFR_RNDN); mpfr_init(r28789); mpfr_init(r28790); mpfr_init(r28791); mpfr_init_set_str(r28792, "1/60", 10, MPFR_RNDN); mpfr_init(r28793); mpfr_init_set_str(r28794, "5", 10, MPFR_RNDN); mpfr_init(r28795); mpfr_init(r28796); mpfr_init_set_str(r28797, "2", 10, MPFR_RNDN); mpfr_init(r28798); mpfr_init_set_str(r28799, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r28800, "3", 10, MPFR_RNDN); mpfr_init(r28801); mpfr_init(r28802); mpfr_init(r28803); mpfr_init(r28804); mpfr_init(r28805); mpfr_init(r28806); } double f_dm(double re, double im) { ; mpfr_set_d(r28789, re, MPFR_RNDN); mpfr_sin(r28790, r28789, MPFR_RNDN); mpfr_mul(r28791, r28788, r28790, MPFR_RNDN); ; mpfr_set_d(r28793, im, MPFR_RNDN); ; mpfr_pow(r28795, r28793, r28794, MPFR_RNDN); mpfr_mul(r28796, r28792, r28795, MPFR_RNDN); ; mpfr_mul(r28798, r28797, r28793, MPFR_RNDN); ; ; mpfr_pow(r28801, r28793, r28800, MPFR_RNDN); mpfr_mul(r28802, r28799, r28801, MPFR_RNDN); mpfr_add(r28803, r28798, r28802, MPFR_RNDN); mpfr_add(r28804, r28796, r28803, MPFR_RNDN); mpfr_neg(r28805, r28804, MPFR_RNDN); mpfr_mul(r28806, r28791, r28805, MPFR_RNDN); return mpfr_get_d(r28806, MPFR_RNDN); }