#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r20807 = b; float r20808 = -r20807; float r20809 = r20807 * r20807; float r20810 = 4; float r20811 = a; float r20812 = r20810 * r20811; float r20813 = c; float r20814 = r20812 * r20813; float r20815 = r20809 - r20814; float r20816 = sqrt(r20815); float r20817 = r20808 + r20816; float r20818 = 2; float r20819 = r20818 * r20811; float r20820 = r20817 / r20819; return r20820; } double f_id(double a, double b, double c) { double r20821 = b; double r20822 = -r20821; double r20823 = r20821 * r20821; double r20824 = 4; double r20825 = a; double r20826 = r20824 * r20825; double r20827 = c; double r20828 = r20826 * r20827; double r20829 = r20823 - r20828; double r20830 = sqrt(r20829); double r20831 = r20822 + r20830; double r20832 = 2; double r20833 = r20832 * r20825; double r20834 = r20831 / r20833; return r20834; } double f_of(float a, float b, float c) { float r20835 = 4; float r20836 = 2; float r20837 = r20835 / r20836; float r20838 = c; float r20839 = b; float r20840 = -r20839; float r20841 = r20839 * r20839; float r20842 = a; float r20843 = r20838 * r20842; float r20844 = r20843 * r20835; float r20845 = r20841 - r20844; float r20846 = sqrt(r20845); float r20847 = r20840 - r20846; float r20848 = r20838 / r20847; float r20849 = r20837 * r20848; return r20849; } double f_od(double a, double b, double c) { double r20850 = 4; double r20851 = 2; double r20852 = r20850 / r20851; double r20853 = c; double r20854 = b; double r20855 = -r20854; double r20856 = r20854 * r20854; double r20857 = a; double r20858 = r20853 * r20857; double r20859 = r20858 * r20850; double r20860 = r20856 - r20859; double r20861 = sqrt(r20860); double r20862 = r20855 - r20861; double r20863 = r20853 / r20862; double r20864 = r20852 * r20863; return r20864; } 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 r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872, r20873, r20874, r20875, r20876, r20877, r20878; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r20865); mpfr_init(r20866); mpfr_init(r20867); mpfr_init_set_str(r20868, "4", 10, MPFR_RNDN); mpfr_init(r20869); mpfr_init(r20870); mpfr_init(r20871); mpfr_init(r20872); mpfr_init(r20873); mpfr_init(r20874); mpfr_init(r20875); mpfr_init_set_str(r20876, "2", 10, MPFR_RNDN); mpfr_init(r20877); mpfr_init(r20878); } double f_im(double a, double b, double c) { mpfr_set_d(r20865, b, MPFR_RNDN); mpfr_neg(r20866, r20865, MPFR_RNDN); mpfr_mul(r20867, r20865, r20865, MPFR_RNDN); ; mpfr_set_d(r20869, a, MPFR_RNDN); mpfr_mul(r20870, r20868, r20869, MPFR_RNDN); mpfr_set_d(r20871, c, MPFR_RNDN); mpfr_mul(r20872, r20870, r20871, MPFR_RNDN); mpfr_sub(r20873, r20867, r20872, MPFR_RNDN); mpfr_sqrt(r20874, r20873, MPFR_RNDN); mpfr_add(r20875, r20866, r20874, MPFR_RNDN); ; mpfr_mul(r20877, r20876, r20869, MPFR_RNDN); mpfr_div(r20878, r20875, r20877, MPFR_RNDN); return mpfr_get_d(r20878, MPFR_RNDN); } static mpfr_t r20879, r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init_set_str(r20879, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20880, "2", 10, MPFR_RNDN); mpfr_init(r20881); mpfr_init(r20882); mpfr_init(r20883); mpfr_init(r20884); mpfr_init(r20885); mpfr_init(r20886); mpfr_init(r20887); mpfr_init(r20888); mpfr_init(r20889); mpfr_init(r20890); mpfr_init(r20891); mpfr_init(r20892); mpfr_init(r20893); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r20881, r20879, r20880, MPFR_RNDN); mpfr_set_d(r20882, c, MPFR_RNDN); mpfr_set_d(r20883, b, MPFR_RNDN); mpfr_neg(r20884, r20883, MPFR_RNDN); mpfr_mul(r20885, r20883, r20883, MPFR_RNDN); mpfr_set_d(r20886, a, MPFR_RNDN); mpfr_mul(r20887, r20882, r20886, MPFR_RNDN); mpfr_mul(r20888, r20887, r20879, MPFR_RNDN); mpfr_sub(r20889, r20885, r20888, MPFR_RNDN); mpfr_sqrt(r20890, r20889, MPFR_RNDN); mpfr_sub(r20891, r20884, r20890, MPFR_RNDN); mpfr_div(r20892, r20882, r20891, MPFR_RNDN); mpfr_mul(r20893, r20881, r20892, MPFR_RNDN); return mpfr_get_d(r20893, MPFR_RNDN); } static mpfr_t r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init_set_str(r20894, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20895, "2", 10, MPFR_RNDN); mpfr_init(r20896); mpfr_init(r20897); mpfr_init(r20898); mpfr_init(r20899); mpfr_init(r20900); mpfr_init(r20901); mpfr_init(r20902); mpfr_init(r20903); mpfr_init(r20904); mpfr_init(r20905); mpfr_init(r20906); mpfr_init(r20907); mpfr_init(r20908); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r20896, r20894, r20895, MPFR_RNDN); mpfr_set_d(r20897, c, MPFR_RNDN); mpfr_set_d(r20898, b, MPFR_RNDN); mpfr_neg(r20899, r20898, MPFR_RNDN); mpfr_mul(r20900, r20898, r20898, MPFR_RNDN); mpfr_set_d(r20901, a, MPFR_RNDN); mpfr_mul(r20902, r20897, r20901, MPFR_RNDN); mpfr_mul(r20903, r20902, r20894, MPFR_RNDN); mpfr_sub(r20904, r20900, r20903, MPFR_RNDN); mpfr_sqrt(r20905, r20904, MPFR_RNDN); mpfr_sub(r20906, r20899, r20905, MPFR_RNDN); mpfr_div(r20907, r20897, r20906, MPFR_RNDN); mpfr_mul(r20908, r20896, r20907, MPFR_RNDN); return mpfr_get_d(r20908, MPFR_RNDN); }