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