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