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