#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20301 = 0.5f; float r20302 = re; float r20303 = sin(r20302); float r20304 = r20301 * r20303; float r20305 = im; float r20306 = -r20305; float r20307 = exp(r20306); float r20308 = exp(r20305); float r20309 = r20307 - r20308; float r20310 = r20304 * r20309; return r20310; } double f_id(double re, double im) { double r20311 = 0.5; double r20312 = re; double r20313 = sin(r20312); double r20314 = r20311 * r20313; double r20315 = im; double r20316 = -r20315; double r20317 = exp(r20316); double r20318 = exp(r20315); double r20319 = r20317 - r20318; double r20320 = r20314 * r20319; return r20320; } double f_of(float re, float im) { float r20321 = 0.5f; float r20322 = re; float r20323 = sin(r20322); float r20324 = r20321 * r20323; float r20325 = 0.01666666753590107f; float r20326 = im; float r20327 = 5.0f; float r20328 = pow(r20326, r20327); float r20329 = r20325 * r20328; float r20330 = 2.0f; float r20331 = r20330 * r20326; float r20332 = 0.3333333432674408f; float r20333 = 3.0f; float r20334 = pow(r20326, r20333); float r20335 = r20332 * r20334; float r20336 = r20331 + r20335; float r20337 = r20329 + r20336; float r20338 = -r20337; float r20339 = r20324 * r20338; return r20339; } double f_od(double re, double im) { double r20340 = 0.5; double r20341 = re; double r20342 = sin(r20341); double r20343 = r20340 * r20342; double r20344 = 0.01666666753590107; double r20345 = im; double r20346 = 5.0; double r20347 = pow(r20345, r20346); double r20348 = r20344 * r20347; double r20349 = 2.0; double r20350 = r20349 * r20345; double r20351 = 0.3333333432674408; double r20352 = 3.0; double r20353 = pow(r20345, r20352); double r20354 = r20351 * r20353; double r20355 = r20350 + r20354; double r20356 = r20348 + r20355; double r20357 = -r20356; double r20358 = r20343 * r20357; return r20358; } 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 r20359, r20360, r20361, r20362, r20363, r20364, r20365, r20366, r20367, r20368; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20359, "0.5", 10, MPFR_RNDN); mpfr_init(r20360); mpfr_init(r20361); mpfr_init(r20362); mpfr_init(r20363); mpfr_init(r20364); mpfr_init(r20365); mpfr_init(r20366); mpfr_init(r20367); mpfr_init(r20368); } double f_im(double re, double im) { ; mpfr_set_d(r20360, re, MPFR_RNDN); mpfr_sin(r20361, r20360, MPFR_RNDN); mpfr_mul(r20362, r20359, r20361, MPFR_RNDN); mpfr_set_d(r20363, im, MPFR_RNDN); mpfr_neg(r20364, r20363, MPFR_RNDN); mpfr_exp(r20365, r20364, MPFR_RNDN); mpfr_exp(r20366, r20363, MPFR_RNDN); mpfr_sub(r20367, r20365, r20366, MPFR_RNDN); mpfr_mul(r20368, r20362, r20367, MPFR_RNDN); return mpfr_get_d(r20368, MPFR_RNDN); } static mpfr_t r20369, r20370, r20371, r20372, r20373, r20374, r20375, r20376, r20377, r20378, r20379, r20380, r20381, r20382, r20383, r20384, r20385, r20386, r20387; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20369, "0.5", 10, MPFR_RNDN); mpfr_init(r20370); mpfr_init(r20371); mpfr_init(r20372); mpfr_init_set_str(r20373, "1/60", 10, MPFR_RNDN); mpfr_init(r20374); mpfr_init_set_str(r20375, "5", 10, MPFR_RNDN); mpfr_init(r20376); mpfr_init(r20377); mpfr_init_set_str(r20378, "2", 10, MPFR_RNDN); mpfr_init(r20379); mpfr_init_set_str(r20380, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20381, "3", 10, MPFR_RNDN); mpfr_init(r20382); mpfr_init(r20383); mpfr_init(r20384); mpfr_init(r20385); mpfr_init(r20386); mpfr_init(r20387); } double f_fm(double re, double im) { ; mpfr_set_d(r20370, re, MPFR_RNDN); mpfr_sin(r20371, r20370, MPFR_RNDN); mpfr_mul(r20372, r20369, r20371, MPFR_RNDN); ; mpfr_set_d(r20374, im, MPFR_RNDN); ; mpfr_pow(r20376, r20374, r20375, MPFR_RNDN); mpfr_mul(r20377, r20373, r20376, MPFR_RNDN); ; mpfr_mul(r20379, r20378, r20374, MPFR_RNDN); ; ; mpfr_pow(r20382, r20374, r20381, MPFR_RNDN); mpfr_mul(r20383, r20380, r20382, MPFR_RNDN); mpfr_add(r20384, r20379, r20383, MPFR_RNDN); mpfr_add(r20385, r20377, r20384, MPFR_RNDN); mpfr_neg(r20386, r20385, MPFR_RNDN); mpfr_mul(r20387, r20372, r20386, MPFR_RNDN); return mpfr_get_d(r20387, MPFR_RNDN); } static mpfr_t r20388, r20389, r20390, r20391, r20392, r20393, r20394, r20395, r20396, r20397, r20398, r20399, r20400, r20401, r20402, r20403, r20404, r20405, r20406; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20388, "0.5", 10, MPFR_RNDN); mpfr_init(r20389); mpfr_init(r20390); mpfr_init(r20391); mpfr_init_set_str(r20392, "1/60", 10, MPFR_RNDN); mpfr_init(r20393); mpfr_init_set_str(r20394, "5", 10, MPFR_RNDN); mpfr_init(r20395); mpfr_init(r20396); mpfr_init_set_str(r20397, "2", 10, MPFR_RNDN); mpfr_init(r20398); mpfr_init_set_str(r20399, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20400, "3", 10, MPFR_RNDN); mpfr_init(r20401); mpfr_init(r20402); mpfr_init(r20403); mpfr_init(r20404); mpfr_init(r20405); mpfr_init(r20406); } double f_dm(double re, double im) { ; mpfr_set_d(r20389, re, MPFR_RNDN); mpfr_sin(r20390, r20389, MPFR_RNDN); mpfr_mul(r20391, r20388, r20390, MPFR_RNDN); ; mpfr_set_d(r20393, im, MPFR_RNDN); ; mpfr_pow(r20395, r20393, r20394, MPFR_RNDN); mpfr_mul(r20396, r20392, r20395, MPFR_RNDN); ; mpfr_mul(r20398, r20397, r20393, MPFR_RNDN); ; ; mpfr_pow(r20401, r20393, r20400, MPFR_RNDN); mpfr_mul(r20402, r20399, r20401, MPFR_RNDN); mpfr_add(r20403, r20398, r20402, MPFR_RNDN); mpfr_add(r20404, r20396, r20403, MPFR_RNDN); mpfr_neg(r20405, r20404, MPFR_RNDN); mpfr_mul(r20406, r20391, r20405, MPFR_RNDN); return mpfr_get_d(r20406, MPFR_RNDN); }