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