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