#include #include #include #include #include char *name = "jeff quadratic root 1"; double f_if(float a, float b, float c) { float r15197 = b; float r15198 = 0.0f; bool r15199 = r15197 >= r15198; float r15200 = -r15197; float r15201 = r15197 * r15197; float r15202 = 4.0f; float r15203 = a; float r15204 = r15202 * r15203; float r15205 = c; float r15206 = r15204 * r15205; float r15207 = r15201 - r15206; float r15208 = sqrt(r15207); float r15209 = r15200 - r15208; float r15210 = 2.0f; float r15211 = r15210 * r15203; float r15212 = r15209 / r15211; float r15213 = r15210 * r15205; float r15214 = r15200 + r15208; float r15215 = r15213 / r15214; float r15216 = r15199 ? r15212 : r15215; return r15216; } double f_id(double a, double b, double c) { double r15217 = b; double r15218 = 0.0; bool r15219 = r15217 >= r15218; double r15220 = -r15217; double r15221 = r15217 * r15217; double r15222 = 4.0; double r15223 = a; double r15224 = r15222 * r15223; double r15225 = c; double r15226 = r15224 * r15225; double r15227 = r15221 - r15226; double r15228 = sqrt(r15227); double r15229 = r15220 - r15228; double r15230 = 2.0; double r15231 = r15230 * r15223; double r15232 = r15229 / r15231; double r15233 = r15230 * r15225; double r15234 = r15220 + r15228; double r15235 = r15233 / r15234; double r15236 = r15219 ? r15232 : r15235; return r15236; } double f_of(float a, float b, float c) { float r15237 = b; float r15238 = 0.0f; bool r15239 = r15237 >= r15238; float r15240 = -r15237; float r15241 = r15237 * r15237; float r15242 = 4.0f; float r15243 = a; float r15244 = r15242 * r15243; float r15245 = c; float r15246 = r15244 * r15245; float r15247 = r15241 - r15246; float r15248 = sqrt(r15247); float r15249 = r15240 - r15248; float r15250 = 2.0f; float r15251 = r15250 * r15243; float r15252 = r15249 / r15251; float r15253 = r15250 * r15245; float r15254 = sqrt(r15248); float r15255 = r15254 * r15254; float r15256 = r15240 + r15255; float r15257 = r15253 / r15256; float r15258 = r15239 ? r15252 : r15257; return r15258; } double f_od(double a, double b, double c) { double r15259 = b; double r15260 = 0.0; bool r15261 = r15259 >= r15260; double r15262 = -r15259; double r15263 = r15259 * r15259; double r15264 = 4.0; double r15265 = a; double r15266 = r15264 * r15265; double r15267 = c; double r15268 = r15266 * r15267; double r15269 = r15263 - r15268; double r15270 = sqrt(r15269); double r15271 = r15262 - r15270; double r15272 = 2.0; double r15273 = r15272 * r15265; double r15274 = r15271 / r15273; double r15275 = r15272 * r15267; double r15276 = sqrt(r15270); double r15277 = r15276 * r15276; double r15278 = r15262 + r15277; double r15279 = r15275 / r15278; double r15280 = r15261 ? r15274 : r15279; return r15280; } 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 r15281, r15282, r15283, r15284, r15285, r15286, r15287, r15288, r15289, r15290, r15291, r15292, r15293, r15294, r15295, r15296, r15297, r15298, r15299, r15300; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15281); mpfr_init_set_str(r15282, "0", 10, MPFR_RNDN); mpfr_init(r15283); mpfr_init(r15284); mpfr_init(r15285); mpfr_init_set_str(r15286, "4", 10, MPFR_RNDN); mpfr_init(r15287); mpfr_init(r15288); mpfr_init(r15289); mpfr_init(r15290); mpfr_init(r15291); mpfr_init(r15292); mpfr_init(r15293); mpfr_init_set_str(r15294, "2", 10, MPFR_RNDN); mpfr_init(r15295); mpfr_init(r15296); mpfr_init(r15297); mpfr_init(r15298); mpfr_init(r15299); mpfr_init(r15300); } double f_im(double a, double b, double c) { mpfr_set_d(r15281, b, MPFR_RNDN); ; mpfr_set_si(r15283, mpfr_cmp(r15281, r15282) >= 0, MPFR_RNDN); mpfr_neg(r15284, r15281, MPFR_RNDN); mpfr_sqr(r15285, r15281, MPFR_RNDN); ; mpfr_set_d(r15287, a, MPFR_RNDN); mpfr_mul(r15288, r15286, r15287, MPFR_RNDN); mpfr_set_d(r15289, c, MPFR_RNDN); mpfr_mul(r15290, r15288, r15289, MPFR_RNDN); mpfr_sub(r15291, r15285, r15290, MPFR_RNDN); mpfr_sqrt(r15292, r15291, MPFR_RNDN); mpfr_sub(r15293, r15284, r15292, MPFR_RNDN); ; mpfr_mul(r15295, r15294, r15287, MPFR_RNDN); mpfr_div(r15296, r15293, r15295, MPFR_RNDN); mpfr_mul(r15297, r15294, r15289, MPFR_RNDN); mpfr_add(r15298, r15284, r15292, MPFR_RNDN); mpfr_div(r15299, r15297, r15298, MPFR_RNDN); if (mpfr_get_si(r15283, MPFR_RNDN)) { mpfr_set(r15300, r15296, MPFR_RNDN); } else { mpfr_set(r15300, r15299, MPFR_RNDN); }; return mpfr_get_d(r15300, MPFR_RNDN); } static mpfr_t r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310, r15311, r15312, r15313, r15314, r15315, r15316, r15317, r15318, r15319, r15320, r15321, r15322; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15301); mpfr_init_set_str(r15302, "0", 10, MPFR_RNDN); mpfr_init(r15303); mpfr_init(r15304); mpfr_init(r15305); mpfr_init_set_str(r15306, "4", 10, MPFR_RNDN); mpfr_init(r15307); mpfr_init(r15308); mpfr_init(r15309); mpfr_init(r15310); mpfr_init(r15311); mpfr_init(r15312); mpfr_init(r15313); mpfr_init_set_str(r15314, "2", 10, MPFR_RNDN); mpfr_init(r15315); mpfr_init(r15316); mpfr_init(r15317); mpfr_init(r15318); mpfr_init(r15319); mpfr_init(r15320); mpfr_init(r15321); mpfr_init(r15322); } double f_fm(double a, double b, double c) { mpfr_set_d(r15301, b, MPFR_RNDN); ; mpfr_set_si(r15303, mpfr_cmp(r15301, r15302) >= 0, MPFR_RNDN); mpfr_neg(r15304, r15301, MPFR_RNDN); mpfr_sqr(r15305, r15301, MPFR_RNDN); ; mpfr_set_d(r15307, a, MPFR_RNDN); mpfr_mul(r15308, r15306, r15307, MPFR_RNDN); mpfr_set_d(r15309, c, MPFR_RNDN); mpfr_mul(r15310, r15308, r15309, MPFR_RNDN); mpfr_sub(r15311, r15305, r15310, MPFR_RNDN); mpfr_sqrt(r15312, r15311, MPFR_RNDN); mpfr_sub(r15313, r15304, r15312, MPFR_RNDN); ; mpfr_mul(r15315, r15314, r15307, MPFR_RNDN); mpfr_div(r15316, r15313, r15315, MPFR_RNDN); mpfr_mul(r15317, r15314, r15309, MPFR_RNDN); mpfr_sqrt(r15318, r15312, MPFR_RNDN); mpfr_sqr(r15319, r15318, MPFR_RNDN); mpfr_add(r15320, r15304, r15319, MPFR_RNDN); mpfr_div(r15321, r15317, r15320, MPFR_RNDN); if (mpfr_get_si(r15303, MPFR_RNDN)) { mpfr_set(r15322, r15316, MPFR_RNDN); } else { mpfr_set(r15322, r15321, MPFR_RNDN); }; return mpfr_get_d(r15322, MPFR_RNDN); } static mpfr_t r15323, r15324, r15325, r15326, r15327, r15328, r15329, r15330, r15331, r15332, r15333, r15334, r15335, r15336, r15337, r15338, r15339, r15340, r15341, r15342, r15343, r15344; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15323); mpfr_init_set_str(r15324, "0", 10, MPFR_RNDN); mpfr_init(r15325); mpfr_init(r15326); mpfr_init(r15327); mpfr_init_set_str(r15328, "4", 10, MPFR_RNDN); mpfr_init(r15329); mpfr_init(r15330); mpfr_init(r15331); mpfr_init(r15332); mpfr_init(r15333); mpfr_init(r15334); mpfr_init(r15335); mpfr_init_set_str(r15336, "2", 10, MPFR_RNDN); mpfr_init(r15337); mpfr_init(r15338); mpfr_init(r15339); mpfr_init(r15340); mpfr_init(r15341); mpfr_init(r15342); mpfr_init(r15343); mpfr_init(r15344); } double f_dm(double a, double b, double c) { mpfr_set_d(r15323, b, MPFR_RNDN); ; mpfr_set_si(r15325, mpfr_cmp(r15323, r15324) >= 0, MPFR_RNDN); mpfr_neg(r15326, r15323, MPFR_RNDN); mpfr_sqr(r15327, r15323, MPFR_RNDN); ; mpfr_set_d(r15329, a, MPFR_RNDN); mpfr_mul(r15330, r15328, r15329, MPFR_RNDN); mpfr_set_d(r15331, c, MPFR_RNDN); mpfr_mul(r15332, r15330, r15331, MPFR_RNDN); mpfr_sub(r15333, r15327, r15332, MPFR_RNDN); mpfr_sqrt(r15334, r15333, MPFR_RNDN); mpfr_sub(r15335, r15326, r15334, MPFR_RNDN); ; mpfr_mul(r15337, r15336, r15329, MPFR_RNDN); mpfr_div(r15338, r15335, r15337, MPFR_RNDN); mpfr_mul(r15339, r15336, r15331, MPFR_RNDN); mpfr_sqrt(r15340, r15334, MPFR_RNDN); mpfr_sqr(r15341, r15340, MPFR_RNDN); mpfr_add(r15342, r15326, r15341, MPFR_RNDN); mpfr_div(r15343, r15339, r15342, MPFR_RNDN); if (mpfr_get_si(r15325, MPFR_RNDN)) { mpfr_set(r15344, r15338, MPFR_RNDN); } else { mpfr_set(r15344, r15343, MPFR_RNDN); }; return mpfr_get_d(r15344, MPFR_RNDN); }