#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r7841 = b; float r7842 = -r7841; float r7843 = r7841 * r7841; float r7844 = 4.0f; float r7845 = a; float r7846 = r7844 * r7845; float r7847 = c; float r7848 = r7846 * r7847; float r7849 = r7843 - r7848; float r7850 = sqrt(r7849); float r7851 = r7842 + r7850; float r7852 = 2.0f; float r7853 = r7852 * r7845; float r7854 = r7851 / r7853; return r7854; } double f_id(double a, double b, double c) { double r7855 = b; double r7856 = -r7855; double r7857 = r7855 * r7855; double r7858 = 4.0; double r7859 = a; double r7860 = r7858 * r7859; double r7861 = c; double r7862 = r7860 * r7861; double r7863 = r7857 - r7862; double r7864 = sqrt(r7863); double r7865 = r7856 + r7864; double r7866 = 2.0; double r7867 = r7866 * r7859; double r7868 = r7865 / r7867; return r7868; } double f_of(float a, float b, float c) { float r7869 = a; float r7870 = c; float r7871 = r7869 * r7870; float r7872 = 4.0f; float r7873 = r7871 * r7872; float r7874 = 2.0f; float r7875 = r7869 * r7874; float r7876 = r7873 / r7875; float r7877 = b; float r7878 = -r7877; float r7879 = r7872 * r7869; float r7880 = -r7870; float r7881 = r7879 * r7880; float r7882 = fma(r7877, r7877, r7881); float r7883 = sqrt(r7882); float r7884 = r7878 - r7883; float r7885 = r7876 / r7884; return r7885; } double f_od(double a, double b, double c) { double r7886 = a; double r7887 = c; double r7888 = r7886 * r7887; double r7889 = 4.0; double r7890 = r7888 * r7889; double r7891 = 2.0; double r7892 = r7886 * r7891; double r7893 = r7890 / r7892; double r7894 = b; double r7895 = -r7894; double r7896 = r7889 * r7886; double r7897 = -r7887; double r7898 = r7896 * r7897; double r7899 = fma(r7894, r7894, r7898); double r7900 = sqrt(r7899); double r7901 = r7895 - r7900; double r7902 = r7893 / r7901; return r7902; } 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 r7903, r7904, r7905, r7906, r7907, r7908, r7909, r7910, r7911, r7912, r7913, r7914, r7915, r7916; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r7903); mpfr_init(r7904); mpfr_init(r7905); mpfr_init_set_str(r7906, "4", 10, MPFR_RNDN); mpfr_init(r7907); mpfr_init(r7908); mpfr_init(r7909); mpfr_init(r7910); mpfr_init(r7911); mpfr_init(r7912); mpfr_init(r7913); mpfr_init_set_str(r7914, "2", 10, MPFR_RNDN); mpfr_init(r7915); mpfr_init(r7916); } double f_im(double a, double b, double c) { mpfr_set_d(r7903, b, MPFR_RNDN); mpfr_neg(r7904, r7903, MPFR_RNDN); mpfr_mul(r7905, r7903, r7903, MPFR_RNDN); ; mpfr_set_d(r7907, a, MPFR_RNDN); mpfr_mul(r7908, r7906, r7907, MPFR_RNDN); mpfr_set_d(r7909, c, MPFR_RNDN); mpfr_mul(r7910, r7908, r7909, MPFR_RNDN); mpfr_sub(r7911, r7905, r7910, MPFR_RNDN); mpfr_sqrt(r7912, r7911, MPFR_RNDN); mpfr_add(r7913, r7904, r7912, MPFR_RNDN); ; mpfr_mul(r7915, r7914, r7907, MPFR_RNDN); mpfr_div(r7916, r7913, r7915, MPFR_RNDN); return mpfr_get_d(r7916, MPFR_RNDN); } static mpfr_t r7917, r7918, r7919, r7920, r7921, r7922, r7923, r7924, r7925, r7926, r7927, r7928, r7929, r7930, r7931, r7932, r7933; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r7917); mpfr_init(r7918); mpfr_init(r7919); mpfr_init_set_str(r7920, "4", 10, MPFR_RNDN); mpfr_init(r7921); mpfr_init_set_str(r7922, "2", 10, MPFR_RNDN); 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); mpfr_init(r7932); mpfr_init(r7933); } double f_fm(double a, double b, double c) { mpfr_set_d(r7917, a, MPFR_RNDN); mpfr_set_d(r7918, c, MPFR_RNDN); mpfr_mul(r7919, r7917, r7918, MPFR_RNDN); ; mpfr_mul(r7921, r7919, r7920, MPFR_RNDN); ; mpfr_mul(r7923, r7917, r7922, MPFR_RNDN); mpfr_div(r7924, r7921, r7923, MPFR_RNDN); mpfr_set_d(r7925, b, MPFR_RNDN); mpfr_neg(r7926, r7925, MPFR_RNDN); mpfr_mul(r7927, r7920, r7917, MPFR_RNDN); mpfr_neg(r7928, r7918, MPFR_RNDN); mpfr_mul(r7929, r7927, r7928, MPFR_RNDN); mpfr_fma(r7930, r7925, r7925, r7929, MPFR_RNDN); mpfr_sqrt(r7931, r7930, MPFR_RNDN); mpfr_sub(r7932, r7926, r7931, MPFR_RNDN); mpfr_div(r7933, r7924, r7932, MPFR_RNDN); return mpfr_get_d(r7933, MPFR_RNDN); } static mpfr_t r7934, r7935, r7936, r7937, r7938, r7939, r7940, r7941, r7942, r7943, r7944, r7945, r7946, r7947, r7948, r7949, r7950; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r7934); mpfr_init(r7935); mpfr_init(r7936); mpfr_init_set_str(r7937, "4", 10, MPFR_RNDN); mpfr_init(r7938); mpfr_init_set_str(r7939, "2", 10, MPFR_RNDN); mpfr_init(r7940); mpfr_init(r7941); mpfr_init(r7942); mpfr_init(r7943); mpfr_init(r7944); mpfr_init(r7945); mpfr_init(r7946); mpfr_init(r7947); mpfr_init(r7948); mpfr_init(r7949); mpfr_init(r7950); } double f_dm(double a, double b, double c) { mpfr_set_d(r7934, a, MPFR_RNDN); mpfr_set_d(r7935, c, MPFR_RNDN); mpfr_mul(r7936, r7934, r7935, MPFR_RNDN); ; mpfr_mul(r7938, r7936, r7937, MPFR_RNDN); ; mpfr_mul(r7940, r7934, r7939, MPFR_RNDN); mpfr_div(r7941, r7938, r7940, MPFR_RNDN); mpfr_set_d(r7942, b, MPFR_RNDN); mpfr_neg(r7943, r7942, MPFR_RNDN); mpfr_mul(r7944, r7937, r7934, MPFR_RNDN); mpfr_neg(r7945, r7935, MPFR_RNDN); mpfr_mul(r7946, r7944, r7945, MPFR_RNDN); mpfr_fma(r7947, r7942, r7942, r7946, MPFR_RNDN); mpfr_sqrt(r7948, r7947, MPFR_RNDN); mpfr_sub(r7949, r7943, r7948, MPFR_RNDN); mpfr_div(r7950, r7941, r7949, MPFR_RNDN); return mpfr_get_d(r7950, MPFR_RNDN); }