#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18239 = 0.5f; float r18240 = re; float r18241 = cos(r18240); float r18242 = r18239 * r18241; float r18243 = 0.0f; float r18244 = im; float r18245 = r18243 - r18244; float r18246 = exp(r18245); float r18247 = exp(r18244); float r18248 = r18246 - r18247; float r18249 = r18242 * r18248; return r18249; } double f_id(double re, double im) { double r18250 = 0.5; double r18251 = re; double r18252 = cos(r18251); double r18253 = r18250 * r18252; double r18254 = 0.0; double r18255 = im; double r18256 = r18254 - r18255; double r18257 = exp(r18256); double r18258 = exp(r18255); double r18259 = r18257 - r18258; double r18260 = r18253 * r18259; return r18260; } double f_of(float re, float im) { float r18261 = 0.01666666753590107f; float r18262 = im; float r18263 = 5.0f; float r18264 = pow(r18262, r18263); float r18265 = r18261 * r18264; float r18266 = 2.0f; float r18267 = r18266 * r18262; float r18268 = 0.3333333432674408f; float r18269 = 3.0f; float r18270 = pow(r18262, r18269); float r18271 = r18268 * r18270; float r18272 = r18267 + r18271; float r18273 = r18265 + r18272; float r18274 = -r18273; float r18275 = re; float r18276 = cos(r18275); float r18277 = 0.5f; float r18278 = r18276 * r18277; float r18279 = r18274 * r18278; return r18279; } double f_od(double re, double im) { double r18280 = 0.01666666753590107; double r18281 = im; double r18282 = 5.0; double r18283 = pow(r18281, r18282); double r18284 = r18280 * r18283; double r18285 = 2.0; double r18286 = r18285 * r18281; double r18287 = 0.3333333432674408; double r18288 = 3.0; double r18289 = pow(r18281, r18288); double r18290 = r18287 * r18289; double r18291 = r18286 + r18290; double r18292 = r18284 + r18291; double r18293 = -r18292; double r18294 = re; double r18295 = cos(r18294); double r18296 = 0.5; double r18297 = r18295 * r18296; double r18298 = r18293 * r18297; return r18298; } 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 r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18299, "0.5", 10, MPFR_RNDN); mpfr_init(r18300); mpfr_init(r18301); mpfr_init(r18302); mpfr_init_set_str(r18303, "0", 10, MPFR_RNDN); mpfr_init(r18304); mpfr_init(r18305); mpfr_init(r18306); mpfr_init(r18307); mpfr_init(r18308); mpfr_init(r18309); } double f_im(double re, double im) { ; mpfr_set_d(r18300, re, MPFR_RNDN); mpfr_cos(r18301, r18300, MPFR_RNDN); mpfr_mul(r18302, r18299, r18301, MPFR_RNDN); ; mpfr_set_d(r18304, im, MPFR_RNDN); mpfr_sub(r18305, r18303, r18304, MPFR_RNDN); mpfr_exp(r18306, r18305, MPFR_RNDN); mpfr_exp(r18307, r18304, MPFR_RNDN); mpfr_sub(r18308, r18306, r18307, MPFR_RNDN); mpfr_mul(r18309, r18302, r18308, MPFR_RNDN); return mpfr_get_d(r18309, MPFR_RNDN); } static mpfr_t r18310, r18311, r18312, r18313, r18314, r18315, r18316, r18317, r18318, r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326, r18327, r18328; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18310, "1/60", 10, MPFR_RNDN); mpfr_init(r18311); mpfr_init_set_str(r18312, "5", 10, MPFR_RNDN); mpfr_init(r18313); mpfr_init(r18314); mpfr_init_set_str(r18315, "2", 10, MPFR_RNDN); mpfr_init(r18316); mpfr_init_set_str(r18317, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18318, "3", 10, MPFR_RNDN); mpfr_init(r18319); mpfr_init(r18320); mpfr_init(r18321); mpfr_init(r18322); mpfr_init(r18323); mpfr_init(r18324); mpfr_init(r18325); mpfr_init_set_str(r18326, "0.5", 10, MPFR_RNDN); mpfr_init(r18327); mpfr_init(r18328); } double f_fm(double re, double im) { ; mpfr_set_d(r18311, im, MPFR_RNDN); ; mpfr_pow(r18313, r18311, r18312, MPFR_RNDN); mpfr_mul(r18314, r18310, r18313, MPFR_RNDN); ; mpfr_mul(r18316, r18315, r18311, MPFR_RNDN); ; ; mpfr_pow(r18319, r18311, r18318, MPFR_RNDN); mpfr_mul(r18320, r18317, r18319, MPFR_RNDN); mpfr_add(r18321, r18316, r18320, MPFR_RNDN); mpfr_add(r18322, r18314, r18321, MPFR_RNDN); mpfr_neg(r18323, r18322, MPFR_RNDN); mpfr_set_d(r18324, re, MPFR_RNDN); mpfr_cos(r18325, r18324, MPFR_RNDN); ; mpfr_mul(r18327, r18325, r18326, MPFR_RNDN); mpfr_mul(r18328, r18323, r18327, MPFR_RNDN); return mpfr_get_d(r18328, MPFR_RNDN); } static mpfr_t r18329, r18330, r18331, r18332, r18333, r18334, r18335, r18336, r18337, r18338, r18339, r18340, r18341, r18342, r18343, r18344, r18345, r18346, r18347; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18329, "1/60", 10, MPFR_RNDN); mpfr_init(r18330); mpfr_init_set_str(r18331, "5", 10, MPFR_RNDN); mpfr_init(r18332); mpfr_init(r18333); mpfr_init_set_str(r18334, "2", 10, MPFR_RNDN); mpfr_init(r18335); mpfr_init_set_str(r18336, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18337, "3", 10, MPFR_RNDN); mpfr_init(r18338); mpfr_init(r18339); mpfr_init(r18340); mpfr_init(r18341); mpfr_init(r18342); mpfr_init(r18343); mpfr_init(r18344); mpfr_init_set_str(r18345, "0.5", 10, MPFR_RNDN); mpfr_init(r18346); mpfr_init(r18347); } double f_dm(double re, double im) { ; mpfr_set_d(r18330, im, MPFR_RNDN); ; mpfr_pow(r18332, r18330, r18331, MPFR_RNDN); mpfr_mul(r18333, r18329, r18332, MPFR_RNDN); ; mpfr_mul(r18335, r18334, r18330, MPFR_RNDN); ; ; mpfr_pow(r18338, r18330, r18337, MPFR_RNDN); mpfr_mul(r18339, r18336, r18338, MPFR_RNDN); mpfr_add(r18340, r18335, r18339, MPFR_RNDN); mpfr_add(r18341, r18333, r18340, MPFR_RNDN); mpfr_neg(r18342, r18341, MPFR_RNDN); mpfr_set_d(r18343, re, MPFR_RNDN); mpfr_cos(r18344, r18343, MPFR_RNDN); ; mpfr_mul(r18346, r18344, r18345, MPFR_RNDN); mpfr_mul(r18347, r18342, r18346, MPFR_RNDN); return mpfr_get_d(r18347, MPFR_RNDN); }