#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19155 = 0.5f; float r19156 = re; float r19157 = cos(r19156); float r19158 = r19155 * r19157; float r19159 = 0.0f; float r19160 = im; float r19161 = r19159 - r19160; float r19162 = exp(r19161); float r19163 = exp(r19160); float r19164 = r19162 - r19163; float r19165 = r19158 * r19164; return r19165; } double f_id(double re, double im) { double r19166 = 0.5; double r19167 = re; double r19168 = cos(r19167); double r19169 = r19166 * r19168; double r19170 = 0.0; double r19171 = im; double r19172 = r19170 - r19171; double r19173 = exp(r19172); double r19174 = exp(r19171); double r19175 = r19173 - r19174; double r19176 = r19169 * r19175; return r19176; } double f_of(float re, float im) { float r19177 = 0.01666666753590107f; float r19178 = im; float r19179 = 5.0f; float r19180 = pow(r19178, r19179); float r19181 = r19177 * r19180; float r19182 = 2.0f; float r19183 = r19182 * r19178; float r19184 = 0.3333333432674408f; float r19185 = 3.0f; float r19186 = pow(r19178, r19185); float r19187 = r19184 * r19186; float r19188 = r19183 + r19187; float r19189 = r19181 + r19188; float r19190 = -r19189; float r19191 = re; float r19192 = cos(r19191); float r19193 = 0.5f; float r19194 = r19192 * r19193; float r19195 = r19190 * r19194; return r19195; } double f_od(double re, double im) { double r19196 = 0.01666666753590107; double r19197 = im; double r19198 = 5.0; double r19199 = pow(r19197, r19198); double r19200 = r19196 * r19199; double r19201 = 2.0; double r19202 = r19201 * r19197; double r19203 = 0.3333333432674408; double r19204 = 3.0; double r19205 = pow(r19197, r19204); double r19206 = r19203 * r19205; double r19207 = r19202 + r19206; double r19208 = r19200 + r19207; double r19209 = -r19208; double r19210 = re; double r19211 = cos(r19210); double r19212 = 0.5; double r19213 = r19211 * r19212; double r19214 = r19209 * r19213; return r19214; } 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 r19215, r19216, r19217, r19218, r19219, r19220, r19221, r19222, r19223, r19224, r19225; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19215, "0.5", 10, MPFR_RNDN); mpfr_init(r19216); mpfr_init(r19217); mpfr_init(r19218); mpfr_init_set_str(r19219, "0", 10, MPFR_RNDN); mpfr_init(r19220); mpfr_init(r19221); mpfr_init(r19222); mpfr_init(r19223); mpfr_init(r19224); mpfr_init(r19225); } double f_im(double re, double im) { ; mpfr_set_d(r19216, re, MPFR_RNDN); mpfr_cos(r19217, r19216, MPFR_RNDN); mpfr_mul(r19218, r19215, r19217, MPFR_RNDN); ; mpfr_set_d(r19220, im, MPFR_RNDN); mpfr_sub(r19221, r19219, r19220, MPFR_RNDN); mpfr_exp(r19222, r19221, MPFR_RNDN); mpfr_exp(r19223, r19220, MPFR_RNDN); mpfr_sub(r19224, r19222, r19223, MPFR_RNDN); mpfr_mul(r19225, r19218, r19224, MPFR_RNDN); return mpfr_get_d(r19225, MPFR_RNDN); } static mpfr_t r19226, r19227, r19228, r19229, r19230, r19231, r19232, r19233, r19234, r19235, r19236, r19237, r19238, r19239, r19240, r19241, r19242, r19243, r19244; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19226, "1/60", 10, MPFR_RNDN); mpfr_init(r19227); mpfr_init_set_str(r19228, "5", 10, MPFR_RNDN); mpfr_init(r19229); mpfr_init(r19230); mpfr_init_set_str(r19231, "2", 10, MPFR_RNDN); mpfr_init(r19232); mpfr_init_set_str(r19233, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19234, "3", 10, MPFR_RNDN); mpfr_init(r19235); mpfr_init(r19236); mpfr_init(r19237); mpfr_init(r19238); mpfr_init(r19239); mpfr_init(r19240); mpfr_init(r19241); mpfr_init_set_str(r19242, "0.5", 10, MPFR_RNDN); mpfr_init(r19243); mpfr_init(r19244); } double f_fm(double re, double im) { ; mpfr_set_d(r19227, im, MPFR_RNDN); ; mpfr_pow(r19229, r19227, r19228, MPFR_RNDN); mpfr_mul(r19230, r19226, r19229, MPFR_RNDN); ; mpfr_mul(r19232, r19231, r19227, MPFR_RNDN); ; ; mpfr_pow(r19235, r19227, r19234, MPFR_RNDN); mpfr_mul(r19236, r19233, r19235, MPFR_RNDN); mpfr_add(r19237, r19232, r19236, MPFR_RNDN); mpfr_add(r19238, r19230, r19237, MPFR_RNDN); mpfr_neg(r19239, r19238, MPFR_RNDN); mpfr_set_d(r19240, re, MPFR_RNDN); mpfr_cos(r19241, r19240, MPFR_RNDN); ; mpfr_mul(r19243, r19241, r19242, MPFR_RNDN); mpfr_mul(r19244, r19239, r19243, MPFR_RNDN); return mpfr_get_d(r19244, MPFR_RNDN); } static mpfr_t r19245, r19246, r19247, r19248, r19249, r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19245, "1/60", 10, MPFR_RNDN); mpfr_init(r19246); mpfr_init_set_str(r19247, "5", 10, MPFR_RNDN); mpfr_init(r19248); mpfr_init(r19249); mpfr_init_set_str(r19250, "2", 10, MPFR_RNDN); mpfr_init(r19251); mpfr_init_set_str(r19252, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19253, "3", 10, MPFR_RNDN); mpfr_init(r19254); mpfr_init(r19255); mpfr_init(r19256); mpfr_init(r19257); mpfr_init(r19258); mpfr_init(r19259); mpfr_init(r19260); mpfr_init_set_str(r19261, "0.5", 10, MPFR_RNDN); mpfr_init(r19262); mpfr_init(r19263); } double f_dm(double re, double im) { ; mpfr_set_d(r19246, im, MPFR_RNDN); ; mpfr_pow(r19248, r19246, r19247, MPFR_RNDN); mpfr_mul(r19249, r19245, r19248, MPFR_RNDN); ; mpfr_mul(r19251, r19250, r19246, MPFR_RNDN); ; ; mpfr_pow(r19254, r19246, r19253, MPFR_RNDN); mpfr_mul(r19255, r19252, r19254, MPFR_RNDN); mpfr_add(r19256, r19251, r19255, MPFR_RNDN); mpfr_add(r19257, r19249, r19256, MPFR_RNDN); mpfr_neg(r19258, r19257, MPFR_RNDN); mpfr_set_d(r19259, re, MPFR_RNDN); mpfr_cos(r19260, r19259, MPFR_RNDN); ; mpfr_mul(r19262, r19260, r19261, MPFR_RNDN); mpfr_mul(r19263, r19258, r19262, MPFR_RNDN); return mpfr_get_d(r19263, MPFR_RNDN); }