#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19033 = 0.5f; float r19034 = re; float r19035 = sin(r19034); float r19036 = r19033 * r19035; float r19037 = im; float r19038 = -r19037; float r19039 = exp(r19038); float r19040 = exp(r19037); float r19041 = r19039 - r19040; float r19042 = r19036 * r19041; return r19042; } double f_id(double re, double im) { double r19043 = 0.5; double r19044 = re; double r19045 = sin(r19044); double r19046 = r19043 * r19045; double r19047 = im; double r19048 = -r19047; double r19049 = exp(r19048); double r19050 = exp(r19047); double r19051 = r19049 - r19050; double r19052 = r19046 * r19051; return r19052; } double f_of(float re, float im) { float r19053 = 0.5f; float r19054 = re; float r19055 = sin(r19054); float r19056 = r19053 * r19055; float r19057 = 0.016666666666666666f; float r19058 = im; float r19059 = 5.0f; float r19060 = pow(r19058, r19059); float r19061 = r19057 * r19060; float r19062 = 2.0f; float r19063 = r19062 * r19058; float r19064 = 0.3333333333333333f; float r19065 = 3.0f; float r19066 = pow(r19058, r19065); float r19067 = r19064 * r19066; float r19068 = r19063 + r19067; float r19069 = r19061 + r19068; float r19070 = -r19069; float r19071 = r19056 * r19070; return r19071; } double f_od(double re, double im) { double r19072 = 0.5; double r19073 = re; double r19074 = sin(r19073); double r19075 = r19072 * r19074; double r19076 = 0.016666666666666666; double r19077 = im; double r19078 = 5.0; double r19079 = pow(r19077, r19078); double r19080 = r19076 * r19079; double r19081 = 2.0; double r19082 = r19081 * r19077; double r19083 = 0.3333333333333333; double r19084 = 3.0; double r19085 = pow(r19077, r19084); double r19086 = r19083 * r19085; double r19087 = r19082 + r19086; double r19088 = r19080 + r19087; double r19089 = -r19088; double r19090 = r19075 * r19089; return r19090; } 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 r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19091, "0.5", 10, MPFR_RNDN); mpfr_init(r19092); mpfr_init(r19093); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); mpfr_init(r19097); mpfr_init(r19098); mpfr_init(r19099); mpfr_init(r19100); } double f_im(double re, double im) { ; mpfr_set_d(r19092, re, MPFR_RNDN); mpfr_sin(r19093, r19092, MPFR_RNDN); mpfr_mul(r19094, r19091, r19093, MPFR_RNDN); mpfr_set_d(r19095, im, MPFR_RNDN); mpfr_neg(r19096, r19095, MPFR_RNDN); mpfr_exp(r19097, r19096, MPFR_RNDN); mpfr_exp(r19098, r19095, MPFR_RNDN); mpfr_sub(r19099, r19097, r19098, MPFR_RNDN); mpfr_mul(r19100, r19094, r19099, MPFR_RNDN); return mpfr_get_d(r19100, MPFR_RNDN); } static mpfr_t r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117, r19118, r19119; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19101, "0.5", 10, MPFR_RNDN); mpfr_init(r19102); mpfr_init(r19103); mpfr_init(r19104); mpfr_init_set_str(r19105, "1/60", 10, MPFR_RNDN); mpfr_init(r19106); mpfr_init_set_str(r19107, "5", 10, MPFR_RNDN); mpfr_init(r19108); mpfr_init(r19109); mpfr_init_set_str(r19110, "2", 10, MPFR_RNDN); mpfr_init(r19111); mpfr_init_set_str(r19112, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19113, "3", 10, MPFR_RNDN); mpfr_init(r19114); mpfr_init(r19115); mpfr_init(r19116); mpfr_init(r19117); mpfr_init(r19118); mpfr_init(r19119); } double f_fm(double re, double im) { ; mpfr_set_d(r19102, re, MPFR_RNDN); mpfr_sin(r19103, r19102, MPFR_RNDN); mpfr_mul(r19104, r19101, r19103, MPFR_RNDN); ; mpfr_set_d(r19106, im, MPFR_RNDN); ; mpfr_pow(r19108, r19106, r19107, MPFR_RNDN); mpfr_mul(r19109, r19105, r19108, MPFR_RNDN); ; mpfr_mul(r19111, r19110, r19106, MPFR_RNDN); ; ; mpfr_pow(r19114, r19106, r19113, MPFR_RNDN); mpfr_mul(r19115, r19112, r19114, MPFR_RNDN); mpfr_add(r19116, r19111, r19115, MPFR_RNDN); mpfr_add(r19117, r19109, r19116, MPFR_RNDN); mpfr_neg(r19118, r19117, MPFR_RNDN); mpfr_mul(r19119, r19104, r19118, MPFR_RNDN); return mpfr_get_d(r19119, MPFR_RNDN); } static mpfr_t r19120, r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19120, "0.5", 10, MPFR_RNDN); mpfr_init(r19121); mpfr_init(r19122); mpfr_init(r19123); mpfr_init_set_str(r19124, "1/60", 10, MPFR_RNDN); mpfr_init(r19125); mpfr_init_set_str(r19126, "5", 10, MPFR_RNDN); mpfr_init(r19127); mpfr_init(r19128); mpfr_init_set_str(r19129, "2", 10, MPFR_RNDN); mpfr_init(r19130); mpfr_init_set_str(r19131, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19132, "3", 10, MPFR_RNDN); mpfr_init(r19133); mpfr_init(r19134); mpfr_init(r19135); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); } double f_dm(double re, double im) { ; mpfr_set_d(r19121, re, MPFR_RNDN); mpfr_sin(r19122, r19121, MPFR_RNDN); mpfr_mul(r19123, r19120, r19122, MPFR_RNDN); ; mpfr_set_d(r19125, im, MPFR_RNDN); ; mpfr_pow(r19127, r19125, r19126, MPFR_RNDN); mpfr_mul(r19128, r19124, r19127, MPFR_RNDN); ; mpfr_mul(r19130, r19129, r19125, MPFR_RNDN); ; ; mpfr_pow(r19133, r19125, r19132, MPFR_RNDN); mpfr_mul(r19134, r19131, r19133, MPFR_RNDN); mpfr_add(r19135, r19130, r19134, MPFR_RNDN); mpfr_add(r19136, r19128, r19135, MPFR_RNDN); mpfr_neg(r19137, r19136, MPFR_RNDN); mpfr_mul(r19138, r19123, r19137, MPFR_RNDN); return mpfr_get_d(r19138, MPFR_RNDN); }