#include #include #include #include #include char *name = "math.cos on complex, real part"; double f_if(float re, float im) { float r22624 = 0.5; float r22625 = re; float r22626 = cos(r22625); float r22627 = r22624 * r22626; float r22628 = im; float r22629 = -r22628; float r22630 = exp(r22629); float r22631 = exp(r22628); float r22632 = r22630 + r22631; float r22633 = r22627 * r22632; return r22633; } double f_id(double re, double im) { double r22634 = 0.5; double r22635 = re; double r22636 = cos(r22635); double r22637 = r22634 * r22636; double r22638 = im; double r22639 = -r22638; double r22640 = exp(r22639); double r22641 = exp(r22638); double r22642 = r22640 + r22641; double r22643 = r22637 * r22642; return r22643; } double f_of(float re, float im) { float r22644 = 0.5; float r22645 = re; float r22646 = cos(r22645); float r22647 = r22644 * r22646; float r22648 = im; float r22649 = -r22648; float r22650 = exp(r22649); float r22651 = 3; float r22652 = pow(r22650, r22651); float r22653 = exp(r22648); float r22654 = pow(r22653, r22651); float r22655 = r22652 + r22654; float r22656 = r22647 * r22655; float r22657 = r22650 * r22650; float r22658 = r22653 * r22653; float r22659 = r22650 * r22653; float r22660 = r22658 - r22659; float r22661 = r22657 + r22660; float r22662 = r22656 / r22661; return r22662; } double f_od(double re, double im) { double r22663 = 0.5; double r22664 = re; double r22665 = cos(r22664); double r22666 = r22663 * r22665; double r22667 = im; double r22668 = -r22667; double r22669 = exp(r22668); double r22670 = 3; double r22671 = pow(r22669, r22670); double r22672 = exp(r22667); double r22673 = pow(r22672, r22670); double r22674 = r22671 + r22673; double r22675 = r22666 * r22674; double r22676 = r22669 * r22669; double r22677 = r22672 * r22672; double r22678 = r22669 * r22672; double r22679 = r22677 - r22678; double r22680 = r22676 + r22679; double r22681 = r22675 / r22680; return r22681; } 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 r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689, r22690, r22691; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22682, "0.5", 10, MPFR_RNDN); mpfr_init(r22683); mpfr_init(r22684); mpfr_init(r22685); mpfr_init(r22686); mpfr_init(r22687); mpfr_init(r22688); mpfr_init(r22689); mpfr_init(r22690); mpfr_init(r22691); } double f_im(double re, double im) { ; mpfr_set_d(r22683, re, MPFR_RNDN); mpfr_cos(r22684, r22683, MPFR_RNDN); mpfr_mul(r22685, r22682, r22684, MPFR_RNDN); mpfr_set_d(r22686, im, MPFR_RNDN); mpfr_neg(r22687, r22686, MPFR_RNDN); mpfr_exp(r22688, r22687, MPFR_RNDN); mpfr_exp(r22689, r22686, MPFR_RNDN); mpfr_add(r22690, r22688, r22689, MPFR_RNDN); mpfr_mul(r22691, r22685, r22690, MPFR_RNDN); return mpfr_get_d(r22691, MPFR_RNDN); } static mpfr_t r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699, r22700, r22701, r22702, r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22692, "0.5", 10, MPFR_RNDN); mpfr_init(r22693); mpfr_init(r22694); mpfr_init(r22695); mpfr_init(r22696); mpfr_init(r22697); mpfr_init(r22698); mpfr_init_set_str(r22699, "3", 10, MPFR_RNDN); mpfr_init(r22700); mpfr_init(r22701); mpfr_init(r22702); mpfr_init(r22703); mpfr_init(r22704); mpfr_init(r22705); mpfr_init(r22706); mpfr_init(r22707); mpfr_init(r22708); mpfr_init(r22709); mpfr_init(r22710); } double f_fm(double re, double im) { ; mpfr_set_d(r22693, re, MPFR_RNDN); mpfr_cos(r22694, r22693, MPFR_RNDN); mpfr_mul(r22695, r22692, r22694, MPFR_RNDN); mpfr_set_d(r22696, im, MPFR_RNDN); mpfr_neg(r22697, r22696, MPFR_RNDN); mpfr_exp(r22698, r22697, MPFR_RNDN); ; mpfr_pow(r22700, r22698, r22699, MPFR_RNDN); mpfr_exp(r22701, r22696, MPFR_RNDN); mpfr_pow(r22702, r22701, r22699, MPFR_RNDN); mpfr_add(r22703, r22700, r22702, MPFR_RNDN); mpfr_mul(r22704, r22695, r22703, MPFR_RNDN); mpfr_mul(r22705, r22698, r22698, MPFR_RNDN); mpfr_mul(r22706, r22701, r22701, MPFR_RNDN); mpfr_mul(r22707, r22698, r22701, MPFR_RNDN); mpfr_sub(r22708, r22706, r22707, MPFR_RNDN); mpfr_add(r22709, r22705, r22708, MPFR_RNDN); mpfr_div(r22710, r22704, r22709, MPFR_RNDN); return mpfr_get_d(r22710, MPFR_RNDN); } static mpfr_t r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22711, "0.5", 10, MPFR_RNDN); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); mpfr_init(r22715); mpfr_init(r22716); mpfr_init(r22717); mpfr_init_set_str(r22718, "3", 10, MPFR_RNDN); mpfr_init(r22719); mpfr_init(r22720); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); mpfr_init(r22725); mpfr_init(r22726); mpfr_init(r22727); mpfr_init(r22728); mpfr_init(r22729); } double f_dm(double re, double im) { ; mpfr_set_d(r22712, re, MPFR_RNDN); mpfr_cos(r22713, r22712, MPFR_RNDN); mpfr_mul(r22714, r22711, r22713, MPFR_RNDN); mpfr_set_d(r22715, im, MPFR_RNDN); mpfr_neg(r22716, r22715, MPFR_RNDN); mpfr_exp(r22717, r22716, MPFR_RNDN); ; mpfr_pow(r22719, r22717, r22718, MPFR_RNDN); mpfr_exp(r22720, r22715, MPFR_RNDN); mpfr_pow(r22721, r22720, r22718, MPFR_RNDN); mpfr_add(r22722, r22719, r22721, MPFR_RNDN); mpfr_mul(r22723, r22714, r22722, MPFR_RNDN); mpfr_mul(r22724, r22717, r22717, MPFR_RNDN); mpfr_mul(r22725, r22720, r22720, MPFR_RNDN); mpfr_mul(r22726, r22717, r22720, MPFR_RNDN); mpfr_sub(r22727, r22725, r22726, MPFR_RNDN); mpfr_add(r22728, r22724, r22727, MPFR_RNDN); mpfr_div(r22729, r22723, r22728, MPFR_RNDN); return mpfr_get_d(r22729, MPFR_RNDN); }