#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17814 = a; float r17815 = r17814 * r17814; float r17816 = b; float r17817 = r17816 * r17816; float r17818 = r17815 + r17817; float r17819 = r17818 * r17818; float r17820 = 4.0f; float r17821 = 1.0f; float r17822 = r17821 + r17814; float r17823 = r17815 * r17822; float r17824 = 3.0f; float r17825 = r17824 * r17814; float r17826 = r17821 - r17825; float r17827 = r17817 * r17826; float r17828 = r17823 + r17827; float r17829 = r17820 * r17828; float r17830 = r17819 + r17829; float r17831 = r17830 - r17821; return r17831; } double f_id(double a, double b) { double r17832 = a; double r17833 = r17832 * r17832; double r17834 = b; double r17835 = r17834 * r17834; double r17836 = r17833 + r17835; double r17837 = r17836 * r17836; double r17838 = 4.0; double r17839 = 1.0; double r17840 = r17839 + r17832; double r17841 = r17833 * r17840; double r17842 = 3.0; double r17843 = r17842 * r17832; double r17844 = r17839 - r17843; double r17845 = r17835 * r17844; double r17846 = r17841 + r17845; double r17847 = r17838 * r17846; double r17848 = r17837 + r17847; double r17849 = r17848 - r17839; return r17849; } double f_of(float a, float b) { float r17850 = 1.0f; float r17851 = a; float r17852 = 3.0f; float r17853 = r17851 * r17852; float r17854 = r17850 - r17853; float r17855 = b; float r17856 = r17855 * r17855; float r17857 = fma(r17851, r17851, r17851); float r17858 = r17857 * r17851; float r17859 = fma(r17854, r17856, r17858); float r17860 = 4.0f; float r17861 = 2.0f; float r17862 = r17855 * r17861; float r17863 = r17862 * r17855; float r17864 = r17851 * r17851; float r17865 = pow(r17851, r17860); float r17866 = pow(r17855, r17860); float r17867 = r17865 + r17866; float r17868 = fma(r17863, r17864, r17867); float r17869 = fma(r17859, r17860, r17868); float r17870 = r17869 - r17850; return r17870; } double f_od(double a, double b) { double r17871 = 1.0; double r17872 = a; double r17873 = 3.0; double r17874 = r17872 * r17873; double r17875 = r17871 - r17874; double r17876 = b; double r17877 = r17876 * r17876; double r17878 = fma(r17872, r17872, r17872); double r17879 = r17878 * r17872; double r17880 = fma(r17875, r17877, r17879); double r17881 = 4.0; double r17882 = 2.0; double r17883 = r17876 * r17882; double r17884 = r17883 * r17876; double r17885 = r17872 * r17872; double r17886 = pow(r17872, r17881); double r17887 = pow(r17876, r17881); double r17888 = r17886 + r17887; double r17889 = fma(r17884, r17885, r17888); double r17890 = fma(r17880, r17881, r17889); double r17891 = r17890 - r17871; return r17891; } 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 r17892, r17893, r17894, r17895, r17896, r17897, r17898, r17899, r17900, r17901, r17902, r17903, r17904, r17905, r17906, r17907, r17908, r17909; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17892); mpfr_init(r17893); mpfr_init(r17894); mpfr_init(r17895); mpfr_init(r17896); mpfr_init(r17897); mpfr_init_set_str(r17898, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17899, "1", 10, MPFR_RNDN); mpfr_init(r17900); mpfr_init(r17901); mpfr_init_set_str(r17902, "3", 10, MPFR_RNDN); mpfr_init(r17903); mpfr_init(r17904); mpfr_init(r17905); mpfr_init(r17906); mpfr_init(r17907); mpfr_init(r17908); mpfr_init(r17909); } double f_im(double a, double b) { mpfr_set_d(r17892, a, MPFR_RNDN); mpfr_sqr(r17893, r17892, MPFR_RNDN); mpfr_set_d(r17894, b, MPFR_RNDN); mpfr_sqr(r17895, r17894, MPFR_RNDN); mpfr_add(r17896, r17893, r17895, MPFR_RNDN); mpfr_sqr(r17897, r17896, MPFR_RNDN); ; ; mpfr_add(r17900, r17899, r17892, MPFR_RNDN); mpfr_mul(r17901, r17893, r17900, MPFR_RNDN); ; mpfr_mul(r17903, r17902, r17892, MPFR_RNDN); mpfr_sub(r17904, r17899, r17903, MPFR_RNDN); mpfr_mul(r17905, r17895, r17904, MPFR_RNDN); mpfr_add(r17906, r17901, r17905, MPFR_RNDN); mpfr_mul(r17907, r17898, r17906, MPFR_RNDN); mpfr_add(r17908, r17897, r17907, MPFR_RNDN); mpfr_sub(r17909, r17908, r17899, MPFR_RNDN); return mpfr_get_d(r17909, MPFR_RNDN); } static mpfr_t r17910, r17911, r17912, r17913, r17914, r17915, r17916, r17917, r17918, r17919, r17920, r17921, r17922, r17923, r17924, r17925, r17926, r17927, r17928, r17929, r17930; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17910, "1", 10, MPFR_RNDN); mpfr_init(r17911); mpfr_init_set_str(r17912, "3", 10, MPFR_RNDN); mpfr_init(r17913); mpfr_init(r17914); mpfr_init(r17915); mpfr_init(r17916); mpfr_init(r17917); mpfr_init(r17918); mpfr_init(r17919); mpfr_init_set_str(r17920, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17921, "2", 10, MPFR_RNDN); mpfr_init(r17922); mpfr_init(r17923); mpfr_init(r17924); mpfr_init(r17925); mpfr_init(r17926); mpfr_init(r17927); mpfr_init(r17928); mpfr_init(r17929); mpfr_init(r17930); } double f_fm(double a, double b) { ; mpfr_set_d(r17911, a, MPFR_RNDN); ; mpfr_mul(r17913, r17911, r17912, MPFR_RNDN); mpfr_sub(r17914, r17910, r17913, MPFR_RNDN); mpfr_set_d(r17915, b, MPFR_RNDN); mpfr_sqr(r17916, r17915, MPFR_RNDN); mpfr_fma(r17917, r17911, r17911, r17911, MPFR_RNDN); mpfr_mul(r17918, r17917, r17911, MPFR_RNDN); mpfr_fma(r17919, r17914, r17916, r17918, MPFR_RNDN); ; ; mpfr_mul(r17922, r17915, r17921, MPFR_RNDN); mpfr_mul(r17923, r17922, r17915, MPFR_RNDN); mpfr_mul(r17924, r17911, r17911, MPFR_RNDN); mpfr_pow(r17925, r17911, r17920, MPFR_RNDN); mpfr_pow(r17926, r17915, r17920, MPFR_RNDN); mpfr_add(r17927, r17925, r17926, MPFR_RNDN); mpfr_fma(r17928, r17923, r17924, r17927, MPFR_RNDN); mpfr_fma(r17929, r17919, r17920, r17928, MPFR_RNDN); mpfr_sub(r17930, r17929, r17910, MPFR_RNDN); return mpfr_get_d(r17930, MPFR_RNDN); } static mpfr_t r17931, r17932, r17933, r17934, r17935, r17936, r17937, r17938, r17939, r17940, r17941, r17942, r17943, r17944, r17945, r17946, r17947, r17948, r17949, r17950, r17951; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17931, "1", 10, MPFR_RNDN); mpfr_init(r17932); mpfr_init_set_str(r17933, "3", 10, MPFR_RNDN); mpfr_init(r17934); mpfr_init(r17935); mpfr_init(r17936); mpfr_init(r17937); mpfr_init(r17938); mpfr_init(r17939); mpfr_init(r17940); mpfr_init_set_str(r17941, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17942, "2", 10, MPFR_RNDN); mpfr_init(r17943); mpfr_init(r17944); mpfr_init(r17945); mpfr_init(r17946); mpfr_init(r17947); mpfr_init(r17948); mpfr_init(r17949); mpfr_init(r17950); mpfr_init(r17951); } double f_dm(double a, double b) { ; mpfr_set_d(r17932, a, MPFR_RNDN); ; mpfr_mul(r17934, r17932, r17933, MPFR_RNDN); mpfr_sub(r17935, r17931, r17934, MPFR_RNDN); mpfr_set_d(r17936, b, MPFR_RNDN); mpfr_sqr(r17937, r17936, MPFR_RNDN); mpfr_fma(r17938, r17932, r17932, r17932, MPFR_RNDN); mpfr_mul(r17939, r17938, r17932, MPFR_RNDN); mpfr_fma(r17940, r17935, r17937, r17939, MPFR_RNDN); ; ; mpfr_mul(r17943, r17936, r17942, MPFR_RNDN); mpfr_mul(r17944, r17943, r17936, MPFR_RNDN); mpfr_mul(r17945, r17932, r17932, MPFR_RNDN); mpfr_pow(r17946, r17932, r17941, MPFR_RNDN); mpfr_pow(r17947, r17936, r17941, MPFR_RNDN); mpfr_add(r17948, r17946, r17947, MPFR_RNDN); mpfr_fma(r17949, r17944, r17945, r17948, MPFR_RNDN); mpfr_fma(r17950, r17940, r17941, r17949, MPFR_RNDN); mpfr_sub(r17951, r17950, r17931, MPFR_RNDN); return mpfr_get_d(r17951, MPFR_RNDN); }