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