#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r18333 = 0.5f; float r18334 = re; float r18335 = sin(r18334); float r18336 = r18333 * r18335; float r18337 = im; float r18338 = -r18337; float r18339 = exp(r18338); float r18340 = exp(r18337); float r18341 = r18339 - r18340; float r18342 = r18336 * r18341; return r18342; } double f_id(double re, double im) { double r18343 = 0.5; double r18344 = re; double r18345 = sin(r18344); double r18346 = r18343 * r18345; double r18347 = im; double r18348 = -r18347; double r18349 = exp(r18348); double r18350 = exp(r18347); double r18351 = r18349 - r18350; double r18352 = r18346 * r18351; return r18352; } double f_of(float re, float im) { float r18353 = im; float r18354 = r18353 * (r18353 * r18353); float r18355 = exp(r18354); float r18356 = log(r18355); float r18357 = 0.3333333432674408f; float r18358 = 5.0f; float r18359 = pow(r18353, r18358); float r18360 = 0.01666666753590107f; float r18361 = 2.0f; float r18362 = r18353 * r18361; float r18363 = fma(r18359, r18360, r18362); float r18364 = fma(r18356, r18357, r18363); float r18365 = re; float r18366 = sin(r18365); float r18367 = 0.5f; float r18368 = -r18367; float r18369 = r18366 * r18368; float r18370 = r18364 * r18369; return r18370; } double f_od(double re, double im) { double r18371 = im; double r18372 = r18371 * (r18371 * r18371); double r18373 = exp(r18372); double r18374 = log(r18373); double r18375 = 0.3333333432674408; double r18376 = 5.0; double r18377 = pow(r18371, r18376); double r18378 = 0.01666666753590107; double r18379 = 2.0; double r18380 = r18371 * r18379; double r18381 = fma(r18377, r18378, r18380); double r18382 = fma(r18374, r18375, r18381); double r18383 = re; double r18384 = sin(r18383); double r18385 = 0.5; double r18386 = -r18385; double r18387 = r18384 * r18386; double r18388 = r18382 * r18387; return r18388; } 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 r18389, r18390, r18391, r18392, r18393, r18394, r18395, r18396, r18397, r18398; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18389, "0.5", 10, MPFR_RNDN); mpfr_init(r18390); mpfr_init(r18391); mpfr_init(r18392); mpfr_init(r18393); mpfr_init(r18394); mpfr_init(r18395); mpfr_init(r18396); mpfr_init(r18397); mpfr_init(r18398); } double f_im(double re, double im) { ; mpfr_set_d(r18390, re, MPFR_RNDN); mpfr_sin(r18391, r18390, MPFR_RNDN); mpfr_mul(r18392, r18389, r18391, MPFR_RNDN); mpfr_set_d(r18393, im, MPFR_RNDN); mpfr_neg(r18394, r18393, MPFR_RNDN); mpfr_exp(r18395, r18394, MPFR_RNDN); mpfr_exp(r18396, r18393, MPFR_RNDN); mpfr_sub(r18397, r18395, r18396, MPFR_RNDN); mpfr_mul(r18398, r18392, r18397, MPFR_RNDN); return mpfr_get_d(r18398, MPFR_RNDN); } static mpfr_t r18399, r18400, r18401, r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410, r18411, r18412, r18413, r18414, r18415, r18416; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18399); mpfr_init(r18400); mpfr_init(r18401); mpfr_init(r18402); mpfr_init_set_str(r18403, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18404, "5", 10, MPFR_RNDN); mpfr_init(r18405); mpfr_init_set_str(r18406, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18407, "2", 10, MPFR_RNDN); mpfr_init(r18408); mpfr_init(r18409); mpfr_init(r18410); mpfr_init(r18411); mpfr_init(r18412); mpfr_init_set_str(r18413, "0.5", 10, MPFR_RNDN); mpfr_init(r18414); mpfr_init(r18415); mpfr_init(r18416); } double f_fm(double re, double im) { mpfr_set_d(r18399, im, MPFR_RNDN); mpfr_mul(r18400, r18399, r18399, MPFR_RNDN); mpfr_mul(r18400, r18400, r18399, MPFR_RNDN); mpfr_exp(r18401, r18400, MPFR_RNDN); mpfr_log(r18402, r18401, MPFR_RNDN); ; ; mpfr_pow(r18405, r18399, r18404, MPFR_RNDN); ; ; mpfr_mul(r18408, r18399, r18407, MPFR_RNDN); mpfr_fma(r18409, r18405, r18406, r18408, MPFR_RNDN); mpfr_fma(r18410, r18402, r18403, r18409, MPFR_RNDN); mpfr_set_d(r18411, re, MPFR_RNDN); mpfr_sin(r18412, r18411, MPFR_RNDN); ; mpfr_neg(r18414, r18413, MPFR_RNDN); mpfr_mul(r18415, r18412, r18414, MPFR_RNDN); mpfr_mul(r18416, r18410, r18415, MPFR_RNDN); return mpfr_get_d(r18416, MPFR_RNDN); } static mpfr_t r18417, r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432, r18433, r18434; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18417); mpfr_init(r18418); mpfr_init(r18419); mpfr_init(r18420); mpfr_init_set_str(r18421, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18422, "5", 10, MPFR_RNDN); mpfr_init(r18423); mpfr_init_set_str(r18424, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18425, "2", 10, MPFR_RNDN); mpfr_init(r18426); mpfr_init(r18427); mpfr_init(r18428); mpfr_init(r18429); mpfr_init(r18430); mpfr_init_set_str(r18431, "0.5", 10, MPFR_RNDN); mpfr_init(r18432); mpfr_init(r18433); mpfr_init(r18434); } double f_dm(double re, double im) { mpfr_set_d(r18417, im, MPFR_RNDN); mpfr_mul(r18418, r18417, r18417, MPFR_RNDN); mpfr_mul(r18418, r18418, r18417, MPFR_RNDN); mpfr_exp(r18419, r18418, MPFR_RNDN); mpfr_log(r18420, r18419, MPFR_RNDN); ; ; mpfr_pow(r18423, r18417, r18422, MPFR_RNDN); ; ; mpfr_mul(r18426, r18417, r18425, MPFR_RNDN); mpfr_fma(r18427, r18423, r18424, r18426, MPFR_RNDN); mpfr_fma(r18428, r18420, r18421, r18427, MPFR_RNDN); mpfr_set_d(r18429, re, MPFR_RNDN); mpfr_sin(r18430, r18429, MPFR_RNDN); ; mpfr_neg(r18432, r18431, MPFR_RNDN); mpfr_mul(r18433, r18430, r18432, MPFR_RNDN); mpfr_mul(r18434, r18428, r18433, MPFR_RNDN); return mpfr_get_d(r18434, MPFR_RNDN); }