#include #include #include #include #include char *name = "math.sqrt on complex, imaginary part, im greater than 0 branch"; double f_if(float re, float im) { float r21922 = 0.5; float r21923 = 2.0; float r21924 = re; float r21925 = r21924 * r21924; float r21926 = im; float r21927 = r21926 * r21926; float r21928 = r21925 - r21927; float r21929 = sqrt(r21928); float r21930 = r21929 + r21924; float r21931 = r21923 * r21930; float r21932 = sqrt(r21931); float r21933 = r21922 * r21932; return r21933; } double f_id(double re, double im) { double r21934 = 0.5; double r21935 = 2.0; double r21936 = re; double r21937 = r21936 * r21936; double r21938 = im; double r21939 = r21938 * r21938; double r21940 = r21937 - r21939; double r21941 = sqrt(r21940); double r21942 = r21941 + r21936; double r21943 = r21935 * r21942; double r21944 = sqrt(r21943); double r21945 = r21934 * r21944; return r21945; } double f_of(float re, float im) { float r21946 = 0.5; float r21947 = 2.0; float r21948 = sqrt(r21947); float r21949 = cbrt(r21948); float r21950 = r21949 * r21949; float r21951 = re; float r21952 = im; float r21953 = r21951 + r21952; float r21954 = sqrt(r21953); float r21955 = r21951 - r21952; float r21956 = sqrt(r21955); float r21957 = fma(r21954, r21956, r21951); float r21958 = sqrt(r21957); float r21959 = r21949 * r21958; float r21960 = r21950 * r21959; float r21961 = r21946 * r21960; return r21961; } double f_od(double re, double im) { double r21962 = 0.5; double r21963 = 2.0; double r21964 = sqrt(r21963); double r21965 = cbrt(r21964); double r21966 = r21965 * r21965; double r21967 = re; double r21968 = im; double r21969 = r21967 + r21968; double r21970 = sqrt(r21969); double r21971 = r21967 - r21968; double r21972 = sqrt(r21971); double r21973 = fma(r21970, r21972, r21967); double r21974 = sqrt(r21973); double r21975 = r21965 * r21974; double r21976 = r21966 * r21975; double r21977 = r21962 * r21976; return r21977; } 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 r21978, r21979, r21980, r21981, r21982, r21983, r21984, r21985, r21986, r21987, r21988, r21989; void setup_mpfr_f_im() { mpfr_set_default_prec(4496); mpfr_init_set_str(r21978, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r21979, "2.0", 10, MPFR_RNDN); 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); mpfr_init(r21989); } double f_im(double re, double im) { ; ; mpfr_set_d(r21980, re, MPFR_RNDN); mpfr_mul(r21981, r21980, r21980, MPFR_RNDN); mpfr_set_d(r21982, im, MPFR_RNDN); mpfr_mul(r21983, r21982, r21982, MPFR_RNDN); mpfr_sub(r21984, r21981, r21983, MPFR_RNDN); mpfr_sqrt(r21985, r21984, MPFR_RNDN); mpfr_add(r21986, r21985, r21980, MPFR_RNDN); mpfr_mul(r21987, r21979, r21986, MPFR_RNDN); mpfr_sqrt(r21988, r21987, MPFR_RNDN); mpfr_mul(r21989, r21978, r21988, MPFR_RNDN); return mpfr_get_d(r21989, MPFR_RNDN); } static mpfr_t r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000, r22001, r22002, r22003, r22004, r22005; void setup_mpfr_f_fm() { mpfr_set_default_prec(4496); mpfr_init_set_str(r21990, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r21991, "2.0", 10, MPFR_RNDN); mpfr_init(r21992); mpfr_init(r21993); mpfr_init(r21994); mpfr_init(r21995); mpfr_init(r21996); mpfr_init(r21997); mpfr_init(r21998); mpfr_init(r21999); mpfr_init(r22000); mpfr_init(r22001); mpfr_init(r22002); mpfr_init(r22003); mpfr_init(r22004); mpfr_init(r22005); } double f_fm(double re, double im) { ; ; mpfr_sqrt(r21992, r21991, MPFR_RNDN); mpfr_cbrt(r21993, r21992, MPFR_RNDN); mpfr_mul(r21994, r21993, r21993, MPFR_RNDN); mpfr_set_d(r21995, re, MPFR_RNDN); mpfr_set_d(r21996, im, MPFR_RNDN); mpfr_add(r21997, r21995, r21996, MPFR_RNDN); mpfr_sqrt(r21998, r21997, MPFR_RNDN); mpfr_sub(r21999, r21995, r21996, MPFR_RNDN); mpfr_sqrt(r22000, r21999, MPFR_RNDN); mpfr_fma(r22001, r21998, r22000, r21995, MPFR_RNDN); mpfr_sqrt(r22002, r22001, MPFR_RNDN); mpfr_mul(r22003, r21993, r22002, MPFR_RNDN); mpfr_mul(r22004, r21994, r22003, MPFR_RNDN); mpfr_mul(r22005, r21990, r22004, MPFR_RNDN); return mpfr_get_d(r22005, MPFR_RNDN); } static mpfr_t r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021; void setup_mpfr_f_dm() { mpfr_set_default_prec(4496); mpfr_init_set_str(r22006, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r22007, "2.0", 10, MPFR_RNDN); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init(r22012); mpfr_init(r22013); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init(r22018); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); } double f_dm(double re, double im) { ; ; mpfr_sqrt(r22008, r22007, MPFR_RNDN); mpfr_cbrt(r22009, r22008, MPFR_RNDN); mpfr_mul(r22010, r22009, r22009, MPFR_RNDN); mpfr_set_d(r22011, re, MPFR_RNDN); mpfr_set_d(r22012, im, MPFR_RNDN); mpfr_add(r22013, r22011, r22012, MPFR_RNDN); mpfr_sqrt(r22014, r22013, MPFR_RNDN); mpfr_sub(r22015, r22011, r22012, MPFR_RNDN); mpfr_sqrt(r22016, r22015, MPFR_RNDN); mpfr_fma(r22017, r22014, r22016, r22011, MPFR_RNDN); mpfr_sqrt(r22018, r22017, MPFR_RNDN); mpfr_mul(r22019, r22009, r22018, MPFR_RNDN); mpfr_mul(r22020, r22010, r22019, MPFR_RNDN); mpfr_mul(r22021, r22006, r22020, MPFR_RNDN); return mpfr_get_d(r22021, MPFR_RNDN); }