#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.5f; float r20014 = re; float r20015 = cos(r20014); float r20016 = r20013 * r20015; float r20017 = 0.016666666666666666f; float r20018 = im; float r20019 = 5.0f; float r20020 = pow(r20018, r20019); float r20021 = r20017 * r20020; float r20022 = 0.3333333333333333f; float r20023 = r20018 * (r20018 * r20018); float r20024 = r20022 * r20023; float r20025 = r20021 + r20024; float r20026 = 2.0f; float r20027 = r20026 * r20018; float r20028 = r20025 + r20027; float r20029 = r20016 * r20028; float r20030 = -r20029; return r20030; } double f_od(double re, double im) { double r20031 = 0.5; double r20032 = re; double r20033 = cos(r20032); double r20034 = r20031 * r20033; double r20035 = 0.016666666666666666; double r20036 = im; double r20037 = 5.0; double r20038 = pow(r20036, r20037); double r20039 = r20035 * r20038; double r20040 = 0.3333333333333333; double r20041 = r20036 * (r20036 * r20036); double r20042 = r20040 * r20041; double r20043 = r20039 + r20042; double r20044 = 2.0; double r20045 = r20044 * r20036; double r20046 = r20043 + r20045; double r20047 = r20034 * r20046; double r20048 = -r20047; return r20048; } 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 r20049, r20050, r20051, r20052, r20053, r20054, r20055, r20056, r20057, r20058, r20059; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20049, "0.5", 10, MPFR_RNDN); mpfr_init(r20050); mpfr_init(r20051); mpfr_init(r20052); mpfr_init_set_str(r20053, "0", 10, MPFR_RNDN); mpfr_init(r20054); mpfr_init(r20055); mpfr_init(r20056); mpfr_init(r20057); mpfr_init(r20058); mpfr_init(r20059); } double f_im(double re, double im) { ; mpfr_set_d(r20050, re, MPFR_RNDN); mpfr_cos(r20051, r20050, MPFR_RNDN); mpfr_mul(r20052, r20049, r20051, MPFR_RNDN); ; mpfr_set_d(r20054, im, MPFR_RNDN); mpfr_sub(r20055, r20053, r20054, MPFR_RNDN); mpfr_exp(r20056, r20055, MPFR_RNDN); mpfr_exp(r20057, r20054, MPFR_RNDN); mpfr_sub(r20058, r20056, r20057, MPFR_RNDN); mpfr_mul(r20059, r20052, r20058, MPFR_RNDN); return mpfr_get_d(r20059, MPFR_RNDN); } static mpfr_t r20060, r20061, r20062, r20063, r20064, r20065, r20066, r20067, r20068, r20069, r20070, r20071, r20072, r20073, r20074, r20075, r20076, r20077; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20060, "0.5", 10, MPFR_RNDN); mpfr_init(r20061); mpfr_init(r20062); mpfr_init(r20063); mpfr_init_set_str(r20064, "1/60", 10, MPFR_RNDN); mpfr_init(r20065); mpfr_init_set_str(r20066, "5", 10, MPFR_RNDN); mpfr_init(r20067); mpfr_init(r20068); mpfr_init_set_str(r20069, "1/3", 10, MPFR_RNDN); mpfr_init(r20070); mpfr_init(r20071); mpfr_init(r20072); mpfr_init_set_str(r20073, "2", 10, MPFR_RNDN); mpfr_init(r20074); mpfr_init(r20075); mpfr_init(r20076); mpfr_init(r20077); } double f_fm(double re, double im) { ; mpfr_set_d(r20061, re, MPFR_RNDN); mpfr_cos(r20062, r20061, MPFR_RNDN); mpfr_mul(r20063, r20060, r20062, MPFR_RNDN); ; mpfr_set_d(r20065, im, MPFR_RNDN); ; mpfr_pow(r20067, r20065, r20066, MPFR_RNDN); mpfr_mul(r20068, r20064, r20067, MPFR_RNDN); ; mpfr_mul(r20070, r20065, r20065, MPFR_RNDN); mpfr_mul(r20070, r20070, r20065, MPFR_RNDN); mpfr_mul(r20071, r20069, r20070, MPFR_RNDN); mpfr_add(r20072, r20068, r20071, MPFR_RNDN); ; mpfr_mul(r20074, r20073, r20065, MPFR_RNDN); mpfr_add(r20075, r20072, r20074, MPFR_RNDN); mpfr_mul(r20076, r20063, r20075, MPFR_RNDN); mpfr_neg(r20077, r20076, MPFR_RNDN); return mpfr_get_d(r20077, MPFR_RNDN); } static mpfr_t 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(r20078, "0.5", 10, MPFR_RNDN); mpfr_init(r20079); mpfr_init(r20080); mpfr_init(r20081); mpfr_init_set_str(r20082, "1/60", 10, MPFR_RNDN); mpfr_init(r20083); mpfr_init_set_str(r20084, "5", 10, MPFR_RNDN); mpfr_init(r20085); mpfr_init(r20086); mpfr_init_set_str(r20087, "1/3", 10, MPFR_RNDN); mpfr_init(r20088); mpfr_init(r20089); mpfr_init(r20090); mpfr_init_set_str(r20091, "2", 10, MPFR_RNDN); mpfr_init(r20092); mpfr_init(r20093); mpfr_init(r20094); mpfr_init(r20095); } double f_dm(double re, double im) { ; mpfr_set_d(r20079, re, MPFR_RNDN); mpfr_cos(r20080, r20079, MPFR_RNDN); mpfr_mul(r20081, r20078, r20080, MPFR_RNDN); ; mpfr_set_d(r20083, im, MPFR_RNDN); ; mpfr_pow(r20085, r20083, r20084, MPFR_RNDN); mpfr_mul(r20086, r20082, r20085, MPFR_RNDN); ; mpfr_mul(r20088, r20083, r20083, MPFR_RNDN); mpfr_mul(r20088, r20088, r20083, MPFR_RNDN); mpfr_mul(r20089, r20087, r20088, MPFR_RNDN); mpfr_add(r20090, r20086, r20089, MPFR_RNDN); ; mpfr_mul(r20092, r20091, r20083, MPFR_RNDN); mpfr_add(r20093, r20090, r20092, MPFR_RNDN); mpfr_mul(r20094, r20081, r20093, MPFR_RNDN); mpfr_neg(r20095, r20094, MPFR_RNDN); return mpfr_get_d(r20095, MPFR_RNDN); }