#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19456 = 0.5f; float r19457 = re; float r19458 = cos(r19457); float r19459 = r19456 * r19458; float r19460 = 0.0f; float r19461 = im; float r19462 = r19460 - r19461; float r19463 = exp(r19462); float r19464 = exp(r19461); float r19465 = r19463 - r19464; float r19466 = r19459 * r19465; return r19466; } double f_id(double re, double im) { double r19467 = 0.5; double r19468 = re; double r19469 = cos(r19468); double r19470 = r19467 * r19469; double r19471 = 0.0; double r19472 = im; double r19473 = r19471 - r19472; double r19474 = exp(r19473); double r19475 = exp(r19472); double r19476 = r19474 - r19475; double r19477 = r19470 * r19476; return r19477; } double f_of(float re, float im) { float r19478 = 0.01666666753590107f; float r19479 = im; float r19480 = 5.0f; float r19481 = pow(r19479, r19480); float r19482 = r19478 * r19481; float r19483 = 2.0f; float r19484 = r19483 * r19479; float r19485 = 0.3333333432674408f; float r19486 = 3.0f; float r19487 = pow(r19479, r19486); float r19488 = r19485 * r19487; float r19489 = r19484 + r19488; float r19490 = r19482 + r19489; float r19491 = -r19490; float r19492 = re; float r19493 = cos(r19492); float r19494 = 0.5f; float r19495 = r19493 * r19494; float r19496 = r19491 * r19495; return r19496; } double f_od(double re, double im) { double r19497 = 0.01666666753590107; double r19498 = im; double r19499 = 5.0; double r19500 = pow(r19498, r19499); double r19501 = r19497 * r19500; double r19502 = 2.0; double r19503 = r19502 * r19498; double r19504 = 0.3333333432674408; double r19505 = 3.0; double r19506 = pow(r19498, r19505); double r19507 = r19504 * r19506; double r19508 = r19503 + r19507; double r19509 = r19501 + r19508; double r19510 = -r19509; double r19511 = re; double r19512 = cos(r19511); double r19513 = 0.5; double r19514 = r19512 * r19513; double r19515 = r19510 * r19514; return r19515; } 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 r19516, r19517, r19518, r19519, r19520, r19521, r19522, r19523, r19524, r19525, r19526; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19516, "0.5", 10, MPFR_RNDN); mpfr_init(r19517); mpfr_init(r19518); mpfr_init(r19519); mpfr_init_set_str(r19520, "0", 10, MPFR_RNDN); mpfr_init(r19521); mpfr_init(r19522); mpfr_init(r19523); mpfr_init(r19524); mpfr_init(r19525); mpfr_init(r19526); } double f_im(double re, double im) { ; mpfr_set_d(r19517, re, MPFR_RNDN); mpfr_cos(r19518, r19517, MPFR_RNDN); mpfr_mul(r19519, r19516, r19518, MPFR_RNDN); ; mpfr_set_d(r19521, im, MPFR_RNDN); mpfr_sub(r19522, r19520, r19521, MPFR_RNDN); mpfr_exp(r19523, r19522, MPFR_RNDN); mpfr_exp(r19524, r19521, MPFR_RNDN); mpfr_sub(r19525, r19523, r19524, MPFR_RNDN); mpfr_mul(r19526, r19519, r19525, MPFR_RNDN); return mpfr_get_d(r19526, MPFR_RNDN); } static mpfr_t r19527, r19528, r19529, r19530, r19531, r19532, r19533, r19534, r19535, r19536, r19537, r19538, r19539, r19540, r19541, r19542, r19543, r19544, r19545; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19527, "1/60", 10, MPFR_RNDN); mpfr_init(r19528); mpfr_init_set_str(r19529, "5", 10, MPFR_RNDN); mpfr_init(r19530); mpfr_init(r19531); mpfr_init_set_str(r19532, "2", 10, MPFR_RNDN); mpfr_init(r19533); mpfr_init_set_str(r19534, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19535, "3", 10, MPFR_RNDN); mpfr_init(r19536); mpfr_init(r19537); mpfr_init(r19538); mpfr_init(r19539); mpfr_init(r19540); mpfr_init(r19541); mpfr_init(r19542); mpfr_init_set_str(r19543, "0.5", 10, MPFR_RNDN); mpfr_init(r19544); mpfr_init(r19545); } double f_fm(double re, double im) { ; mpfr_set_d(r19528, im, MPFR_RNDN); ; mpfr_pow(r19530, r19528, r19529, MPFR_RNDN); mpfr_mul(r19531, r19527, r19530, MPFR_RNDN); ; mpfr_mul(r19533, r19532, r19528, MPFR_RNDN); ; ; mpfr_pow(r19536, r19528, r19535, MPFR_RNDN); mpfr_mul(r19537, r19534, r19536, MPFR_RNDN); mpfr_add(r19538, r19533, r19537, MPFR_RNDN); mpfr_add(r19539, r19531, r19538, MPFR_RNDN); mpfr_neg(r19540, r19539, MPFR_RNDN); mpfr_set_d(r19541, re, MPFR_RNDN); mpfr_cos(r19542, r19541, MPFR_RNDN); ; mpfr_mul(r19544, r19542, r19543, MPFR_RNDN); mpfr_mul(r19545, r19540, r19544, MPFR_RNDN); return mpfr_get_d(r19545, MPFR_RNDN); } static mpfr_t r19546, r19547, r19548, r19549, r19550, r19551, r19552, r19553, r19554, r19555, r19556, r19557, r19558, r19559, r19560, r19561, r19562, r19563, r19564; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19546, "1/60", 10, MPFR_RNDN); mpfr_init(r19547); mpfr_init_set_str(r19548, "5", 10, MPFR_RNDN); mpfr_init(r19549); mpfr_init(r19550); mpfr_init_set_str(r19551, "2", 10, MPFR_RNDN); mpfr_init(r19552); mpfr_init_set_str(r19553, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19554, "3", 10, MPFR_RNDN); mpfr_init(r19555); mpfr_init(r19556); mpfr_init(r19557); mpfr_init(r19558); mpfr_init(r19559); mpfr_init(r19560); mpfr_init(r19561); mpfr_init_set_str(r19562, "0.5", 10, MPFR_RNDN); mpfr_init(r19563); mpfr_init(r19564); } double f_dm(double re, double im) { ; mpfr_set_d(r19547, im, MPFR_RNDN); ; mpfr_pow(r19549, r19547, r19548, MPFR_RNDN); mpfr_mul(r19550, r19546, r19549, MPFR_RNDN); ; mpfr_mul(r19552, r19551, r19547, MPFR_RNDN); ; ; mpfr_pow(r19555, r19547, r19554, MPFR_RNDN); mpfr_mul(r19556, r19553, r19555, MPFR_RNDN); mpfr_add(r19557, r19552, r19556, MPFR_RNDN); mpfr_add(r19558, r19550, r19557, MPFR_RNDN); mpfr_neg(r19559, r19558, MPFR_RNDN); mpfr_set_d(r19560, re, MPFR_RNDN); mpfr_cos(r19561, r19560, MPFR_RNDN); ; mpfr_mul(r19563, r19561, r19562, MPFR_RNDN); mpfr_mul(r19564, r19559, r19563, MPFR_RNDN); return mpfr_get_d(r19564, MPFR_RNDN); }