#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25344 = 0.5; float r25345 = re; float r25346 = cos(r25345); float r25347 = r25344 * r25346; float r25348 = 0; float r25349 = im; float r25350 = r25348 - r25349; float r25351 = exp(r25350); float r25352 = exp(r25349); float r25353 = r25351 - r25352; float r25354 = r25347 * r25353; return r25354; } double f_id(double re, double im) { double r25355 = 0.5; double r25356 = re; double r25357 = cos(r25356); double r25358 = r25355 * r25357; double r25359 = 0; double r25360 = im; double r25361 = r25359 - r25360; double r25362 = exp(r25361); double r25363 = exp(r25360); double r25364 = r25362 - r25363; double r25365 = r25358 * r25364; return r25365; } double f_of(float re, float im) { float r25366 = 0.5; float r25367 = re; float r25368 = cos(r25367); float r25369 = r25366 * r25368; float r25370 = 1/60; float r25371 = im; float r25372 = 5; float r25373 = pow(r25371, r25372); float r25374 = r25370 * r25373; float r25375 = 2; float r25376 = r25375 * r25371; float r25377 = 1/3; float r25378 = 3; float r25379 = pow(r25371, r25378); float r25380 = r25377 * r25379; float r25381 = exp(r25380); float r25382 = log(r25381); float r25383 = r25376 + r25382; float r25384 = r25374 + r25383; float r25385 = -r25384; float r25386 = r25369 * r25385; return r25386; } double f_od(double re, double im) { double r25387 = 0.5; double r25388 = re; double r25389 = cos(r25388); double r25390 = r25387 * r25389; double r25391 = 1/60; double r25392 = im; double r25393 = 5; double r25394 = pow(r25392, r25393); double r25395 = r25391 * r25394; double r25396 = 2; double r25397 = r25396 * r25392; double r25398 = 1/3; double r25399 = 3; double r25400 = pow(r25392, r25399); double r25401 = r25398 * r25400; double r25402 = exp(r25401); double r25403 = log(r25402); double r25404 = r25397 + r25403; double r25405 = r25395 + r25404; double r25406 = -r25405; double r25407 = r25390 * r25406; return r25407; } 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 r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25408, "0.5", 10, MPFR_RNDN); mpfr_init(r25409); mpfr_init(r25410); mpfr_init(r25411); mpfr_init_set_str(r25412, "0", 10, MPFR_RNDN); mpfr_init(r25413); mpfr_init(r25414); mpfr_init(r25415); mpfr_init(r25416); mpfr_init(r25417); mpfr_init(r25418); } double f_im(double re, double im) { ; mpfr_set_d(r25409, re, MPFR_RNDN); mpfr_cos(r25410, r25409, MPFR_RNDN); mpfr_mul(r25411, r25408, r25410, MPFR_RNDN); ; mpfr_set_d(r25413, im, MPFR_RNDN); mpfr_sub(r25414, r25412, r25413, MPFR_RNDN); mpfr_exp(r25415, r25414, MPFR_RNDN); mpfr_exp(r25416, r25413, MPFR_RNDN); mpfr_sub(r25417, r25415, r25416, MPFR_RNDN); mpfr_mul(r25418, r25411, r25417, MPFR_RNDN); return mpfr_get_d(r25418, MPFR_RNDN); } static mpfr_t r25419, r25420, r25421, r25422, r25423, r25424, r25425, r25426, r25427, r25428, r25429, r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25419, "0.5", 10, MPFR_RNDN); mpfr_init(r25420); mpfr_init(r25421); mpfr_init(r25422); mpfr_init_set_str(r25423, "1/60", 10, MPFR_RNDN); mpfr_init(r25424); mpfr_init_set_str(r25425, "5", 10, MPFR_RNDN); mpfr_init(r25426); mpfr_init(r25427); mpfr_init_set_str(r25428, "2", 10, MPFR_RNDN); mpfr_init(r25429); mpfr_init_set_str(r25430, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25431, "3", 10, MPFR_RNDN); mpfr_init(r25432); mpfr_init(r25433); mpfr_init(r25434); mpfr_init(r25435); mpfr_init(r25436); mpfr_init(r25437); mpfr_init(r25438); mpfr_init(r25439); } double f_fm(double re, double im) { ; mpfr_set_d(r25420, re, MPFR_RNDN); mpfr_cos(r25421, r25420, MPFR_RNDN); mpfr_mul(r25422, r25419, r25421, MPFR_RNDN); ; mpfr_set_d(r25424, im, MPFR_RNDN); ; mpfr_pow(r25426, r25424, r25425, MPFR_RNDN); mpfr_mul(r25427, r25423, r25426, MPFR_RNDN); ; mpfr_mul(r25429, r25428, r25424, MPFR_RNDN); ; ; mpfr_pow(r25432, r25424, r25431, MPFR_RNDN); mpfr_mul(r25433, r25430, r25432, MPFR_RNDN); mpfr_exp(r25434, r25433, MPFR_RNDN); mpfr_log(r25435, r25434, MPFR_RNDN); mpfr_add(r25436, r25429, r25435, MPFR_RNDN); mpfr_add(r25437, r25427, r25436, MPFR_RNDN); mpfr_neg(r25438, r25437, MPFR_RNDN); mpfr_mul(r25439, r25422, r25438, MPFR_RNDN); return mpfr_get_d(r25439, MPFR_RNDN); } static mpfr_t r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25440, "0.5", 10, MPFR_RNDN); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init_set_str(r25444, "1/60", 10, MPFR_RNDN); mpfr_init(r25445); mpfr_init_set_str(r25446, "5", 10, MPFR_RNDN); mpfr_init(r25447); mpfr_init(r25448); mpfr_init_set_str(r25449, "2", 10, MPFR_RNDN); mpfr_init(r25450); mpfr_init_set_str(r25451, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25452, "3", 10, MPFR_RNDN); mpfr_init(r25453); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init(r25457); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); } double f_dm(double re, double im) { ; mpfr_set_d(r25441, re, MPFR_RNDN); mpfr_cos(r25442, r25441, MPFR_RNDN); mpfr_mul(r25443, r25440, r25442, MPFR_RNDN); ; mpfr_set_d(r25445, im, MPFR_RNDN); ; mpfr_pow(r25447, r25445, r25446, MPFR_RNDN); mpfr_mul(r25448, r25444, r25447, MPFR_RNDN); ; mpfr_mul(r25450, r25449, r25445, MPFR_RNDN); ; ; mpfr_pow(r25453, r25445, r25452, MPFR_RNDN); mpfr_mul(r25454, r25451, r25453, MPFR_RNDN); mpfr_exp(r25455, r25454, MPFR_RNDN); mpfr_log(r25456, r25455, MPFR_RNDN); mpfr_add(r25457, r25450, r25456, MPFR_RNDN); mpfr_add(r25458, r25448, r25457, MPFR_RNDN); mpfr_neg(r25459, r25458, MPFR_RNDN); mpfr_mul(r25460, r25443, r25459, MPFR_RNDN); return mpfr_get_d(r25460, MPFR_RNDN); }