#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18931 = 0.5f; float r18932 = re; float r18933 = cos(r18932); float r18934 = r18931 * r18933; float r18935 = 0.0f; float r18936 = im; float r18937 = r18935 - r18936; float r18938 = exp(r18937); float r18939 = exp(r18936); float r18940 = r18938 - r18939; float r18941 = r18934 * r18940; return r18941; } double f_id(double re, double im) { double r18942 = 0.5; double r18943 = re; double r18944 = cos(r18943); double r18945 = r18942 * r18944; double r18946 = 0.0; double r18947 = im; double r18948 = r18946 - r18947; double r18949 = exp(r18948); double r18950 = exp(r18947); double r18951 = r18949 - r18950; double r18952 = r18945 * r18951; return r18952; } double f_of(float re, float im) { float r18953 = 0.016666666666666666f; float r18954 = im; float r18955 = 5.0f; float r18956 = pow(r18954, r18955); float r18957 = r18953 * r18956; float r18958 = 2.0f; float r18959 = r18958 * r18954; float r18960 = 0.3333333333333333f; float r18961 = 3.0f; float r18962 = pow(r18954, r18961); float r18963 = r18960 * r18962; float r18964 = r18959 + r18963; float r18965 = r18957 + r18964; float r18966 = -r18965; float r18967 = re; float r18968 = cos(r18967); float r18969 = 0.5f; float r18970 = r18968 * r18969; float r18971 = r18966 * r18970; return r18971; } double f_od(double re, double im) { double r18972 = 0.016666666666666666; double r18973 = im; double r18974 = 5.0; double r18975 = pow(r18973, r18974); double r18976 = r18972 * r18975; double r18977 = 2.0; double r18978 = r18977 * r18973; double r18979 = 0.3333333333333333; double r18980 = 3.0; double r18981 = pow(r18973, r18980); double r18982 = r18979 * r18981; double r18983 = r18978 + r18982; double r18984 = r18976 + r18983; double r18985 = -r18984; double r18986 = re; double r18987 = cos(r18986); double r18988 = 0.5; double r18989 = r18987 * r18988; double r18990 = r18985 * r18989; return r18990; } 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 r18991, r18992, r18993, r18994, r18995, r18996, r18997, r18998, r18999, r19000, r19001; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18991, "0.5", 10, MPFR_RNDN); mpfr_init(r18992); mpfr_init(r18993); mpfr_init(r18994); mpfr_init_set_str(r18995, "0", 10, MPFR_RNDN); mpfr_init(r18996); mpfr_init(r18997); mpfr_init(r18998); mpfr_init(r18999); mpfr_init(r19000); mpfr_init(r19001); } double f_im(double re, double im) { ; mpfr_set_d(r18992, re, MPFR_RNDN); mpfr_cos(r18993, r18992, MPFR_RNDN); mpfr_mul(r18994, r18991, r18993, MPFR_RNDN); ; mpfr_set_d(r18996, im, MPFR_RNDN); mpfr_sub(r18997, r18995, r18996, MPFR_RNDN); mpfr_exp(r18998, r18997, MPFR_RNDN); mpfr_exp(r18999, r18996, MPFR_RNDN); mpfr_sub(r19000, r18998, r18999, MPFR_RNDN); mpfr_mul(r19001, r18994, r19000, MPFR_RNDN); return mpfr_get_d(r19001, MPFR_RNDN); } static mpfr_t r19002, r19003, r19004, r19005, r19006, r19007, r19008, r19009, r19010, r19011, r19012, r19013, r19014, r19015, r19016, r19017, r19018, r19019, r19020; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19002, "1/60", 10, MPFR_RNDN); mpfr_init(r19003); mpfr_init_set_str(r19004, "5", 10, MPFR_RNDN); mpfr_init(r19005); mpfr_init(r19006); mpfr_init_set_str(r19007, "2", 10, MPFR_RNDN); mpfr_init(r19008); mpfr_init_set_str(r19009, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19010, "3", 10, MPFR_RNDN); mpfr_init(r19011); mpfr_init(r19012); mpfr_init(r19013); mpfr_init(r19014); mpfr_init(r19015); mpfr_init(r19016); mpfr_init(r19017); mpfr_init_set_str(r19018, "0.5", 10, MPFR_RNDN); mpfr_init(r19019); mpfr_init(r19020); } double f_fm(double re, double im) { ; mpfr_set_d(r19003, im, MPFR_RNDN); ; mpfr_pow(r19005, r19003, r19004, MPFR_RNDN); mpfr_mul(r19006, r19002, r19005, MPFR_RNDN); ; mpfr_mul(r19008, r19007, r19003, MPFR_RNDN); ; ; mpfr_pow(r19011, r19003, r19010, MPFR_RNDN); mpfr_mul(r19012, r19009, r19011, MPFR_RNDN); mpfr_add(r19013, r19008, r19012, MPFR_RNDN); mpfr_add(r19014, r19006, r19013, MPFR_RNDN); mpfr_neg(r19015, r19014, MPFR_RNDN); mpfr_set_d(r19016, re, MPFR_RNDN); mpfr_cos(r19017, r19016, MPFR_RNDN); ; mpfr_mul(r19019, r19017, r19018, MPFR_RNDN); mpfr_mul(r19020, r19015, r19019, MPFR_RNDN); return mpfr_get_d(r19020, MPFR_RNDN); } static mpfr_t r19021, r19022, r19023, r19024, r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033, r19034, r19035, r19036, r19037, r19038, r19039; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19021, "1/60", 10, MPFR_RNDN); mpfr_init(r19022); mpfr_init_set_str(r19023, "5", 10, MPFR_RNDN); mpfr_init(r19024); mpfr_init(r19025); mpfr_init_set_str(r19026, "2", 10, MPFR_RNDN); mpfr_init(r19027); mpfr_init_set_str(r19028, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19029, "3", 10, MPFR_RNDN); mpfr_init(r19030); mpfr_init(r19031); mpfr_init(r19032); mpfr_init(r19033); mpfr_init(r19034); mpfr_init(r19035); mpfr_init(r19036); mpfr_init_set_str(r19037, "0.5", 10, MPFR_RNDN); mpfr_init(r19038); mpfr_init(r19039); } double f_dm(double re, double im) { ; mpfr_set_d(r19022, im, MPFR_RNDN); ; mpfr_pow(r19024, r19022, r19023, MPFR_RNDN); mpfr_mul(r19025, r19021, r19024, MPFR_RNDN); ; mpfr_mul(r19027, r19026, r19022, MPFR_RNDN); ; ; mpfr_pow(r19030, r19022, r19029, MPFR_RNDN); mpfr_mul(r19031, r19028, r19030, MPFR_RNDN); mpfr_add(r19032, r19027, r19031, MPFR_RNDN); mpfr_add(r19033, r19025, r19032, MPFR_RNDN); mpfr_neg(r19034, r19033, MPFR_RNDN); mpfr_set_d(r19035, re, MPFR_RNDN); mpfr_cos(r19036, r19035, MPFR_RNDN); ; mpfr_mul(r19038, r19036, r19037, MPFR_RNDN); mpfr_mul(r19039, r19034, r19038, MPFR_RNDN); return mpfr_get_d(r19039, MPFR_RNDN); }