#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.01666666753590107f; float r19126 = im; float r19127 = 5.0f; float r19128 = pow(r19126, r19127); float r19129 = r19125 * r19128; float r19130 = 2.0f; float r19131 = r19130 * r19126; float r19132 = 0.3333333432674408f; float r19133 = 3.0f; float r19134 = pow(r19126, r19133); float r19135 = r19132 * r19134; float r19136 = r19131 + r19135; float r19137 = r19129 + r19136; float r19138 = -r19137; float r19139 = re; float r19140 = cos(r19139); float r19141 = 0.5f; float r19142 = r19140 * r19141; float r19143 = r19138 * r19142; return r19143; } double f_od(double re, double im) { double r19144 = 0.01666666753590107; double r19145 = im; double r19146 = 5.0; double r19147 = pow(r19145, r19146); double r19148 = r19144 * r19147; double r19149 = 2.0; double r19150 = r19149 * r19145; double r19151 = 0.3333333432674408; double r19152 = 3.0; double r19153 = pow(r19145, r19152); double r19154 = r19151 * r19153; double r19155 = r19150 + r19154; double r19156 = r19148 + r19155; double r19157 = -r19156; double r19158 = re; double r19159 = cos(r19158); double r19160 = 0.5; double r19161 = r19159 * r19160; double r19162 = r19157 * r19161; return r19162; } 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 r19163, r19164, r19165, r19166, r19167, r19168, r19169, r19170, r19171, r19172, r19173; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19163, "0.5", 10, MPFR_RNDN); mpfr_init(r19164); mpfr_init(r19165); mpfr_init(r19166); mpfr_init_set_str(r19167, "0", 10, MPFR_RNDN); mpfr_init(r19168); mpfr_init(r19169); mpfr_init(r19170); mpfr_init(r19171); mpfr_init(r19172); mpfr_init(r19173); } double f_im(double re, double im) { ; mpfr_set_d(r19164, re, MPFR_RNDN); mpfr_cos(r19165, r19164, MPFR_RNDN); mpfr_mul(r19166, r19163, r19165, MPFR_RNDN); ; mpfr_set_d(r19168, im, MPFR_RNDN); mpfr_sub(r19169, r19167, r19168, MPFR_RNDN); mpfr_exp(r19170, r19169, MPFR_RNDN); mpfr_exp(r19171, r19168, MPFR_RNDN); mpfr_sub(r19172, r19170, r19171, MPFR_RNDN); mpfr_mul(r19173, r19166, r19172, MPFR_RNDN); return mpfr_get_d(r19173, MPFR_RNDN); } static mpfr_t r19174, r19175, r19176, r19177, r19178, r19179, r19180, r19181, r19182, r19183, r19184, r19185, r19186, r19187, r19188, r19189, r19190, r19191, r19192; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19174, "1/60", 10, MPFR_RNDN); mpfr_init(r19175); mpfr_init_set_str(r19176, "5", 10, MPFR_RNDN); mpfr_init(r19177); mpfr_init(r19178); mpfr_init_set_str(r19179, "2", 10, MPFR_RNDN); mpfr_init(r19180); mpfr_init_set_str(r19181, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19182, "3", 10, MPFR_RNDN); mpfr_init(r19183); mpfr_init(r19184); mpfr_init(r19185); mpfr_init(r19186); mpfr_init(r19187); mpfr_init(r19188); mpfr_init(r19189); mpfr_init_set_str(r19190, "0.5", 10, MPFR_RNDN); mpfr_init(r19191); mpfr_init(r19192); } double f_fm(double re, double im) { ; mpfr_set_d(r19175, im, MPFR_RNDN); ; mpfr_pow(r19177, r19175, r19176, MPFR_RNDN); mpfr_mul(r19178, r19174, r19177, MPFR_RNDN); ; mpfr_mul(r19180, r19179, r19175, MPFR_RNDN); ; ; mpfr_pow(r19183, r19175, r19182, MPFR_RNDN); mpfr_mul(r19184, r19181, r19183, MPFR_RNDN); mpfr_add(r19185, r19180, r19184, MPFR_RNDN); mpfr_add(r19186, r19178, r19185, MPFR_RNDN); mpfr_neg(r19187, r19186, MPFR_RNDN); mpfr_set_d(r19188, re, MPFR_RNDN); mpfr_cos(r19189, r19188, MPFR_RNDN); ; mpfr_mul(r19191, r19189, r19190, MPFR_RNDN); mpfr_mul(r19192, r19187, r19191, MPFR_RNDN); return mpfr_get_d(r19192, MPFR_RNDN); } static mpfr_t r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201, r19202, r19203, r19204, r19205, r19206, r19207, r19208, r19209, r19210, r19211; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19193, "1/60", 10, MPFR_RNDN); mpfr_init(r19194); mpfr_init_set_str(r19195, "5", 10, MPFR_RNDN); mpfr_init(r19196); mpfr_init(r19197); mpfr_init_set_str(r19198, "2", 10, MPFR_RNDN); mpfr_init(r19199); mpfr_init_set_str(r19200, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19201, "3", 10, MPFR_RNDN); mpfr_init(r19202); mpfr_init(r19203); mpfr_init(r19204); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); mpfr_init(r19208); mpfr_init_set_str(r19209, "0.5", 10, MPFR_RNDN); mpfr_init(r19210); mpfr_init(r19211); } double f_dm(double re, double im) { ; mpfr_set_d(r19194, im, MPFR_RNDN); ; mpfr_pow(r19196, r19194, r19195, MPFR_RNDN); mpfr_mul(r19197, r19193, r19196, MPFR_RNDN); ; mpfr_mul(r19199, r19198, r19194, MPFR_RNDN); ; ; mpfr_pow(r19202, r19194, r19201, MPFR_RNDN); mpfr_mul(r19203, r19200, r19202, MPFR_RNDN); mpfr_add(r19204, r19199, r19203, MPFR_RNDN); mpfr_add(r19205, r19197, r19204, MPFR_RNDN); mpfr_neg(r19206, r19205, MPFR_RNDN); mpfr_set_d(r19207, re, MPFR_RNDN); mpfr_cos(r19208, r19207, MPFR_RNDN); ; mpfr_mul(r19210, r19208, r19209, MPFR_RNDN); mpfr_mul(r19211, r19206, r19210, MPFR_RNDN); return mpfr_get_d(r19211, MPFR_RNDN); }