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