#include #include #include #include #include char *name = "cos2 (problem 3.4.1)"; double f_if(float x) { float r21691 = 1; float r21692 = x; float r21693 = cos(r21692); float r21694 = r21691 - r21693; float r21695 = r21692 * r21692; float r21696 = r21694 / r21695; return r21696; } double f_id(double x) { double r21697 = 1; double r21698 = x; double r21699 = cos(r21698); double r21700 = r21697 - r21699; double r21701 = r21698 * r21698; double r21702 = r21700 / r21701; return r21702; } double f_of(float x) { float r21703 = x; float r21704 = -0.03021091695231464; bool r21705 = r21703 <= r21704; float r21706 = 1; float r21707 = cos(r21703); float r21708 = r21706 - r21707; float r21709 = r21708 / r21703; float r21710 = r21709 / r21703; float r21711 = 0.03040127665186386; bool r21712 = r21703 <= r21711; float r21713 = 4; float r21714 = pow(r21703, r21713); float r21715 = 1/720; float r21716 = 1/2; float r21717 = fma(r21714, r21715, r21716); float r21718 = 1/24; float r21719 = r21718 * r21703; float r21720 = r21703 * r21719; float r21721 = r21717 - r21720; float r21722 = r21712 ? r21721 : r21710; float r21723 = r21705 ? r21710 : r21722; return r21723; } double f_od(double x) { double r21724 = x; double r21725 = -0.03021091695231464; bool r21726 = r21724 <= r21725; double r21727 = 1; double r21728 = cos(r21724); double r21729 = r21727 - r21728; double r21730 = r21729 / r21724; double r21731 = r21730 / r21724; double r21732 = 0.03040127665186386; bool r21733 = r21724 <= r21732; double r21734 = 4; double r21735 = pow(r21724, r21734); double r21736 = 1/720; double r21737 = 1/2; double r21738 = fma(r21735, r21736, r21737); double r21739 = 1/24; double r21740 = r21739 * r21724; double r21741 = r21724 * r21740; double r21742 = r21738 - r21741; double r21743 = r21733 ? r21742 : r21731; double r21744 = r21726 ? r21731 : r21743; return r21744; } 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 r21745, r21746, r21747, r21748, r21749, r21750; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init_set_str(r21745, "1", 10, MPFR_RNDN); mpfr_init(r21746); mpfr_init(r21747); mpfr_init(r21748); mpfr_init(r21749); mpfr_init(r21750); } double f_im(double x) { ; mpfr_set_d(r21746, x, MPFR_RNDN); mpfr_cos(r21747, r21746, MPFR_RNDN); mpfr_sub(r21748, r21745, r21747, MPFR_RNDN); mpfr_mul(r21749, r21746, r21746, MPFR_RNDN); mpfr_div(r21750, r21748, r21749, MPFR_RNDN); return mpfr_get_d(r21750, MPFR_RNDN); } static mpfr_t r21751, r21752, r21753, r21754, r21755, r21756, r21757, r21758, r21759, r21760, r21761, r21762, r21763, r21764, r21765, r21766, r21767, r21768, r21769, r21770, r21771; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r21751); mpfr_init_set_str(r21752, "-0.03021091695231464", 10, MPFR_RNDN); mpfr_init(r21753); mpfr_init_set_str(r21754, "1", 10, MPFR_RNDN); mpfr_init(r21755); mpfr_init(r21756); mpfr_init(r21757); mpfr_init(r21758); mpfr_init_set_str(r21759, "0.03040127665186386", 10, MPFR_RNDN); mpfr_init(r21760); mpfr_init_set_str(r21761, "4", 10, MPFR_RNDN); mpfr_init(r21762); mpfr_init_set_str(r21763, "1/720", 10, MPFR_RNDN); mpfr_init_set_str(r21764, "1/2", 10, MPFR_RNDN); mpfr_init(r21765); mpfr_init_set_str(r21766, "1/24", 10, MPFR_RNDN); mpfr_init(r21767); mpfr_init(r21768); mpfr_init(r21769); mpfr_init(r21770); mpfr_init(r21771); } double f_fm(double x) { mpfr_set_d(r21751, x, MPFR_RNDN); ; mpfr_set_si(r21753, mpfr_cmp(r21751, r21752) <= 0, MPFR_RNDN); ; mpfr_cos(r21755, r21751, MPFR_RNDN); mpfr_sub(r21756, r21754, r21755, MPFR_RNDN); mpfr_div(r21757, r21756, r21751, MPFR_RNDN); mpfr_div(r21758, r21757, r21751, MPFR_RNDN); ; mpfr_set_si(r21760, mpfr_cmp(r21751, r21759) <= 0, MPFR_RNDN); ; mpfr_pow(r21762, r21751, r21761, MPFR_RNDN); ; ; mpfr_fma(r21765, r21762, r21763, r21764, MPFR_RNDN); ; mpfr_mul(r21767, r21766, r21751, MPFR_RNDN); mpfr_mul(r21768, r21751, r21767, MPFR_RNDN); mpfr_sub(r21769, r21765, r21768, MPFR_RNDN); if (mpfr_get_si(r21760, MPFR_RNDN)) { mpfr_set(r21770, r21769, MPFR_RNDN); } else { mpfr_set(r21770, r21758, MPFR_RNDN); }; if (mpfr_get_si(r21753, MPFR_RNDN)) { mpfr_set(r21771, r21758, MPFR_RNDN); } else { mpfr_set(r21771, r21770, MPFR_RNDN); }; return mpfr_get_d(r21771, MPFR_RNDN); } static mpfr_t r21772, r21773, r21774, r21775, r21776, r21777, r21778, r21779, r21780, r21781, r21782, r21783, r21784, r21785, r21786, r21787, r21788, r21789, r21790, r21791, r21792; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r21772); mpfr_init_set_str(r21773, "-0.03021091695231464", 10, MPFR_RNDN); mpfr_init(r21774); mpfr_init_set_str(r21775, "1", 10, MPFR_RNDN); mpfr_init(r21776); mpfr_init(r21777); mpfr_init(r21778); mpfr_init(r21779); mpfr_init_set_str(r21780, "0.03040127665186386", 10, MPFR_RNDN); mpfr_init(r21781); mpfr_init_set_str(r21782, "4", 10, MPFR_RNDN); mpfr_init(r21783); mpfr_init_set_str(r21784, "1/720", 10, MPFR_RNDN); mpfr_init_set_str(r21785, "1/2", 10, MPFR_RNDN); mpfr_init(r21786); mpfr_init_set_str(r21787, "1/24", 10, MPFR_RNDN); mpfr_init(r21788); mpfr_init(r21789); mpfr_init(r21790); mpfr_init(r21791); mpfr_init(r21792); } double f_dm(double x) { mpfr_set_d(r21772, x, MPFR_RNDN); ; mpfr_set_si(r21774, mpfr_cmp(r21772, r21773) <= 0, MPFR_RNDN); ; mpfr_cos(r21776, r21772, MPFR_RNDN); mpfr_sub(r21777, r21775, r21776, MPFR_RNDN); mpfr_div(r21778, r21777, r21772, MPFR_RNDN); mpfr_div(r21779, r21778, r21772, MPFR_RNDN); ; mpfr_set_si(r21781, mpfr_cmp(r21772, r21780) <= 0, MPFR_RNDN); ; mpfr_pow(r21783, r21772, r21782, MPFR_RNDN); ; ; mpfr_fma(r21786, r21783, r21784, r21785, MPFR_RNDN); ; mpfr_mul(r21788, r21787, r21772, MPFR_RNDN); mpfr_mul(r21789, r21772, r21788, MPFR_RNDN); mpfr_sub(r21790, r21786, r21789, MPFR_RNDN); if (mpfr_get_si(r21781, MPFR_RNDN)) { mpfr_set(r21791, r21790, MPFR_RNDN); } else { mpfr_set(r21791, r21779, MPFR_RNDN); }; if (mpfr_get_si(r21774, MPFR_RNDN)) { mpfr_set(r21792, r21779, MPFR_RNDN); } else { mpfr_set(r21792, r21791, MPFR_RNDN); }; return mpfr_get_d(r21792, MPFR_RNDN); }