#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r27342 = b; float r27343 = -r27342; float r27344 = r27342 * r27342; float r27345 = 4; float r27346 = a; float r27347 = r27345 * r27346; float r27348 = c; float r27349 = r27347 * r27348; float r27350 = r27344 - r27349; float r27351 = sqrt(r27350); float r27352 = r27343 + r27351; float r27353 = 2; float r27354 = r27353 * r27346; float r27355 = r27352 / r27354; return r27355; } double f_id(double a, double b, double c) { double r27356 = b; double r27357 = -r27356; double r27358 = r27356 * r27356; double r27359 = 4; double r27360 = a; double r27361 = r27359 * r27360; double r27362 = c; double r27363 = r27361 * r27362; double r27364 = r27358 - r27363; double r27365 = sqrt(r27364); double r27366 = r27357 + r27365; double r27367 = 2; double r27368 = r27367 * r27360; double r27369 = r27366 / r27368; return r27369; } double f_of(float a, float b, float c) { float r27370 = b; float r27371 = -r27370; float r27372 = -5.7583530328199956e-65; bool r27373 = r27371 <= r27372; float r27374 = c; float r27375 = -r27374; float r27376 = r27375 / r27370; float r27377 = 7.455938382452604e+91; bool r27378 = r27371 <= r27377; float r27379 = r27370 * r27370; float r27380 = a; float r27381 = 4; float r27382 = r27380 * r27381; float r27383 = r27382 * r27374; float r27384 = r27379 - r27383; float r27385 = sqrt(r27384); float r27386 = r27385 + r27371; float r27387 = 2; float r27388 = r27387 * r27380; float r27389 = r27386 / r27388; float r27390 = r27374 / r27370; float r27391 = r27370 + r27370; float r27392 = r27391 / r27388; float r27393 = r27390 - r27392; float r27394 = r27378 ? r27389 : r27393; float r27395 = r27373 ? r27376 : r27394; return r27395; } double f_od(double a, double b, double c) { double r27396 = b; double r27397 = -r27396; double r27398 = -5.7583530328199956e-65; bool r27399 = r27397 <= r27398; double r27400 = c; double r27401 = -r27400; double r27402 = r27401 / r27396; double r27403 = 7.455938382452604e+91; bool r27404 = r27397 <= r27403; double r27405 = r27396 * r27396; double r27406 = a; double r27407 = 4; double r27408 = r27406 * r27407; double r27409 = r27408 * r27400; double r27410 = r27405 - r27409; double r27411 = sqrt(r27410); double r27412 = r27411 + r27397; double r27413 = 2; double r27414 = r27413 * r27406; double r27415 = r27412 / r27414; double r27416 = r27400 / r27396; double r27417 = r27396 + r27396; double r27418 = r27417 / r27414; double r27419 = r27416 - r27418; double r27420 = r27404 ? r27415 : r27419; double r27421 = r27399 ? r27402 : r27420; return r27421; } 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 r27422, r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430, r27431, r27432, r27433, r27434, r27435; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init(r27422); mpfr_init(r27423); mpfr_init(r27424); mpfr_init_set_str(r27425, "4", 10, MPFR_RNDN); mpfr_init(r27426); mpfr_init(r27427); mpfr_init(r27428); mpfr_init(r27429); mpfr_init(r27430); mpfr_init(r27431); mpfr_init(r27432); mpfr_init_set_str(r27433, "2", 10, MPFR_RNDN); mpfr_init(r27434); mpfr_init(r27435); } double f_im(double a, double b, double c) { mpfr_set_d(r27422, b, MPFR_RNDN); mpfr_neg(r27423, r27422, MPFR_RNDN); mpfr_mul(r27424, r27422, r27422, MPFR_RNDN); ; mpfr_set_d(r27426, a, MPFR_RNDN); mpfr_mul(r27427, r27425, r27426, MPFR_RNDN); mpfr_set_d(r27428, c, MPFR_RNDN); mpfr_mul(r27429, r27427, r27428, MPFR_RNDN); mpfr_sub(r27430, r27424, r27429, MPFR_RNDN); mpfr_sqrt(r27431, r27430, MPFR_RNDN); mpfr_add(r27432, r27423, r27431, MPFR_RNDN); ; mpfr_mul(r27434, r27433, r27426, MPFR_RNDN); mpfr_div(r27435, r27432, r27434, MPFR_RNDN); return mpfr_get_d(r27435, MPFR_RNDN); } static mpfr_t r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init(r27436); mpfr_init(r27437); mpfr_init_set_str(r27438, "-5.7583530328199956e-65", 10, MPFR_RNDN); mpfr_init(r27439); mpfr_init(r27440); mpfr_init(r27441); mpfr_init(r27442); mpfr_init_set_str(r27443, "7.455938382452604e+91", 10, MPFR_RNDN); mpfr_init(r27444); mpfr_init(r27445); mpfr_init(r27446); mpfr_init_set_str(r27447, "4", 10, MPFR_RNDN); mpfr_init(r27448); mpfr_init(r27449); mpfr_init(r27450); mpfr_init(r27451); mpfr_init(r27452); mpfr_init_set_str(r27453, "2", 10, MPFR_RNDN); mpfr_init(r27454); mpfr_init(r27455); mpfr_init(r27456); mpfr_init(r27457); mpfr_init(r27458); mpfr_init(r27459); mpfr_init(r27460); mpfr_init(r27461); } double f_fm(double a, double b, double c) { mpfr_set_d(r27436, b, MPFR_RNDN); mpfr_neg(r27437, r27436, MPFR_RNDN); ; mpfr_set_si(r27439, mpfr_cmp(r27437, r27438) <= 0, MPFR_RNDN); mpfr_set_d(r27440, c, MPFR_RNDN); mpfr_neg(r27441, r27440, MPFR_RNDN); mpfr_div(r27442, r27441, r27436, MPFR_RNDN); ; mpfr_set_si(r27444, mpfr_cmp(r27437, r27443) <= 0, MPFR_RNDN); mpfr_mul(r27445, r27436, r27436, MPFR_RNDN); mpfr_set_d(r27446, a, MPFR_RNDN); ; mpfr_mul(r27448, r27446, r27447, MPFR_RNDN); mpfr_mul(r27449, r27448, r27440, MPFR_RNDN); mpfr_sub(r27450, r27445, r27449, MPFR_RNDN); mpfr_sqrt(r27451, r27450, MPFR_RNDN); mpfr_add(r27452, r27451, r27437, MPFR_RNDN); ; mpfr_mul(r27454, r27453, r27446, MPFR_RNDN); mpfr_div(r27455, r27452, r27454, MPFR_RNDN); mpfr_div(r27456, r27440, r27436, MPFR_RNDN); mpfr_add(r27457, r27436, r27436, MPFR_RNDN); mpfr_div(r27458, r27457, r27454, MPFR_RNDN); mpfr_sub(r27459, r27456, r27458, MPFR_RNDN); if (mpfr_get_si(r27444, MPFR_RNDN)) { mpfr_set(r27460, r27455, MPFR_RNDN); } else { mpfr_set(r27460, r27459, MPFR_RNDN); }; if (mpfr_get_si(r27439, MPFR_RNDN)) { mpfr_set(r27461, r27442, MPFR_RNDN); } else { mpfr_set(r27461, r27460, MPFR_RNDN); }; return mpfr_get_d(r27461, MPFR_RNDN); } static mpfr_t r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483, r27484, r27485, r27486, r27487; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init(r27462); mpfr_init(r27463); mpfr_init_set_str(r27464, "-5.7583530328199956e-65", 10, MPFR_RNDN); mpfr_init(r27465); mpfr_init(r27466); mpfr_init(r27467); mpfr_init(r27468); mpfr_init_set_str(r27469, "7.455938382452604e+91", 10, MPFR_RNDN); mpfr_init(r27470); mpfr_init(r27471); mpfr_init(r27472); mpfr_init_set_str(r27473, "4", 10, MPFR_RNDN); mpfr_init(r27474); mpfr_init(r27475); mpfr_init(r27476); mpfr_init(r27477); mpfr_init(r27478); mpfr_init_set_str(r27479, "2", 10, MPFR_RNDN); mpfr_init(r27480); mpfr_init(r27481); mpfr_init(r27482); mpfr_init(r27483); mpfr_init(r27484); mpfr_init(r27485); mpfr_init(r27486); mpfr_init(r27487); } double f_dm(double a, double b, double c) { mpfr_set_d(r27462, b, MPFR_RNDN); mpfr_neg(r27463, r27462, MPFR_RNDN); ; mpfr_set_si(r27465, mpfr_cmp(r27463, r27464) <= 0, MPFR_RNDN); mpfr_set_d(r27466, c, MPFR_RNDN); mpfr_neg(r27467, r27466, MPFR_RNDN); mpfr_div(r27468, r27467, r27462, MPFR_RNDN); ; mpfr_set_si(r27470, mpfr_cmp(r27463, r27469) <= 0, MPFR_RNDN); mpfr_mul(r27471, r27462, r27462, MPFR_RNDN); mpfr_set_d(r27472, a, MPFR_RNDN); ; mpfr_mul(r27474, r27472, r27473, MPFR_RNDN); mpfr_mul(r27475, r27474, r27466, MPFR_RNDN); mpfr_sub(r27476, r27471, r27475, MPFR_RNDN); mpfr_sqrt(r27477, r27476, MPFR_RNDN); mpfr_add(r27478, r27477, r27463, MPFR_RNDN); ; mpfr_mul(r27480, r27479, r27472, MPFR_RNDN); mpfr_div(r27481, r27478, r27480, MPFR_RNDN); mpfr_div(r27482, r27466, r27462, MPFR_RNDN); mpfr_add(r27483, r27462, r27462, MPFR_RNDN); mpfr_div(r27484, r27483, r27480, MPFR_RNDN); mpfr_sub(r27485, r27482, r27484, MPFR_RNDN); if (mpfr_get_si(r27470, MPFR_RNDN)) { mpfr_set(r27486, r27481, MPFR_RNDN); } else { mpfr_set(r27486, r27485, MPFR_RNDN); }; if (mpfr_get_si(r27465, MPFR_RNDN)) { mpfr_set(r27487, r27468, MPFR_RNDN); } else { mpfr_set(r27487, r27486, MPFR_RNDN); }; return mpfr_get_d(r27487, MPFR_RNDN); }