#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19761 = 0.5f; float r19762 = re; float r19763 = sin(r19762); float r19764 = r19761 * r19763; float r19765 = im; float r19766 = -r19765; float r19767 = exp(r19766); float r19768 = exp(r19765); float r19769 = r19767 - r19768; float r19770 = r19764 * r19769; return r19770; } double f_id(double re, double im) { double r19771 = 0.5; double r19772 = re; double r19773 = sin(r19772); double r19774 = r19771 * r19773; double r19775 = im; double r19776 = -r19775; double r19777 = exp(r19776); double r19778 = exp(r19775); double r19779 = r19777 - r19778; double r19780 = r19774 * r19779; return r19780; } double f_of(float re, float im) { float r19781 = 0.5f; float r19782 = re; float r19783 = sin(r19782); float r19784 = r19781 * r19783; float r19785 = 0.016666666666666666f; float r19786 = im; float r19787 = 5.0f; float r19788 = pow(r19786, r19787); float r19789 = r19785 * r19788; float r19790 = 2.0f; float r19791 = r19790 * r19786; float r19792 = 0.3333333333333333f; float r19793 = 3.0f; float r19794 = pow(r19786, r19793); float r19795 = r19792 * r19794; float r19796 = r19791 + r19795; float r19797 = r19789 + r19796; float r19798 = -r19797; float r19799 = r19784 * r19798; return r19799; } double f_od(double re, double im) { double r19800 = 0.5; double r19801 = re; double r19802 = sin(r19801); double r19803 = r19800 * r19802; double r19804 = 0.016666666666666666; double r19805 = im; double r19806 = 5.0; double r19807 = pow(r19805, r19806); double r19808 = r19804 * r19807; double r19809 = 2.0; double r19810 = r19809 * r19805; double r19811 = 0.3333333333333333; double r19812 = 3.0; double r19813 = pow(r19805, r19812); double r19814 = r19811 * r19813; double r19815 = r19810 + r19814; double r19816 = r19808 + r19815; double r19817 = -r19816; double r19818 = r19803 * r19817; return r19818; } 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 r19819, r19820, r19821, r19822, r19823, r19824, r19825, r19826, r19827, r19828; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19819, "0.5", 10, MPFR_RNDN); mpfr_init(r19820); mpfr_init(r19821); mpfr_init(r19822); mpfr_init(r19823); mpfr_init(r19824); mpfr_init(r19825); mpfr_init(r19826); mpfr_init(r19827); mpfr_init(r19828); } double f_im(double re, double im) { ; mpfr_set_d(r19820, re, MPFR_RNDN); mpfr_sin(r19821, r19820, MPFR_RNDN); mpfr_mul(r19822, r19819, r19821, MPFR_RNDN); mpfr_set_d(r19823, im, MPFR_RNDN); mpfr_neg(r19824, r19823, MPFR_RNDN); mpfr_exp(r19825, r19824, MPFR_RNDN); mpfr_exp(r19826, r19823, MPFR_RNDN); mpfr_sub(r19827, r19825, r19826, MPFR_RNDN); mpfr_mul(r19828, r19822, r19827, MPFR_RNDN); return mpfr_get_d(r19828, MPFR_RNDN); } static mpfr_t r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846, r19847; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19829, "0.5", 10, MPFR_RNDN); mpfr_init(r19830); mpfr_init(r19831); mpfr_init(r19832); mpfr_init_set_str(r19833, "1/60", 10, MPFR_RNDN); mpfr_init(r19834); mpfr_init_set_str(r19835, "5", 10, MPFR_RNDN); mpfr_init(r19836); mpfr_init(r19837); mpfr_init_set_str(r19838, "2", 10, MPFR_RNDN); mpfr_init(r19839); mpfr_init_set_str(r19840, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19841, "3", 10, MPFR_RNDN); mpfr_init(r19842); mpfr_init(r19843); mpfr_init(r19844); mpfr_init(r19845); mpfr_init(r19846); mpfr_init(r19847); } double f_fm(double re, double im) { ; mpfr_set_d(r19830, re, MPFR_RNDN); mpfr_sin(r19831, r19830, MPFR_RNDN); mpfr_mul(r19832, r19829, r19831, MPFR_RNDN); ; mpfr_set_d(r19834, im, MPFR_RNDN); ; mpfr_pow(r19836, r19834, r19835, MPFR_RNDN); mpfr_mul(r19837, r19833, r19836, MPFR_RNDN); ; mpfr_mul(r19839, r19838, r19834, MPFR_RNDN); ; ; mpfr_pow(r19842, r19834, r19841, MPFR_RNDN); mpfr_mul(r19843, r19840, r19842, MPFR_RNDN); mpfr_add(r19844, r19839, r19843, MPFR_RNDN); mpfr_add(r19845, r19837, r19844, MPFR_RNDN); mpfr_neg(r19846, r19845, MPFR_RNDN); mpfr_mul(r19847, r19832, r19846, MPFR_RNDN); return mpfr_get_d(r19847, MPFR_RNDN); } static mpfr_t r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19848, "0.5", 10, MPFR_RNDN); mpfr_init(r19849); mpfr_init(r19850); mpfr_init(r19851); mpfr_init_set_str(r19852, "1/60", 10, MPFR_RNDN); mpfr_init(r19853); mpfr_init_set_str(r19854, "5", 10, MPFR_RNDN); mpfr_init(r19855); mpfr_init(r19856); mpfr_init_set_str(r19857, "2", 10, MPFR_RNDN); mpfr_init(r19858); mpfr_init_set_str(r19859, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19860, "3", 10, MPFR_RNDN); mpfr_init(r19861); mpfr_init(r19862); mpfr_init(r19863); mpfr_init(r19864); mpfr_init(r19865); mpfr_init(r19866); } double f_dm(double re, double im) { ; mpfr_set_d(r19849, re, MPFR_RNDN); mpfr_sin(r19850, r19849, MPFR_RNDN); mpfr_mul(r19851, r19848, r19850, MPFR_RNDN); ; mpfr_set_d(r19853, im, MPFR_RNDN); ; mpfr_pow(r19855, r19853, r19854, MPFR_RNDN); mpfr_mul(r19856, r19852, r19855, MPFR_RNDN); ; mpfr_mul(r19858, r19857, r19853, MPFR_RNDN); ; ; mpfr_pow(r19861, r19853, r19860, MPFR_RNDN); mpfr_mul(r19862, r19859, r19861, MPFR_RNDN); mpfr_add(r19863, r19858, r19862, MPFR_RNDN); mpfr_add(r19864, r19856, r19863, MPFR_RNDN); mpfr_neg(r19865, r19864, MPFR_RNDN); mpfr_mul(r19866, r19851, r19865, MPFR_RNDN); return mpfr_get_d(r19866, MPFR_RNDN); }