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