#include #include #include #include #include char *name = "The quadratic formula (r2)"; double f_if(float a, float b, float c) { float r15382 = b; float r15383 = -r15382; float r15384 = r15382 * r15382; float r15385 = 4.0f; float r15386 = a; float r15387 = c; float r15388 = r15386 * r15387; float r15389 = r15385 * r15388; float r15390 = r15384 - r15389; float r15391 = sqrt(r15390); float r15392 = r15383 - r15391; float r15393 = 2.0f; float r15394 = r15393 * r15386; float r15395 = r15392 / r15394; return r15395; } double f_id(double a, double b, double c) { double r15396 = b; double r15397 = -r15396; double r15398 = r15396 * r15396; double r15399 = 4.0; double r15400 = a; double r15401 = c; double r15402 = r15400 * r15401; double r15403 = r15399 * r15402; double r15404 = r15398 - r15403; double r15405 = sqrt(r15404); double r15406 = r15397 - r15405; double r15407 = 2.0; double r15408 = r15407 * r15400; double r15409 = r15406 / r15408; return r15409; } double f_of(float a, float b, float c) { float r15410 = c; float r15411 = 4.0f; float r15412 = a; float r15413 = r15411 * r15412; float r15414 = r15410 * r15413; float r15415 = b; float r15416 = -r15415; float r15417 = r15415 * r15415; float r15418 = r15412 * r15410; float r15419 = r15411 * r15418; float r15420 = r15417 - r15419; float r15421 = sqrt(r15420); float r15422 = r15416 + r15421; float r15423 = r15414 / r15422; float r15424 = 2.0f; float r15425 = r15424 * r15412; float r15426 = r15423 / r15425; return r15426; } double f_od(double a, double b, double c) { double r15427 = c; double r15428 = 4.0; double r15429 = a; double r15430 = r15428 * r15429; double r15431 = r15427 * r15430; double r15432 = b; double r15433 = -r15432; double r15434 = r15432 * r15432; double r15435 = r15429 * r15427; double r15436 = r15428 * r15435; double r15437 = r15434 - r15436; double r15438 = sqrt(r15437); double r15439 = r15433 + r15438; double r15440 = r15431 / r15439; double r15441 = 2.0; double r15442 = r15441 * r15429; double r15443 = r15440 / r15442; return r15443; } 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 r15444, r15445, r15446, r15447, r15448, r15449, r15450, r15451, r15452, r15453, r15454, r15455, r15456, r15457; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15444); mpfr_init(r15445); mpfr_init(r15446); mpfr_init_set_str(r15447, "4", 10, MPFR_RNDN); mpfr_init(r15448); mpfr_init(r15449); mpfr_init(r15450); mpfr_init(r15451); mpfr_init(r15452); mpfr_init(r15453); mpfr_init(r15454); mpfr_init_set_str(r15455, "2", 10, MPFR_RNDN); mpfr_init(r15456); mpfr_init(r15457); } double f_im(double a, double b, double c) { mpfr_set_d(r15444, b, MPFR_RNDN); mpfr_neg(r15445, r15444, MPFR_RNDN); mpfr_sqr(r15446, r15444, MPFR_RNDN); ; mpfr_set_d(r15448, a, MPFR_RNDN); mpfr_set_d(r15449, c, MPFR_RNDN); mpfr_mul(r15450, r15448, r15449, MPFR_RNDN); mpfr_mul(r15451, r15447, r15450, MPFR_RNDN); mpfr_sub(r15452, r15446, r15451, MPFR_RNDN); mpfr_sqrt(r15453, r15452, MPFR_RNDN); mpfr_sub(r15454, r15445, r15453, MPFR_RNDN); ; mpfr_mul(r15456, r15455, r15448, MPFR_RNDN); mpfr_div(r15457, r15454, r15456, MPFR_RNDN); return mpfr_get_d(r15457, MPFR_RNDN); } static mpfr_t r15458, r15459, r15460, r15461, r15462, r15463, r15464, r15465, r15466, r15467, r15468, r15469, r15470, r15471, r15472, r15473, r15474; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15458); mpfr_init_set_str(r15459, "4", 10, MPFR_RNDN); mpfr_init(r15460); mpfr_init(r15461); mpfr_init(r15462); mpfr_init(r15463); mpfr_init(r15464); mpfr_init(r15465); mpfr_init(r15466); mpfr_init(r15467); mpfr_init(r15468); mpfr_init(r15469); mpfr_init(r15470); mpfr_init(r15471); mpfr_init_set_str(r15472, "2", 10, MPFR_RNDN); mpfr_init(r15473); mpfr_init(r15474); } double f_fm(double a, double b, double c) { mpfr_set_d(r15458, c, MPFR_RNDN); ; mpfr_set_d(r15460, a, MPFR_RNDN); mpfr_mul(r15461, r15459, r15460, MPFR_RNDN); mpfr_mul(r15462, r15458, r15461, MPFR_RNDN); mpfr_set_d(r15463, b, MPFR_RNDN); mpfr_neg(r15464, r15463, MPFR_RNDN); mpfr_sqr(r15465, r15463, MPFR_RNDN); mpfr_mul(r15466, r15460, r15458, MPFR_RNDN); mpfr_mul(r15467, r15459, r15466, MPFR_RNDN); mpfr_sub(r15468, r15465, r15467, MPFR_RNDN); mpfr_sqrt(r15469, r15468, MPFR_RNDN); mpfr_add(r15470, r15464, r15469, MPFR_RNDN); mpfr_div(r15471, r15462, r15470, MPFR_RNDN); ; mpfr_mul(r15473, r15472, r15460, MPFR_RNDN); mpfr_div(r15474, r15471, r15473, MPFR_RNDN); return mpfr_get_d(r15474, MPFR_RNDN); } static mpfr_t r15475, r15476, r15477, r15478, r15479, r15480, r15481, r15482, r15483, r15484, r15485, r15486, r15487, r15488, r15489, r15490, r15491; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15475); mpfr_init_set_str(r15476, "4", 10, MPFR_RNDN); mpfr_init(r15477); mpfr_init(r15478); mpfr_init(r15479); mpfr_init(r15480); mpfr_init(r15481); mpfr_init(r15482); mpfr_init(r15483); mpfr_init(r15484); mpfr_init(r15485); mpfr_init(r15486); mpfr_init(r15487); mpfr_init(r15488); mpfr_init_set_str(r15489, "2", 10, MPFR_RNDN); mpfr_init(r15490); mpfr_init(r15491); } double f_dm(double a, double b, double c) { mpfr_set_d(r15475, c, MPFR_RNDN); ; mpfr_set_d(r15477, a, MPFR_RNDN); mpfr_mul(r15478, r15476, r15477, MPFR_RNDN); mpfr_mul(r15479, r15475, r15478, MPFR_RNDN); mpfr_set_d(r15480, b, MPFR_RNDN); mpfr_neg(r15481, r15480, MPFR_RNDN); mpfr_sqr(r15482, r15480, MPFR_RNDN); mpfr_mul(r15483, r15477, r15475, MPFR_RNDN); mpfr_mul(r15484, r15476, r15483, MPFR_RNDN); mpfr_sub(r15485, r15482, r15484, MPFR_RNDN); mpfr_sqrt(r15486, r15485, MPFR_RNDN); mpfr_add(r15487, r15481, r15486, MPFR_RNDN); mpfr_div(r15488, r15479, r15487, MPFR_RNDN); ; mpfr_mul(r15490, r15489, r15477, MPFR_RNDN); mpfr_div(r15491, r15488, r15490, MPFR_RNDN); return mpfr_get_d(r15491, MPFR_RNDN); }