#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r21916 = atan2(1.0, 0.0); float r21917 = 2; float r21918 = r21916 / r21917; float r21919 = 1; float r21920 = b; float r21921 = r21920 * r21920; float r21922 = a; float r21923 = r21922 * r21922; float r21924 = r21921 - r21923; float r21925 = r21919 / r21924; float r21926 = r21918 * r21925; float r21927 = r21919 / r21922; float r21928 = r21919 / r21920; float r21929 = r21927 - r21928; float r21930 = r21926 * r21929; return r21930; } double f_id(double a, double b) { double r21931 = atan2(1.0, 0.0); double r21932 = 2; double r21933 = r21931 / r21932; double r21934 = 1; double r21935 = b; double r21936 = r21935 * r21935; double r21937 = a; double r21938 = r21937 * r21937; double r21939 = r21936 - r21938; double r21940 = r21934 / r21939; double r21941 = r21933 * r21940; double r21942 = r21934 / r21937; double r21943 = r21934 / r21935; double r21944 = r21942 - r21943; double r21945 = r21941 * r21944; return r21945; } double f_of(float a, float b) { float r21946 = 1; float r21947 = b; float r21948 = a; float r21949 = r21947 + r21948; float r21950 = r21946 / r21949; float r21951 = atan2(1.0, 0.0); float r21952 = r21951 / r21948; float r21953 = r21951 / r21947; float r21954 = r21952 - r21953; float r21955 = r21947 - r21948; float r21956 = 2; float r21957 = r21955 * r21956; float r21958 = r21954 / r21957; float r21959 = r21950 * r21958; return r21959; } double f_od(double a, double b) { double r21960 = 1; double r21961 = b; double r21962 = a; double r21963 = r21961 + r21962; double r21964 = r21960 / r21963; double r21965 = atan2(1.0, 0.0); double r21966 = r21965 / r21962; double r21967 = r21965 / r21961; double r21968 = r21966 - r21967; double r21969 = r21961 - r21962; double r21970 = 2; double r21971 = r21969 * r21970; double r21972 = r21968 / r21971; double r21973 = r21964 * r21972; return r21973; } 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 r21974, r21975, r21976, r21977, r21978, r21979, r21980, r21981, r21982, r21983, r21984, r21985, r21986, r21987, r21988; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21974); mpfr_init_set_str(r21975, "2", 10, MPFR_RNDN); mpfr_init(r21976); mpfr_init_set_str(r21977, "1", 10, MPFR_RNDN); mpfr_init(r21978); mpfr_init(r21979); mpfr_init(r21980); mpfr_init(r21981); mpfr_init(r21982); mpfr_init(r21983); mpfr_init(r21984); mpfr_init(r21985); mpfr_init(r21986); mpfr_init(r21987); mpfr_init(r21988); } double f_im(double a, double b) { mpfr_const_pi(r21974, MPFR_RNDN); ; mpfr_div(r21976, r21974, r21975, MPFR_RNDN); ; mpfr_set_d(r21978, b, MPFR_RNDN); mpfr_mul(r21979, r21978, r21978, MPFR_RNDN); mpfr_set_d(r21980, a, MPFR_RNDN); mpfr_mul(r21981, r21980, r21980, MPFR_RNDN); mpfr_sub(r21982, r21979, r21981, MPFR_RNDN); mpfr_div(r21983, r21977, r21982, MPFR_RNDN); mpfr_mul(r21984, r21976, r21983, MPFR_RNDN); mpfr_div(r21985, r21977, r21980, MPFR_RNDN); mpfr_div(r21986, r21977, r21978, MPFR_RNDN); mpfr_sub(r21987, r21985, r21986, MPFR_RNDN); mpfr_mul(r21988, r21984, r21987, MPFR_RNDN); return mpfr_get_d(r21988, MPFR_RNDN); } static mpfr_t r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000, r22001, r22002; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21989, "1", 10, MPFR_RNDN); mpfr_init(r21990); mpfr_init(r21991); mpfr_init(r21992); mpfr_init(r21993); mpfr_init(r21994); mpfr_init(r21995); mpfr_init(r21996); mpfr_init(r21997); mpfr_init(r21998); mpfr_init_set_str(r21999, "2", 10, MPFR_RNDN); mpfr_init(r22000); mpfr_init(r22001); mpfr_init(r22002); } double f_fm(double a, double b) { ; mpfr_set_d(r21990, b, MPFR_RNDN); mpfr_set_d(r21991, a, MPFR_RNDN); mpfr_add(r21992, r21990, r21991, MPFR_RNDN); mpfr_div(r21993, r21989, r21992, MPFR_RNDN); mpfr_const_pi(r21994, MPFR_RNDN); mpfr_div(r21995, r21994, r21991, MPFR_RNDN); mpfr_div(r21996, r21994, r21990, MPFR_RNDN); mpfr_sub(r21997, r21995, r21996, MPFR_RNDN); mpfr_sub(r21998, r21990, r21991, MPFR_RNDN); ; mpfr_mul(r22000, r21998, r21999, MPFR_RNDN); mpfr_div(r22001, r21997, r22000, MPFR_RNDN); mpfr_mul(r22002, r21993, r22001, MPFR_RNDN); return mpfr_get_d(r22002, MPFR_RNDN); } static mpfr_t r22003, r22004, r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22003, "1", 10, MPFR_RNDN); mpfr_init(r22004); mpfr_init(r22005); mpfr_init(r22006); mpfr_init(r22007); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init(r22012); mpfr_init_set_str(r22013, "2", 10, MPFR_RNDN); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); } double f_dm(double a, double b) { ; mpfr_set_d(r22004, b, MPFR_RNDN); mpfr_set_d(r22005, a, MPFR_RNDN); mpfr_add(r22006, r22004, r22005, MPFR_RNDN); mpfr_div(r22007, r22003, r22006, MPFR_RNDN); mpfr_const_pi(r22008, MPFR_RNDN); mpfr_div(r22009, r22008, r22005, MPFR_RNDN); mpfr_div(r22010, r22008, r22004, MPFR_RNDN); mpfr_sub(r22011, r22009, r22010, MPFR_RNDN); mpfr_sub(r22012, r22004, r22005, MPFR_RNDN); ; mpfr_mul(r22014, r22012, r22013, MPFR_RNDN); mpfr_div(r22015, r22011, r22014, MPFR_RNDN); mpfr_mul(r22016, r22007, r22015, MPFR_RNDN); return mpfr_get_d(r22016, MPFR_RNDN); }