#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25426 = 0.5; float r25427 = re; float r25428 = cos(r25427); float r25429 = r25426 * r25428; float r25430 = 0; float r25431 = im; float r25432 = r25430 - r25431; float r25433 = exp(r25432); float r25434 = exp(r25431); float r25435 = r25433 - r25434; float r25436 = r25429 * r25435; return r25436; } double f_id(double re, double im) { double r25437 = 0.5; double r25438 = re; double r25439 = cos(r25438); double r25440 = r25437 * r25439; double r25441 = 0; double r25442 = im; double r25443 = r25441 - r25442; double r25444 = exp(r25443); double r25445 = exp(r25442); double r25446 = r25444 - r25445; double r25447 = r25440 * r25446; return r25447; } double f_of(float re, float im) { float r25448 = 0.5; float r25449 = re; float r25450 = cos(r25449); float r25451 = r25448 * r25450; float r25452 = 1/60; float r25453 = im; float r25454 = 5; float r25455 = pow(r25453, r25454); float r25456 = r25452 * r25455; float r25457 = 2; float r25458 = r25457 * r25453; float r25459 = 1/3; float r25460 = 3; float r25461 = pow(r25453, r25460); float r25462 = r25459 * r25461; float r25463 = r25458 + r25462; float r25464 = r25456 + r25463; float r25465 = -r25464; float r25466 = r25451 * r25465; return r25466; } double f_od(double re, double im) { double r25467 = 0.5; double r25468 = re; double r25469 = cos(r25468); double r25470 = r25467 * r25469; double r25471 = 1/60; double r25472 = im; double r25473 = 5; double r25474 = pow(r25472, r25473); double r25475 = r25471 * r25474; double r25476 = 2; double r25477 = r25476 * r25472; double r25478 = 1/3; double r25479 = 3; double r25480 = pow(r25472, r25479); double r25481 = r25478 * r25480; double r25482 = r25477 + r25481; double r25483 = r25475 + r25482; double r25484 = -r25483; double r25485 = r25470 * r25484; return r25485; } 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 r25486, r25487, r25488, r25489, r25490, r25491, r25492, r25493, r25494, r25495, r25496; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25486, "0.5", 10, MPFR_RNDN); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init_set_str(r25490, "0", 10, MPFR_RNDN); mpfr_init(r25491); mpfr_init(r25492); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); } double f_im(double re, double im) { ; mpfr_set_d(r25487, re, MPFR_RNDN); mpfr_cos(r25488, r25487, MPFR_RNDN); mpfr_mul(r25489, r25486, r25488, MPFR_RNDN); ; mpfr_set_d(r25491, im, MPFR_RNDN); mpfr_sub(r25492, r25490, r25491, MPFR_RNDN); mpfr_exp(r25493, r25492, MPFR_RNDN); mpfr_exp(r25494, r25491, MPFR_RNDN); mpfr_sub(r25495, r25493, r25494, MPFR_RNDN); mpfr_mul(r25496, r25489, r25495, MPFR_RNDN); return mpfr_get_d(r25496, MPFR_RNDN); } static mpfr_t r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25497, "0.5", 10, MPFR_RNDN); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init_set_str(r25501, "1/60", 10, MPFR_RNDN); mpfr_init(r25502); mpfr_init_set_str(r25503, "5", 10, MPFR_RNDN); mpfr_init(r25504); mpfr_init(r25505); mpfr_init_set_str(r25506, "2", 10, MPFR_RNDN); mpfr_init(r25507); mpfr_init_set_str(r25508, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25509, "3", 10, MPFR_RNDN); mpfr_init(r25510); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); } double f_fm(double re, double im) { ; mpfr_set_d(r25498, re, MPFR_RNDN); mpfr_cos(r25499, r25498, MPFR_RNDN); mpfr_mul(r25500, r25497, r25499, MPFR_RNDN); ; mpfr_set_d(r25502, im, MPFR_RNDN); ; mpfr_pow(r25504, r25502, r25503, MPFR_RNDN); mpfr_mul(r25505, r25501, r25504, MPFR_RNDN); ; mpfr_mul(r25507, r25506, r25502, MPFR_RNDN); ; ; mpfr_pow(r25510, r25502, r25509, MPFR_RNDN); mpfr_mul(r25511, r25508, r25510, MPFR_RNDN); mpfr_add(r25512, r25507, r25511, MPFR_RNDN); mpfr_add(r25513, r25505, r25512, MPFR_RNDN); mpfr_neg(r25514, r25513, MPFR_RNDN); mpfr_mul(r25515, r25500, r25514, MPFR_RNDN); return mpfr_get_d(r25515, MPFR_RNDN); } static mpfr_t r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25516, "0.5", 10, MPFR_RNDN); mpfr_init(r25517); mpfr_init(r25518); mpfr_init(r25519); mpfr_init_set_str(r25520, "1/60", 10, MPFR_RNDN); mpfr_init(r25521); mpfr_init_set_str(r25522, "5", 10, MPFR_RNDN); mpfr_init(r25523); mpfr_init(r25524); mpfr_init_set_str(r25525, "2", 10, MPFR_RNDN); mpfr_init(r25526); mpfr_init_set_str(r25527, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25528, "3", 10, MPFR_RNDN); mpfr_init(r25529); mpfr_init(r25530); mpfr_init(r25531); mpfr_init(r25532); mpfr_init(r25533); mpfr_init(r25534); } double f_dm(double re, double im) { ; mpfr_set_d(r25517, re, MPFR_RNDN); mpfr_cos(r25518, r25517, MPFR_RNDN); mpfr_mul(r25519, r25516, r25518, MPFR_RNDN); ; mpfr_set_d(r25521, im, MPFR_RNDN); ; mpfr_pow(r25523, r25521, r25522, MPFR_RNDN); mpfr_mul(r25524, r25520, r25523, MPFR_RNDN); ; mpfr_mul(r25526, r25525, r25521, MPFR_RNDN); ; ; mpfr_pow(r25529, r25521, r25528, MPFR_RNDN); mpfr_mul(r25530, r25527, r25529, MPFR_RNDN); mpfr_add(r25531, r25526, r25530, MPFR_RNDN); mpfr_add(r25532, r25524, r25531, MPFR_RNDN); mpfr_neg(r25533, r25532, MPFR_RNDN); mpfr_mul(r25534, r25519, r25533, MPFR_RNDN); return mpfr_get_d(r25534, MPFR_RNDN); }