#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18595 = 0.5f; float r18596 = re; float r18597 = cos(r18596); float r18598 = r18595 * r18597; float r18599 = 0.0f; float r18600 = im; float r18601 = r18599 - r18600; float r18602 = exp(r18601); float r18603 = exp(r18600); float r18604 = r18602 - r18603; float r18605 = r18598 * r18604; return r18605; } double f_id(double re, double im) { double r18606 = 0.5; double r18607 = re; double r18608 = cos(r18607); double r18609 = r18606 * r18608; double r18610 = 0.0; double r18611 = im; double r18612 = r18610 - r18611; double r18613 = exp(r18612); double r18614 = exp(r18611); double r18615 = r18613 - r18614; double r18616 = r18609 * r18615; return r18616; } double f_of(float re, float im) { float r18617 = 0.016666666666666666f; float r18618 = im; float r18619 = 5.0f; float r18620 = pow(r18618, r18619); float r18621 = r18618 * (r18618 * r18618); float r18622 = 0.3333333333333333f; float r18623 = 2.0f; float r18624 = r18623 * r18618; float r18625 = fma(r18621, r18622, r18624); float r18626 = fma(r18617, r18620, r18625); float r18627 = 1.0f; float r18628 = r18626 / r18627; float r18629 = re; float r18630 = cos(r18629); float r18631 = 0.5f; float r18632 = -r18631; float r18633 = r18630 * r18632; float r18634 = r18628 * r18633; return r18634; } double f_od(double re, double im) { double r18635 = 0.016666666666666666; double r18636 = im; double r18637 = 5.0; double r18638 = pow(r18636, r18637); double r18639 = r18636 * (r18636 * r18636); double r18640 = 0.3333333333333333; double r18641 = 2.0; double r18642 = r18641 * r18636; double r18643 = fma(r18639, r18640, r18642); double r18644 = fma(r18635, r18638, r18643); double r18645 = 1.0; double r18646 = r18644 / r18645; double r18647 = re; double r18648 = cos(r18647); double r18649 = 0.5; double r18650 = -r18649; double r18651 = r18648 * r18650; double r18652 = r18646 * r18651; return r18652; } 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 r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18653, "0.5", 10, MPFR_RNDN); mpfr_init(r18654); mpfr_init(r18655); mpfr_init(r18656); mpfr_init_set_str(r18657, "0", 10, MPFR_RNDN); mpfr_init(r18658); mpfr_init(r18659); mpfr_init(r18660); mpfr_init(r18661); mpfr_init(r18662); mpfr_init(r18663); } double f_im(double re, double im) { ; mpfr_set_d(r18654, re, MPFR_RNDN); mpfr_cos(r18655, r18654, MPFR_RNDN); mpfr_mul(r18656, r18653, r18655, MPFR_RNDN); ; mpfr_set_d(r18658, im, MPFR_RNDN); mpfr_sub(r18659, r18657, r18658, MPFR_RNDN); mpfr_exp(r18660, r18659, MPFR_RNDN); mpfr_exp(r18661, r18658, MPFR_RNDN); mpfr_sub(r18662, r18660, r18661, MPFR_RNDN); mpfr_mul(r18663, r18656, r18662, MPFR_RNDN); return mpfr_get_d(r18663, MPFR_RNDN); } static mpfr_t r18664, r18665, r18666, r18667, r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18664, "1/60", 10, MPFR_RNDN); mpfr_init(r18665); mpfr_init_set_str(r18666, "5", 10, MPFR_RNDN); mpfr_init(r18667); mpfr_init(r18668); mpfr_init_set_str(r18669, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18670, "2", 10, MPFR_RNDN); mpfr_init(r18671); mpfr_init(r18672); mpfr_init(r18673); mpfr_init_set_str(r18674, "1", 10, MPFR_RNDN); mpfr_init(r18675); mpfr_init(r18676); mpfr_init(r18677); mpfr_init_set_str(r18678, "0.5", 10, MPFR_RNDN); mpfr_init(r18679); mpfr_init(r18680); mpfr_init(r18681); } double f_fm(double re, double im) { ; mpfr_set_d(r18665, im, MPFR_RNDN); ; mpfr_pow(r18667, r18665, r18666, MPFR_RNDN); mpfr_mul(r18668, r18665, r18665, MPFR_RNDN); mpfr_mul(r18668, r18668, r18665, MPFR_RNDN); ; ; mpfr_mul(r18671, r18670, r18665, MPFR_RNDN); mpfr_fma(r18672, r18668, r18669, r18671, MPFR_RNDN); mpfr_fma(r18673, r18664, r18667, r18672, MPFR_RNDN); ; mpfr_div(r18675, r18673, r18674, MPFR_RNDN); mpfr_set_d(r18676, re, MPFR_RNDN); mpfr_cos(r18677, r18676, MPFR_RNDN); ; mpfr_neg(r18679, r18678, MPFR_RNDN); mpfr_mul(r18680, r18677, r18679, MPFR_RNDN); mpfr_mul(r18681, r18675, r18680, MPFR_RNDN); return mpfr_get_d(r18681, MPFR_RNDN); } static mpfr_t r18682, r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18682, "1/60", 10, MPFR_RNDN); mpfr_init(r18683); mpfr_init_set_str(r18684, "5", 10, MPFR_RNDN); mpfr_init(r18685); mpfr_init(r18686); mpfr_init_set_str(r18687, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18688, "2", 10, MPFR_RNDN); mpfr_init(r18689); mpfr_init(r18690); mpfr_init(r18691); mpfr_init_set_str(r18692, "1", 10, MPFR_RNDN); mpfr_init(r18693); mpfr_init(r18694); mpfr_init(r18695); mpfr_init_set_str(r18696, "0.5", 10, MPFR_RNDN); mpfr_init(r18697); mpfr_init(r18698); mpfr_init(r18699); } double f_dm(double re, double im) { ; mpfr_set_d(r18683, im, MPFR_RNDN); ; mpfr_pow(r18685, r18683, r18684, MPFR_RNDN); mpfr_mul(r18686, r18683, r18683, MPFR_RNDN); mpfr_mul(r18686, r18686, r18683, MPFR_RNDN); ; ; mpfr_mul(r18689, r18688, r18683, MPFR_RNDN); mpfr_fma(r18690, r18686, r18687, r18689, MPFR_RNDN); mpfr_fma(r18691, r18682, r18685, r18690, MPFR_RNDN); ; mpfr_div(r18693, r18691, r18692, MPFR_RNDN); mpfr_set_d(r18694, re, MPFR_RNDN); mpfr_cos(r18695, r18694, MPFR_RNDN); ; mpfr_neg(r18697, r18696, MPFR_RNDN); mpfr_mul(r18698, r18695, r18697, MPFR_RNDN); mpfr_mul(r18699, r18693, r18698, MPFR_RNDN); return mpfr_get_d(r18699, MPFR_RNDN); }