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