#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r24323 = b; float r24324 = -r24323; float r24325 = r24323 * r24323; float r24326 = 4; float r24327 = a; float r24328 = r24326 * r24327; float r24329 = c; float r24330 = r24328 * r24329; float r24331 = r24325 - r24330; float r24332 = sqrt(r24331); float r24333 = r24324 + r24332; float r24334 = 2; float r24335 = r24334 * r24327; float r24336 = r24333 / r24335; return r24336; } double f_id(double a, double b, double c) { double r24337 = b; double r24338 = -r24337; double r24339 = r24337 * r24337; double r24340 = 4; double r24341 = a; double r24342 = r24340 * r24341; double r24343 = c; double r24344 = r24342 * r24343; double r24345 = r24339 - r24344; double r24346 = sqrt(r24345); double r24347 = r24338 + r24346; double r24348 = 2; double r24349 = r24348 * r24341; double r24350 = r24347 / r24349; return r24350; } double f_of(float a, float b, float c) { float r24351 = 1; float r24352 = 2; float r24353 = r24351 / r24352; float r24354 = 4; float r24355 = c; float r24356 = r24354 * r24355; float r24357 = b; float r24358 = -r24357; float r24359 = -r24355; float r24360 = a; float r24361 = r24354 * r24360; float r24362 = r24357 * r24357; float r24363 = fma(r24359, r24361, r24362); float r24364 = sqrt(r24363); float r24365 = r24358 - r24364; float r24366 = r24356 / r24365; float r24367 = r24353 * r24366; return r24367; } double f_od(double a, double b, double c) { double r24368 = 1; double r24369 = 2; double r24370 = r24368 / r24369; double r24371 = 4; double r24372 = c; double r24373 = r24371 * r24372; double r24374 = b; double r24375 = -r24374; double r24376 = -r24372; double r24377 = a; double r24378 = r24371 * r24377; double r24379 = r24374 * r24374; double r24380 = fma(r24376, r24378, r24379); double r24381 = sqrt(r24380); double r24382 = r24375 - r24381; double r24383 = r24373 / r24382; double r24384 = r24370 * r24383; return r24384; } 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 r24385, r24386, r24387, r24388, r24389, r24390, r24391, r24392, r24393, r24394, r24395, r24396, r24397, r24398; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r24385); mpfr_init(r24386); mpfr_init(r24387); mpfr_init_set_str(r24388, "4", 10, MPFR_RNDN); mpfr_init(r24389); mpfr_init(r24390); mpfr_init(r24391); mpfr_init(r24392); mpfr_init(r24393); mpfr_init(r24394); mpfr_init(r24395); mpfr_init_set_str(r24396, "2", 10, MPFR_RNDN); mpfr_init(r24397); mpfr_init(r24398); } double f_im(double a, double b, double c) { mpfr_set_d(r24385, b, MPFR_RNDN); mpfr_neg(r24386, r24385, MPFR_RNDN); mpfr_mul(r24387, r24385, r24385, MPFR_RNDN); ; mpfr_set_d(r24389, a, MPFR_RNDN); mpfr_mul(r24390, r24388, r24389, MPFR_RNDN); mpfr_set_d(r24391, c, MPFR_RNDN); mpfr_mul(r24392, r24390, r24391, MPFR_RNDN); mpfr_sub(r24393, r24387, r24392, MPFR_RNDN); mpfr_sqrt(r24394, r24393, MPFR_RNDN); mpfr_add(r24395, r24386, r24394, MPFR_RNDN); ; mpfr_mul(r24397, r24396, r24389, MPFR_RNDN); mpfr_div(r24398, r24395, r24397, MPFR_RNDN); return mpfr_get_d(r24398, MPFR_RNDN); } static mpfr_t r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414, r24415; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init_set_str(r24399, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24400, "2", 10, MPFR_RNDN); mpfr_init(r24401); mpfr_init_set_str(r24402, "4", 10, MPFR_RNDN); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init(r24406); mpfr_init(r24407); mpfr_init(r24408); mpfr_init(r24409); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init(r24414); mpfr_init(r24415); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r24401, r24399, r24400, MPFR_RNDN); ; mpfr_set_d(r24403, c, MPFR_RNDN); mpfr_mul(r24404, r24402, r24403, MPFR_RNDN); mpfr_set_d(r24405, b, MPFR_RNDN); mpfr_neg(r24406, r24405, MPFR_RNDN); mpfr_neg(r24407, r24403, MPFR_RNDN); mpfr_set_d(r24408, a, MPFR_RNDN); mpfr_mul(r24409, r24402, r24408, MPFR_RNDN); mpfr_mul(r24410, r24405, r24405, MPFR_RNDN); mpfr_fma(r24411, r24407, r24409, r24410, MPFR_RNDN); mpfr_sqrt(r24412, r24411, MPFR_RNDN); mpfr_sub(r24413, r24406, r24412, MPFR_RNDN); mpfr_div(r24414, r24404, r24413, MPFR_RNDN); mpfr_mul(r24415, r24401, r24414, MPFR_RNDN); return mpfr_get_d(r24415, MPFR_RNDN); } static mpfr_t r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init_set_str(r24416, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24417, "2", 10, MPFR_RNDN); mpfr_init(r24418); mpfr_init_set_str(r24419, "4", 10, MPFR_RNDN); mpfr_init(r24420); mpfr_init(r24421); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r24418, r24416, r24417, MPFR_RNDN); ; mpfr_set_d(r24420, c, MPFR_RNDN); mpfr_mul(r24421, r24419, r24420, MPFR_RNDN); mpfr_set_d(r24422, b, MPFR_RNDN); mpfr_neg(r24423, r24422, MPFR_RNDN); mpfr_neg(r24424, r24420, MPFR_RNDN); mpfr_set_d(r24425, a, MPFR_RNDN); mpfr_mul(r24426, r24419, r24425, MPFR_RNDN); mpfr_mul(r24427, r24422, r24422, MPFR_RNDN); mpfr_fma(r24428, r24424, r24426, r24427, MPFR_RNDN); mpfr_sqrt(r24429, r24428, MPFR_RNDN); mpfr_sub(r24430, r24423, r24429, MPFR_RNDN); mpfr_div(r24431, r24421, r24430, MPFR_RNDN); mpfr_mul(r24432, r24418, r24431, MPFR_RNDN); return mpfr_get_d(r24432, MPFR_RNDN); }