#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19259 = 0.5f; float r19260 = re; float r19261 = cos(r19260); float r19262 = r19259 * r19261; float r19263 = 0.0f; float r19264 = im; float r19265 = r19263 - r19264; float r19266 = exp(r19265); float r19267 = exp(r19264); float r19268 = r19266 - r19267; float r19269 = r19262 * r19268; return r19269; } double f_id(double re, double im) { double r19270 = 0.5; double r19271 = re; double r19272 = cos(r19271); double r19273 = r19270 * r19272; double r19274 = 0.0; double r19275 = im; double r19276 = r19274 - r19275; double r19277 = exp(r19276); double r19278 = exp(r19275); double r19279 = r19277 - r19278; double r19280 = r19273 * r19279; return r19280; } double f_of(float re, float im) { float r19281 = 0.016666666666666666f; float r19282 = im; float r19283 = 5.0f; float r19284 = pow(r19282, r19283); float r19285 = r19282 * (r19282 * r19282); float r19286 = 0.3333333333333333f; float r19287 = 2.0f; float r19288 = r19287 * r19282; float r19289 = fma(r19285, r19286, r19288); float r19290 = fma(r19281, r19284, r19289); float r19291 = 1.0f; float r19292 = r19290 / r19291; float r19293 = re; float r19294 = cos(r19293); float r19295 = 0.5f; float r19296 = -r19295; float r19297 = r19294 * r19296; float r19298 = r19292 * r19297; return r19298; } double f_od(double re, double im) { double r19299 = 0.016666666666666666; double r19300 = im; double r19301 = 5.0; double r19302 = pow(r19300, r19301); double r19303 = r19300 * (r19300 * r19300); double r19304 = 0.3333333333333333; double r19305 = 2.0; double r19306 = r19305 * r19300; double r19307 = fma(r19303, r19304, r19306); double r19308 = fma(r19299, r19302, r19307); double r19309 = 1.0; double r19310 = r19308 / r19309; double r19311 = re; double r19312 = cos(r19311); double r19313 = 0.5; double r19314 = -r19313; double r19315 = r19312 * r19314; double r19316 = r19310 * r19315; return r19316; } 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 r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19317, "0.5", 10, MPFR_RNDN); mpfr_init(r19318); mpfr_init(r19319); mpfr_init(r19320); mpfr_init_set_str(r19321, "0", 10, MPFR_RNDN); mpfr_init(r19322); mpfr_init(r19323); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init(r19327); } double f_im(double re, double im) { ; mpfr_set_d(r19318, re, MPFR_RNDN); mpfr_cos(r19319, r19318, MPFR_RNDN); mpfr_mul(r19320, r19317, r19319, MPFR_RNDN); ; mpfr_set_d(r19322, im, MPFR_RNDN); mpfr_sub(r19323, r19321, r19322, MPFR_RNDN); mpfr_exp(r19324, r19323, MPFR_RNDN); mpfr_exp(r19325, r19322, MPFR_RNDN); mpfr_sub(r19326, r19324, r19325, MPFR_RNDN); mpfr_mul(r19327, r19320, r19326, MPFR_RNDN); return mpfr_get_d(r19327, MPFR_RNDN); } static mpfr_t r19328, r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336, r19337, r19338, r19339, r19340, r19341, r19342, r19343, r19344, r19345; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19328, "1/60", 10, MPFR_RNDN); mpfr_init(r19329); mpfr_init_set_str(r19330, "5", 10, MPFR_RNDN); mpfr_init(r19331); mpfr_init(r19332); mpfr_init_set_str(r19333, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19334, "2", 10, MPFR_RNDN); mpfr_init(r19335); mpfr_init(r19336); mpfr_init(r19337); mpfr_init_set_str(r19338, "1", 10, MPFR_RNDN); mpfr_init(r19339); mpfr_init(r19340); mpfr_init(r19341); mpfr_init_set_str(r19342, "0.5", 10, MPFR_RNDN); mpfr_init(r19343); mpfr_init(r19344); mpfr_init(r19345); } double f_fm(double re, double im) { ; mpfr_set_d(r19329, im, MPFR_RNDN); ; mpfr_pow(r19331, r19329, r19330, MPFR_RNDN); mpfr_mul(r19332, r19329, r19329, MPFR_RNDN); mpfr_mul(r19332, r19332, r19329, MPFR_RNDN); ; ; mpfr_mul(r19335, r19334, r19329, MPFR_RNDN); mpfr_fma(r19336, r19332, r19333, r19335, MPFR_RNDN); mpfr_fma(r19337, r19328, r19331, r19336, MPFR_RNDN); ; mpfr_div(r19339, r19337, r19338, MPFR_RNDN); mpfr_set_d(r19340, re, MPFR_RNDN); mpfr_cos(r19341, r19340, MPFR_RNDN); ; mpfr_neg(r19343, r19342, MPFR_RNDN); mpfr_mul(r19344, r19341, r19343, MPFR_RNDN); mpfr_mul(r19345, r19339, r19344, MPFR_RNDN); return mpfr_get_d(r19345, MPFR_RNDN); } static mpfr_t r19346, r19347, r19348, r19349, r19350, r19351, r19352, r19353, r19354, r19355, r19356, r19357, r19358, r19359, r19360, r19361, r19362, r19363; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19346, "1/60", 10, MPFR_RNDN); mpfr_init(r19347); mpfr_init_set_str(r19348, "5", 10, MPFR_RNDN); mpfr_init(r19349); mpfr_init(r19350); mpfr_init_set_str(r19351, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19352, "2", 10, MPFR_RNDN); mpfr_init(r19353); mpfr_init(r19354); mpfr_init(r19355); mpfr_init_set_str(r19356, "1", 10, MPFR_RNDN); mpfr_init(r19357); mpfr_init(r19358); mpfr_init(r19359); mpfr_init_set_str(r19360, "0.5", 10, MPFR_RNDN); mpfr_init(r19361); mpfr_init(r19362); mpfr_init(r19363); } double f_dm(double re, double im) { ; mpfr_set_d(r19347, im, MPFR_RNDN); ; mpfr_pow(r19349, r19347, r19348, MPFR_RNDN); mpfr_mul(r19350, r19347, r19347, MPFR_RNDN); mpfr_mul(r19350, r19350, r19347, MPFR_RNDN); ; ; mpfr_mul(r19353, r19352, r19347, MPFR_RNDN); mpfr_fma(r19354, r19350, r19351, r19353, MPFR_RNDN); mpfr_fma(r19355, r19346, r19349, r19354, MPFR_RNDN); ; mpfr_div(r19357, r19355, r19356, MPFR_RNDN); mpfr_set_d(r19358, re, MPFR_RNDN); mpfr_cos(r19359, r19358, MPFR_RNDN); ; mpfr_neg(r19361, r19360, MPFR_RNDN); mpfr_mul(r19362, r19359, r19361, MPFR_RNDN); mpfr_mul(r19363, r19357, r19362, MPFR_RNDN); return mpfr_get_d(r19363, MPFR_RNDN); }