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