#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r24441 = b; float r24442 = -r24441; float r24443 = r24441 * r24441; float r24444 = 4; float r24445 = a; float r24446 = r24444 * r24445; float r24447 = c; float r24448 = r24446 * r24447; float r24449 = r24443 - r24448; float r24450 = sqrt(r24449); float r24451 = r24442 + r24450; float r24452 = 2; float r24453 = r24452 * r24445; float r24454 = r24451 / r24453; return r24454; } double f_id(double a, double b, double c) { double r24455 = b; double r24456 = -r24455; double r24457 = r24455 * r24455; double r24458 = 4; double r24459 = a; double r24460 = r24458 * r24459; double r24461 = c; double r24462 = r24460 * r24461; double r24463 = r24457 - r24462; double r24464 = sqrt(r24463); double r24465 = r24456 + r24464; double r24466 = 2; double r24467 = r24466 * r24459; double r24468 = r24465 / r24467; return r24468; } double f_of(float a, float b, float c) { float r24469 = c; float r24470 = 2; float r24471 = r24469 / r24470; float r24472 = 4; float r24473 = b; float r24474 = -r24473; float r24475 = r24473 * r24473; float r24476 = a; float r24477 = r24472 * r24476; float r24478 = r24469 * r24477; float r24479 = r24475 - r24478; float r24480 = sqrt(r24479); float r24481 = r24474 - r24480; float r24482 = r24472 / r24481; float r24483 = r24471 * r24482; return r24483; } double f_od(double a, double b, double c) { double r24484 = c; double r24485 = 2; double r24486 = r24484 / r24485; double r24487 = 4; double r24488 = b; double r24489 = -r24488; double r24490 = r24488 * r24488; double r24491 = a; double r24492 = r24487 * r24491; double r24493 = r24484 * r24492; double r24494 = r24490 - r24493; double r24495 = sqrt(r24494); double r24496 = r24489 - r24495; double r24497 = r24487 / r24496; double r24498 = r24486 * r24497; return r24498; } 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 r24499, r24500, r24501, r24502, r24503, r24504, r24505, r24506, r24507, r24508, r24509, r24510, r24511, r24512; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24499); mpfr_init(r24500); mpfr_init(r24501); mpfr_init_set_str(r24502, "4", 10, MPFR_RNDN); mpfr_init(r24503); mpfr_init(r24504); mpfr_init(r24505); mpfr_init(r24506); mpfr_init(r24507); mpfr_init(r24508); mpfr_init(r24509); mpfr_init_set_str(r24510, "2", 10, MPFR_RNDN); mpfr_init(r24511); mpfr_init(r24512); } double f_im(double a, double b, double c) { mpfr_set_d(r24499, b, MPFR_RNDN); mpfr_neg(r24500, r24499, MPFR_RNDN); mpfr_mul(r24501, r24499, r24499, MPFR_RNDN); ; mpfr_set_d(r24503, a, MPFR_RNDN); mpfr_mul(r24504, r24502, r24503, MPFR_RNDN); mpfr_set_d(r24505, c, MPFR_RNDN); mpfr_mul(r24506, r24504, r24505, MPFR_RNDN); mpfr_sub(r24507, r24501, r24506, MPFR_RNDN); mpfr_sqrt(r24508, r24507, MPFR_RNDN); mpfr_add(r24509, r24500, r24508, MPFR_RNDN); ; mpfr_mul(r24511, r24510, r24503, MPFR_RNDN); mpfr_div(r24512, r24509, r24511, MPFR_RNDN); return mpfr_get_d(r24512, MPFR_RNDN); } static mpfr_t r24513, r24514, r24515, r24516, r24517, r24518, r24519, r24520, r24521, r24522, r24523, r24524, r24525, r24526, r24527; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24513); mpfr_init_set_str(r24514, "2", 10, MPFR_RNDN); mpfr_init(r24515); mpfr_init_set_str(r24516, "4", 10, MPFR_RNDN); mpfr_init(r24517); mpfr_init(r24518); mpfr_init(r24519); mpfr_init(r24520); mpfr_init(r24521); mpfr_init(r24522); mpfr_init(r24523); mpfr_init(r24524); mpfr_init(r24525); mpfr_init(r24526); mpfr_init(r24527); } double f_fm(double a, double b, double c) { mpfr_set_d(r24513, c, MPFR_RNDN); ; mpfr_div(r24515, r24513, r24514, MPFR_RNDN); ; mpfr_set_d(r24517, b, MPFR_RNDN); mpfr_neg(r24518, r24517, MPFR_RNDN); mpfr_mul(r24519, r24517, r24517, MPFR_RNDN); mpfr_set_d(r24520, a, MPFR_RNDN); mpfr_mul(r24521, r24516, r24520, MPFR_RNDN); mpfr_mul(r24522, r24513, r24521, MPFR_RNDN); mpfr_sub(r24523, r24519, r24522, MPFR_RNDN); mpfr_sqrt(r24524, r24523, MPFR_RNDN); mpfr_sub(r24525, r24518, r24524, MPFR_RNDN); mpfr_div(r24526, r24516, r24525, MPFR_RNDN); mpfr_mul(r24527, r24515, r24526, MPFR_RNDN); return mpfr_get_d(r24527, MPFR_RNDN); } static mpfr_t r24528, r24529, r24530, r24531, r24532, r24533, r24534, r24535, r24536, r24537, r24538, r24539, r24540, r24541, r24542; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24528); mpfr_init_set_str(r24529, "2", 10, MPFR_RNDN); mpfr_init(r24530); mpfr_init_set_str(r24531, "4", 10, MPFR_RNDN); mpfr_init(r24532); mpfr_init(r24533); mpfr_init(r24534); mpfr_init(r24535); mpfr_init(r24536); mpfr_init(r24537); mpfr_init(r24538); mpfr_init(r24539); mpfr_init(r24540); mpfr_init(r24541); mpfr_init(r24542); } double f_dm(double a, double b, double c) { mpfr_set_d(r24528, c, MPFR_RNDN); ; mpfr_div(r24530, r24528, r24529, MPFR_RNDN); ; mpfr_set_d(r24532, b, MPFR_RNDN); mpfr_neg(r24533, r24532, MPFR_RNDN); mpfr_mul(r24534, r24532, r24532, MPFR_RNDN); mpfr_set_d(r24535, a, MPFR_RNDN); mpfr_mul(r24536, r24531, r24535, MPFR_RNDN); mpfr_mul(r24537, r24528, r24536, MPFR_RNDN); mpfr_sub(r24538, r24534, r24537, MPFR_RNDN); mpfr_sqrt(r24539, r24538, MPFR_RNDN); mpfr_sub(r24540, r24533, r24539, MPFR_RNDN); mpfr_div(r24541, r24531, r24540, MPFR_RNDN); mpfr_mul(r24542, r24530, r24541, MPFR_RNDN); return mpfr_get_d(r24542, MPFR_RNDN); }