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