#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r18477 = 0.5f; float r18478 = re; float r18479 = sin(r18478); float r18480 = r18477 * r18479; float r18481 = im; float r18482 = -r18481; float r18483 = exp(r18482); float r18484 = exp(r18481); float r18485 = r18483 - r18484; float r18486 = r18480 * r18485; return r18486; } double f_id(double re, double im) { double r18487 = 0.5; double r18488 = re; double r18489 = sin(r18488); double r18490 = r18487 * r18489; double r18491 = im; double r18492 = -r18491; double r18493 = exp(r18492); double r18494 = exp(r18491); double r18495 = r18493 - r18494; double r18496 = r18490 * r18495; return r18496; } double f_of(float re, float im) { float r18497 = 0.5f; float r18498 = re; float r18499 = sin(r18498); float r18500 = r18497 * r18499; float r18501 = 0.01666666753590107f; float r18502 = im; float r18503 = 5.0f; float r18504 = pow(r18502, r18503); float r18505 = r18501 * r18504; float r18506 = 2.0f; float r18507 = r18506 * r18502; float r18508 = 0.3333333432674408f; float r18509 = 3.0f; float r18510 = pow(r18502, r18509); float r18511 = r18508 * r18510; float r18512 = r18507 + r18511; float r18513 = r18505 + r18512; float r18514 = -r18513; float r18515 = r18500 * r18514; return r18515; } double f_od(double re, double im) { double r18516 = 0.5; double r18517 = re; double r18518 = sin(r18517); double r18519 = r18516 * r18518; double r18520 = 0.01666666753590107; double r18521 = im; double r18522 = 5.0; double r18523 = pow(r18521, r18522); double r18524 = r18520 * r18523; double r18525 = 2.0; double r18526 = r18525 * r18521; double r18527 = 0.3333333432674408; double r18528 = 3.0; double r18529 = pow(r18521, r18528); double r18530 = r18527 * r18529; double r18531 = r18526 + r18530; double r18532 = r18524 + r18531; double r18533 = -r18532; double r18534 = r18519 * r18533; return r18534; } 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 r18535, r18536, r18537, r18538, r18539, r18540, r18541, r18542, r18543, r18544; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18535, "0.5", 10, MPFR_RNDN); mpfr_init(r18536); mpfr_init(r18537); mpfr_init(r18538); mpfr_init(r18539); mpfr_init(r18540); mpfr_init(r18541); mpfr_init(r18542); mpfr_init(r18543); mpfr_init(r18544); } double f_im(double re, double im) { ; mpfr_set_d(r18536, re, MPFR_RNDN); mpfr_sin(r18537, r18536, MPFR_RNDN); mpfr_mul(r18538, r18535, r18537, MPFR_RNDN); mpfr_set_d(r18539, im, MPFR_RNDN); mpfr_neg(r18540, r18539, MPFR_RNDN); mpfr_exp(r18541, r18540, MPFR_RNDN); mpfr_exp(r18542, r18539, MPFR_RNDN); mpfr_sub(r18543, r18541, r18542, MPFR_RNDN); mpfr_mul(r18544, r18538, r18543, MPFR_RNDN); return mpfr_get_d(r18544, MPFR_RNDN); } static mpfr_t r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18545, "0.5", 10, MPFR_RNDN); mpfr_init(r18546); mpfr_init(r18547); mpfr_init(r18548); mpfr_init_set_str(r18549, "1/60", 10, MPFR_RNDN); mpfr_init(r18550); mpfr_init_set_str(r18551, "5", 10, MPFR_RNDN); mpfr_init(r18552); mpfr_init(r18553); mpfr_init_set_str(r18554, "2", 10, MPFR_RNDN); mpfr_init(r18555); mpfr_init_set_str(r18556, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18557, "3", 10, MPFR_RNDN); mpfr_init(r18558); mpfr_init(r18559); mpfr_init(r18560); mpfr_init(r18561); mpfr_init(r18562); mpfr_init(r18563); } double f_fm(double re, double im) { ; mpfr_set_d(r18546, re, MPFR_RNDN); mpfr_sin(r18547, r18546, MPFR_RNDN); mpfr_mul(r18548, r18545, r18547, MPFR_RNDN); ; mpfr_set_d(r18550, im, MPFR_RNDN); ; mpfr_pow(r18552, r18550, r18551, MPFR_RNDN); mpfr_mul(r18553, r18549, r18552, MPFR_RNDN); ; mpfr_mul(r18555, r18554, r18550, MPFR_RNDN); ; ; mpfr_pow(r18558, r18550, r18557, MPFR_RNDN); mpfr_mul(r18559, r18556, r18558, MPFR_RNDN); mpfr_add(r18560, r18555, r18559, MPFR_RNDN); mpfr_add(r18561, r18553, r18560, MPFR_RNDN); mpfr_neg(r18562, r18561, MPFR_RNDN); mpfr_mul(r18563, r18548, r18562, MPFR_RNDN); return mpfr_get_d(r18563, MPFR_RNDN); } static mpfr_t r18564, r18565, r18566, r18567, r18568, r18569, r18570, r18571, r18572, r18573, r18574, r18575, r18576, r18577, r18578, r18579, r18580, r18581, r18582; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18564, "0.5", 10, MPFR_RNDN); mpfr_init(r18565); mpfr_init(r18566); mpfr_init(r18567); mpfr_init_set_str(r18568, "1/60", 10, MPFR_RNDN); mpfr_init(r18569); mpfr_init_set_str(r18570, "5", 10, MPFR_RNDN); mpfr_init(r18571); mpfr_init(r18572); mpfr_init_set_str(r18573, "2", 10, MPFR_RNDN); mpfr_init(r18574); mpfr_init_set_str(r18575, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18576, "3", 10, MPFR_RNDN); mpfr_init(r18577); mpfr_init(r18578); mpfr_init(r18579); mpfr_init(r18580); mpfr_init(r18581); mpfr_init(r18582); } double f_dm(double re, double im) { ; mpfr_set_d(r18565, re, MPFR_RNDN); mpfr_sin(r18566, r18565, MPFR_RNDN); mpfr_mul(r18567, r18564, r18566, MPFR_RNDN); ; mpfr_set_d(r18569, im, MPFR_RNDN); ; mpfr_pow(r18571, r18569, r18570, MPFR_RNDN); mpfr_mul(r18572, r18568, r18571, MPFR_RNDN); ; mpfr_mul(r18574, r18573, r18569, MPFR_RNDN); ; ; mpfr_pow(r18577, r18569, r18576, MPFR_RNDN); mpfr_mul(r18578, r18575, r18577, MPFR_RNDN); mpfr_add(r18579, r18574, r18578, MPFR_RNDN); mpfr_add(r18580, r18572, r18579, MPFR_RNDN); mpfr_neg(r18581, r18580, MPFR_RNDN); mpfr_mul(r18582, r18567, r18581, MPFR_RNDN); return mpfr_get_d(r18582, MPFR_RNDN); }