#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19207 = 0.5f; float r19208 = re; float r19209 = cos(r19208); float r19210 = r19207 * r19209; float r19211 = 0.0f; float r19212 = im; float r19213 = r19211 - r19212; float r19214 = exp(r19213); float r19215 = exp(r19212); float r19216 = r19214 - r19215; float r19217 = r19210 * r19216; return r19217; } double f_id(double re, double im) { double r19218 = 0.5; double r19219 = re; double r19220 = cos(r19219); double r19221 = r19218 * r19220; double r19222 = 0.0; double r19223 = im; double r19224 = r19222 - r19223; double r19225 = exp(r19224); double r19226 = exp(r19223); double r19227 = r19225 - r19226; double r19228 = r19221 * r19227; return r19228; } double f_of(float re, float im) { float r19229 = 0.016666666666666666f; float r19230 = im; float r19231 = 5.0f; float r19232 = pow(r19230, r19231); float r19233 = r19230 * (r19230 * r19230); float r19234 = 0.3333333333333333f; float r19235 = 2.0f; float r19236 = r19235 * r19230; float r19237 = fma(r19233, r19234, r19236); float r19238 = fma(r19229, r19232, r19237); float r19239 = 1.0f; float r19240 = r19238 / r19239; float r19241 = re; float r19242 = cos(r19241); float r19243 = 0.5f; float r19244 = -r19243; float r19245 = r19242 * r19244; float r19246 = r19240 * r19245; return r19246; } double f_od(double re, double im) { double r19247 = 0.016666666666666666; double r19248 = im; double r19249 = 5.0; double r19250 = pow(r19248, r19249); double r19251 = r19248 * (r19248 * r19248); double r19252 = 0.3333333333333333; double r19253 = 2.0; double r19254 = r19253 * r19248; double r19255 = fma(r19251, r19252, r19254); double r19256 = fma(r19247, r19250, r19255); double r19257 = 1.0; double r19258 = r19256 / r19257; double r19259 = re; double r19260 = cos(r19259); double r19261 = 0.5; double r19262 = -r19261; double r19263 = r19260 * r19262; double r19264 = r19258 * r19263; return r19264; } 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 r19265, r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19265, "0.5", 10, MPFR_RNDN); mpfr_init(r19266); mpfr_init(r19267); mpfr_init(r19268); mpfr_init_set_str(r19269, "0", 10, MPFR_RNDN); mpfr_init(r19270); mpfr_init(r19271); mpfr_init(r19272); mpfr_init(r19273); mpfr_init(r19274); mpfr_init(r19275); } double f_im(double re, double im) { ; mpfr_set_d(r19266, re, MPFR_RNDN); mpfr_cos(r19267, r19266, MPFR_RNDN); mpfr_mul(r19268, r19265, r19267, MPFR_RNDN); ; mpfr_set_d(r19270, im, MPFR_RNDN); mpfr_sub(r19271, r19269, r19270, MPFR_RNDN); mpfr_exp(r19272, r19271, MPFR_RNDN); mpfr_exp(r19273, r19270, MPFR_RNDN); mpfr_sub(r19274, r19272, r19273, MPFR_RNDN); mpfr_mul(r19275, r19268, r19274, MPFR_RNDN); return mpfr_get_d(r19275, MPFR_RNDN); } static mpfr_t r19276, r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284, r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19276, "1/60", 10, MPFR_RNDN); mpfr_init(r19277); mpfr_init_set_str(r19278, "5", 10, MPFR_RNDN); mpfr_init(r19279); mpfr_init(r19280); mpfr_init_set_str(r19281, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19282, "2", 10, MPFR_RNDN); mpfr_init(r19283); mpfr_init(r19284); mpfr_init(r19285); mpfr_init_set_str(r19286, "1", 10, MPFR_RNDN); mpfr_init(r19287); mpfr_init(r19288); mpfr_init(r19289); mpfr_init_set_str(r19290, "0.5", 10, MPFR_RNDN); mpfr_init(r19291); mpfr_init(r19292); mpfr_init(r19293); } double f_fm(double re, double im) { ; mpfr_set_d(r19277, im, MPFR_RNDN); ; mpfr_pow(r19279, r19277, r19278, MPFR_RNDN); mpfr_mul(r19280, r19277, r19277, MPFR_RNDN); mpfr_mul(r19280, r19280, r19277, MPFR_RNDN); ; ; mpfr_mul(r19283, r19282, r19277, MPFR_RNDN); mpfr_fma(r19284, r19280, r19281, r19283, MPFR_RNDN); mpfr_fma(r19285, r19276, r19279, r19284, MPFR_RNDN); ; mpfr_div(r19287, r19285, r19286, MPFR_RNDN); mpfr_set_d(r19288, re, MPFR_RNDN); mpfr_cos(r19289, r19288, MPFR_RNDN); ; mpfr_neg(r19291, r19290, MPFR_RNDN); mpfr_mul(r19292, r19289, r19291, MPFR_RNDN); mpfr_mul(r19293, r19287, r19292, MPFR_RNDN); return mpfr_get_d(r19293, MPFR_RNDN); } static mpfr_t r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19294, "1/60", 10, MPFR_RNDN); mpfr_init(r19295); mpfr_init_set_str(r19296, "5", 10, MPFR_RNDN); mpfr_init(r19297); mpfr_init(r19298); mpfr_init_set_str(r19299, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19300, "2", 10, MPFR_RNDN); mpfr_init(r19301); mpfr_init(r19302); mpfr_init(r19303); mpfr_init_set_str(r19304, "1", 10, MPFR_RNDN); mpfr_init(r19305); mpfr_init(r19306); mpfr_init(r19307); mpfr_init_set_str(r19308, "0.5", 10, MPFR_RNDN); mpfr_init(r19309); mpfr_init(r19310); mpfr_init(r19311); } double f_dm(double re, double im) { ; mpfr_set_d(r19295, im, MPFR_RNDN); ; mpfr_pow(r19297, r19295, r19296, MPFR_RNDN); mpfr_mul(r19298, r19295, r19295, MPFR_RNDN); mpfr_mul(r19298, r19298, r19295, MPFR_RNDN); ; ; mpfr_mul(r19301, r19300, r19295, MPFR_RNDN); mpfr_fma(r19302, r19298, r19299, r19301, MPFR_RNDN); mpfr_fma(r19303, r19294, r19297, r19302, MPFR_RNDN); ; mpfr_div(r19305, r19303, r19304, MPFR_RNDN); mpfr_set_d(r19306, re, MPFR_RNDN); mpfr_cos(r19307, r19306, MPFR_RNDN); ; mpfr_neg(r19309, r19308, MPFR_RNDN); mpfr_mul(r19310, r19307, r19309, MPFR_RNDN); mpfr_mul(r19311, r19305, r19310, MPFR_RNDN); return mpfr_get_d(r19311, MPFR_RNDN); }