#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18395 = 0.5f; float r18396 = re; float r18397 = cos(r18396); float r18398 = r18395 * r18397; float r18399 = 0.0f; float r18400 = im; float r18401 = r18399 - r18400; float r18402 = exp(r18401); float r18403 = exp(r18400); float r18404 = r18402 - r18403; float r18405 = r18398 * r18404; return r18405; } double f_id(double re, double im) { double r18406 = 0.5; double r18407 = re; double r18408 = cos(r18407); double r18409 = r18406 * r18408; double r18410 = 0.0; double r18411 = im; double r18412 = r18410 - r18411; double r18413 = exp(r18412); double r18414 = exp(r18411); double r18415 = r18413 - r18414; double r18416 = r18409 * r18415; return r18416; } double f_of(float re, float im) { float r18417 = im; float r18418 = -0.12694145739078522f; bool r18419 = r18417 <= r18418; float r18420 = -r18417; float r18421 = exp(r18420); float r18422 = exp(r18417); float r18423 = r18421 - r18422; float r18424 = log(r18423); float r18425 = exp(r18424); float r18426 = re; float r18427 = cos(r18426); float r18428 = 0.5f; float r18429 = r18427 * r18428; float r18430 = r18425 * r18429; float r18431 = 0.01666666753590107f; float r18432 = 5.0f; float r18433 = pow(r18417, r18432); float r18434 = r18431 * r18433; float r18435 = 2.0f; float r18436 = r18435 * r18417; float r18437 = 0.3333333432674408f; float r18438 = 3.0f; float r18439 = pow(r18417, r18438); float r18440 = r18437 * r18439; float r18441 = r18436 + r18440; float r18442 = r18434 + r18441; float r18443 = -r18442; float r18444 = r18443 * r18429; float r18445 = r18419 ? r18430 : r18444; return r18445; } double f_od(double re, double im) { double r18446 = im; double r18447 = -0.12694145739078522; bool r18448 = r18446 <= r18447; double r18449 = -r18446; double r18450 = exp(r18449); double r18451 = exp(r18446); double r18452 = r18450 - r18451; double r18453 = log(r18452); double r18454 = exp(r18453); double r18455 = re; double r18456 = cos(r18455); double r18457 = 0.5; double r18458 = r18456 * r18457; double r18459 = r18454 * r18458; double r18460 = 0.01666666753590107; double r18461 = 5.0; double r18462 = pow(r18446, r18461); double r18463 = r18460 * r18462; double r18464 = 2.0; double r18465 = r18464 * r18446; double r18466 = 0.3333333432674408; double r18467 = 3.0; double r18468 = pow(r18446, r18467); double r18469 = r18466 * r18468; double r18470 = r18465 + r18469; double r18471 = r18463 + r18470; double r18472 = -r18471; double r18473 = r18472 * r18458; double r18474 = r18448 ? r18459 : r18473; return r18474; } 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 r18475, r18476, r18477, r18478, r18479, r18480, r18481, r18482, r18483, r18484, r18485; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18475, "0.5", 10, MPFR_RNDN); mpfr_init(r18476); mpfr_init(r18477); mpfr_init(r18478); mpfr_init_set_str(r18479, "0", 10, MPFR_RNDN); mpfr_init(r18480); mpfr_init(r18481); mpfr_init(r18482); mpfr_init(r18483); mpfr_init(r18484); mpfr_init(r18485); } double f_im(double re, double im) { ; mpfr_set_d(r18476, re, MPFR_RNDN); mpfr_cos(r18477, r18476, MPFR_RNDN); mpfr_mul(r18478, r18475, r18477, MPFR_RNDN); ; mpfr_set_d(r18480, im, MPFR_RNDN); mpfr_sub(r18481, r18479, r18480, MPFR_RNDN); mpfr_exp(r18482, r18481, MPFR_RNDN); mpfr_exp(r18483, r18480, MPFR_RNDN); mpfr_sub(r18484, r18482, r18483, MPFR_RNDN); mpfr_mul(r18485, r18478, r18484, MPFR_RNDN); return mpfr_get_d(r18485, MPFR_RNDN); } static mpfr_t r18486, r18487, r18488, r18489, r18490, r18491, r18492, r18493, r18494, r18495, r18496, r18497, r18498, r18499, r18500, r18501, r18502, r18503, r18504, r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513, r18514; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18486); mpfr_init_set_str(r18487, "-0.12694146f0", 10, MPFR_RNDN); mpfr_init(r18488); mpfr_init(r18489); mpfr_init(r18490); mpfr_init(r18491); mpfr_init(r18492); mpfr_init(r18493); mpfr_init(r18494); mpfr_init(r18495); mpfr_init(r18496); mpfr_init_set_str(r18497, "0.5", 10, MPFR_RNDN); mpfr_init(r18498); mpfr_init(r18499); mpfr_init_set_str(r18500, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18501, "5", 10, MPFR_RNDN); mpfr_init(r18502); mpfr_init(r18503); mpfr_init_set_str(r18504, "2", 10, MPFR_RNDN); mpfr_init(r18505); mpfr_init_set_str(r18506, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18507, "3", 10, MPFR_RNDN); mpfr_init(r18508); mpfr_init(r18509); mpfr_init(r18510); mpfr_init(r18511); mpfr_init(r18512); mpfr_init(r18513); mpfr_init(r18514); } double f_fm(double re, double im) { mpfr_set_d(r18486, im, MPFR_RNDN); ; mpfr_set_si(r18488, mpfr_cmp(r18486, r18487) <= 0, MPFR_RNDN); mpfr_neg(r18489, r18486, MPFR_RNDN); mpfr_exp(r18490, r18489, MPFR_RNDN); mpfr_exp(r18491, r18486, MPFR_RNDN); mpfr_sub(r18492, r18490, r18491, MPFR_RNDN); mpfr_log(r18493, r18492, MPFR_RNDN); mpfr_exp(r18494, r18493, MPFR_RNDN); mpfr_set_d(r18495, re, MPFR_RNDN); mpfr_cos(r18496, r18495, MPFR_RNDN); ; mpfr_mul(r18498, r18496, r18497, MPFR_RNDN); mpfr_mul(r18499, r18494, r18498, MPFR_RNDN); ; ; mpfr_pow(r18502, r18486, r18501, MPFR_RNDN); mpfr_mul(r18503, r18500, r18502, MPFR_RNDN); ; mpfr_mul(r18505, r18504, r18486, MPFR_RNDN); ; ; mpfr_pow(r18508, r18486, r18507, MPFR_RNDN); mpfr_mul(r18509, r18506, r18508, MPFR_RNDN); mpfr_add(r18510, r18505, r18509, MPFR_RNDN); mpfr_add(r18511, r18503, r18510, MPFR_RNDN); mpfr_neg(r18512, r18511, MPFR_RNDN); mpfr_mul(r18513, r18512, r18498, MPFR_RNDN); if (mpfr_get_si(r18488, MPFR_RNDN)) { mpfr_set(r18514, r18499, MPFR_RNDN); } else { mpfr_set(r18514, r18513, MPFR_RNDN); }; return mpfr_get_d(r18514, MPFR_RNDN); } static mpfr_t r18515, r18516, r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, 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(r18515); mpfr_init_set_str(r18516, "-0.12694146f0", 10, MPFR_RNDN); mpfr_init(r18517); mpfr_init(r18518); 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); mpfr_init_set_str(r18529, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18530, "5", 10, MPFR_RNDN); mpfr_init(r18531); mpfr_init(r18532); mpfr_init_set_str(r18533, "2", 10, MPFR_RNDN); mpfr_init(r18534); mpfr_init_set_str(r18535, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18536, "3", 10, MPFR_RNDN); mpfr_init(r18537); mpfr_init(r18538); mpfr_init(r18539); mpfr_init(r18540); mpfr_init(r18541); mpfr_init(r18542); mpfr_init(r18543); } double f_dm(double re, double im) { mpfr_set_d(r18515, im, MPFR_RNDN); ; mpfr_set_si(r18517, mpfr_cmp(r18515, r18516) <= 0, MPFR_RNDN); mpfr_neg(r18518, r18515, MPFR_RNDN); mpfr_exp(r18519, r18518, MPFR_RNDN); mpfr_exp(r18520, r18515, MPFR_RNDN); mpfr_sub(r18521, r18519, r18520, MPFR_RNDN); mpfr_log(r18522, r18521, MPFR_RNDN); mpfr_exp(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); ; ; mpfr_pow(r18531, r18515, r18530, MPFR_RNDN); mpfr_mul(r18532, r18529, r18531, MPFR_RNDN); ; mpfr_mul(r18534, r18533, r18515, MPFR_RNDN); ; ; mpfr_pow(r18537, r18515, r18536, MPFR_RNDN); mpfr_mul(r18538, r18535, r18537, MPFR_RNDN); mpfr_add(r18539, r18534, r18538, MPFR_RNDN); mpfr_add(r18540, r18532, r18539, MPFR_RNDN); mpfr_neg(r18541, r18540, MPFR_RNDN); mpfr_mul(r18542, r18541, r18527, MPFR_RNDN); if (mpfr_get_si(r18517, MPFR_RNDN)) { mpfr_set(r18543, r18528, MPFR_RNDN); } else { mpfr_set(r18543, r18542, MPFR_RNDN); }; return mpfr_get_d(r18543, MPFR_RNDN); }