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