#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r16852 = b; float r16853 = -r16852; float r16854 = r16852 * r16852; float r16855 = 4.0f; float r16856 = a; float r16857 = r16855 * r16856; float r16858 = c; float r16859 = r16857 * r16858; float r16860 = r16854 - r16859; float r16861 = sqrt(r16860); float r16862 = r16853 + r16861; float r16863 = 2.0f; float r16864 = r16863 * r16856; float r16865 = r16862 / r16864; return r16865; } double f_id(double a, double b, double c) { double r16866 = b; double r16867 = -r16866; double r16868 = r16866 * r16866; double r16869 = 4.0; double r16870 = a; double r16871 = r16869 * r16870; double r16872 = c; double r16873 = r16871 * r16872; double r16874 = r16868 - r16873; double r16875 = sqrt(r16874); double r16876 = r16867 + r16875; double r16877 = 2.0; double r16878 = r16877 * r16870; double r16879 = r16876 / r16878; return r16879; } double f_of(float a, float b, float c) { float r16880 = 4.0f; float r16881 = a; float r16882 = r16880 * r16881; float r16883 = c; float r16884 = r16882 * r16883; float r16885 = b; float r16886 = -r16885; float r16887 = r16885 * r16885; float r16888 = r16887 - r16884; float r16889 = sqrt(r16888); float r16890 = r16886 - r16889; float r16891 = r16884 / r16890; float r16892 = 2.0f; float r16893 = r16892 * r16881; float r16894 = r16891 / r16893; return r16894; } double f_od(double a, double b, double c) { double r16895 = 4.0; double r16896 = a; double r16897 = r16895 * r16896; double r16898 = c; double r16899 = r16897 * r16898; double r16900 = b; double r16901 = -r16900; double r16902 = r16900 * r16900; double r16903 = r16902 - r16899; double r16904 = sqrt(r16903); double r16905 = r16901 - r16904; double r16906 = r16899 / r16905; double r16907 = 2.0; double r16908 = r16907 * r16896; double r16909 = r16906 / r16908; return r16909; } 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 r16910, r16911, r16912, r16913, r16914, r16915, r16916, r16917, r16918, r16919, r16920, r16921, r16922, r16923; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16910); mpfr_init(r16911); mpfr_init(r16912); mpfr_init_set_str(r16913, "4", 10, MPFR_RNDN); mpfr_init(r16914); mpfr_init(r16915); mpfr_init(r16916); mpfr_init(r16917); mpfr_init(r16918); mpfr_init(r16919); mpfr_init(r16920); mpfr_init_set_str(r16921, "2", 10, MPFR_RNDN); mpfr_init(r16922); mpfr_init(r16923); } double f_im(double a, double b, double c) { mpfr_set_d(r16910, b, MPFR_RNDN); mpfr_neg(r16911, r16910, MPFR_RNDN); mpfr_sqr(r16912, r16910, MPFR_RNDN); ; mpfr_set_d(r16914, a, MPFR_RNDN); mpfr_mul(r16915, r16913, r16914, MPFR_RNDN); mpfr_set_d(r16916, c, MPFR_RNDN); mpfr_mul(r16917, r16915, r16916, MPFR_RNDN); mpfr_sub(r16918, r16912, r16917, MPFR_RNDN); mpfr_sqrt(r16919, r16918, MPFR_RNDN); mpfr_add(r16920, r16911, r16919, MPFR_RNDN); ; mpfr_mul(r16922, r16921, r16914, MPFR_RNDN); mpfr_div(r16923, r16920, r16922, MPFR_RNDN); return mpfr_get_d(r16923, MPFR_RNDN); } static mpfr_t r16924, r16925, r16926, r16927, r16928, r16929, r16930, r16931, r16932, r16933, r16934, r16935, r16936, r16937, r16938; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r16924, "4", 10, MPFR_RNDN); mpfr_init(r16925); mpfr_init(r16926); mpfr_init(r16927); mpfr_init(r16928); mpfr_init(r16929); mpfr_init(r16930); mpfr_init(r16931); mpfr_init(r16932); mpfr_init(r16933); mpfr_init(r16934); mpfr_init(r16935); mpfr_init_set_str(r16936, "2", 10, MPFR_RNDN); mpfr_init(r16937); mpfr_init(r16938); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r16925, a, MPFR_RNDN); mpfr_mul(r16926, r16924, r16925, MPFR_RNDN); mpfr_set_d(r16927, c, MPFR_RNDN); mpfr_mul(r16928, r16926, r16927, MPFR_RNDN); mpfr_set_d(r16929, b, MPFR_RNDN); mpfr_neg(r16930, r16929, MPFR_RNDN); mpfr_sqr(r16931, r16929, MPFR_RNDN); mpfr_sub(r16932, r16931, r16928, MPFR_RNDN); mpfr_sqrt(r16933, r16932, MPFR_RNDN); mpfr_sub(r16934, r16930, r16933, MPFR_RNDN); mpfr_div(r16935, r16928, r16934, MPFR_RNDN); ; mpfr_mul(r16937, r16936, r16925, MPFR_RNDN); mpfr_div(r16938, r16935, r16937, MPFR_RNDN); return mpfr_get_d(r16938, MPFR_RNDN); } static mpfr_t r16939, r16940, r16941, r16942, r16943, r16944, r16945, r16946, r16947, r16948, r16949, r16950, r16951, r16952, r16953; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r16939, "4", 10, MPFR_RNDN); mpfr_init(r16940); mpfr_init(r16941); mpfr_init(r16942); mpfr_init(r16943); mpfr_init(r16944); mpfr_init(r16945); mpfr_init(r16946); mpfr_init(r16947); mpfr_init(r16948); mpfr_init(r16949); mpfr_init(r16950); mpfr_init_set_str(r16951, "2", 10, MPFR_RNDN); mpfr_init(r16952); mpfr_init(r16953); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r16940, a, MPFR_RNDN); mpfr_mul(r16941, r16939, r16940, MPFR_RNDN); mpfr_set_d(r16942, c, MPFR_RNDN); mpfr_mul(r16943, r16941, r16942, MPFR_RNDN); mpfr_set_d(r16944, b, MPFR_RNDN); mpfr_neg(r16945, r16944, MPFR_RNDN); mpfr_sqr(r16946, r16944, MPFR_RNDN); mpfr_sub(r16947, r16946, r16943, MPFR_RNDN); mpfr_sqrt(r16948, r16947, MPFR_RNDN); mpfr_sub(r16949, r16945, r16948, MPFR_RNDN); mpfr_div(r16950, r16943, r16949, MPFR_RNDN); ; mpfr_mul(r16952, r16951, r16940, MPFR_RNDN); mpfr_div(r16953, r16950, r16952, MPFR_RNDN); return mpfr_get_d(r16953, MPFR_RNDN); }