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