#include #include #include #include #include char *name = "cos(2*x)/(cos^2(x)*sin^2(x))"; double f_if(float x, float cos, float sin) { float r22623 = 2; float r22624 = x; float r22625 = r22623 * r22624; float r22626 = cos(r22625); float r22627 = cos; float r22628 = pow(r22627, r22623); float r22629 = sin; float r22630 = pow(r22629, r22623); float r22631 = r22624 * r22630; float r22632 = r22631 * r22624; float r22633 = r22628 * r22632; float r22634 = r22626 / r22633; return r22634; } double f_id(double x, double cos, double sin) { double r22635 = 2; double r22636 = x; double r22637 = r22635 * r22636; double r22638 = cos(r22637); double r22639 = cos; double r22640 = pow(r22639, r22635); double r22641 = sin; double r22642 = pow(r22641, r22635); double r22643 = r22636 * r22642; double r22644 = r22643 * r22636; double r22645 = r22640 * r22644; double r22646 = r22638 / r22645; return r22646; } double f_of(float x, float cos, float sin) { float r22647 = 2; float r22648 = x; float r22649 = r22647 * r22648; float r22650 = cos(r22649); float r22651 = cos; float r22652 = pow(r22651, r22647); float r22653 = sqrt(r22652); float r22654 = sin; float r22655 = r22654 * r22648; float r22656 = r22655 * r22655; float r22657 = fabs(r22651); float r22658 = r22656 * r22657; float r22659 = r22653 * r22658; float r22660 = r22650 / r22659; float r22661 = 4.082477361157426e+244; bool r22662 = r22660 <= r22661; float r22663 = r22648 * r22654; float r22664 = r22651 * r22663; float r22665 = fabs(r22664); float r22666 = r22650 / r22665; float r22667 = r22666 / r22665; float r22668 = r22654 * r22651; float r22669 = r22668 * r22648; float r22670 = fabs(r22669); float r22671 = fabs(r22670); float r22672 = r22671 * r22671; float r22673 = r22650 / r22672; float r22674 = r22662 ? r22667 : r22673; return r22674; } double f_od(double x, double cos, double sin) { double r22675 = 2; double r22676 = x; double r22677 = r22675 * r22676; double r22678 = cos(r22677); double r22679 = cos; double r22680 = pow(r22679, r22675); double r22681 = sqrt(r22680); double r22682 = sin; double r22683 = r22682 * r22676; double r22684 = r22683 * r22683; double r22685 = fabs(r22679); double r22686 = r22684 * r22685; double r22687 = r22681 * r22686; double r22688 = r22678 / r22687; double r22689 = 4.082477361157426e+244; bool r22690 = r22688 <= r22689; double r22691 = r22676 * r22682; double r22692 = r22679 * r22691; double r22693 = fabs(r22692); double r22694 = r22678 / r22693; double r22695 = r22694 / r22693; double r22696 = r22682 * r22679; double r22697 = r22696 * r22676; double r22698 = fabs(r22697); double r22699 = fabs(r22698); double r22700 = r22699 * r22699; double r22701 = r22678 / r22700; double r22702 = r22690 ? r22695 : r22701; return r22702; } 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 r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r22703, "2", 10, MPFR_RNDN); mpfr_init(r22704); mpfr_init(r22705); mpfr_init(r22706); mpfr_init(r22707); mpfr_init(r22708); mpfr_init(r22709); mpfr_init(r22710); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r22704, x, MPFR_RNDN); mpfr_mul(r22705, r22703, r22704, MPFR_RNDN); mpfr_cos(r22706, r22705, MPFR_RNDN); mpfr_set_d(r22707, cos, MPFR_RNDN); mpfr_pow(r22708, r22707, r22703, MPFR_RNDN); mpfr_set_d(r22709, sin, MPFR_RNDN); mpfr_pow(r22710, r22709, r22703, MPFR_RNDN); mpfr_mul(r22711, r22704, r22710, MPFR_RNDN); mpfr_mul(r22712, r22711, r22704, MPFR_RNDN); mpfr_mul(r22713, r22708, r22712, MPFR_RNDN); mpfr_div(r22714, r22706, r22713, MPFR_RNDN); return mpfr_get_d(r22714, MPFR_RNDN); } static mpfr_t r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22715, "2", 10, MPFR_RNDN); mpfr_init(r22716); mpfr_init(r22717); mpfr_init(r22718); 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_set_str(r22729, "4.082477361157426e+244", 10, MPFR_RNDN); mpfr_init(r22730); mpfr_init(r22731); mpfr_init(r22732); mpfr_init(r22733); mpfr_init(r22734); mpfr_init(r22735); mpfr_init(r22736); mpfr_init(r22737); mpfr_init(r22738); mpfr_init(r22739); mpfr_init(r22740); mpfr_init(r22741); mpfr_init(r22742); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r22716, x, MPFR_RNDN); mpfr_mul(r22717, r22715, r22716, MPFR_RNDN); mpfr_cos(r22718, r22717, MPFR_RNDN); mpfr_set_d(r22719, cos, MPFR_RNDN); mpfr_pow(r22720, r22719, r22715, MPFR_RNDN); mpfr_sqrt(r22721, r22720, MPFR_RNDN); mpfr_set_d(r22722, sin, MPFR_RNDN); mpfr_mul(r22723, r22722, r22716, MPFR_RNDN); mpfr_mul(r22724, r22723, r22723, MPFR_RNDN); mpfr_abs(r22725, r22719, MPFR_RNDN); mpfr_mul(r22726, r22724, r22725, MPFR_RNDN); mpfr_mul(r22727, r22721, r22726, MPFR_RNDN); mpfr_div(r22728, r22718, r22727, MPFR_RNDN); ; mpfr_set_si(r22730, mpfr_cmp(r22728, r22729) <= 0, MPFR_RNDN); mpfr_mul(r22731, r22716, r22722, MPFR_RNDN); mpfr_mul(r22732, r22719, r22731, MPFR_RNDN); mpfr_abs(r22733, r22732, MPFR_RNDN); mpfr_div(r22734, r22718, r22733, MPFR_RNDN); mpfr_div(r22735, r22734, r22733, MPFR_RNDN); mpfr_mul(r22736, r22722, r22719, MPFR_RNDN); mpfr_mul(r22737, r22736, r22716, MPFR_RNDN); mpfr_abs(r22738, r22737, MPFR_RNDN); mpfr_abs(r22739, r22738, MPFR_RNDN); mpfr_mul(r22740, r22739, r22739, MPFR_RNDN); mpfr_div(r22741, r22718, r22740, MPFR_RNDN); if (mpfr_get_si(r22730, MPFR_RNDN)) { mpfr_set(r22742, r22735, MPFR_RNDN); } else { mpfr_set(r22742, r22741, MPFR_RNDN); }; return mpfr_get_d(r22742, MPFR_RNDN); } static mpfr_t r22743, r22744, r22745, r22746, r22747, r22748, r22749, r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764, r22765, r22766, r22767, r22768, r22769, r22770; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22743, "2", 10, MPFR_RNDN); mpfr_init(r22744); mpfr_init(r22745); mpfr_init(r22746); mpfr_init(r22747); mpfr_init(r22748); mpfr_init(r22749); mpfr_init(r22750); mpfr_init(r22751); mpfr_init(r22752); mpfr_init(r22753); mpfr_init(r22754); mpfr_init(r22755); mpfr_init(r22756); mpfr_init_set_str(r22757, "4.082477361157426e+244", 10, MPFR_RNDN); mpfr_init(r22758); mpfr_init(r22759); mpfr_init(r22760); mpfr_init(r22761); mpfr_init(r22762); mpfr_init(r22763); mpfr_init(r22764); mpfr_init(r22765); mpfr_init(r22766); mpfr_init(r22767); mpfr_init(r22768); mpfr_init(r22769); mpfr_init(r22770); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r22744, x, MPFR_RNDN); mpfr_mul(r22745, r22743, r22744, MPFR_RNDN); mpfr_cos(r22746, r22745, MPFR_RNDN); mpfr_set_d(r22747, cos, MPFR_RNDN); mpfr_pow(r22748, r22747, r22743, MPFR_RNDN); mpfr_sqrt(r22749, r22748, MPFR_RNDN); mpfr_set_d(r22750, sin, MPFR_RNDN); mpfr_mul(r22751, r22750, r22744, MPFR_RNDN); mpfr_mul(r22752, r22751, r22751, MPFR_RNDN); mpfr_abs(r22753, r22747, MPFR_RNDN); mpfr_mul(r22754, r22752, r22753, MPFR_RNDN); mpfr_mul(r22755, r22749, r22754, MPFR_RNDN); mpfr_div(r22756, r22746, r22755, MPFR_RNDN); ; mpfr_set_si(r22758, mpfr_cmp(r22756, r22757) <= 0, MPFR_RNDN); mpfr_mul(r22759, r22744, r22750, MPFR_RNDN); mpfr_mul(r22760, r22747, r22759, MPFR_RNDN); mpfr_abs(r22761, r22760, MPFR_RNDN); mpfr_div(r22762, r22746, r22761, MPFR_RNDN); mpfr_div(r22763, r22762, r22761, MPFR_RNDN); mpfr_mul(r22764, r22750, r22747, MPFR_RNDN); mpfr_mul(r22765, r22764, r22744, MPFR_RNDN); mpfr_abs(r22766, r22765, MPFR_RNDN); mpfr_abs(r22767, r22766, MPFR_RNDN); mpfr_mul(r22768, r22767, r22767, MPFR_RNDN); mpfr_div(r22769, r22746, r22768, MPFR_RNDN); if (mpfr_get_si(r22758, MPFR_RNDN)) { mpfr_set(r22770, r22763, MPFR_RNDN); } else { mpfr_set(r22770, r22769, MPFR_RNDN); }; return mpfr_get_d(r22770, MPFR_RNDN); }