#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r24344 = b; float r24345 = -r24344; float r24346 = r24344 * r24344; float r24347 = 4; float r24348 = a; float r24349 = r24347 * r24348; float r24350 = c; float r24351 = r24349 * r24350; float r24352 = r24346 - r24351; float r24353 = sqrt(r24352); float r24354 = r24345 + r24353; float r24355 = 2; float r24356 = r24355 * r24348; float r24357 = r24354 / r24356; return r24357; } double f_id(double a, double b, double c) { double r24358 = b; double r24359 = -r24358; double r24360 = r24358 * r24358; double r24361 = 4; double r24362 = a; double r24363 = r24361 * r24362; double r24364 = c; double r24365 = r24363 * r24364; double r24366 = r24360 - r24365; double r24367 = sqrt(r24366); double r24368 = r24359 + r24367; double r24369 = 2; double r24370 = r24369 * r24362; double r24371 = r24368 / r24370; return r24371; } double f_of(float a, float b, float c) { float r24372 = a; float r24373 = 2; float r24374 = r24372 / r24373; float r24375 = c; float r24376 = r24374 * r24375; float r24377 = 4; float r24378 = -r24377; float r24379 = r24377 * r24372; float r24380 = -r24375; float r24381 = b; float r24382 = r24381 * r24381; float r24383 = fma(r24379, r24380, r24382); float r24384 = sqrt(r24383); float r24385 = r24384 + r24381; float r24386 = r24378 / r24385; float r24387 = r24386 / r24372; float r24388 = r24376 * r24387; return r24388; } double f_od(double a, double b, double c) { double r24389 = a; double r24390 = 2; double r24391 = r24389 / r24390; double r24392 = c; double r24393 = r24391 * r24392; double r24394 = 4; double r24395 = -r24394; double r24396 = r24394 * r24389; double r24397 = -r24392; double r24398 = b; double r24399 = r24398 * r24398; double r24400 = fma(r24396, r24397, r24399); double r24401 = sqrt(r24400); double r24402 = r24401 + r24398; double r24403 = r24395 / r24402; double r24404 = r24403 / r24389; double r24405 = r24393 * r24404; return r24405; } 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 r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24406); mpfr_init(r24407); mpfr_init(r24408); mpfr_init_set_str(r24409, "4", 10, MPFR_RNDN); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init(r24414); mpfr_init(r24415); mpfr_init(r24416); mpfr_init_set_str(r24417, "2", 10, MPFR_RNDN); mpfr_init(r24418); mpfr_init(r24419); } double f_im(double a, double b, double c) { mpfr_set_d(r24406, b, MPFR_RNDN); mpfr_neg(r24407, r24406, MPFR_RNDN); mpfr_mul(r24408, r24406, r24406, MPFR_RNDN); ; mpfr_set_d(r24410, a, MPFR_RNDN); mpfr_mul(r24411, r24409, r24410, MPFR_RNDN); mpfr_set_d(r24412, c, MPFR_RNDN); mpfr_mul(r24413, r24411, r24412, MPFR_RNDN); mpfr_sub(r24414, r24408, r24413, MPFR_RNDN); mpfr_sqrt(r24415, r24414, MPFR_RNDN); mpfr_add(r24416, r24407, r24415, MPFR_RNDN); ; mpfr_mul(r24418, r24417, r24410, MPFR_RNDN); mpfr_div(r24419, r24416, r24418, MPFR_RNDN); return mpfr_get_d(r24419, MPFR_RNDN); } static mpfr_t r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24420); mpfr_init_set_str(r24421, "2", 10, MPFR_RNDN); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init_set_str(r24425, "4", 10, MPFR_RNDN); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); } double f_fm(double a, double b, double c) { mpfr_set_d(r24420, a, MPFR_RNDN); ; mpfr_div(r24422, r24420, r24421, MPFR_RNDN); mpfr_set_d(r24423, c, MPFR_RNDN); mpfr_mul(r24424, r24422, r24423, MPFR_RNDN); ; mpfr_neg(r24426, r24425, MPFR_RNDN); mpfr_mul(r24427, r24425, r24420, MPFR_RNDN); mpfr_neg(r24428, r24423, MPFR_RNDN); mpfr_set_d(r24429, b, MPFR_RNDN); mpfr_mul(r24430, r24429, r24429, MPFR_RNDN); mpfr_fma(r24431, r24427, r24428, r24430, MPFR_RNDN); mpfr_sqrt(r24432, r24431, MPFR_RNDN); mpfr_add(r24433, r24432, r24429, MPFR_RNDN); mpfr_div(r24434, r24426, r24433, MPFR_RNDN); mpfr_div(r24435, r24434, r24420, MPFR_RNDN); mpfr_mul(r24436, r24424, r24435, MPFR_RNDN); return mpfr_get_d(r24436, MPFR_RNDN); } static mpfr_t r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24437); mpfr_init_set_str(r24438, "2", 10, MPFR_RNDN); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init_set_str(r24442, "4", 10, MPFR_RNDN); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); } double f_dm(double a, double b, double c) { mpfr_set_d(r24437, a, MPFR_RNDN); ; mpfr_div(r24439, r24437, r24438, MPFR_RNDN); mpfr_set_d(r24440, c, MPFR_RNDN); mpfr_mul(r24441, r24439, r24440, MPFR_RNDN); ; mpfr_neg(r24443, r24442, MPFR_RNDN); mpfr_mul(r24444, r24442, r24437, MPFR_RNDN); mpfr_neg(r24445, r24440, MPFR_RNDN); mpfr_set_d(r24446, b, MPFR_RNDN); mpfr_mul(r24447, r24446, r24446, MPFR_RNDN); mpfr_fma(r24448, r24444, r24445, r24447, MPFR_RNDN); mpfr_sqrt(r24449, r24448, MPFR_RNDN); mpfr_add(r24450, r24449, r24446, MPFR_RNDN); mpfr_div(r24451, r24443, r24450, MPFR_RNDN); mpfr_div(r24452, r24451, r24437, MPFR_RNDN); mpfr_mul(r24453, r24441, r24452, MPFR_RNDN); return mpfr_get_d(r24453, MPFR_RNDN); }