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