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