#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17771 = b; float r17772 = -r17771; float r17773 = r17771 * r17771; float r17774 = 4.0f; float r17775 = a; float r17776 = r17774 * r17775; float r17777 = c; float r17778 = r17776 * r17777; float r17779 = r17773 - r17778; float r17780 = sqrt(r17779); float r17781 = r17772 + r17780; float r17782 = 2.0f; float r17783 = r17782 * r17775; float r17784 = r17781 / r17783; return r17784; } double f_id(double a, double b, double c) { double r17785 = b; double r17786 = -r17785; double r17787 = r17785 * r17785; double r17788 = 4.0; double r17789 = a; double r17790 = r17788 * r17789; double r17791 = c; double r17792 = r17790 * r17791; double r17793 = r17787 - r17792; double r17794 = sqrt(r17793); double r17795 = r17786 + r17794; double r17796 = 2.0; double r17797 = r17796 * r17789; double r17798 = r17795 / r17797; return r17798; } double f_of(float a, float b, float c) { float r17799 = 1.0f; float r17800 = 2.0f; float r17801 = r17799 / r17800; float r17802 = 4.0f; float r17803 = c; float r17804 = r17802 * r17803; float r17805 = b; float r17806 = -r17805; float r17807 = r17805 * r17805; float r17808 = a; float r17809 = r17803 * r17808; float r17810 = r17809 * r17802; float r17811 = r17807 - r17810; float r17812 = sqrt(r17811); float r17813 = r17806 - r17812; float r17814 = r17804 / r17813; float r17815 = r17801 * r17814; return r17815; } double f_od(double a, double b, double c) { double r17816 = 1.0; double r17817 = 2.0; double r17818 = r17816 / r17817; double r17819 = 4.0; double r17820 = c; double r17821 = r17819 * r17820; double r17822 = b; double r17823 = -r17822; double r17824 = r17822 * r17822; double r17825 = a; double r17826 = r17820 * r17825; double r17827 = r17826 * r17819; double r17828 = r17824 - r17827; double r17829 = sqrt(r17828); double r17830 = r17823 - r17829; double r17831 = r17821 / r17830; double r17832 = r17818 * r17831; return r17832; } 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 r17833, r17834, r17835, r17836, r17837, r17838, r17839, r17840, r17841, r17842, r17843, r17844, r17845, r17846; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17833); mpfr_init(r17834); mpfr_init(r17835); mpfr_init_set_str(r17836, "4", 10, MPFR_RNDN); mpfr_init(r17837); mpfr_init(r17838); mpfr_init(r17839); mpfr_init(r17840); mpfr_init(r17841); mpfr_init(r17842); mpfr_init(r17843); mpfr_init_set_str(r17844, "2", 10, MPFR_RNDN); mpfr_init(r17845); mpfr_init(r17846); } double f_im(double a, double b, double c) { mpfr_set_d(r17833, b, MPFR_RNDN); mpfr_neg(r17834, r17833, MPFR_RNDN); mpfr_sqr(r17835, r17833, MPFR_RNDN); ; mpfr_set_d(r17837, a, MPFR_RNDN); mpfr_mul(r17838, r17836, r17837, MPFR_RNDN); mpfr_set_d(r17839, c, MPFR_RNDN); mpfr_mul(r17840, r17838, r17839, MPFR_RNDN); mpfr_sub(r17841, r17835, r17840, MPFR_RNDN); mpfr_sqrt(r17842, r17841, MPFR_RNDN); mpfr_add(r17843, r17834, r17842, MPFR_RNDN); ; mpfr_mul(r17845, r17844, r17837, MPFR_RNDN); mpfr_div(r17846, r17843, r17845, MPFR_RNDN); return mpfr_get_d(r17846, MPFR_RNDN); } static mpfr_t r17847, r17848, r17849, r17850, r17851, r17852, r17853, r17854, r17855, r17856, r17857, r17858, r17859, r17860, r17861, r17862, r17863; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17847, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17848, "2", 10, MPFR_RNDN); mpfr_init(r17849); mpfr_init_set_str(r17850, "4", 10, MPFR_RNDN); mpfr_init(r17851); mpfr_init(r17852); mpfr_init(r17853); mpfr_init(r17854); mpfr_init(r17855); mpfr_init(r17856); mpfr_init(r17857); mpfr_init(r17858); mpfr_init(r17859); mpfr_init(r17860); mpfr_init(r17861); mpfr_init(r17862); mpfr_init(r17863); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r17849, r17847, r17848, MPFR_RNDN); ; mpfr_set_d(r17851, c, MPFR_RNDN); mpfr_mul(r17852, r17850, r17851, MPFR_RNDN); mpfr_set_d(r17853, b, MPFR_RNDN); mpfr_neg(r17854, r17853, MPFR_RNDN); mpfr_sqr(r17855, r17853, MPFR_RNDN); mpfr_set_d(r17856, a, MPFR_RNDN); mpfr_mul(r17857, r17851, r17856, MPFR_RNDN); mpfr_mul(r17858, r17857, r17850, MPFR_RNDN); mpfr_sub(r17859, r17855, r17858, MPFR_RNDN); mpfr_sqrt(r17860, r17859, MPFR_RNDN); mpfr_sub(r17861, r17854, r17860, MPFR_RNDN); mpfr_div(r17862, r17852, r17861, MPFR_RNDN); mpfr_mul(r17863, r17849, r17862, MPFR_RNDN); return mpfr_get_d(r17863, MPFR_RNDN); } static mpfr_t r17864, r17865, r17866, r17867, r17868, r17869, r17870, r17871, r17872, r17873, r17874, r17875, r17876, r17877, r17878, r17879, r17880; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17864, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17865, "2", 10, MPFR_RNDN); mpfr_init(r17866); mpfr_init_set_str(r17867, "4", 10, MPFR_RNDN); mpfr_init(r17868); mpfr_init(r17869); mpfr_init(r17870); mpfr_init(r17871); mpfr_init(r17872); mpfr_init(r17873); mpfr_init(r17874); mpfr_init(r17875); mpfr_init(r17876); mpfr_init(r17877); mpfr_init(r17878); mpfr_init(r17879); mpfr_init(r17880); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r17866, r17864, r17865, MPFR_RNDN); ; mpfr_set_d(r17868, c, MPFR_RNDN); mpfr_mul(r17869, r17867, r17868, MPFR_RNDN); mpfr_set_d(r17870, b, MPFR_RNDN); mpfr_neg(r17871, r17870, MPFR_RNDN); mpfr_sqr(r17872, r17870, MPFR_RNDN); mpfr_set_d(r17873, a, MPFR_RNDN); mpfr_mul(r17874, r17868, r17873, MPFR_RNDN); mpfr_mul(r17875, r17874, r17867, MPFR_RNDN); mpfr_sub(r17876, r17872, r17875, MPFR_RNDN); mpfr_sqrt(r17877, r17876, MPFR_RNDN); mpfr_sub(r17878, r17871, r17877, MPFR_RNDN); mpfr_div(r17879, r17869, r17878, MPFR_RNDN); mpfr_mul(r17880, r17866, r17879, MPFR_RNDN); return mpfr_get_d(r17880, MPFR_RNDN); }