#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r3911286 = b; float r3911287 = -r3911286; float r3911288 = r3911286 * r3911286; float r3911289 = 4; float r3911290 = a; float r3911291 = r3911289 * r3911290; float r3911292 = c; float r3911293 = r3911291 * r3911292; float r3911294 = r3911288 - r3911293; float r3911295 = sqrt(r3911294); float r3911296 = r3911287 + r3911295; float r3911297 = 2; float r3911298 = r3911297 * r3911290; float r3911299 = r3911296 / r3911298; return r3911299; } double f_id(double a, double b, double c) { double r3911300 = b; double r3911301 = -r3911300; double r3911302 = r3911300 * r3911300; double r3911303 = 4; double r3911304 = a; double r3911305 = r3911303 * r3911304; double r3911306 = c; double r3911307 = r3911305 * r3911306; double r3911308 = r3911302 - r3911307; double r3911309 = sqrt(r3911308); double r3911310 = r3911301 + r3911309; double r3911311 = 2; double r3911312 = r3911311 * r3911304; double r3911313 = r3911310 / r3911312; return r3911313; } double f_of(float a, float b, float c) { float r3911314 = 4; float r3911315 = 2; float r3911316 = r3911314 / r3911315; float r3911317 = c; float r3911318 = b; float r3911319 = -r3911318; float r3911320 = r3911318 * r3911318; float r3911321 = a; float r3911322 = r3911317 * r3911321; float r3911323 = r3911322 * r3911314; float r3911324 = r3911320 - r3911323; float r3911325 = sqrt(r3911324); float r3911326 = r3911319 - r3911325; float r3911327 = r3911317 / r3911326; float r3911328 = r3911316 * r3911327; return r3911328; } double f_od(double a, double b, double c) { double r3911329 = 4; double r3911330 = 2; double r3911331 = r3911329 / r3911330; double r3911332 = c; double r3911333 = b; double r3911334 = -r3911333; double r3911335 = r3911333 * r3911333; double r3911336 = a; double r3911337 = r3911332 * r3911336; double r3911338 = r3911337 * r3911329; double r3911339 = r3911335 - r3911338; double r3911340 = sqrt(r3911339); double r3911341 = r3911334 - r3911340; double r3911342 = r3911332 / r3911341; double r3911343 = r3911331 * r3911342; return r3911343; } 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 r3911344, r3911345, r3911346, r3911347, r3911348, r3911349, r3911350, r3911351, r3911352, r3911353, r3911354, r3911355, r3911356, r3911357; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r3911344); mpfr_init(r3911345); mpfr_init(r3911346); mpfr_init_set_str(r3911347, "4", 10, MPFR_RNDN); mpfr_init(r3911348); mpfr_init(r3911349); mpfr_init(r3911350); mpfr_init(r3911351); mpfr_init(r3911352); mpfr_init(r3911353); mpfr_init(r3911354); mpfr_init_set_str(r3911355, "2", 10, MPFR_RNDN); mpfr_init(r3911356); mpfr_init(r3911357); } double f_im(double a, double b, double c) { mpfr_set_d(r3911344, b, MPFR_RNDN); mpfr_neg(r3911345, r3911344, MPFR_RNDN); mpfr_mul(r3911346, r3911344, r3911344, MPFR_RNDN); ; mpfr_set_d(r3911348, a, MPFR_RNDN); mpfr_mul(r3911349, r3911347, r3911348, MPFR_RNDN); mpfr_set_d(r3911350, c, MPFR_RNDN); mpfr_mul(r3911351, r3911349, r3911350, MPFR_RNDN); mpfr_sub(r3911352, r3911346, r3911351, MPFR_RNDN); mpfr_sqrt(r3911353, r3911352, MPFR_RNDN); mpfr_add(r3911354, r3911345, r3911353, MPFR_RNDN); ; mpfr_mul(r3911356, r3911355, r3911348, MPFR_RNDN); mpfr_div(r3911357, r3911354, r3911356, MPFR_RNDN); return mpfr_get_d(r3911357, MPFR_RNDN); } static mpfr_t r3911358, r3911359, r3911360, r3911361, r3911362, r3911363, r3911364, r3911365, r3911366, r3911367, r3911368, r3911369, r3911370, r3911371, r3911372; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init_set_str(r3911358, "4", 10, MPFR_RNDN); mpfr_init_set_str(r3911359, "2", 10, MPFR_RNDN); mpfr_init(r3911360); mpfr_init(r3911361); mpfr_init(r3911362); mpfr_init(r3911363); mpfr_init(r3911364); mpfr_init(r3911365); mpfr_init(r3911366); mpfr_init(r3911367); mpfr_init(r3911368); mpfr_init(r3911369); mpfr_init(r3911370); mpfr_init(r3911371); mpfr_init(r3911372); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r3911360, r3911358, r3911359, MPFR_RNDN); mpfr_set_d(r3911361, c, MPFR_RNDN); mpfr_set_d(r3911362, b, MPFR_RNDN); mpfr_neg(r3911363, r3911362, MPFR_RNDN); mpfr_mul(r3911364, r3911362, r3911362, MPFR_RNDN); mpfr_set_d(r3911365, a, MPFR_RNDN); mpfr_mul(r3911366, r3911361, r3911365, MPFR_RNDN); mpfr_mul(r3911367, r3911366, r3911358, MPFR_RNDN); mpfr_sub(r3911368, r3911364, r3911367, MPFR_RNDN); mpfr_sqrt(r3911369, r3911368, MPFR_RNDN); mpfr_sub(r3911370, r3911363, r3911369, MPFR_RNDN); mpfr_div(r3911371, r3911361, r3911370, MPFR_RNDN); mpfr_mul(r3911372, r3911360, r3911371, MPFR_RNDN); return mpfr_get_d(r3911372, MPFR_RNDN); } static mpfr_t r3911373, r3911374, r3911375, r3911376, r3911377, r3911378, r3911379, r3911380, r3911381, r3911382, r3911383, r3911384, r3911385, r3911386, r3911387; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init_set_str(r3911373, "4", 10, MPFR_RNDN); mpfr_init_set_str(r3911374, "2", 10, MPFR_RNDN); mpfr_init(r3911375); mpfr_init(r3911376); mpfr_init(r3911377); mpfr_init(r3911378); mpfr_init(r3911379); mpfr_init(r3911380); mpfr_init(r3911381); mpfr_init(r3911382); mpfr_init(r3911383); mpfr_init(r3911384); mpfr_init(r3911385); mpfr_init(r3911386); mpfr_init(r3911387); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r3911375, r3911373, r3911374, MPFR_RNDN); mpfr_set_d(r3911376, c, MPFR_RNDN); mpfr_set_d(r3911377, b, MPFR_RNDN); mpfr_neg(r3911378, r3911377, MPFR_RNDN); mpfr_mul(r3911379, r3911377, r3911377, MPFR_RNDN); mpfr_set_d(r3911380, a, MPFR_RNDN); mpfr_mul(r3911381, r3911376, r3911380, MPFR_RNDN); mpfr_mul(r3911382, r3911381, r3911373, MPFR_RNDN); mpfr_sub(r3911383, r3911379, r3911382, MPFR_RNDN); mpfr_sqrt(r3911384, r3911383, MPFR_RNDN); mpfr_sub(r3911385, r3911378, r3911384, MPFR_RNDN); mpfr_div(r3911386, r3911376, r3911385, MPFR_RNDN); mpfr_mul(r3911387, r3911375, r3911386, MPFR_RNDN); return mpfr_get_d(r3911387, MPFR_RNDN); }