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