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