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