#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19937 = 0.5f; float r19938 = re; float r19939 = sin(r19938); float r19940 = r19937 * r19939; float r19941 = im; float r19942 = -r19941; float r19943 = exp(r19942); float r19944 = exp(r19941); float r19945 = r19943 - r19944; float r19946 = r19940 * r19945; return r19946; } double f_id(double re, double im) { double r19947 = 0.5; double r19948 = re; double r19949 = sin(r19948); double r19950 = r19947 * r19949; double r19951 = im; double r19952 = -r19951; double r19953 = exp(r19952); double r19954 = exp(r19951); double r19955 = r19953 - r19954; double r19956 = r19950 * r19955; return r19956; } double f_of(float re, float im) { float r19957 = 0.5f; float r19958 = re; float r19959 = sin(r19958); float r19960 = r19957 * r19959; float r19961 = 0.016666666666666666f; float r19962 = im; float r19963 = 5.0f; float r19964 = pow(r19962, r19963); float r19965 = r19961 * r19964; float r19966 = 2.0f; float r19967 = r19966 * r19962; float r19968 = 0.3333333333333333f; float r19969 = 3.0f; float r19970 = pow(r19962, r19969); float r19971 = r19968 * r19970; float r19972 = r19967 + r19971; float r19973 = r19965 + r19972; float r19974 = -r19973; float r19975 = r19960 * r19974; return r19975; } double f_od(double re, double im) { double r19976 = 0.5; double r19977 = re; double r19978 = sin(r19977); double r19979 = r19976 * r19978; double r19980 = 0.016666666666666666; double r19981 = im; double r19982 = 5.0; double r19983 = pow(r19981, r19982); double r19984 = r19980 * r19983; double r19985 = 2.0; double r19986 = r19985 * r19981; double r19987 = 0.3333333333333333; double r19988 = 3.0; double r19989 = pow(r19981, r19988); double r19990 = r19987 * r19989; double r19991 = r19986 + r19990; double r19992 = r19984 + r19991; double r19993 = -r19992; double r19994 = r19979 * r19993; return r19994; } 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 r19995, r19996, r19997, r19998, r19999, r20000, r20001, r20002, r20003, r20004; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19995, "0.5", 10, MPFR_RNDN); mpfr_init(r19996); mpfr_init(r19997); mpfr_init(r19998); mpfr_init(r19999); mpfr_init(r20000); mpfr_init(r20001); mpfr_init(r20002); mpfr_init(r20003); mpfr_init(r20004); } double f_im(double re, double im) { ; mpfr_set_d(r19996, re, MPFR_RNDN); mpfr_sin(r19997, r19996, MPFR_RNDN); mpfr_mul(r19998, r19995, r19997, MPFR_RNDN); mpfr_set_d(r19999, im, MPFR_RNDN); mpfr_neg(r20000, r19999, MPFR_RNDN); mpfr_exp(r20001, r20000, MPFR_RNDN); mpfr_exp(r20002, r19999, MPFR_RNDN); mpfr_sub(r20003, r20001, r20002, MPFR_RNDN); mpfr_mul(r20004, r19998, r20003, MPFR_RNDN); return mpfr_get_d(r20004, MPFR_RNDN); } static mpfr_t r20005, r20006, r20007, r20008, r20009, r20010, r20011, r20012, r20013, r20014, r20015, r20016, r20017, r20018, r20019, r20020, r20021, r20022, r20023; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20005, "0.5", 10, MPFR_RNDN); mpfr_init(r20006); mpfr_init(r20007); mpfr_init(r20008); mpfr_init_set_str(r20009, "1/60", 10, MPFR_RNDN); mpfr_init(r20010); mpfr_init_set_str(r20011, "5", 10, MPFR_RNDN); mpfr_init(r20012); mpfr_init(r20013); mpfr_init_set_str(r20014, "2", 10, MPFR_RNDN); mpfr_init(r20015); mpfr_init_set_str(r20016, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20017, "3", 10, MPFR_RNDN); mpfr_init(r20018); mpfr_init(r20019); mpfr_init(r20020); mpfr_init(r20021); mpfr_init(r20022); mpfr_init(r20023); } double f_fm(double re, double im) { ; mpfr_set_d(r20006, re, MPFR_RNDN); mpfr_sin(r20007, r20006, MPFR_RNDN); mpfr_mul(r20008, r20005, r20007, MPFR_RNDN); ; mpfr_set_d(r20010, im, MPFR_RNDN); ; mpfr_pow(r20012, r20010, r20011, MPFR_RNDN); mpfr_mul(r20013, r20009, r20012, MPFR_RNDN); ; mpfr_mul(r20015, r20014, r20010, MPFR_RNDN); ; ; mpfr_pow(r20018, r20010, r20017, MPFR_RNDN); mpfr_mul(r20019, r20016, r20018, MPFR_RNDN); mpfr_add(r20020, r20015, r20019, MPFR_RNDN); mpfr_add(r20021, r20013, r20020, MPFR_RNDN); mpfr_neg(r20022, r20021, MPFR_RNDN); mpfr_mul(r20023, r20008, r20022, MPFR_RNDN); return mpfr_get_d(r20023, MPFR_RNDN); } static mpfr_t r20024, r20025, r20026, r20027, r20028, r20029, r20030, r20031, r20032, r20033, r20034, r20035, r20036, r20037, r20038, r20039, r20040, r20041, r20042; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20024, "0.5", 10, MPFR_RNDN); mpfr_init(r20025); mpfr_init(r20026); mpfr_init(r20027); mpfr_init_set_str(r20028, "1/60", 10, MPFR_RNDN); mpfr_init(r20029); mpfr_init_set_str(r20030, "5", 10, MPFR_RNDN); mpfr_init(r20031); mpfr_init(r20032); mpfr_init_set_str(r20033, "2", 10, MPFR_RNDN); mpfr_init(r20034); mpfr_init_set_str(r20035, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20036, "3", 10, MPFR_RNDN); mpfr_init(r20037); mpfr_init(r20038); mpfr_init(r20039); mpfr_init(r20040); mpfr_init(r20041); mpfr_init(r20042); } double f_dm(double re, double im) { ; mpfr_set_d(r20025, re, MPFR_RNDN); mpfr_sin(r20026, r20025, MPFR_RNDN); mpfr_mul(r20027, r20024, r20026, MPFR_RNDN); ; mpfr_set_d(r20029, im, MPFR_RNDN); ; mpfr_pow(r20031, r20029, r20030, MPFR_RNDN); mpfr_mul(r20032, r20028, r20031, MPFR_RNDN); ; mpfr_mul(r20034, r20033, r20029, MPFR_RNDN); ; ; mpfr_pow(r20037, r20029, r20036, MPFR_RNDN); mpfr_mul(r20038, r20035, r20037, MPFR_RNDN); mpfr_add(r20039, r20034, r20038, MPFR_RNDN); mpfr_add(r20040, r20032, r20039, MPFR_RNDN); mpfr_neg(r20041, r20040, MPFR_RNDN); mpfr_mul(r20042, r20027, r20041, MPFR_RNDN); return mpfr_get_d(r20042, MPFR_RNDN); }