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