#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19103 = 0.5f; float r19104 = re; float r19105 = cos(r19104); float r19106 = r19103 * r19105; float r19107 = 0.0f; float r19108 = im; float r19109 = r19107 - r19108; float r19110 = exp(r19109); float r19111 = exp(r19108); float r19112 = r19110 - r19111; float r19113 = r19106 * r19112; return r19113; } double f_id(double re, double im) { double r19114 = 0.5; double r19115 = re; double r19116 = cos(r19115); double r19117 = r19114 * r19116; double r19118 = 0.0; double r19119 = im; double r19120 = r19118 - r19119; double r19121 = exp(r19120); double r19122 = exp(r19119); double r19123 = r19121 - r19122; double r19124 = r19117 * r19123; return r19124; } double f_of(float re, float im) { float r19125 = 0.016666666666666666f; float r19126 = im; float r19127 = 5.0f; float r19128 = pow(r19126, r19127); float r19129 = r19126 * (r19126 * r19126); float r19130 = 0.3333333333333333f; float r19131 = 2.0f; float r19132 = r19131 * r19126; float r19133 = fma(r19129, r19130, r19132); float r19134 = fma(r19125, r19128, r19133); float r19135 = 1.0f; float r19136 = r19134 / r19135; float r19137 = re; float r19138 = cos(r19137); float r19139 = 0.5f; float r19140 = -r19139; float r19141 = r19138 * r19140; float r19142 = r19136 * r19141; return r19142; } double f_od(double re, double im) { double r19143 = 0.016666666666666666; double r19144 = im; double r19145 = 5.0; double r19146 = pow(r19144, r19145); double r19147 = r19144 * (r19144 * r19144); double r19148 = 0.3333333333333333; double r19149 = 2.0; double r19150 = r19149 * r19144; double r19151 = fma(r19147, r19148, r19150); double r19152 = fma(r19143, r19146, r19151); double r19153 = 1.0; double r19154 = r19152 / r19153; double r19155 = re; double r19156 = cos(r19155); double r19157 = 0.5; double r19158 = -r19157; double r19159 = r19156 * r19158; double r19160 = r19154 * r19159; return r19160; } 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 r19161, r19162, r19163, r19164, r19165, r19166, r19167, r19168, r19169, r19170, r19171; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19161, "0.5", 10, MPFR_RNDN); mpfr_init(r19162); mpfr_init(r19163); mpfr_init(r19164); mpfr_init_set_str(r19165, "0", 10, MPFR_RNDN); mpfr_init(r19166); mpfr_init(r19167); mpfr_init(r19168); mpfr_init(r19169); mpfr_init(r19170); mpfr_init(r19171); } double f_im(double re, double im) { ; mpfr_set_d(r19162, re, MPFR_RNDN); mpfr_cos(r19163, r19162, MPFR_RNDN); mpfr_mul(r19164, r19161, r19163, MPFR_RNDN); ; mpfr_set_d(r19166, im, MPFR_RNDN); mpfr_sub(r19167, r19165, r19166, MPFR_RNDN); mpfr_exp(r19168, r19167, MPFR_RNDN); mpfr_exp(r19169, r19166, MPFR_RNDN); mpfr_sub(r19170, r19168, r19169, MPFR_RNDN); mpfr_mul(r19171, r19164, r19170, MPFR_RNDN); return mpfr_get_d(r19171, MPFR_RNDN); } static mpfr_t r19172, r19173, r19174, r19175, r19176, r19177, r19178, r19179, r19180, r19181, r19182, r19183, r19184, r19185, r19186, r19187, r19188, r19189; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19172, "1/60", 10, MPFR_RNDN); mpfr_init(r19173); mpfr_init_set_str(r19174, "5", 10, MPFR_RNDN); mpfr_init(r19175); mpfr_init(r19176); mpfr_init_set_str(r19177, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19178, "2", 10, MPFR_RNDN); mpfr_init(r19179); mpfr_init(r19180); mpfr_init(r19181); mpfr_init_set_str(r19182, "1", 10, MPFR_RNDN); mpfr_init(r19183); mpfr_init(r19184); mpfr_init(r19185); mpfr_init_set_str(r19186, "0.5", 10, MPFR_RNDN); mpfr_init(r19187); mpfr_init(r19188); mpfr_init(r19189); } double f_fm(double re, double im) { ; mpfr_set_d(r19173, im, MPFR_RNDN); ; mpfr_pow(r19175, r19173, r19174, MPFR_RNDN); mpfr_mul(r19176, r19173, r19173, MPFR_RNDN); mpfr_mul(r19176, r19176, r19173, MPFR_RNDN); ; ; mpfr_mul(r19179, r19178, r19173, MPFR_RNDN); mpfr_fma(r19180, r19176, r19177, r19179, MPFR_RNDN); mpfr_fma(r19181, r19172, r19175, r19180, MPFR_RNDN); ; mpfr_div(r19183, r19181, r19182, MPFR_RNDN); mpfr_set_d(r19184, re, MPFR_RNDN); mpfr_cos(r19185, r19184, MPFR_RNDN); ; mpfr_neg(r19187, r19186, MPFR_RNDN); mpfr_mul(r19188, r19185, r19187, MPFR_RNDN); mpfr_mul(r19189, r19183, r19188, MPFR_RNDN); return mpfr_get_d(r19189, MPFR_RNDN); } static mpfr_t r19190, r19191, r19192, r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201, r19202, r19203, r19204, r19205, r19206, r19207; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19190, "1/60", 10, MPFR_RNDN); mpfr_init(r19191); mpfr_init_set_str(r19192, "5", 10, MPFR_RNDN); mpfr_init(r19193); mpfr_init(r19194); mpfr_init_set_str(r19195, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19196, "2", 10, MPFR_RNDN); mpfr_init(r19197); mpfr_init(r19198); mpfr_init(r19199); mpfr_init_set_str(r19200, "1", 10, MPFR_RNDN); mpfr_init(r19201); mpfr_init(r19202); mpfr_init(r19203); mpfr_init_set_str(r19204, "0.5", 10, MPFR_RNDN); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); } double f_dm(double re, double im) { ; mpfr_set_d(r19191, im, MPFR_RNDN); ; mpfr_pow(r19193, r19191, r19192, MPFR_RNDN); mpfr_mul(r19194, r19191, r19191, MPFR_RNDN); mpfr_mul(r19194, r19194, r19191, MPFR_RNDN); ; ; mpfr_mul(r19197, r19196, r19191, MPFR_RNDN); mpfr_fma(r19198, r19194, r19195, r19197, MPFR_RNDN); mpfr_fma(r19199, r19190, r19193, r19198, MPFR_RNDN); ; mpfr_div(r19201, r19199, r19200, MPFR_RNDN); mpfr_set_d(r19202, re, MPFR_RNDN); mpfr_cos(r19203, r19202, MPFR_RNDN); ; mpfr_neg(r19205, r19204, MPFR_RNDN); mpfr_mul(r19206, r19203, r19205, MPFR_RNDN); mpfr_mul(r19207, r19201, r19206, MPFR_RNDN); return mpfr_get_d(r19207, MPFR_RNDN); }