#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r20836 = b; float r20837 = -r20836; float r20838 = r20836 * r20836; float r20839 = 4; float r20840 = a; float r20841 = r20839 * r20840; float r20842 = c; float r20843 = r20841 * r20842; float r20844 = r20838 - r20843; float r20845 = sqrt(r20844); float r20846 = r20837 + r20845; float r20847 = 2; float r20848 = r20847 * r20840; float r20849 = r20846 / r20848; return r20849; } double f_id(double a, double b, double c) { double r20850 = b; double r20851 = -r20850; double r20852 = r20850 * r20850; double r20853 = 4; double r20854 = a; double r20855 = r20853 * r20854; double r20856 = c; double r20857 = r20855 * r20856; double r20858 = r20852 - r20857; double r20859 = sqrt(r20858); double r20860 = r20851 + r20859; double r20861 = 2; double r20862 = r20861 * r20854; double r20863 = r20860 / r20862; return r20863; } double f_of(float a, float b, float c) { float r20864 = 4; float r20865 = 2; float r20866 = r20864 / r20865; float r20867 = c; float r20868 = b; float r20869 = -r20868; float r20870 = r20868 * r20868; float r20871 = a; float r20872 = r20867 * r20871; float r20873 = r20872 * r20864; float r20874 = r20870 - r20873; float r20875 = sqrt(r20874); float r20876 = r20869 - r20875; float r20877 = r20867 / r20876; float r20878 = r20866 * r20877; return r20878; } double f_od(double a, double b, double c) { double r20879 = 4; double r20880 = 2; double r20881 = r20879 / r20880; double r20882 = c; double r20883 = b; double r20884 = -r20883; double r20885 = r20883 * r20883; double r20886 = a; double r20887 = r20882 * r20886; double r20888 = r20887 * r20879; double r20889 = r20885 - r20888; double r20890 = sqrt(r20889); double r20891 = r20884 - r20890; double r20892 = r20882 / r20891; double r20893 = r20881 * r20892; return r20893; } 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 r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r20894); mpfr_init(r20895); mpfr_init(r20896); mpfr_init_set_str(r20897, "4", 10, MPFR_RNDN); mpfr_init(r20898); mpfr_init(r20899); mpfr_init(r20900); mpfr_init(r20901); mpfr_init(r20902); mpfr_init(r20903); mpfr_init(r20904); mpfr_init_set_str(r20905, "2", 10, MPFR_RNDN); mpfr_init(r20906); mpfr_init(r20907); } double f_im(double a, double b, double c) { mpfr_set_d(r20894, b, MPFR_RNDN); mpfr_neg(r20895, r20894, MPFR_RNDN); mpfr_mul(r20896, r20894, r20894, MPFR_RNDN); ; mpfr_set_d(r20898, a, MPFR_RNDN); mpfr_mul(r20899, r20897, r20898, MPFR_RNDN); mpfr_set_d(r20900, c, MPFR_RNDN); mpfr_mul(r20901, r20899, r20900, MPFR_RNDN); mpfr_sub(r20902, r20896, r20901, MPFR_RNDN); mpfr_sqrt(r20903, r20902, MPFR_RNDN); mpfr_add(r20904, r20895, r20903, MPFR_RNDN); ; mpfr_mul(r20906, r20905, r20898, MPFR_RNDN); mpfr_div(r20907, r20904, r20906, MPFR_RNDN); return mpfr_get_d(r20907, MPFR_RNDN); } static mpfr_t r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init_set_str(r20908, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20909, "2", 10, MPFR_RNDN); mpfr_init(r20910); mpfr_init(r20911); mpfr_init(r20912); mpfr_init(r20913); mpfr_init(r20914); mpfr_init(r20915); mpfr_init(r20916); mpfr_init(r20917); mpfr_init(r20918); mpfr_init(r20919); mpfr_init(r20920); mpfr_init(r20921); mpfr_init(r20922); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r20910, r20908, r20909, MPFR_RNDN); mpfr_set_d(r20911, c, MPFR_RNDN); mpfr_set_d(r20912, b, MPFR_RNDN); mpfr_neg(r20913, r20912, MPFR_RNDN); mpfr_mul(r20914, r20912, r20912, MPFR_RNDN); mpfr_set_d(r20915, a, MPFR_RNDN); mpfr_mul(r20916, r20911, r20915, MPFR_RNDN); mpfr_mul(r20917, r20916, r20908, MPFR_RNDN); mpfr_sub(r20918, r20914, r20917, MPFR_RNDN); mpfr_sqrt(r20919, r20918, MPFR_RNDN); mpfr_sub(r20920, r20913, r20919, MPFR_RNDN); mpfr_div(r20921, r20911, r20920, MPFR_RNDN); mpfr_mul(r20922, r20910, r20921, MPFR_RNDN); return mpfr_get_d(r20922, MPFR_RNDN); } static mpfr_t r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init_set_str(r20923, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20924, "2", 10, MPFR_RNDN); mpfr_init(r20925); mpfr_init(r20926); mpfr_init(r20927); mpfr_init(r20928); mpfr_init(r20929); mpfr_init(r20930); mpfr_init(r20931); mpfr_init(r20932); mpfr_init(r20933); mpfr_init(r20934); mpfr_init(r20935); mpfr_init(r20936); mpfr_init(r20937); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r20925, r20923, r20924, MPFR_RNDN); mpfr_set_d(r20926, c, MPFR_RNDN); mpfr_set_d(r20927, b, MPFR_RNDN); mpfr_neg(r20928, r20927, MPFR_RNDN); mpfr_mul(r20929, r20927, r20927, MPFR_RNDN); mpfr_set_d(r20930, a, MPFR_RNDN); mpfr_mul(r20931, r20926, r20930, MPFR_RNDN); mpfr_mul(r20932, r20931, r20923, MPFR_RNDN); mpfr_sub(r20933, r20929, r20932, MPFR_RNDN); mpfr_sqrt(r20934, r20933, MPFR_RNDN); mpfr_sub(r20935, r20928, r20934, MPFR_RNDN); mpfr_div(r20936, r20926, r20935, MPFR_RNDN); mpfr_mul(r20937, r20925, r20936, MPFR_RNDN); return mpfr_get_d(r20937, MPFR_RNDN); }