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