#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r20842 = atan2(1.0, 0.0); float r20843 = 2; float r20844 = r20842 / r20843; float r20845 = 1; float r20846 = b; float r20847 = r20846 * r20846; float r20848 = a; float r20849 = r20848 * r20848; float r20850 = r20847 - r20849; float r20851 = r20845 / r20850; float r20852 = r20844 * r20851; float r20853 = r20845 / r20848; float r20854 = r20845 / r20846; float r20855 = r20853 - r20854; float r20856 = r20852 * r20855; return r20856; } double f_id(double a, double b) { double r20857 = atan2(1.0, 0.0); double r20858 = 2; double r20859 = r20857 / r20858; double r20860 = 1; double r20861 = b; double r20862 = r20861 * r20861; double r20863 = a; double r20864 = r20863 * r20863; double r20865 = r20862 - r20864; double r20866 = r20860 / r20865; double r20867 = r20859 * r20866; double r20868 = r20860 / r20863; double r20869 = r20860 / r20861; double r20870 = r20868 - r20869; double r20871 = r20867 * r20870; return r20871; } double f_of(float a, float b) { float r20872 = 1; float r20873 = a; float r20874 = r20872 / r20873; float r20875 = b; float r20876 = r20872 / r20875; float r20877 = r20874 - r20876; float r20878 = r20875 + r20873; float r20879 = r20877 / r20878; float r20880 = atan2(1.0, 0.0); float r20881 = r20879 * r20880; float r20882 = r20875 - r20873; float r20883 = 2; float r20884 = r20882 * r20883; float r20885 = r20881 / r20884; return r20885; } double f_od(double a, double b) { double r20886 = 1; double r20887 = a; double r20888 = r20886 / r20887; double r20889 = b; double r20890 = r20886 / r20889; double r20891 = r20888 - r20890; double r20892 = r20889 + r20887; double r20893 = r20891 / r20892; double r20894 = atan2(1.0, 0.0); double r20895 = r20893 * r20894; double r20896 = r20889 - r20887; double r20897 = 2; double r20898 = r20896 * r20897; double r20899 = r20895 / r20898; return r20899; } 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 r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20900); mpfr_init_set_str(r20901, "2", 10, MPFR_RNDN); mpfr_init(r20902); mpfr_init_set_str(r20903, "1", 10, MPFR_RNDN); mpfr_init(r20904); mpfr_init(r20905); mpfr_init(r20906); mpfr_init(r20907); mpfr_init(r20908); mpfr_init(r20909); mpfr_init(r20910); mpfr_init(r20911); mpfr_init(r20912); mpfr_init(r20913); mpfr_init(r20914); } double f_im(double a, double b) { mpfr_const_pi(r20900, MPFR_RNDN); ; mpfr_div(r20902, r20900, r20901, MPFR_RNDN); ; mpfr_set_d(r20904, b, MPFR_RNDN); mpfr_mul(r20905, r20904, r20904, MPFR_RNDN); mpfr_set_d(r20906, a, MPFR_RNDN); mpfr_mul(r20907, r20906, r20906, MPFR_RNDN); mpfr_sub(r20908, r20905, r20907, MPFR_RNDN); mpfr_div(r20909, r20903, r20908, MPFR_RNDN); mpfr_mul(r20910, r20902, r20909, MPFR_RNDN); mpfr_div(r20911, r20903, r20906, MPFR_RNDN); mpfr_div(r20912, r20903, r20904, MPFR_RNDN); mpfr_sub(r20913, r20911, r20912, MPFR_RNDN); mpfr_mul(r20914, r20910, r20913, MPFR_RNDN); return mpfr_get_d(r20914, MPFR_RNDN); } static mpfr_t r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922, r20923, r20924, r20925, r20926, r20927, r20928; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20915, "1", 10, MPFR_RNDN); mpfr_init(r20916); mpfr_init(r20917); mpfr_init(r20918); mpfr_init(r20919); mpfr_init(r20920); mpfr_init(r20921); mpfr_init(r20922); mpfr_init(r20923); mpfr_init(r20924); mpfr_init(r20925); mpfr_init_set_str(r20926, "2", 10, MPFR_RNDN); mpfr_init(r20927); mpfr_init(r20928); } double f_fm(double a, double b) { ; mpfr_set_d(r20916, a, MPFR_RNDN); mpfr_div(r20917, r20915, r20916, MPFR_RNDN); mpfr_set_d(r20918, b, MPFR_RNDN); mpfr_div(r20919, r20915, r20918, MPFR_RNDN); mpfr_sub(r20920, r20917, r20919, MPFR_RNDN); mpfr_add(r20921, r20918, r20916, MPFR_RNDN); mpfr_div(r20922, r20920, r20921, MPFR_RNDN); mpfr_const_pi(r20923, MPFR_RNDN); mpfr_mul(r20924, r20922, r20923, MPFR_RNDN); mpfr_sub(r20925, r20918, r20916, MPFR_RNDN); ; mpfr_mul(r20927, r20925, r20926, MPFR_RNDN); mpfr_div(r20928, r20924, r20927, MPFR_RNDN); return mpfr_get_d(r20928, MPFR_RNDN); } static mpfr_t r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937, r20938, r20939, r20940, r20941, r20942; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20929, "1", 10, MPFR_RNDN); mpfr_init(r20930); mpfr_init(r20931); mpfr_init(r20932); mpfr_init(r20933); mpfr_init(r20934); mpfr_init(r20935); mpfr_init(r20936); mpfr_init(r20937); mpfr_init(r20938); mpfr_init(r20939); mpfr_init_set_str(r20940, "2", 10, MPFR_RNDN); mpfr_init(r20941); mpfr_init(r20942); } double f_dm(double a, double b) { ; mpfr_set_d(r20930, a, MPFR_RNDN); mpfr_div(r20931, r20929, r20930, MPFR_RNDN); mpfr_set_d(r20932, b, MPFR_RNDN); mpfr_div(r20933, r20929, r20932, MPFR_RNDN); mpfr_sub(r20934, r20931, r20933, MPFR_RNDN); mpfr_add(r20935, r20932, r20930, MPFR_RNDN); mpfr_div(r20936, r20934, r20935, MPFR_RNDN); mpfr_const_pi(r20937, MPFR_RNDN); mpfr_mul(r20938, r20936, r20937, MPFR_RNDN); mpfr_sub(r20939, r20932, r20930, MPFR_RNDN); ; mpfr_mul(r20941, r20939, r20940, MPFR_RNDN); mpfr_div(r20942, r20938, r20941, MPFR_RNDN); return mpfr_get_d(r20942, MPFR_RNDN); }