#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19985 = 0.5f; float r19986 = re; float r19987 = sin(r19986); float r19988 = r19985 * r19987; float r19989 = im; float r19990 = -r19989; float r19991 = exp(r19990); float r19992 = exp(r19989); float r19993 = r19991 - r19992; float r19994 = r19988 * r19993; return r19994; } double f_id(double re, double im) { double r19995 = 0.5; double r19996 = re; double r19997 = sin(r19996); double r19998 = r19995 * r19997; double r19999 = im; double r20000 = -r19999; double r20001 = exp(r20000); double r20002 = exp(r19999); double r20003 = r20001 - r20002; double r20004 = r19998 * r20003; return r20004; } double f_of(float re, float im) { float r20005 = 0.5f; float r20006 = re; float r20007 = sin(r20006); float r20008 = r20005 * r20007; float r20009 = 0.016666666666666666f; 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.3333333333333333f; 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 = r20008 * r20022; return r20023; } double f_od(double re, double im) { double r20024 = 0.5; double r20025 = re; double r20026 = sin(r20025); double r20027 = r20024 * r20026; double r20028 = 0.016666666666666666; 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.3333333333333333; 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 = r20027 * r20041; return r20042; } 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 r20043, r20044, r20045, r20046, r20047, r20048, r20049, r20050, r20051, r20052; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20043, "0.5", 10, MPFR_RNDN); mpfr_init(r20044); mpfr_init(r20045); mpfr_init(r20046); mpfr_init(r20047); mpfr_init(r20048); mpfr_init(r20049); mpfr_init(r20050); mpfr_init(r20051); mpfr_init(r20052); } double f_im(double re, double im) { ; mpfr_set_d(r20044, re, MPFR_RNDN); mpfr_sin(r20045, r20044, MPFR_RNDN); mpfr_mul(r20046, r20043, r20045, MPFR_RNDN); mpfr_set_d(r20047, im, MPFR_RNDN); mpfr_neg(r20048, r20047, MPFR_RNDN); mpfr_exp(r20049, r20048, MPFR_RNDN); mpfr_exp(r20050, r20047, MPFR_RNDN); mpfr_sub(r20051, r20049, r20050, MPFR_RNDN); mpfr_mul(r20052, r20046, 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_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20053, "0.5", 10, MPFR_RNDN); mpfr_init(r20054); mpfr_init(r20055); mpfr_init(r20056); mpfr_init_set_str(r20057, "1/60", 10, MPFR_RNDN); mpfr_init(r20058); mpfr_init_set_str(r20059, "5", 10, MPFR_RNDN); mpfr_init(r20060); mpfr_init(r20061); mpfr_init_set_str(r20062, "2", 10, MPFR_RNDN); mpfr_init(r20063); mpfr_init_set_str(r20064, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20065, "3", 10, MPFR_RNDN); mpfr_init(r20066); mpfr_init(r20067); mpfr_init(r20068); mpfr_init(r20069); mpfr_init(r20070); mpfr_init(r20071); } double f_fm(double re, double im) { ; mpfr_set_d(r20054, re, MPFR_RNDN); mpfr_sin(r20055, r20054, MPFR_RNDN); mpfr_mul(r20056, r20053, r20055, MPFR_RNDN); ; mpfr_set_d(r20058, im, MPFR_RNDN); ; mpfr_pow(r20060, r20058, r20059, MPFR_RNDN); mpfr_mul(r20061, r20057, r20060, MPFR_RNDN); ; mpfr_mul(r20063, r20062, r20058, MPFR_RNDN); ; ; mpfr_pow(r20066, r20058, r20065, MPFR_RNDN); mpfr_mul(r20067, r20064, r20066, MPFR_RNDN); mpfr_add(r20068, r20063, r20067, MPFR_RNDN); mpfr_add(r20069, r20061, r20068, MPFR_RNDN); mpfr_neg(r20070, r20069, MPFR_RNDN); mpfr_mul(r20071, r20056, r20070, MPFR_RNDN); return mpfr_get_d(r20071, MPFR_RNDN); } static mpfr_t r20072, r20073, r20074, r20075, r20076, r20077, r20078, r20079, r20080, r20081, r20082, r20083, r20084, r20085, r20086, r20087, r20088, r20089, r20090; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20072, "0.5", 10, MPFR_RNDN); mpfr_init(r20073); mpfr_init(r20074); mpfr_init(r20075); mpfr_init_set_str(r20076, "1/60", 10, MPFR_RNDN); mpfr_init(r20077); mpfr_init_set_str(r20078, "5", 10, MPFR_RNDN); mpfr_init(r20079); mpfr_init(r20080); mpfr_init_set_str(r20081, "2", 10, MPFR_RNDN); mpfr_init(r20082); mpfr_init_set_str(r20083, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20084, "3", 10, MPFR_RNDN); mpfr_init(r20085); mpfr_init(r20086); mpfr_init(r20087); mpfr_init(r20088); mpfr_init(r20089); mpfr_init(r20090); } double f_dm(double re, double im) { ; mpfr_set_d(r20073, re, MPFR_RNDN); mpfr_sin(r20074, r20073, MPFR_RNDN); mpfr_mul(r20075, r20072, r20074, MPFR_RNDN); ; mpfr_set_d(r20077, im, MPFR_RNDN); ; mpfr_pow(r20079, r20077, r20078, MPFR_RNDN); mpfr_mul(r20080, r20076, r20079, MPFR_RNDN); ; mpfr_mul(r20082, r20081, r20077, MPFR_RNDN); ; ; mpfr_pow(r20085, r20077, r20084, MPFR_RNDN); mpfr_mul(r20086, r20083, r20085, MPFR_RNDN); mpfr_add(r20087, r20082, r20086, MPFR_RNDN); mpfr_add(r20088, r20080, r20087, MPFR_RNDN); mpfr_neg(r20089, r20088, MPFR_RNDN); mpfr_mul(r20090, r20075, r20089, MPFR_RNDN); return mpfr_get_d(r20090, MPFR_RNDN); }