#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17158 = b; float r17159 = -r17158; float r17160 = r17158 * r17158; float r17161 = 4.0f; float r17162 = a; float r17163 = r17161 * r17162; float r17164 = c; float r17165 = r17163 * r17164; float r17166 = r17160 - r17165; float r17167 = sqrt(r17166); float r17168 = r17159 + r17167; float r17169 = 2.0f; float r17170 = r17169 * r17162; float r17171 = r17168 / r17170; return r17171; } double f_id(double a, double b, double c) { double r17172 = b; double r17173 = -r17172; double r17174 = r17172 * r17172; double r17175 = 4.0; double r17176 = a; double r17177 = r17175 * r17176; double r17178 = c; double r17179 = r17177 * r17178; double r17180 = r17174 - r17179; double r17181 = sqrt(r17180); double r17182 = r17173 + r17181; double r17183 = 2.0; double r17184 = r17183 * r17176; double r17185 = r17182 / r17184; return r17185; } double f_of(float a, float b, float c) { float r17186 = 4.0f; float r17187 = a; float r17188 = r17186 * r17187; float r17189 = c; float r17190 = r17188 * r17189; float r17191 = b; float r17192 = -r17191; float r17193 = r17191 * r17191; float r17194 = r17193 - r17190; float r17195 = sqrt(r17194); float r17196 = r17192 - r17195; float r17197 = r17190 / r17196; float r17198 = 2.0f; float r17199 = r17198 * r17187; float r17200 = r17197 / r17199; return r17200; } double f_od(double a, double b, double c) { double r17201 = 4.0; double r17202 = a; double r17203 = r17201 * r17202; double r17204 = c; double r17205 = r17203 * r17204; double r17206 = b; double r17207 = -r17206; double r17208 = r17206 * r17206; double r17209 = r17208 - r17205; double r17210 = sqrt(r17209); double r17211 = r17207 - r17210; double r17212 = r17205 / r17211; double r17213 = 2.0; double r17214 = r17213 * r17202; double r17215 = r17212 / r17214; return r17215; } 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 r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17216); mpfr_init(r17217); mpfr_init(r17218); mpfr_init_set_str(r17219, "4", 10, MPFR_RNDN); mpfr_init(r17220); mpfr_init(r17221); mpfr_init(r17222); mpfr_init(r17223); mpfr_init(r17224); mpfr_init(r17225); mpfr_init(r17226); mpfr_init_set_str(r17227, "2", 10, MPFR_RNDN); mpfr_init(r17228); mpfr_init(r17229); } double f_im(double a, double b, double c) { mpfr_set_d(r17216, b, MPFR_RNDN); mpfr_neg(r17217, r17216, MPFR_RNDN); mpfr_sqr(r17218, r17216, MPFR_RNDN); ; mpfr_set_d(r17220, a, MPFR_RNDN); mpfr_mul(r17221, r17219, r17220, MPFR_RNDN); mpfr_set_d(r17222, c, MPFR_RNDN); mpfr_mul(r17223, r17221, r17222, MPFR_RNDN); mpfr_sub(r17224, r17218, r17223, MPFR_RNDN); mpfr_sqrt(r17225, r17224, MPFR_RNDN); mpfr_add(r17226, r17217, r17225, MPFR_RNDN); ; mpfr_mul(r17228, r17227, r17220, MPFR_RNDN); mpfr_div(r17229, r17226, r17228, MPFR_RNDN); return mpfr_get_d(r17229, MPFR_RNDN); } static mpfr_t r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240, r17241, r17242, r17243, r17244; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17230, "4", 10, MPFR_RNDN); mpfr_init(r17231); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); mpfr_init(r17235); mpfr_init(r17236); mpfr_init(r17237); mpfr_init(r17238); mpfr_init(r17239); mpfr_init(r17240); mpfr_init(r17241); mpfr_init_set_str(r17242, "2", 10, MPFR_RNDN); mpfr_init(r17243); mpfr_init(r17244); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17231, a, MPFR_RNDN); mpfr_mul(r17232, r17230, r17231, MPFR_RNDN); mpfr_set_d(r17233, c, MPFR_RNDN); mpfr_mul(r17234, r17232, r17233, MPFR_RNDN); mpfr_set_d(r17235, b, MPFR_RNDN); mpfr_neg(r17236, r17235, MPFR_RNDN); mpfr_sqr(r17237, r17235, MPFR_RNDN); mpfr_sub(r17238, r17237, r17234, MPFR_RNDN); mpfr_sqrt(r17239, r17238, MPFR_RNDN); mpfr_sub(r17240, r17236, r17239, MPFR_RNDN); mpfr_div(r17241, r17234, r17240, MPFR_RNDN); ; mpfr_mul(r17243, r17242, r17231, MPFR_RNDN); mpfr_div(r17244, r17241, r17243, MPFR_RNDN); return mpfr_get_d(r17244, MPFR_RNDN); } static mpfr_t r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257, r17258, r17259; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17245, "4", 10, MPFR_RNDN); mpfr_init(r17246); mpfr_init(r17247); mpfr_init(r17248); mpfr_init(r17249); mpfr_init(r17250); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); mpfr_init(r17256); mpfr_init_set_str(r17257, "2", 10, MPFR_RNDN); mpfr_init(r17258); mpfr_init(r17259); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17246, a, MPFR_RNDN); mpfr_mul(r17247, r17245, r17246, MPFR_RNDN); mpfr_set_d(r17248, c, MPFR_RNDN); mpfr_mul(r17249, r17247, r17248, MPFR_RNDN); mpfr_set_d(r17250, b, MPFR_RNDN); mpfr_neg(r17251, r17250, MPFR_RNDN); mpfr_sqr(r17252, r17250, MPFR_RNDN); mpfr_sub(r17253, r17252, r17249, MPFR_RNDN); mpfr_sqrt(r17254, r17253, MPFR_RNDN); mpfr_sub(r17255, r17251, r17254, MPFR_RNDN); mpfr_div(r17256, r17249, r17255, MPFR_RNDN); ; mpfr_mul(r17258, r17257, r17246, MPFR_RNDN); mpfr_div(r17259, r17256, r17258, MPFR_RNDN); return mpfr_get_d(r17259, MPFR_RNDN); }