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