#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r18352 = 0.5f; float r18353 = re; float r18354 = sin(r18353); float r18355 = r18352 * r18354; float r18356 = im; float r18357 = -r18356; float r18358 = exp(r18357); float r18359 = exp(r18356); float r18360 = r18358 - r18359; float r18361 = r18355 * r18360; return r18361; } double f_id(double re, double im) { double r18362 = 0.5; double r18363 = re; double r18364 = sin(r18363); double r18365 = r18362 * r18364; double r18366 = im; double r18367 = -r18366; double r18368 = exp(r18367); double r18369 = exp(r18366); double r18370 = r18368 - r18369; double r18371 = r18365 * r18370; return r18371; } double f_of(float re, float im) { float r18372 = 0.5f; float r18373 = re; float r18374 = sin(r18373); float r18375 = r18372 * r18374; float r18376 = 0.01666666753590107f; float r18377 = im; float r18378 = 5.0f; float r18379 = pow(r18377, r18378); float r18380 = r18376 * r18379; float r18381 = 2.0f; float r18382 = r18381 * r18377; float r18383 = 0.3333333432674408f; float r18384 = 3.0f; float r18385 = pow(r18377, r18384); float r18386 = r18383 * r18385; float r18387 = r18382 + r18386; float r18388 = r18380 + r18387; float r18389 = -r18388; float r18390 = r18375 * r18389; return r18390; } double f_od(double re, double im) { double r18391 = 0.5; double r18392 = re; double r18393 = sin(r18392); double r18394 = r18391 * r18393; double r18395 = 0.01666666753590107; double r18396 = im; double r18397 = 5.0; double r18398 = pow(r18396, r18397); double r18399 = r18395 * r18398; double r18400 = 2.0; double r18401 = r18400 * r18396; double r18402 = 0.3333333432674408; double r18403 = 3.0; double r18404 = pow(r18396, r18403); double r18405 = r18402 * r18404; double r18406 = r18401 + r18405; double r18407 = r18399 + r18406; double r18408 = -r18407; double r18409 = r18394 * r18408; return r18409; } 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 r18410, r18411, r18412, r18413, r18414, r18415, r18416, r18417, r18418, r18419; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18410, "0.5", 10, MPFR_RNDN); mpfr_init(r18411); mpfr_init(r18412); mpfr_init(r18413); mpfr_init(r18414); mpfr_init(r18415); mpfr_init(r18416); mpfr_init(r18417); mpfr_init(r18418); mpfr_init(r18419); } double f_im(double re, double im) { ; mpfr_set_d(r18411, re, MPFR_RNDN); mpfr_sin(r18412, r18411, MPFR_RNDN); mpfr_mul(r18413, r18410, r18412, MPFR_RNDN); mpfr_set_d(r18414, im, MPFR_RNDN); mpfr_neg(r18415, r18414, MPFR_RNDN); mpfr_exp(r18416, r18415, MPFR_RNDN); mpfr_exp(r18417, r18414, MPFR_RNDN); mpfr_sub(r18418, r18416, r18417, MPFR_RNDN); mpfr_mul(r18419, r18413, r18418, MPFR_RNDN); return mpfr_get_d(r18419, MPFR_RNDN); } static mpfr_t r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432, r18433, r18434, r18435, r18436, r18437, r18438; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18420, "0.5", 10, MPFR_RNDN); mpfr_init(r18421); mpfr_init(r18422); mpfr_init(r18423); mpfr_init_set_str(r18424, "1/60", 10, MPFR_RNDN); mpfr_init(r18425); mpfr_init_set_str(r18426, "5", 10, MPFR_RNDN); mpfr_init(r18427); mpfr_init(r18428); mpfr_init_set_str(r18429, "2", 10, MPFR_RNDN); mpfr_init(r18430); mpfr_init_set_str(r18431, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18432, "3", 10, MPFR_RNDN); mpfr_init(r18433); mpfr_init(r18434); mpfr_init(r18435); mpfr_init(r18436); mpfr_init(r18437); mpfr_init(r18438); } double f_fm(double re, double im) { ; mpfr_set_d(r18421, re, MPFR_RNDN); mpfr_sin(r18422, r18421, MPFR_RNDN); mpfr_mul(r18423, r18420, r18422, MPFR_RNDN); ; mpfr_set_d(r18425, im, MPFR_RNDN); ; mpfr_pow(r18427, r18425, r18426, MPFR_RNDN); mpfr_mul(r18428, r18424, r18427, MPFR_RNDN); ; mpfr_mul(r18430, r18429, r18425, MPFR_RNDN); ; ; mpfr_pow(r18433, r18425, r18432, MPFR_RNDN); mpfr_mul(r18434, r18431, r18433, MPFR_RNDN); mpfr_add(r18435, r18430, r18434, MPFR_RNDN); mpfr_add(r18436, r18428, r18435, MPFR_RNDN); mpfr_neg(r18437, r18436, MPFR_RNDN); mpfr_mul(r18438, r18423, r18437, MPFR_RNDN); return mpfr_get_d(r18438, MPFR_RNDN); } static mpfr_t r18439, r18440, r18441, r18442, r18443, r18444, r18445, r18446, r18447, r18448, r18449, r18450, r18451, r18452, r18453, r18454, r18455, r18456, r18457; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18439, "0.5", 10, MPFR_RNDN); mpfr_init(r18440); mpfr_init(r18441); mpfr_init(r18442); mpfr_init_set_str(r18443, "1/60", 10, MPFR_RNDN); mpfr_init(r18444); mpfr_init_set_str(r18445, "5", 10, MPFR_RNDN); mpfr_init(r18446); mpfr_init(r18447); mpfr_init_set_str(r18448, "2", 10, MPFR_RNDN); mpfr_init(r18449); mpfr_init_set_str(r18450, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18451, "3", 10, MPFR_RNDN); mpfr_init(r18452); mpfr_init(r18453); mpfr_init(r18454); mpfr_init(r18455); mpfr_init(r18456); mpfr_init(r18457); } double f_dm(double re, double im) { ; mpfr_set_d(r18440, re, MPFR_RNDN); mpfr_sin(r18441, r18440, MPFR_RNDN); mpfr_mul(r18442, r18439, r18441, MPFR_RNDN); ; mpfr_set_d(r18444, im, MPFR_RNDN); ; mpfr_pow(r18446, r18444, r18445, MPFR_RNDN); mpfr_mul(r18447, r18443, r18446, MPFR_RNDN); ; mpfr_mul(r18449, r18448, r18444, MPFR_RNDN); ; ; mpfr_pow(r18452, r18444, r18451, MPFR_RNDN); mpfr_mul(r18453, r18450, r18452, MPFR_RNDN); mpfr_add(r18454, r18449, r18453, MPFR_RNDN); mpfr_add(r18455, r18447, r18454, MPFR_RNDN); mpfr_neg(r18456, r18455, MPFR_RNDN); mpfr_mul(r18457, r18442, r18456, MPFR_RNDN); return mpfr_get_d(r18457, MPFR_RNDN); }