#include #include #include #include #include char *name = "jeff quadratic root 1"; double f_if(float a, float b, float c) { float r15101 = b; float r15102 = 0.0f; bool r15103 = r15101 >= r15102; float r15104 = -r15101; float r15105 = r15101 * r15101; float r15106 = 4.0f; float r15107 = a; float r15108 = r15106 * r15107; float r15109 = c; float r15110 = r15108 * r15109; float r15111 = r15105 - r15110; float r15112 = sqrt(r15111); float r15113 = r15104 - r15112; float r15114 = 2.0f; float r15115 = r15114 * r15107; float r15116 = r15113 / r15115; float r15117 = r15114 * r15109; float r15118 = r15104 + r15112; float r15119 = r15117 / r15118; float r15120 = r15103 ? r15116 : r15119; return r15120; } double f_id(double a, double b, double c) { double r15121 = b; double r15122 = 0.0; bool r15123 = r15121 >= r15122; double r15124 = -r15121; double r15125 = r15121 * r15121; double r15126 = 4.0; double r15127 = a; double r15128 = r15126 * r15127; double r15129 = c; double r15130 = r15128 * r15129; double r15131 = r15125 - r15130; double r15132 = sqrt(r15131); double r15133 = r15124 - r15132; double r15134 = 2.0; double r15135 = r15134 * r15127; double r15136 = r15133 / r15135; double r15137 = r15134 * r15129; double r15138 = r15124 + r15132; double r15139 = r15137 / r15138; double r15140 = r15123 ? r15136 : r15139; return r15140; } double f_of(float a, float b, float c) { float r15141 = b; float r15142 = 0.0f; bool r15143 = r15141 >= r15142; float r15144 = a; float r15145 = r15141 / r15144; float r15146 = -2.0f; float r15147 = 2.0f; float r15148 = r15146 / r15147; float r15149 = r15145 * r15148; float r15150 = c; float r15151 = r15147 * r15150; float r15152 = -r15141; float r15153 = r15141 * r15141; float r15154 = 4.0f; float r15155 = r15144 * r15150; float r15156 = r15154 * r15155; float r15157 = r15153 - r15156; float r15158 = sqrt(r15157); float r15159 = r15152 + r15158; float r15160 = r15151 / r15159; float r15161 = r15143 ? r15149 : r15160; return r15161; } double f_od(double a, double b, double c) { double r15162 = b; double r15163 = 0.0; bool r15164 = r15162 >= r15163; double r15165 = a; double r15166 = r15162 / r15165; double r15167 = -2.0; double r15168 = 2.0; double r15169 = r15167 / r15168; double r15170 = r15166 * r15169; double r15171 = c; double r15172 = r15168 * r15171; double r15173 = -r15162; double r15174 = r15162 * r15162; double r15175 = 4.0; double r15176 = r15165 * r15171; double r15177 = r15175 * r15176; double r15178 = r15174 - r15177; double r15179 = sqrt(r15178); double r15180 = r15173 + r15179; double r15181 = r15172 / r15180; double r15182 = r15164 ? r15170 : r15181; return r15182; } 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 r15183, r15184, r15185, r15186, r15187, r15188, r15189, r15190, r15191, r15192, r15193, r15194, r15195, r15196, r15197, r15198, r15199, r15200, r15201, r15202; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15183); mpfr_init_set_str(r15184, "0", 10, MPFR_RNDN); mpfr_init(r15185); mpfr_init(r15186); mpfr_init(r15187); mpfr_init_set_str(r15188, "4", 10, MPFR_RNDN); mpfr_init(r15189); mpfr_init(r15190); mpfr_init(r15191); mpfr_init(r15192); mpfr_init(r15193); mpfr_init(r15194); mpfr_init(r15195); mpfr_init_set_str(r15196, "2", 10, MPFR_RNDN); mpfr_init(r15197); mpfr_init(r15198); mpfr_init(r15199); mpfr_init(r15200); mpfr_init(r15201); mpfr_init(r15202); } double f_im(double a, double b, double c) { mpfr_set_d(r15183, b, MPFR_RNDN); ; mpfr_set_si(r15185, mpfr_cmp(r15183, r15184) >= 0, MPFR_RNDN); mpfr_neg(r15186, r15183, MPFR_RNDN); mpfr_sqr(r15187, r15183, MPFR_RNDN); ; mpfr_set_d(r15189, a, MPFR_RNDN); mpfr_mul(r15190, r15188, r15189, MPFR_RNDN); mpfr_set_d(r15191, c, MPFR_RNDN); mpfr_mul(r15192, r15190, r15191, MPFR_RNDN); mpfr_sub(r15193, r15187, r15192, MPFR_RNDN); mpfr_sqrt(r15194, r15193, MPFR_RNDN); mpfr_sub(r15195, r15186, r15194, MPFR_RNDN); ; mpfr_mul(r15197, r15196, r15189, MPFR_RNDN); mpfr_div(r15198, r15195, r15197, MPFR_RNDN); mpfr_mul(r15199, r15196, r15191, MPFR_RNDN); mpfr_add(r15200, r15186, r15194, MPFR_RNDN); mpfr_div(r15201, r15199, r15200, MPFR_RNDN); if (mpfr_get_si(r15185, MPFR_RNDN)) { mpfr_set(r15202, r15198, MPFR_RNDN); } else { mpfr_set(r15202, r15201, MPFR_RNDN); }; return mpfr_get_d(r15202, MPFR_RNDN); } static mpfr_t r15203, r15204, r15205, r15206, r15207, r15208, r15209, r15210, r15211, r15212, r15213, r15214, r15215, r15216, r15217, r15218, r15219, r15220, r15221, r15222, r15223; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15203); mpfr_init_set_str(r15204, "0", 10, MPFR_RNDN); mpfr_init(r15205); mpfr_init(r15206); mpfr_init(r15207); mpfr_init_set_str(r15208, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r15209, "2", 10, MPFR_RNDN); mpfr_init(r15210); mpfr_init(r15211); mpfr_init(r15212); mpfr_init(r15213); mpfr_init(r15214); mpfr_init(r15215); mpfr_init_set_str(r15216, "4", 10, MPFR_RNDN); mpfr_init(r15217); mpfr_init(r15218); mpfr_init(r15219); mpfr_init(r15220); mpfr_init(r15221); mpfr_init(r15222); mpfr_init(r15223); } double f_fm(double a, double b, double c) { mpfr_set_d(r15203, b, MPFR_RNDN); ; mpfr_set_si(r15205, mpfr_cmp(r15203, r15204) >= 0, MPFR_RNDN); mpfr_set_d(r15206, a, MPFR_RNDN); mpfr_div(r15207, r15203, r15206, MPFR_RNDN); ; ; mpfr_div(r15210, r15208, r15209, MPFR_RNDN); mpfr_mul(r15211, r15207, r15210, MPFR_RNDN); mpfr_set_d(r15212, c, MPFR_RNDN); mpfr_mul(r15213, r15209, r15212, MPFR_RNDN); mpfr_neg(r15214, r15203, MPFR_RNDN); mpfr_sqr(r15215, r15203, MPFR_RNDN); ; mpfr_mul(r15217, r15206, r15212, MPFR_RNDN); mpfr_mul(r15218, r15216, r15217, MPFR_RNDN); mpfr_sub(r15219, r15215, r15218, MPFR_RNDN); mpfr_sqrt(r15220, r15219, MPFR_RNDN); mpfr_add(r15221, r15214, r15220, MPFR_RNDN); mpfr_div(r15222, r15213, r15221, MPFR_RNDN); if (mpfr_get_si(r15205, MPFR_RNDN)) { mpfr_set(r15223, r15211, MPFR_RNDN); } else { mpfr_set(r15223, r15222, MPFR_RNDN); }; return mpfr_get_d(r15223, MPFR_RNDN); } static mpfr_t r15224, r15225, r15226, r15227, r15228, r15229, r15230, r15231, r15232, r15233, r15234, r15235, r15236, r15237, r15238, r15239, r15240, r15241, r15242, r15243, r15244; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15224); mpfr_init_set_str(r15225, "0", 10, MPFR_RNDN); mpfr_init(r15226); mpfr_init(r15227); mpfr_init(r15228); mpfr_init_set_str(r15229, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r15230, "2", 10, MPFR_RNDN); mpfr_init(r15231); mpfr_init(r15232); mpfr_init(r15233); mpfr_init(r15234); mpfr_init(r15235); mpfr_init(r15236); mpfr_init_set_str(r15237, "4", 10, MPFR_RNDN); mpfr_init(r15238); mpfr_init(r15239); mpfr_init(r15240); mpfr_init(r15241); mpfr_init(r15242); mpfr_init(r15243); mpfr_init(r15244); } double f_dm(double a, double b, double c) { mpfr_set_d(r15224, b, MPFR_RNDN); ; mpfr_set_si(r15226, mpfr_cmp(r15224, r15225) >= 0, MPFR_RNDN); mpfr_set_d(r15227, a, MPFR_RNDN); mpfr_div(r15228, r15224, r15227, MPFR_RNDN); ; ; mpfr_div(r15231, r15229, r15230, MPFR_RNDN); mpfr_mul(r15232, r15228, r15231, MPFR_RNDN); mpfr_set_d(r15233, c, MPFR_RNDN); mpfr_mul(r15234, r15230, r15233, MPFR_RNDN); mpfr_neg(r15235, r15224, MPFR_RNDN); mpfr_sqr(r15236, r15224, MPFR_RNDN); ; mpfr_mul(r15238, r15227, r15233, MPFR_RNDN); mpfr_mul(r15239, r15237, r15238, MPFR_RNDN); mpfr_sub(r15240, r15236, r15239, MPFR_RNDN); mpfr_sqrt(r15241, r15240, MPFR_RNDN); mpfr_add(r15242, r15235, r15241, MPFR_RNDN); mpfr_div(r15243, r15234, r15242, MPFR_RNDN); if (mpfr_get_si(r15226, MPFR_RNDN)) { mpfr_set(r15244, r15232, MPFR_RNDN); } else { mpfr_set(r15244, r15243, MPFR_RNDN); }; return mpfr_get_d(r15244, MPFR_RNDN); }