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