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