#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18439 = 0.5f; float r18440 = re; float r18441 = cos(r18440); float r18442 = r18439 * r18441; float r18443 = 0.0f; float r18444 = im; float r18445 = r18443 - r18444; float r18446 = exp(r18445); float r18447 = exp(r18444); float r18448 = r18446 - r18447; float r18449 = r18442 * r18448; return r18449; } double f_id(double re, double im) { double r18450 = 0.5; double r18451 = re; double r18452 = cos(r18451); double r18453 = r18450 * r18452; double r18454 = 0.0; double r18455 = im; double r18456 = r18454 - r18455; double r18457 = exp(r18456); double r18458 = exp(r18455); double r18459 = r18457 - r18458; double r18460 = r18453 * r18459; return r18460; } double f_of(float re, float im) { float r18461 = 0.01666666753590107f; float r18462 = im; float r18463 = 5.0f; float r18464 = pow(r18462, r18463); float r18465 = r18461 * r18464; float r18466 = 2.0f; float r18467 = r18466 * r18462; float r18468 = 0.3333333432674408f; float r18469 = 3.0f; float r18470 = pow(r18462, r18469); float r18471 = r18468 * r18470; float r18472 = r18467 + r18471; float r18473 = r18465 + r18472; float r18474 = -r18473; float r18475 = re; float r18476 = cos(r18475); float r18477 = 0.5f; float r18478 = r18476 * r18477; float r18479 = r18474 * r18478; return r18479; } double f_od(double re, double im) { double r18480 = 0.01666666753590107; double r18481 = im; double r18482 = 5.0; double r18483 = pow(r18481, r18482); double r18484 = r18480 * r18483; double r18485 = 2.0; double r18486 = r18485 * r18481; double r18487 = 0.3333333432674408; double r18488 = 3.0; double r18489 = pow(r18481, r18488); double r18490 = r18487 * r18489; double r18491 = r18486 + r18490; double r18492 = r18484 + r18491; double r18493 = -r18492; double r18494 = re; double r18495 = cos(r18494); double r18496 = 0.5; double r18497 = r18495 * r18496; double r18498 = r18493 * r18497; return r18498; } 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 r18499, r18500, r18501, r18502, r18503, r18504, r18505, r18506, r18507, r18508, r18509; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18499, "0.5", 10, MPFR_RNDN); mpfr_init(r18500); mpfr_init(r18501); mpfr_init(r18502); mpfr_init_set_str(r18503, "0", 10, MPFR_RNDN); mpfr_init(r18504); mpfr_init(r18505); mpfr_init(r18506); mpfr_init(r18507); mpfr_init(r18508); mpfr_init(r18509); } double f_im(double re, double im) { ; mpfr_set_d(r18500, re, MPFR_RNDN); mpfr_cos(r18501, r18500, MPFR_RNDN); mpfr_mul(r18502, r18499, r18501, MPFR_RNDN); ; mpfr_set_d(r18504, im, MPFR_RNDN); mpfr_sub(r18505, r18503, r18504, MPFR_RNDN); mpfr_exp(r18506, r18505, MPFR_RNDN); mpfr_exp(r18507, r18504, MPFR_RNDN); mpfr_sub(r18508, r18506, r18507, MPFR_RNDN); mpfr_mul(r18509, r18502, r18508, MPFR_RNDN); return mpfr_get_d(r18509, MPFR_RNDN); } static mpfr_t r18510, r18511, r18512, r18513, r18514, r18515, r18516, r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525, r18526, r18527, r18528; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18510, "1/60", 10, MPFR_RNDN); mpfr_init(r18511); mpfr_init_set_str(r18512, "5", 10, MPFR_RNDN); mpfr_init(r18513); mpfr_init(r18514); mpfr_init_set_str(r18515, "2", 10, MPFR_RNDN); mpfr_init(r18516); mpfr_init_set_str(r18517, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18518, "3", 10, MPFR_RNDN); mpfr_init(r18519); mpfr_init(r18520); mpfr_init(r18521); mpfr_init(r18522); mpfr_init(r18523); mpfr_init(r18524); mpfr_init(r18525); mpfr_init_set_str(r18526, "0.5", 10, MPFR_RNDN); mpfr_init(r18527); mpfr_init(r18528); } double f_fm(double re, double im) { ; mpfr_set_d(r18511, im, MPFR_RNDN); ; mpfr_pow(r18513, r18511, r18512, MPFR_RNDN); mpfr_mul(r18514, r18510, r18513, MPFR_RNDN); ; mpfr_mul(r18516, r18515, r18511, MPFR_RNDN); ; ; mpfr_pow(r18519, r18511, r18518, MPFR_RNDN); mpfr_mul(r18520, r18517, r18519, MPFR_RNDN); mpfr_add(r18521, r18516, r18520, MPFR_RNDN); mpfr_add(r18522, r18514, r18521, MPFR_RNDN); mpfr_neg(r18523, r18522, MPFR_RNDN); mpfr_set_d(r18524, re, MPFR_RNDN); mpfr_cos(r18525, r18524, MPFR_RNDN); ; mpfr_mul(r18527, r18525, r18526, MPFR_RNDN); mpfr_mul(r18528, r18523, r18527, MPFR_RNDN); return mpfr_get_d(r18528, MPFR_RNDN); } static mpfr_t r18529, r18530, r18531, r18532, r18533, r18534, r18535, r18536, r18537, r18538, r18539, r18540, r18541, r18542, r18543, r18544, r18545, r18546, r18547; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18529, "1/60", 10, MPFR_RNDN); mpfr_init(r18530); mpfr_init_set_str(r18531, "5", 10, MPFR_RNDN); mpfr_init(r18532); mpfr_init(r18533); mpfr_init_set_str(r18534, "2", 10, MPFR_RNDN); mpfr_init(r18535); mpfr_init_set_str(r18536, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18537, "3", 10, MPFR_RNDN); mpfr_init(r18538); mpfr_init(r18539); mpfr_init(r18540); mpfr_init(r18541); mpfr_init(r18542); mpfr_init(r18543); mpfr_init(r18544); mpfr_init_set_str(r18545, "0.5", 10, MPFR_RNDN); mpfr_init(r18546); mpfr_init(r18547); } double f_dm(double re, double im) { ; mpfr_set_d(r18530, im, MPFR_RNDN); ; mpfr_pow(r18532, r18530, r18531, MPFR_RNDN); mpfr_mul(r18533, r18529, r18532, MPFR_RNDN); ; mpfr_mul(r18535, r18534, r18530, MPFR_RNDN); ; ; mpfr_pow(r18538, r18530, r18537, MPFR_RNDN); mpfr_mul(r18539, r18536, r18538, MPFR_RNDN); mpfr_add(r18540, r18535, r18539, MPFR_RNDN); mpfr_add(r18541, r18533, r18540, MPFR_RNDN); mpfr_neg(r18542, r18541, MPFR_RNDN); mpfr_set_d(r18543, re, MPFR_RNDN); mpfr_cos(r18544, r18543, MPFR_RNDN); ; mpfr_mul(r18546, r18544, r18545, MPFR_RNDN); mpfr_mul(r18547, r18542, r18546, MPFR_RNDN); return mpfr_get_d(r18547, MPFR_RNDN); }