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