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