#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20253 = 0.5f; float r20254 = re; float r20255 = sin(r20254); float r20256 = r20253 * r20255; float r20257 = im; float r20258 = -r20257; float r20259 = exp(r20258); float r20260 = exp(r20257); float r20261 = r20259 - r20260; float r20262 = r20256 * r20261; return r20262; } double f_id(double re, double im) { double r20263 = 0.5; double r20264 = re; double r20265 = sin(r20264); double r20266 = r20263 * r20265; double r20267 = im; double r20268 = -r20267; double r20269 = exp(r20268); double r20270 = exp(r20267); double r20271 = r20269 - r20270; double r20272 = r20266 * r20271; return r20272; } double f_of(float re, float im) { float r20273 = 0.5f; float r20274 = re; float r20275 = sin(r20274); float r20276 = r20273 * r20275; float r20277 = 0.016666666666666666f; float r20278 = im; float r20279 = 5.0f; float r20280 = pow(r20278, r20279); float r20281 = r20277 * r20280; float r20282 = 2.0f; float r20283 = r20282 * r20278; float r20284 = 0.3333333333333333f; float r20285 = 3.0f; float r20286 = pow(r20278, r20285); float r20287 = r20284 * r20286; float r20288 = r20283 + r20287; float r20289 = r20281 + r20288; float r20290 = -r20289; float r20291 = r20276 * r20290; return r20291; } double f_od(double re, double im) { double r20292 = 0.5; double r20293 = re; double r20294 = sin(r20293); double r20295 = r20292 * r20294; double r20296 = 0.016666666666666666; double r20297 = im; double r20298 = 5.0; double r20299 = pow(r20297, r20298); double r20300 = r20296 * r20299; double r20301 = 2.0; double r20302 = r20301 * r20297; double r20303 = 0.3333333333333333; double r20304 = 3.0; double r20305 = pow(r20297, r20304); double r20306 = r20303 * r20305; double r20307 = r20302 + r20306; double r20308 = r20300 + r20307; double r20309 = -r20308; double r20310 = r20295 * r20309; return r20310; } 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 r20311, r20312, r20313, r20314, r20315, r20316, r20317, r20318, r20319, r20320; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20311, "0.5", 10, MPFR_RNDN); mpfr_init(r20312); mpfr_init(r20313); mpfr_init(r20314); mpfr_init(r20315); mpfr_init(r20316); mpfr_init(r20317); mpfr_init(r20318); mpfr_init(r20319); mpfr_init(r20320); } double f_im(double re, double im) { ; mpfr_set_d(r20312, re, MPFR_RNDN); mpfr_sin(r20313, r20312, MPFR_RNDN); mpfr_mul(r20314, r20311, r20313, MPFR_RNDN); mpfr_set_d(r20315, im, MPFR_RNDN); mpfr_neg(r20316, r20315, MPFR_RNDN); mpfr_exp(r20317, r20316, MPFR_RNDN); mpfr_exp(r20318, r20315, MPFR_RNDN); mpfr_sub(r20319, r20317, r20318, MPFR_RNDN); mpfr_mul(r20320, r20314, r20319, MPFR_RNDN); return mpfr_get_d(r20320, MPFR_RNDN); } static mpfr_t r20321, r20322, r20323, r20324, r20325, r20326, r20327, r20328, r20329, r20330, r20331, r20332, r20333, r20334, r20335, r20336, r20337, r20338, r20339; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20321, "0.5", 10, MPFR_RNDN); mpfr_init(r20322); mpfr_init(r20323); mpfr_init(r20324); mpfr_init_set_str(r20325, "1/60", 10, MPFR_RNDN); mpfr_init(r20326); mpfr_init_set_str(r20327, "5", 10, MPFR_RNDN); mpfr_init(r20328); mpfr_init(r20329); mpfr_init_set_str(r20330, "2", 10, MPFR_RNDN); mpfr_init(r20331); mpfr_init_set_str(r20332, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20333, "3", 10, MPFR_RNDN); mpfr_init(r20334); mpfr_init(r20335); mpfr_init(r20336); mpfr_init(r20337); mpfr_init(r20338); mpfr_init(r20339); } double f_fm(double re, double im) { ; mpfr_set_d(r20322, re, MPFR_RNDN); mpfr_sin(r20323, r20322, MPFR_RNDN); mpfr_mul(r20324, r20321, r20323, MPFR_RNDN); ; mpfr_set_d(r20326, im, MPFR_RNDN); ; mpfr_pow(r20328, r20326, r20327, MPFR_RNDN); mpfr_mul(r20329, r20325, r20328, MPFR_RNDN); ; mpfr_mul(r20331, r20330, r20326, MPFR_RNDN); ; ; mpfr_pow(r20334, r20326, r20333, MPFR_RNDN); mpfr_mul(r20335, r20332, r20334, MPFR_RNDN); mpfr_add(r20336, r20331, r20335, MPFR_RNDN); mpfr_add(r20337, r20329, r20336, MPFR_RNDN); mpfr_neg(r20338, r20337, MPFR_RNDN); mpfr_mul(r20339, r20324, r20338, MPFR_RNDN); return mpfr_get_d(r20339, MPFR_RNDN); } static mpfr_t r20340, r20341, r20342, r20343, r20344, r20345, r20346, r20347, r20348, r20349, r20350, r20351, r20352, r20353, r20354, r20355, r20356, r20357, r20358; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20340, "0.5", 10, MPFR_RNDN); mpfr_init(r20341); mpfr_init(r20342); mpfr_init(r20343); mpfr_init_set_str(r20344, "1/60", 10, MPFR_RNDN); mpfr_init(r20345); mpfr_init_set_str(r20346, "5", 10, MPFR_RNDN); mpfr_init(r20347); mpfr_init(r20348); mpfr_init_set_str(r20349, "2", 10, MPFR_RNDN); mpfr_init(r20350); mpfr_init_set_str(r20351, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20352, "3", 10, MPFR_RNDN); mpfr_init(r20353); mpfr_init(r20354); mpfr_init(r20355); mpfr_init(r20356); mpfr_init(r20357); mpfr_init(r20358); } double f_dm(double re, double im) { ; mpfr_set_d(r20341, re, MPFR_RNDN); mpfr_sin(r20342, r20341, MPFR_RNDN); mpfr_mul(r20343, r20340, r20342, MPFR_RNDN); ; mpfr_set_d(r20345, im, MPFR_RNDN); ; mpfr_pow(r20347, r20345, r20346, MPFR_RNDN); mpfr_mul(r20348, r20344, r20347, MPFR_RNDN); ; mpfr_mul(r20350, r20349, r20345, MPFR_RNDN); ; ; mpfr_pow(r20353, r20345, r20352, MPFR_RNDN); mpfr_mul(r20354, r20351, r20353, MPFR_RNDN); mpfr_add(r20355, r20350, r20354, MPFR_RNDN); mpfr_add(r20356, r20348, r20355, MPFR_RNDN); mpfr_neg(r20357, r20356, MPFR_RNDN); mpfr_mul(r20358, r20343, r20357, MPFR_RNDN); return mpfr_get_d(r20358, MPFR_RNDN); }