#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r22152 = b; float r22153 = -r22152; float r22154 = r22152 * r22152; float r22155 = 4; float r22156 = a; float r22157 = r22155 * r22156; float r22158 = c; float r22159 = r22157 * r22158; float r22160 = r22154 - r22159; float r22161 = sqrt(r22160); float r22162 = r22153 + r22161; float r22163 = 2; float r22164 = r22163 * r22156; float r22165 = r22162 / r22164; return r22165; } double f_id(double a, double b, double c) { double r22166 = b; double r22167 = -r22166; double r22168 = r22166 * r22166; double r22169 = 4; double r22170 = a; double r22171 = r22169 * r22170; double r22172 = c; double r22173 = r22171 * r22172; double r22174 = r22168 - r22173; double r22175 = sqrt(r22174); double r22176 = r22167 + r22175; double r22177 = 2; double r22178 = r22177 * r22170; double r22179 = r22176 / r22178; return r22179; } double f_of(float a, float b, float c) { float r22180 = c; float r22181 = 2; float r22182 = r22180 / r22181; float r22183 = 4; float r22184 = b; float r22185 = -r22184; float r22186 = -r22183; float r22187 = a; float r22188 = r22180 * r22187; float r22189 = r22184 * r22184; float r22190 = fma(r22186, r22188, r22189); float r22191 = sqrt(r22190); float r22192 = r22185 - r22191; float r22193 = r22183 / r22192; float r22194 = r22182 * r22193; return r22194; } double f_od(double a, double b, double c) { double r22195 = c; double r22196 = 2; double r22197 = r22195 / r22196; double r22198 = 4; double r22199 = b; double r22200 = -r22199; double r22201 = -r22198; double r22202 = a; double r22203 = r22195 * r22202; double r22204 = r22199 * r22199; double r22205 = fma(r22201, r22203, r22204); double r22206 = sqrt(r22205); double r22207 = r22200 - r22206; double r22208 = r22198 / r22207; double r22209 = r22197 * r22208; return r22209; } 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 r22210, r22211, r22212, r22213, r22214, r22215, r22216, r22217, r22218, r22219, r22220, r22221, r22222, r22223; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22210); mpfr_init(r22211); mpfr_init(r22212); mpfr_init_set_str(r22213, "4", 10, MPFR_RNDN); mpfr_init(r22214); mpfr_init(r22215); mpfr_init(r22216); mpfr_init(r22217); mpfr_init(r22218); mpfr_init(r22219); mpfr_init(r22220); mpfr_init_set_str(r22221, "2", 10, MPFR_RNDN); mpfr_init(r22222); mpfr_init(r22223); } double f_im(double a, double b, double c) { mpfr_set_d(r22210, b, MPFR_RNDN); mpfr_neg(r22211, r22210, MPFR_RNDN); mpfr_mul(r22212, r22210, r22210, MPFR_RNDN); ; mpfr_set_d(r22214, a, MPFR_RNDN); mpfr_mul(r22215, r22213, r22214, MPFR_RNDN); mpfr_set_d(r22216, c, MPFR_RNDN); mpfr_mul(r22217, r22215, r22216, MPFR_RNDN); mpfr_sub(r22218, r22212, r22217, MPFR_RNDN); mpfr_sqrt(r22219, r22218, MPFR_RNDN); mpfr_add(r22220, r22211, r22219, MPFR_RNDN); ; mpfr_mul(r22222, r22221, r22214, MPFR_RNDN); mpfr_div(r22223, r22220, r22222, MPFR_RNDN); return mpfr_get_d(r22223, MPFR_RNDN); } static mpfr_t r22224, r22225, r22226, r22227, r22228, r22229, r22230, r22231, r22232, r22233, r22234, r22235, r22236, r22237, r22238; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r22224); mpfr_init_set_str(r22225, "2", 10, MPFR_RNDN); mpfr_init(r22226); mpfr_init_set_str(r22227, "4", 10, MPFR_RNDN); mpfr_init(r22228); mpfr_init(r22229); mpfr_init(r22230); mpfr_init(r22231); mpfr_init(r22232); mpfr_init(r22233); mpfr_init(r22234); mpfr_init(r22235); mpfr_init(r22236); mpfr_init(r22237); mpfr_init(r22238); } double f_fm(double a, double b, double c) { mpfr_set_d(r22224, c, MPFR_RNDN); ; mpfr_div(r22226, r22224, r22225, MPFR_RNDN); ; mpfr_set_d(r22228, b, MPFR_RNDN); mpfr_neg(r22229, r22228, MPFR_RNDN); mpfr_neg(r22230, r22227, MPFR_RNDN); mpfr_set_d(r22231, a, MPFR_RNDN); mpfr_mul(r22232, r22224, r22231, MPFR_RNDN); mpfr_mul(r22233, r22228, r22228, MPFR_RNDN); mpfr_fma(r22234, r22230, r22232, r22233, MPFR_RNDN); mpfr_sqrt(r22235, r22234, MPFR_RNDN); mpfr_sub(r22236, r22229, r22235, MPFR_RNDN); mpfr_div(r22237, r22227, r22236, MPFR_RNDN); mpfr_mul(r22238, r22226, r22237, MPFR_RNDN); return mpfr_get_d(r22238, MPFR_RNDN); } static mpfr_t r22239, r22240, r22241, r22242, r22243, r22244, r22245, r22246, r22247, r22248, r22249, r22250, r22251, r22252, r22253; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22239); mpfr_init_set_str(r22240, "2", 10, MPFR_RNDN); mpfr_init(r22241); mpfr_init_set_str(r22242, "4", 10, MPFR_RNDN); mpfr_init(r22243); mpfr_init(r22244); mpfr_init(r22245); mpfr_init(r22246); mpfr_init(r22247); mpfr_init(r22248); mpfr_init(r22249); mpfr_init(r22250); mpfr_init(r22251); mpfr_init(r22252); mpfr_init(r22253); } double f_dm(double a, double b, double c) { mpfr_set_d(r22239, c, MPFR_RNDN); ; mpfr_div(r22241, r22239, r22240, MPFR_RNDN); ; mpfr_set_d(r22243, b, MPFR_RNDN); mpfr_neg(r22244, r22243, MPFR_RNDN); mpfr_neg(r22245, r22242, MPFR_RNDN); mpfr_set_d(r22246, a, MPFR_RNDN); mpfr_mul(r22247, r22239, r22246, MPFR_RNDN); mpfr_mul(r22248, r22243, r22243, MPFR_RNDN); mpfr_fma(r22249, r22245, r22247, r22248, MPFR_RNDN); mpfr_sqrt(r22250, r22249, MPFR_RNDN); mpfr_sub(r22251, r22244, r22250, MPFR_RNDN); mpfr_div(r22252, r22242, r22251, MPFR_RNDN); mpfr_mul(r22253, r22241, r22252, MPFR_RNDN); return mpfr_get_d(r22253, MPFR_RNDN); }