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