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