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