#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17409 = b; float r17410 = -r17409; float r17411 = r17409 * r17409; float r17412 = 4.0f; float r17413 = a; float r17414 = r17412 * r17413; float r17415 = c; float r17416 = r17414 * r17415; float r17417 = r17411 - r17416; float r17418 = sqrt(r17417); float r17419 = r17410 + r17418; float r17420 = 2.0f; float r17421 = r17420 * r17413; float r17422 = r17419 / r17421; return r17422; } double f_id(double a, double b, double c) { double r17423 = b; double r17424 = -r17423; double r17425 = r17423 * r17423; double r17426 = 4.0; double r17427 = a; double r17428 = r17426 * r17427; double r17429 = c; double r17430 = r17428 * r17429; double r17431 = r17425 - r17430; double r17432 = sqrt(r17431); double r17433 = r17424 + r17432; double r17434 = 2.0; double r17435 = r17434 * r17427; double r17436 = r17433 / r17435; return r17436; } double f_of(float a, float b, float c) { float r17437 = 4.0f; float r17438 = a; float r17439 = r17437 * r17438; float r17440 = c; float r17441 = r17439 * r17440; float r17442 = b; float r17443 = -r17442; float r17444 = r17442 * r17442; float r17445 = r17444 - r17441; float r17446 = sqrt(r17445); float r17447 = r17443 - r17446; float r17448 = r17441 / r17447; float r17449 = 2.0f; float r17450 = r17449 * r17438; float r17451 = r17448 / r17450; return r17451; } double f_od(double a, double b, double c) { double r17452 = 4.0; double r17453 = a; double r17454 = r17452 * r17453; double r17455 = c; double r17456 = r17454 * r17455; double r17457 = b; double r17458 = -r17457; double r17459 = r17457 * r17457; double r17460 = r17459 - r17456; double r17461 = sqrt(r17460); double r17462 = r17458 - r17461; double r17463 = r17456 / r17462; double r17464 = 2.0; double r17465 = r17464 * r17453; double r17466 = r17463 / r17465; return r17466; } 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 r17467, r17468, r17469, r17470, r17471, r17472, r17473, r17474, r17475, r17476, r17477, r17478, r17479, r17480; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17467); mpfr_init(r17468); mpfr_init(r17469); mpfr_init_set_str(r17470, "4", 10, MPFR_RNDN); mpfr_init(r17471); mpfr_init(r17472); mpfr_init(r17473); mpfr_init(r17474); mpfr_init(r17475); mpfr_init(r17476); mpfr_init(r17477); mpfr_init_set_str(r17478, "2", 10, MPFR_RNDN); mpfr_init(r17479); mpfr_init(r17480); } double f_im(double a, double b, double c) { mpfr_set_d(r17467, b, MPFR_RNDN); mpfr_neg(r17468, r17467, MPFR_RNDN); mpfr_sqr(r17469, r17467, MPFR_RNDN); ; mpfr_set_d(r17471, a, MPFR_RNDN); mpfr_mul(r17472, r17470, r17471, MPFR_RNDN); mpfr_set_d(r17473, c, MPFR_RNDN); mpfr_mul(r17474, r17472, r17473, MPFR_RNDN); mpfr_sub(r17475, r17469, r17474, MPFR_RNDN); mpfr_sqrt(r17476, r17475, MPFR_RNDN); mpfr_add(r17477, r17468, r17476, MPFR_RNDN); ; mpfr_mul(r17479, r17478, r17471, MPFR_RNDN); mpfr_div(r17480, r17477, r17479, MPFR_RNDN); return mpfr_get_d(r17480, MPFR_RNDN); } static mpfr_t r17481, r17482, r17483, r17484, r17485, r17486, r17487, r17488, r17489, r17490, r17491, r17492, r17493, r17494, r17495; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17481, "4", 10, MPFR_RNDN); mpfr_init(r17482); mpfr_init(r17483); mpfr_init(r17484); mpfr_init(r17485); mpfr_init(r17486); mpfr_init(r17487); mpfr_init(r17488); mpfr_init(r17489); mpfr_init(r17490); mpfr_init(r17491); mpfr_init(r17492); mpfr_init_set_str(r17493, "2", 10, MPFR_RNDN); mpfr_init(r17494); mpfr_init(r17495); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17482, a, MPFR_RNDN); mpfr_mul(r17483, r17481, r17482, MPFR_RNDN); mpfr_set_d(r17484, c, MPFR_RNDN); mpfr_mul(r17485, r17483, r17484, MPFR_RNDN); mpfr_set_d(r17486, b, MPFR_RNDN); mpfr_neg(r17487, r17486, MPFR_RNDN); mpfr_sqr(r17488, r17486, MPFR_RNDN); mpfr_sub(r17489, r17488, r17485, MPFR_RNDN); mpfr_sqrt(r17490, r17489, MPFR_RNDN); mpfr_sub(r17491, r17487, r17490, MPFR_RNDN); mpfr_div(r17492, r17485, r17491, MPFR_RNDN); ; mpfr_mul(r17494, r17493, r17482, MPFR_RNDN); mpfr_div(r17495, r17492, r17494, MPFR_RNDN); return mpfr_get_d(r17495, MPFR_RNDN); } static mpfr_t r17496, r17497, r17498, r17499, r17500, r17501, r17502, r17503, r17504, r17505, r17506, r17507, r17508, r17509, r17510; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17496, "4", 10, MPFR_RNDN); mpfr_init(r17497); mpfr_init(r17498); mpfr_init(r17499); mpfr_init(r17500); mpfr_init(r17501); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init(r17505); mpfr_init(r17506); mpfr_init(r17507); mpfr_init_set_str(r17508, "2", 10, MPFR_RNDN); mpfr_init(r17509); mpfr_init(r17510); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17497, a, MPFR_RNDN); mpfr_mul(r17498, r17496, r17497, MPFR_RNDN); mpfr_set_d(r17499, c, MPFR_RNDN); mpfr_mul(r17500, r17498, r17499, MPFR_RNDN); mpfr_set_d(r17501, b, MPFR_RNDN); mpfr_neg(r17502, r17501, MPFR_RNDN); mpfr_sqr(r17503, r17501, MPFR_RNDN); mpfr_sub(r17504, r17503, r17500, MPFR_RNDN); mpfr_sqrt(r17505, r17504, MPFR_RNDN); mpfr_sub(r17506, r17502, r17505, MPFR_RNDN); mpfr_div(r17507, r17500, r17506, MPFR_RNDN); ; mpfr_mul(r17509, r17508, r17497, MPFR_RNDN); mpfr_div(r17510, r17507, r17509, MPFR_RNDN); return mpfr_get_d(r17510, MPFR_RNDN); }