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