#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17429 = b; float r17430 = -r17429; float r17431 = r17429 * r17429; float r17432 = 4.0f; float r17433 = a; float r17434 = r17432 * r17433; float r17435 = c; float r17436 = r17434 * r17435; float r17437 = r17431 - r17436; float r17438 = sqrt(r17437); float r17439 = r17430 + r17438; float r17440 = 2.0f; float r17441 = r17440 * r17433; float r17442 = r17439 / r17441; return r17442; } double f_id(double a, double b, double c) { double r17443 = b; double r17444 = -r17443; double r17445 = r17443 * r17443; double r17446 = 4.0; double r17447 = a; double r17448 = r17446 * r17447; double r17449 = c; double r17450 = r17448 * r17449; double r17451 = r17445 - r17450; double r17452 = sqrt(r17451); double r17453 = r17444 + r17452; double r17454 = 2.0; double r17455 = r17454 * r17447; double r17456 = r17453 / r17455; return r17456; } double f_of(float a, float b, float c) { float r17457 = 4.0f; float r17458 = a; float r17459 = r17457 * r17458; float r17460 = c; float r17461 = r17459 * r17460; float r17462 = b; float r17463 = -r17462; float r17464 = r17462 * r17462; float r17465 = r17464 - r17461; float r17466 = sqrt(r17465); float r17467 = r17463 - r17466; float r17468 = r17461 / r17467; float r17469 = 2.0f; float r17470 = r17469 * r17458; float r17471 = r17468 / r17470; return r17471; } double f_od(double a, double b, double c) { double r17472 = 4.0; double r17473 = a; double r17474 = r17472 * r17473; double r17475 = c; double r17476 = r17474 * r17475; double r17477 = b; double r17478 = -r17477; double r17479 = r17477 * r17477; double r17480 = r17479 - r17476; double r17481 = sqrt(r17480); double r17482 = r17478 - r17481; double r17483 = r17476 / r17482; double r17484 = 2.0; double r17485 = r17484 * r17473; double r17486 = r17483 / r17485; return r17486; } 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 r17487, r17488, r17489, r17490, r17491, r17492, r17493, r17494, r17495, r17496, r17497, r17498, r17499, r17500; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17487); mpfr_init(r17488); mpfr_init(r17489); mpfr_init_set_str(r17490, "4", 10, MPFR_RNDN); mpfr_init(r17491); mpfr_init(r17492); mpfr_init(r17493); mpfr_init(r17494); mpfr_init(r17495); mpfr_init(r17496); mpfr_init(r17497); mpfr_init_set_str(r17498, "2", 10, MPFR_RNDN); mpfr_init(r17499); mpfr_init(r17500); } double f_im(double a, double b, double c) { mpfr_set_d(r17487, b, MPFR_RNDN); mpfr_neg(r17488, r17487, MPFR_RNDN); mpfr_sqr(r17489, r17487, MPFR_RNDN); ; mpfr_set_d(r17491, a, MPFR_RNDN); mpfr_mul(r17492, r17490, r17491, MPFR_RNDN); mpfr_set_d(r17493, c, MPFR_RNDN); mpfr_mul(r17494, r17492, r17493, MPFR_RNDN); mpfr_sub(r17495, r17489, r17494, MPFR_RNDN); mpfr_sqrt(r17496, r17495, MPFR_RNDN); mpfr_add(r17497, r17488, r17496, MPFR_RNDN); ; mpfr_mul(r17499, r17498, r17491, MPFR_RNDN); mpfr_div(r17500, r17497, r17499, MPFR_RNDN); return mpfr_get_d(r17500, MPFR_RNDN); } static mpfr_t r17501, r17502, r17503, r17504, r17505, r17506, r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17501, "4", 10, MPFR_RNDN); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init(r17505); mpfr_init(r17506); mpfr_init(r17507); mpfr_init(r17508); mpfr_init(r17509); mpfr_init(r17510); mpfr_init(r17511); mpfr_init(r17512); mpfr_init_set_str(r17513, "2", 10, MPFR_RNDN); mpfr_init(r17514); mpfr_init(r17515); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17502, a, MPFR_RNDN); mpfr_mul(r17503, r17501, r17502, MPFR_RNDN); mpfr_set_d(r17504, c, MPFR_RNDN); mpfr_mul(r17505, r17503, r17504, MPFR_RNDN); mpfr_set_d(r17506, b, MPFR_RNDN); mpfr_neg(r17507, r17506, MPFR_RNDN); mpfr_sqr(r17508, r17506, MPFR_RNDN); mpfr_sub(r17509, r17508, r17505, MPFR_RNDN); mpfr_sqrt(r17510, r17509, MPFR_RNDN); mpfr_sub(r17511, r17507, r17510, MPFR_RNDN); mpfr_div(r17512, r17505, r17511, MPFR_RNDN); ; mpfr_mul(r17514, r17513, r17502, MPFR_RNDN); mpfr_div(r17515, r17512, r17514, MPFR_RNDN); return mpfr_get_d(r17515, MPFR_RNDN); } static mpfr_t r17516, r17517, r17518, r17519, r17520, r17521, r17522, r17523, r17524, r17525, r17526, r17527, r17528, r17529, r17530; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17516, "4", 10, MPFR_RNDN); mpfr_init(r17517); mpfr_init(r17518); mpfr_init(r17519); mpfr_init(r17520); mpfr_init(r17521); mpfr_init(r17522); mpfr_init(r17523); mpfr_init(r17524); mpfr_init(r17525); mpfr_init(r17526); mpfr_init(r17527); mpfr_init_set_str(r17528, "2", 10, MPFR_RNDN); mpfr_init(r17529); mpfr_init(r17530); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17517, a, MPFR_RNDN); mpfr_mul(r17518, r17516, r17517, MPFR_RNDN); mpfr_set_d(r17519, c, MPFR_RNDN); mpfr_mul(r17520, r17518, r17519, MPFR_RNDN); mpfr_set_d(r17521, b, MPFR_RNDN); mpfr_neg(r17522, r17521, MPFR_RNDN); mpfr_sqr(r17523, r17521, MPFR_RNDN); mpfr_sub(r17524, r17523, r17520, MPFR_RNDN); mpfr_sqrt(r17525, r17524, MPFR_RNDN); mpfr_sub(r17526, r17522, r17525, MPFR_RNDN); mpfr_div(r17527, r17520, r17526, MPFR_RNDN); ; mpfr_mul(r17529, r17528, r17517, MPFR_RNDN); mpfr_div(r17530, r17527, r17529, MPFR_RNDN); return mpfr_get_d(r17530, MPFR_RNDN); }