#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17341 = b; float r17342 = -r17341; float r17343 = r17341 * r17341; float r17344 = 4.0f; float r17345 = a; float r17346 = r17344 * r17345; float r17347 = c; float r17348 = r17346 * r17347; float r17349 = r17343 - r17348; float r17350 = sqrt(r17349); float r17351 = r17342 + r17350; float r17352 = 2.0f; float r17353 = r17352 * r17345; float r17354 = r17351 / r17353; return r17354; } double f_id(double a, double b, double c) { double r17355 = b; double r17356 = -r17355; double r17357 = r17355 * r17355; double r17358 = 4.0; double r17359 = a; double r17360 = r17358 * r17359; double r17361 = c; double r17362 = r17360 * r17361; double r17363 = r17357 - r17362; double r17364 = sqrt(r17363); double r17365 = r17356 + r17364; double r17366 = 2.0; double r17367 = r17366 * r17359; double r17368 = r17365 / r17367; return r17368; } double f_of(float a, float b, float c) { float r17369 = c; float r17370 = 1.0f; float r17371 = r17369 / r17370; float r17372 = 4.0f; float r17373 = 2.0f; float r17374 = r17372 / r17373; float r17375 = r17371 * r17374; float r17376 = b; float r17377 = -r17376; float r17378 = r17376 * r17376; float r17379 = a; float r17380 = r17372 * r17379; float r17381 = r17380 * r17369; float r17382 = r17378 - r17381; float r17383 = sqrt(r17382); float r17384 = r17377 - r17383; float r17385 = r17375 / r17384; return r17385; } double f_od(double a, double b, double c) { double r17386 = c; double r17387 = 1.0; double r17388 = r17386 / r17387; double r17389 = 4.0; double r17390 = 2.0; double r17391 = r17389 / r17390; double r17392 = r17388 * r17391; double r17393 = b; double r17394 = -r17393; double r17395 = r17393 * r17393; double r17396 = a; double r17397 = r17389 * r17396; double r17398 = r17397 * r17386; double r17399 = r17395 - r17398; double r17400 = sqrt(r17399); double r17401 = r17394 - r17400; double r17402 = r17392 / r17401; return r17402; } 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 r17403, r17404, r17405, r17406, r17407, r17408, r17409, r17410, r17411, r17412, r17413, r17414, r17415, r17416; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17403); mpfr_init(r17404); mpfr_init(r17405); mpfr_init_set_str(r17406, "4", 10, MPFR_RNDN); mpfr_init(r17407); mpfr_init(r17408); mpfr_init(r17409); mpfr_init(r17410); mpfr_init(r17411); mpfr_init(r17412); mpfr_init(r17413); mpfr_init_set_str(r17414, "2", 10, MPFR_RNDN); mpfr_init(r17415); mpfr_init(r17416); } double f_im(double a, double b, double c) { mpfr_set_d(r17403, b, MPFR_RNDN); mpfr_neg(r17404, r17403, MPFR_RNDN); mpfr_sqr(r17405, r17403, MPFR_RNDN); ; mpfr_set_d(r17407, a, MPFR_RNDN); mpfr_mul(r17408, r17406, r17407, MPFR_RNDN); mpfr_set_d(r17409, c, MPFR_RNDN); mpfr_mul(r17410, r17408, r17409, MPFR_RNDN); mpfr_sub(r17411, r17405, r17410, MPFR_RNDN); mpfr_sqrt(r17412, r17411, MPFR_RNDN); mpfr_add(r17413, r17404, r17412, MPFR_RNDN); ; mpfr_mul(r17415, r17414, r17407, MPFR_RNDN); mpfr_div(r17416, r17413, r17415, MPFR_RNDN); return mpfr_get_d(r17416, MPFR_RNDN); } static mpfr_t r17417, r17418, r17419, r17420, r17421, r17422, r17423, r17424, r17425, r17426, r17427, r17428, r17429, r17430, r17431, r17432, r17433; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17417); mpfr_init_set_str(r17418, "1", 10, MPFR_RNDN); mpfr_init(r17419); mpfr_init_set_str(r17420, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17421, "2", 10, MPFR_RNDN); mpfr_init(r17422); mpfr_init(r17423); mpfr_init(r17424); mpfr_init(r17425); mpfr_init(r17426); mpfr_init(r17427); mpfr_init(r17428); mpfr_init(r17429); mpfr_init(r17430); mpfr_init(r17431); mpfr_init(r17432); mpfr_init(r17433); } double f_fm(double a, double b, double c) { mpfr_set_d(r17417, c, MPFR_RNDN); ; mpfr_div(r17419, r17417, r17418, MPFR_RNDN); ; ; mpfr_div(r17422, r17420, r17421, MPFR_RNDN); mpfr_mul(r17423, r17419, r17422, MPFR_RNDN); mpfr_set_d(r17424, b, MPFR_RNDN); mpfr_neg(r17425, r17424, MPFR_RNDN); mpfr_sqr(r17426, r17424, MPFR_RNDN); mpfr_set_d(r17427, a, MPFR_RNDN); mpfr_mul(r17428, r17420, r17427, MPFR_RNDN); mpfr_mul(r17429, r17428, r17417, MPFR_RNDN); mpfr_sub(r17430, r17426, r17429, MPFR_RNDN); mpfr_sqrt(r17431, r17430, MPFR_RNDN); mpfr_sub(r17432, r17425, r17431, MPFR_RNDN); mpfr_div(r17433, r17423, r17432, MPFR_RNDN); return mpfr_get_d(r17433, MPFR_RNDN); } static mpfr_t r17434, r17435, r17436, r17437, r17438, r17439, r17440, r17441, r17442, r17443, r17444, r17445, r17446, r17447, r17448, r17449, r17450; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17434); mpfr_init_set_str(r17435, "1", 10, MPFR_RNDN); mpfr_init(r17436); mpfr_init_set_str(r17437, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17438, "2", 10, MPFR_RNDN); mpfr_init(r17439); mpfr_init(r17440); mpfr_init(r17441); mpfr_init(r17442); mpfr_init(r17443); mpfr_init(r17444); mpfr_init(r17445); mpfr_init(r17446); mpfr_init(r17447); mpfr_init(r17448); mpfr_init(r17449); mpfr_init(r17450); } double f_dm(double a, double b, double c) { mpfr_set_d(r17434, c, MPFR_RNDN); ; mpfr_div(r17436, r17434, r17435, MPFR_RNDN); ; ; mpfr_div(r17439, r17437, r17438, MPFR_RNDN); mpfr_mul(r17440, r17436, r17439, MPFR_RNDN); mpfr_set_d(r17441, b, MPFR_RNDN); mpfr_neg(r17442, r17441, MPFR_RNDN); mpfr_sqr(r17443, r17441, MPFR_RNDN); mpfr_set_d(r17444, a, MPFR_RNDN); mpfr_mul(r17445, r17437, r17444, MPFR_RNDN); mpfr_mul(r17446, r17445, r17434, MPFR_RNDN); mpfr_sub(r17447, r17443, r17446, MPFR_RNDN); mpfr_sqrt(r17448, r17447, MPFR_RNDN); mpfr_sub(r17449, r17442, r17448, MPFR_RNDN); mpfr_div(r17450, r17440, r17449, MPFR_RNDN); return mpfr_get_d(r17450, MPFR_RNDN); }