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