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