#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r25461 = 0.5; float r25462 = re; float r25463 = sin(r25462); float r25464 = r25461 * r25463; float r25465 = im; float r25466 = -r25465; float r25467 = exp(r25466); float r25468 = exp(r25465); float r25469 = r25467 - r25468; float r25470 = r25464 * r25469; return r25470; } double f_id(double re, double im) { double r25471 = 0.5; double r25472 = re; double r25473 = sin(r25472); double r25474 = r25471 * r25473; double r25475 = im; double r25476 = -r25475; double r25477 = exp(r25476); double r25478 = exp(r25475); double r25479 = r25477 - r25478; double r25480 = r25474 * r25479; return r25480; } double f_of(float re, float im) { float r25481 = 1/3; float r25482 = im; float r25483 = r25482 * r25482; float r25484 = r25481 * r25483; float r25485 = 2; float r25486 = r25484 + r25485; float r25487 = re; float r25488 = sin(r25487); float r25489 = 0.5; float r25490 = -r25489; float r25491 = r25488 * r25490; float r25492 = r25482 * r25491; float r25493 = r25486 * r25492; float r25494 = 5; float r25495 = pow(r25482, r25494); float r25496 = 1/60; float r25497 = r25495 * r25496; float r25498 = r25497 * r25491; float r25499 = r25493 + r25498; float r25500 = 1; float r25501 = pow(r25499, r25500); return r25501; } double f_od(double re, double im) { double r25502 = 1/3; double r25503 = im; double r25504 = r25503 * r25503; double r25505 = r25502 * r25504; double r25506 = 2; double r25507 = r25505 + r25506; double r25508 = re; double r25509 = sin(r25508); double r25510 = 0.5; double r25511 = -r25510; double r25512 = r25509 * r25511; double r25513 = r25503 * r25512; double r25514 = r25507 * r25513; double r25515 = 5; double r25516 = pow(r25503, r25515); double r25517 = 1/60; double r25518 = r25516 * r25517; double r25519 = r25518 * r25512; double r25520 = r25514 + r25519; double r25521 = 1; double r25522 = pow(r25520, r25521); return r25522; } 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 r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25523, "0.5", 10, MPFR_RNDN); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init(r25530); mpfr_init(r25531); mpfr_init(r25532); } double f_im(double re, double im) { ; mpfr_set_d(r25524, re, MPFR_RNDN); mpfr_sin(r25525, r25524, MPFR_RNDN); mpfr_mul(r25526, r25523, r25525, MPFR_RNDN); mpfr_set_d(r25527, im, MPFR_RNDN); mpfr_neg(r25528, r25527, MPFR_RNDN); mpfr_exp(r25529, r25528, MPFR_RNDN); mpfr_exp(r25530, r25527, MPFR_RNDN); mpfr_sub(r25531, r25529, r25530, MPFR_RNDN); mpfr_mul(r25532, r25526, r25531, MPFR_RNDN); return mpfr_get_d(r25532, MPFR_RNDN); } static mpfr_t r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25533, "1/3", 10, MPFR_RNDN); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init_set_str(r25537, "2", 10, MPFR_RNDN); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init_set_str(r25541, "0.5", 10, MPFR_RNDN); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init_set_str(r25546, "5", 10, MPFR_RNDN); mpfr_init(r25547); mpfr_init_set_str(r25548, "1/60", 10, MPFR_RNDN); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); mpfr_init_set_str(r25552, "1", 10, MPFR_RNDN); mpfr_init(r25553); } double f_fm(double re, double im) { ; mpfr_set_d(r25534, im, MPFR_RNDN); mpfr_mul(r25535, r25534, r25534, MPFR_RNDN); mpfr_mul(r25536, r25533, r25535, MPFR_RNDN); ; mpfr_add(r25538, r25536, r25537, MPFR_RNDN); mpfr_set_d(r25539, re, MPFR_RNDN); mpfr_sin(r25540, r25539, MPFR_RNDN); ; mpfr_neg(r25542, r25541, MPFR_RNDN); mpfr_mul(r25543, r25540, r25542, MPFR_RNDN); mpfr_mul(r25544, r25534, r25543, MPFR_RNDN); mpfr_mul(r25545, r25538, r25544, MPFR_RNDN); ; mpfr_pow(r25547, r25534, r25546, MPFR_RNDN); ; mpfr_mul(r25549, r25547, r25548, MPFR_RNDN); mpfr_mul(r25550, r25549, r25543, MPFR_RNDN); mpfr_add(r25551, r25545, r25550, MPFR_RNDN); ; mpfr_pow(r25553, r25551, r25552, MPFR_RNDN); return mpfr_get_d(r25553, MPFR_RNDN); } static mpfr_t r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25554, "1/3", 10, MPFR_RNDN); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init_set_str(r25558, "2", 10, MPFR_RNDN); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init_set_str(r25562, "0.5", 10, MPFR_RNDN); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init_set_str(r25567, "5", 10, MPFR_RNDN); mpfr_init(r25568); mpfr_init_set_str(r25569, "1/60", 10, MPFR_RNDN); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init_set_str(r25573, "1", 10, MPFR_RNDN); mpfr_init(r25574); } double f_dm(double re, double im) { ; mpfr_set_d(r25555, im, MPFR_RNDN); mpfr_mul(r25556, r25555, r25555, MPFR_RNDN); mpfr_mul(r25557, r25554, r25556, MPFR_RNDN); ; mpfr_add(r25559, r25557, r25558, MPFR_RNDN); mpfr_set_d(r25560, re, MPFR_RNDN); mpfr_sin(r25561, r25560, MPFR_RNDN); ; mpfr_neg(r25563, r25562, MPFR_RNDN); mpfr_mul(r25564, r25561, r25563, MPFR_RNDN); mpfr_mul(r25565, r25555, r25564, MPFR_RNDN); mpfr_mul(r25566, r25559, r25565, MPFR_RNDN); ; mpfr_pow(r25568, r25555, r25567, MPFR_RNDN); ; mpfr_mul(r25570, r25568, r25569, MPFR_RNDN); mpfr_mul(r25571, r25570, r25564, MPFR_RNDN); mpfr_add(r25572, r25566, r25571, MPFR_RNDN); ; mpfr_pow(r25574, r25572, r25573, MPFR_RNDN); return mpfr_get_d(r25574, MPFR_RNDN); }