#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17671 = b; float r17672 = -r17671; float r17673 = r17671 * r17671; float r17674 = 4.0f; float r17675 = a; float r17676 = r17674 * r17675; float r17677 = c; float r17678 = r17676 * r17677; float r17679 = r17673 - r17678; float r17680 = sqrt(r17679); float r17681 = r17672 + r17680; float r17682 = 2.0f; float r17683 = r17682 * r17675; float r17684 = r17681 / r17683; return r17684; } double f_id(double a, double b, double c) { double r17685 = b; double r17686 = -r17685; double r17687 = r17685 * r17685; double r17688 = 4.0; double r17689 = a; double r17690 = r17688 * r17689; double r17691 = c; double r17692 = r17690 * r17691; double r17693 = r17687 - r17692; double r17694 = sqrt(r17693); double r17695 = r17686 + r17694; double r17696 = 2.0; double r17697 = r17696 * r17689; double r17698 = r17695 / r17697; return r17698; } double f_of(float a, float b, float c) { float r17699 = 4.0f; float r17700 = a; float r17701 = r17699 * r17700; float r17702 = 1.0f; float r17703 = r17701 / r17702; float r17704 = c; float r17705 = b; float r17706 = -r17705; float r17707 = r17705 * r17705; float r17708 = r17701 * r17704; float r17709 = r17707 - r17708; float r17710 = sqrt(r17709); float r17711 = r17706 - r17710; float r17712 = r17704 / r17711; float r17713 = r17703 * r17712; float r17714 = 2.0f; float r17715 = r17714 * r17700; float r17716 = r17713 / r17715; return r17716; } double f_od(double a, double b, double c) { double r17717 = 4.0; double r17718 = a; double r17719 = r17717 * r17718; double r17720 = 1.0; double r17721 = r17719 / r17720; double r17722 = c; double r17723 = b; double r17724 = -r17723; double r17725 = r17723 * r17723; double r17726 = r17719 * r17722; double r17727 = r17725 - r17726; double r17728 = sqrt(r17727); double r17729 = r17724 - r17728; double r17730 = r17722 / r17729; double r17731 = r17721 * r17730; double r17732 = 2.0; double r17733 = r17732 * r17718; double r17734 = r17731 / r17733; return r17734; } 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 r17735, r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743, r17744, r17745, r17746, r17747, r17748; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17735); mpfr_init(r17736); mpfr_init(r17737); mpfr_init_set_str(r17738, "4", 10, MPFR_RNDN); mpfr_init(r17739); mpfr_init(r17740); mpfr_init(r17741); mpfr_init(r17742); mpfr_init(r17743); mpfr_init(r17744); mpfr_init(r17745); mpfr_init_set_str(r17746, "2", 10, MPFR_RNDN); mpfr_init(r17747); mpfr_init(r17748); } double f_im(double a, double b, double c) { mpfr_set_d(r17735, b, MPFR_RNDN); mpfr_neg(r17736, r17735, MPFR_RNDN); mpfr_sqr(r17737, r17735, MPFR_RNDN); ; mpfr_set_d(r17739, a, MPFR_RNDN); mpfr_mul(r17740, r17738, r17739, MPFR_RNDN); mpfr_set_d(r17741, c, MPFR_RNDN); mpfr_mul(r17742, r17740, r17741, MPFR_RNDN); mpfr_sub(r17743, r17737, r17742, MPFR_RNDN); mpfr_sqrt(r17744, r17743, MPFR_RNDN); mpfr_add(r17745, r17736, r17744, MPFR_RNDN); ; mpfr_mul(r17747, r17746, r17739, MPFR_RNDN); mpfr_div(r17748, r17745, r17747, MPFR_RNDN); return mpfr_get_d(r17748, MPFR_RNDN); } static mpfr_t r17749, r17750, r17751, r17752, r17753, r17754, r17755, r17756, r17757, r17758, r17759, r17760, r17761, r17762, r17763, r17764, r17765, r17766; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17749, "4", 10, MPFR_RNDN); mpfr_init(r17750); mpfr_init(r17751); mpfr_init_set_str(r17752, "1", 10, MPFR_RNDN); mpfr_init(r17753); mpfr_init(r17754); mpfr_init(r17755); mpfr_init(r17756); mpfr_init(r17757); mpfr_init(r17758); mpfr_init(r17759); mpfr_init(r17760); mpfr_init(r17761); mpfr_init(r17762); mpfr_init(r17763); mpfr_init_set_str(r17764, "2", 10, MPFR_RNDN); mpfr_init(r17765); mpfr_init(r17766); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17750, a, MPFR_RNDN); mpfr_mul(r17751, r17749, r17750, MPFR_RNDN); ; mpfr_div(r17753, r17751, r17752, MPFR_RNDN); mpfr_set_d(r17754, c, MPFR_RNDN); mpfr_set_d(r17755, b, MPFR_RNDN); mpfr_neg(r17756, r17755, MPFR_RNDN); mpfr_sqr(r17757, r17755, MPFR_RNDN); mpfr_mul(r17758, r17751, r17754, MPFR_RNDN); mpfr_sub(r17759, r17757, r17758, MPFR_RNDN); mpfr_sqrt(r17760, r17759, MPFR_RNDN); mpfr_sub(r17761, r17756, r17760, MPFR_RNDN); mpfr_div(r17762, r17754, r17761, MPFR_RNDN); mpfr_mul(r17763, r17753, r17762, MPFR_RNDN); ; mpfr_mul(r17765, r17764, r17750, MPFR_RNDN); mpfr_div(r17766, r17763, r17765, MPFR_RNDN); return mpfr_get_d(r17766, MPFR_RNDN); } static mpfr_t r17767, r17768, r17769, r17770, r17771, r17772, r17773, r17774, r17775, r17776, r17777, r17778, r17779, r17780, r17781, r17782, r17783, r17784; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17767, "4", 10, MPFR_RNDN); mpfr_init(r17768); mpfr_init(r17769); mpfr_init_set_str(r17770, "1", 10, MPFR_RNDN); mpfr_init(r17771); mpfr_init(r17772); mpfr_init(r17773); mpfr_init(r17774); mpfr_init(r17775); mpfr_init(r17776); mpfr_init(r17777); mpfr_init(r17778); mpfr_init(r17779); mpfr_init(r17780); mpfr_init(r17781); mpfr_init_set_str(r17782, "2", 10, MPFR_RNDN); mpfr_init(r17783); mpfr_init(r17784); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17768, a, MPFR_RNDN); mpfr_mul(r17769, r17767, r17768, MPFR_RNDN); ; mpfr_div(r17771, r17769, r17770, MPFR_RNDN); mpfr_set_d(r17772, c, MPFR_RNDN); mpfr_set_d(r17773, b, MPFR_RNDN); mpfr_neg(r17774, r17773, MPFR_RNDN); mpfr_sqr(r17775, r17773, MPFR_RNDN); mpfr_mul(r17776, r17769, r17772, MPFR_RNDN); mpfr_sub(r17777, r17775, r17776, MPFR_RNDN); mpfr_sqrt(r17778, r17777, MPFR_RNDN); mpfr_sub(r17779, r17774, r17778, MPFR_RNDN); mpfr_div(r17780, r17772, r17779, MPFR_RNDN); mpfr_mul(r17781, r17771, r17780, MPFR_RNDN); ; mpfr_mul(r17783, r17782, r17768, MPFR_RNDN); mpfr_div(r17784, r17781, r17783, MPFR_RNDN); return mpfr_get_d(r17784, MPFR_RNDN); }