#include #include #include #include #include char *name = "Cubic critical, medium range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r26136 = b; float r26137 = -r26136; float r26138 = r26136 * r26136; float r26139 = 3; float r26140 = a; float r26141 = r26139 * r26140; float r26142 = c; float r26143 = r26141 * r26142; float r26144 = r26138 - r26143; float r26145 = sqrt(r26144); float r26146 = r26137 + r26145; float r26147 = r26146 / r26141; return r26147; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26148 = b; double r26149 = -r26148; double r26150 = r26148 * r26148; double r26151 = 3; double r26152 = a; double r26153 = r26151 * r26152; double r26154 = c; double r26155 = r26153 * r26154; double r26156 = r26150 - r26155; double r26157 = sqrt(r26156); double r26158 = r26149 + r26157; double r26159 = r26158 / r26153; return r26159; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26160 = c; float r26161 = -r26160; float r26162 = a; float r26163 = r26161 * r26162; float r26164 = 1; float r26165 = 3; float r26166 = r26165 * r26162; float r26167 = b; float r26168 = r26167 * r26167; float r26169 = fma(r26166, r26161, r26168); float r26170 = sqrt(r26169); float r26171 = r26170 + r26167; float r26172 = r26164 / r26171; float r26173 = r26172 / r26162; float r26174 = r26163 * r26173; return r26174; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26175 = c; double r26176 = -r26175; double r26177 = a; double r26178 = r26176 * r26177; double r26179 = 1; double r26180 = 3; double r26181 = r26180 * r26177; double r26182 = b; double r26183 = r26182 * r26182; double r26184 = fma(r26181, r26176, r26183); double r26185 = sqrt(r26184); double r26186 = r26185 + r26182; double r26187 = r26179 / r26186; double r26188 = r26187 / r26177; double r26189 = r26178 * r26188; return r26189; } 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 r26190, r26191, r26192, r26193, r26194, r26195, r26196, r26197, r26198, r26199, r26200, r26201; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26190); mpfr_init(r26191); mpfr_init(r26192); mpfr_init_set_str(r26193, "3", 10, MPFR_RNDN); mpfr_init(r26194); mpfr_init(r26195); mpfr_init(r26196); mpfr_init(r26197); mpfr_init(r26198); mpfr_init(r26199); mpfr_init(r26200); mpfr_init(r26201); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26190, b, MPFR_RNDN); mpfr_neg(r26191, r26190, MPFR_RNDN); mpfr_mul(r26192, r26190, r26190, MPFR_RNDN); ; mpfr_set_d(r26194, a, MPFR_RNDN); mpfr_mul(r26195, r26193, r26194, MPFR_RNDN); mpfr_set_d(r26196, c, MPFR_RNDN); mpfr_mul(r26197, r26195, r26196, MPFR_RNDN); mpfr_sub(r26198, r26192, r26197, MPFR_RNDN); mpfr_sqrt(r26199, r26198, MPFR_RNDN); mpfr_add(r26200, r26191, r26199, MPFR_RNDN); mpfr_div(r26201, r26200, r26195, MPFR_RNDN); return mpfr_get_d(r26201, MPFR_RNDN); } static mpfr_t r26202, r26203, r26204, r26205, r26206, r26207, r26208, r26209, r26210, r26211, r26212, r26213, r26214, r26215, r26216; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26202); mpfr_init(r26203); mpfr_init(r26204); mpfr_init(r26205); mpfr_init_set_str(r26206, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26207, "3", 10, MPFR_RNDN); mpfr_init(r26208); mpfr_init(r26209); mpfr_init(r26210); mpfr_init(r26211); mpfr_init(r26212); mpfr_init(r26213); mpfr_init(r26214); mpfr_init(r26215); mpfr_init(r26216); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26202, c, MPFR_RNDN); mpfr_neg(r26203, r26202, MPFR_RNDN); mpfr_set_d(r26204, a, MPFR_RNDN); mpfr_mul(r26205, r26203, r26204, MPFR_RNDN); ; ; mpfr_mul(r26208, r26207, r26204, MPFR_RNDN); mpfr_set_d(r26209, b, MPFR_RNDN); mpfr_mul(r26210, r26209, r26209, MPFR_RNDN); mpfr_fma(r26211, r26208, r26203, r26210, MPFR_RNDN); mpfr_sqrt(r26212, r26211, MPFR_RNDN); mpfr_add(r26213, r26212, r26209, MPFR_RNDN); mpfr_div(r26214, r26206, r26213, MPFR_RNDN); mpfr_div(r26215, r26214, r26204, MPFR_RNDN); mpfr_mul(r26216, r26205, r26215, MPFR_RNDN); return mpfr_get_d(r26216, MPFR_RNDN); } static mpfr_t r26217, r26218, r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229, r26230, r26231; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26217); mpfr_init(r26218); mpfr_init(r26219); mpfr_init(r26220); mpfr_init_set_str(r26221, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26222, "3", 10, MPFR_RNDN); mpfr_init(r26223); mpfr_init(r26224); mpfr_init(r26225); mpfr_init(r26226); mpfr_init(r26227); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init(r26231); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26217, c, MPFR_RNDN); mpfr_neg(r26218, r26217, MPFR_RNDN); mpfr_set_d(r26219, a, MPFR_RNDN); mpfr_mul(r26220, r26218, r26219, MPFR_RNDN); ; ; mpfr_mul(r26223, r26222, r26219, MPFR_RNDN); mpfr_set_d(r26224, b, MPFR_RNDN); mpfr_mul(r26225, r26224, r26224, MPFR_RNDN); mpfr_fma(r26226, r26223, r26218, r26225, MPFR_RNDN); mpfr_sqrt(r26227, r26226, MPFR_RNDN); mpfr_add(r26228, r26227, r26224, MPFR_RNDN); mpfr_div(r26229, r26221, r26228, MPFR_RNDN); mpfr_div(r26230, r26229, r26219, MPFR_RNDN); mpfr_mul(r26231, r26220, r26230, MPFR_RNDN); return mpfr_get_d(r26231, MPFR_RNDN); }