#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r7822 = b; float r7823 = -r7822; float r7824 = r7822 * r7822; float r7825 = 4.0f; float r7826 = a; float r7827 = r7825 * r7826; float r7828 = c; float r7829 = r7827 * r7828; float r7830 = r7824 - r7829; float r7831 = sqrt(r7830); float r7832 = r7823 + r7831; float r7833 = 2.0f; float r7834 = r7833 * r7826; float r7835 = r7832 / r7834; return r7835; } double f_id(double a, double b, double c) { double r7836 = b; double r7837 = -r7836; double r7838 = r7836 * r7836; double r7839 = 4.0; double r7840 = a; double r7841 = r7839 * r7840; double r7842 = c; double r7843 = r7841 * r7842; double r7844 = r7838 - r7843; double r7845 = sqrt(r7844); double r7846 = r7837 + r7845; double r7847 = 2.0; double r7848 = r7847 * r7840; double r7849 = r7846 / r7848; return r7849; } double f_of(float a, float b, float c) { float r7850 = a; float r7851 = c; float r7852 = r7850 * r7851; float r7853 = 4.0f; float r7854 = r7852 * r7853; float r7855 = 2.0f; float r7856 = r7850 * r7855; float r7857 = r7854 / r7856; float r7858 = b; float r7859 = -r7858; float r7860 = r7858 * r7858; float r7861 = r7853 * r7850; float r7862 = r7861 * r7851; float r7863 = r7860 - r7862; float r7864 = sqrt(r7863); float r7865 = r7859 - r7864; float r7866 = r7857 / r7865; return r7866; } double f_od(double a, double b, double c) { double r7867 = a; double r7868 = c; double r7869 = r7867 * r7868; double r7870 = 4.0; double r7871 = r7869 * r7870; double r7872 = 2.0; double r7873 = r7867 * r7872; double r7874 = r7871 / r7873; double r7875 = b; double r7876 = -r7875; double r7877 = r7875 * r7875; double r7878 = r7870 * r7867; double r7879 = r7878 * r7868; double r7880 = r7877 - r7879; double r7881 = sqrt(r7880); double r7882 = r7876 - r7881; double r7883 = r7874 / r7882; return r7883; } 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 r7884, r7885, r7886, r7887, r7888, r7889, r7890, r7891, r7892, r7893, r7894, r7895, r7896, r7897; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r7884); mpfr_init(r7885); mpfr_init(r7886); mpfr_init_set_str(r7887, "4", 10, MPFR_RNDN); mpfr_init(r7888); mpfr_init(r7889); mpfr_init(r7890); mpfr_init(r7891); mpfr_init(r7892); mpfr_init(r7893); mpfr_init(r7894); mpfr_init_set_str(r7895, "2", 10, MPFR_RNDN); mpfr_init(r7896); mpfr_init(r7897); } double f_im(double a, double b, double c) { mpfr_set_d(r7884, b, MPFR_RNDN); mpfr_neg(r7885, r7884, MPFR_RNDN); mpfr_mul(r7886, r7884, r7884, MPFR_RNDN); ; mpfr_set_d(r7888, a, MPFR_RNDN); mpfr_mul(r7889, r7887, r7888, MPFR_RNDN); mpfr_set_d(r7890, c, MPFR_RNDN); mpfr_mul(r7891, r7889, r7890, MPFR_RNDN); mpfr_sub(r7892, r7886, r7891, MPFR_RNDN); mpfr_sqrt(r7893, r7892, MPFR_RNDN); mpfr_add(r7894, r7885, r7893, MPFR_RNDN); ; mpfr_mul(r7896, r7895, r7888, MPFR_RNDN); mpfr_div(r7897, r7894, r7896, MPFR_RNDN); return mpfr_get_d(r7897, MPFR_RNDN); } static mpfr_t r7898, r7899, r7900, r7901, r7902, r7903, r7904, r7905, r7906, r7907, r7908, r7909, r7910, r7911, r7912, r7913, r7914; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r7898); mpfr_init(r7899); mpfr_init(r7900); mpfr_init_set_str(r7901, "4", 10, MPFR_RNDN); mpfr_init(r7902); mpfr_init_set_str(r7903, "2", 10, MPFR_RNDN); mpfr_init(r7904); mpfr_init(r7905); mpfr_init(r7906); mpfr_init(r7907); mpfr_init(r7908); mpfr_init(r7909); mpfr_init(r7910); mpfr_init(r7911); mpfr_init(r7912); mpfr_init(r7913); mpfr_init(r7914); } double f_fm(double a, double b, double c) { mpfr_set_d(r7898, a, MPFR_RNDN); mpfr_set_d(r7899, c, MPFR_RNDN); mpfr_mul(r7900, r7898, r7899, MPFR_RNDN); ; mpfr_mul(r7902, r7900, r7901, MPFR_RNDN); ; mpfr_mul(r7904, r7898, r7903, MPFR_RNDN); mpfr_div(r7905, r7902, r7904, MPFR_RNDN); mpfr_set_d(r7906, b, MPFR_RNDN); mpfr_neg(r7907, r7906, MPFR_RNDN); mpfr_mul(r7908, r7906, r7906, MPFR_RNDN); mpfr_mul(r7909, r7901, r7898, MPFR_RNDN); mpfr_mul(r7910, r7909, r7899, MPFR_RNDN); mpfr_sub(r7911, r7908, r7910, MPFR_RNDN); mpfr_sqrt(r7912, r7911, MPFR_RNDN); mpfr_sub(r7913, r7907, r7912, MPFR_RNDN); mpfr_div(r7914, r7905, r7913, MPFR_RNDN); return mpfr_get_d(r7914, MPFR_RNDN); } static mpfr_t r7915, r7916, r7917, r7918, r7919, r7920, r7921, r7922, r7923, r7924, r7925, r7926, r7927, r7928, r7929, r7930, r7931; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r7915); mpfr_init(r7916); mpfr_init(r7917); mpfr_init_set_str(r7918, "4", 10, MPFR_RNDN); mpfr_init(r7919); mpfr_init_set_str(r7920, "2", 10, MPFR_RNDN); mpfr_init(r7921); mpfr_init(r7922); mpfr_init(r7923); mpfr_init(r7924); mpfr_init(r7925); mpfr_init(r7926); mpfr_init(r7927); mpfr_init(r7928); mpfr_init(r7929); mpfr_init(r7930); mpfr_init(r7931); } double f_dm(double a, double b, double c) { mpfr_set_d(r7915, a, MPFR_RNDN); mpfr_set_d(r7916, c, MPFR_RNDN); mpfr_mul(r7917, r7915, r7916, MPFR_RNDN); ; mpfr_mul(r7919, r7917, r7918, MPFR_RNDN); ; mpfr_mul(r7921, r7915, r7920, MPFR_RNDN); mpfr_div(r7922, r7919, r7921, MPFR_RNDN); mpfr_set_d(r7923, b, MPFR_RNDN); mpfr_neg(r7924, r7923, MPFR_RNDN); mpfr_mul(r7925, r7923, r7923, MPFR_RNDN); mpfr_mul(r7926, r7918, r7915, MPFR_RNDN); mpfr_mul(r7927, r7926, r7916, MPFR_RNDN); mpfr_sub(r7928, r7925, r7927, MPFR_RNDN); mpfr_sqrt(r7929, r7928, MPFR_RNDN); mpfr_sub(r7930, r7924, r7929, MPFR_RNDN); mpfr_div(r7931, r7922, r7930, MPFR_RNDN); return mpfr_get_d(r7931, MPFR_RNDN); }