#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 = r18617 * r18620; float r18622 = 2.0f; float r18623 = r18622 * r18618; float r18624 = 0.3333333333333333f; float r18625 = 3.0f; float r18626 = pow(r18618, r18625); float r18627 = r18624 * r18626; float r18628 = r18623 + r18627; float r18629 = r18621 + r18628; float r18630 = -r18629; float r18631 = re; float r18632 = cos(r18631); float r18633 = 0.5f; float r18634 = r18632 * r18633; float r18635 = r18630 * r18634; return r18635; } double f_od(double re, double im) { double r18636 = 0.016666666666666666; double r18637 = im; double r18638 = 5.0; double r18639 = pow(r18637, r18638); double r18640 = r18636 * r18639; double r18641 = 2.0; double r18642 = r18641 * r18637; double r18643 = 0.3333333333333333; double r18644 = 3.0; double r18645 = pow(r18637, r18644); double r18646 = r18643 * r18645; double r18647 = r18642 + r18646; double r18648 = r18640 + r18647; double r18649 = -r18648; double r18650 = re; double r18651 = cos(r18650); double r18652 = 0.5; double r18653 = r18651 * r18652; double r18654 = r18649 * r18653; return r18654; } 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 r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18655, "0.5", 10, MPFR_RNDN); mpfr_init(r18656); mpfr_init(r18657); mpfr_init(r18658); mpfr_init_set_str(r18659, "0", 10, MPFR_RNDN); mpfr_init(r18660); mpfr_init(r18661); mpfr_init(r18662); mpfr_init(r18663); mpfr_init(r18664); mpfr_init(r18665); } double f_im(double re, double im) { ; mpfr_set_d(r18656, re, MPFR_RNDN); mpfr_cos(r18657, r18656, MPFR_RNDN); mpfr_mul(r18658, r18655, r18657, MPFR_RNDN); ; mpfr_set_d(r18660, im, MPFR_RNDN); mpfr_sub(r18661, r18659, r18660, MPFR_RNDN); mpfr_exp(r18662, r18661, MPFR_RNDN); mpfr_exp(r18663, r18660, MPFR_RNDN); mpfr_sub(r18664, r18662, r18663, MPFR_RNDN); mpfr_mul(r18665, r18658, r18664, MPFR_RNDN); return mpfr_get_d(r18665, MPFR_RNDN); } static mpfr_t r18666, r18667, r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681, r18682, r18683, r18684; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18666, "1/60", 10, MPFR_RNDN); mpfr_init(r18667); mpfr_init_set_str(r18668, "5", 10, MPFR_RNDN); mpfr_init(r18669); mpfr_init(r18670); mpfr_init_set_str(r18671, "2", 10, MPFR_RNDN); mpfr_init(r18672); mpfr_init_set_str(r18673, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18674, "3", 10, MPFR_RNDN); mpfr_init(r18675); mpfr_init(r18676); mpfr_init(r18677); mpfr_init(r18678); mpfr_init(r18679); mpfr_init(r18680); mpfr_init(r18681); mpfr_init_set_str(r18682, "0.5", 10, MPFR_RNDN); mpfr_init(r18683); mpfr_init(r18684); } double f_fm(double re, double im) { ; mpfr_set_d(r18667, im, MPFR_RNDN); ; mpfr_pow(r18669, r18667, r18668, MPFR_RNDN); mpfr_mul(r18670, r18666, r18669, MPFR_RNDN); ; mpfr_mul(r18672, r18671, r18667, MPFR_RNDN); ; ; mpfr_pow(r18675, r18667, r18674, MPFR_RNDN); mpfr_mul(r18676, r18673, r18675, MPFR_RNDN); mpfr_add(r18677, r18672, r18676, MPFR_RNDN); mpfr_add(r18678, r18670, r18677, MPFR_RNDN); mpfr_neg(r18679, r18678, MPFR_RNDN); mpfr_set_d(r18680, re, MPFR_RNDN); mpfr_cos(r18681, r18680, MPFR_RNDN); ; mpfr_mul(r18683, r18681, r18682, MPFR_RNDN); mpfr_mul(r18684, r18679, r18683, MPFR_RNDN); return mpfr_get_d(r18684, MPFR_RNDN); } static mpfr_t r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18685, "1/60", 10, MPFR_RNDN); mpfr_init(r18686); mpfr_init_set_str(r18687, "5", 10, MPFR_RNDN); mpfr_init(r18688); mpfr_init(r18689); mpfr_init_set_str(r18690, "2", 10, MPFR_RNDN); mpfr_init(r18691); mpfr_init_set_str(r18692, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18693, "3", 10, MPFR_RNDN); mpfr_init(r18694); mpfr_init(r18695); mpfr_init(r18696); mpfr_init(r18697); mpfr_init(r18698); mpfr_init(r18699); mpfr_init(r18700); mpfr_init_set_str(r18701, "0.5", 10, MPFR_RNDN); mpfr_init(r18702); mpfr_init(r18703); } double f_dm(double re, double im) { ; mpfr_set_d(r18686, im, MPFR_RNDN); ; mpfr_pow(r18688, r18686, r18687, MPFR_RNDN); mpfr_mul(r18689, r18685, r18688, MPFR_RNDN); ; mpfr_mul(r18691, r18690, r18686, MPFR_RNDN); ; ; mpfr_pow(r18694, r18686, r18693, MPFR_RNDN); mpfr_mul(r18695, r18692, r18694, MPFR_RNDN); mpfr_add(r18696, r18691, r18695, MPFR_RNDN); mpfr_add(r18697, r18689, r18696, MPFR_RNDN); mpfr_neg(r18698, r18697, MPFR_RNDN); mpfr_set_d(r18699, re, MPFR_RNDN); mpfr_cos(r18700, r18699, MPFR_RNDN); ; mpfr_mul(r18702, r18700, r18701, MPFR_RNDN); mpfr_mul(r18703, r18698, r18702, MPFR_RNDN); return mpfr_get_d(r18703, MPFR_RNDN); }