#include #include #include #include #include char *name = "_divideComplex, real part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r22903 = x_re; float r22904 = y_re; float r22905 = r22903 * r22904; float r22906 = x_im; float r22907 = y_im; float r22908 = r22906 * r22907; float r22909 = r22905 + r22908; float r22910 = r22904 * r22904; float r22911 = r22907 * r22907; float r22912 = r22910 + r22911; float r22913 = r22909 / r22912; return r22913; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r22914 = x_re; double r22915 = y_re; double r22916 = r22914 * r22915; double r22917 = x_im; double r22918 = y_im; double r22919 = r22917 * r22918; double r22920 = r22916 + r22919; double r22921 = r22915 * r22915; double r22922 = r22918 * r22918; double r22923 = r22921 + r22922; double r22924 = r22920 / r22923; return r22924; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r22925 = y_re; float r22926 = 3.637924346615467e+157; bool r22927 = r22925 <= r22926; float r22928 = 1; float r22929 = y_im; float r22930 = hypot(r22929, r22925); float r22931 = r22928 / r22930; float r22932 = x_im; float r22933 = x_re; float r22934 = r22925 * r22933; float r22935 = fma(r22932, r22929, r22934); float r22936 = r22935 / r22930; float r22937 = r22931 * r22936; float r22938 = r22933 / r22930; float r22939 = r22927 ? r22937 : r22938; return r22939; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r22940 = y_re; double r22941 = 3.637924346615467e+157; bool r22942 = r22940 <= r22941; double r22943 = 1; double r22944 = y_im; double r22945 = hypot(r22944, r22940); double r22946 = r22943 / r22945; double r22947 = x_im; double r22948 = x_re; double r22949 = r22940 * r22948; double r22950 = fma(r22947, r22944, r22949); double r22951 = r22950 / r22945; double r22952 = r22946 * r22951; double r22953 = r22948 / r22945; double r22954 = r22942 ? r22952 : r22953; return r22954; } 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 r22955, r22956, r22957, r22958, r22959, r22960, r22961, r22962, r22963, r22964, r22965; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22955); mpfr_init(r22956); mpfr_init(r22957); mpfr_init(r22958); mpfr_init(r22959); mpfr_init(r22960); mpfr_init(r22961); mpfr_init(r22962); mpfr_init(r22963); mpfr_init(r22964); mpfr_init(r22965); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22955, x_re, MPFR_RNDN); mpfr_set_d(r22956, y_re, MPFR_RNDN); mpfr_mul(r22957, r22955, r22956, MPFR_RNDN); mpfr_set_d(r22958, x_im, MPFR_RNDN); mpfr_set_d(r22959, y_im, MPFR_RNDN); mpfr_mul(r22960, r22958, r22959, MPFR_RNDN); mpfr_add(r22961, r22957, r22960, MPFR_RNDN); mpfr_mul(r22962, r22956, r22956, MPFR_RNDN); mpfr_mul(r22963, r22959, r22959, MPFR_RNDN); mpfr_add(r22964, r22962, r22963, MPFR_RNDN); mpfr_div(r22965, r22961, r22964, MPFR_RNDN); return mpfr_get_d(r22965, MPFR_RNDN); } static mpfr_t r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r22966); mpfr_init_set_str(r22967, "3.637924346615467e+157", 10, MPFR_RNDN); mpfr_init(r22968); mpfr_init_set_str(r22969, "1", 10, MPFR_RNDN); mpfr_init(r22970); mpfr_init(r22971); mpfr_init(r22972); mpfr_init(r22973); mpfr_init(r22974); mpfr_init(r22975); mpfr_init(r22976); mpfr_init(r22977); mpfr_init(r22978); mpfr_init(r22979); mpfr_init(r22980); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22966, y_re, MPFR_RNDN); ; mpfr_set_si(r22968, mpfr_cmp(r22966, r22967) <= 0, MPFR_RNDN); ; mpfr_set_d(r22970, y_im, MPFR_RNDN); mpfr_hypot(r22971, r22970, r22966, MPFR_RNDN); mpfr_div(r22972, r22969, r22971, MPFR_RNDN); mpfr_set_d(r22973, x_im, MPFR_RNDN); mpfr_set_d(r22974, x_re, MPFR_RNDN); mpfr_mul(r22975, r22966, r22974, MPFR_RNDN); mpfr_fma(r22976, r22973, r22970, r22975, MPFR_RNDN); mpfr_div(r22977, r22976, r22971, MPFR_RNDN); mpfr_mul(r22978, r22972, r22977, MPFR_RNDN); mpfr_div(r22979, r22974, r22971, MPFR_RNDN); if (mpfr_get_si(r22968, MPFR_RNDN)) { mpfr_set(r22980, r22978, MPFR_RNDN); } else { mpfr_set(r22980, r22979, MPFR_RNDN); }; return mpfr_get_d(r22980, MPFR_RNDN); } static mpfr_t r22981, r22982, r22983, r22984, r22985, r22986, r22987, r22988, r22989, r22990, r22991, r22992, r22993, r22994, r22995; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22981); mpfr_init_set_str(r22982, "3.637924346615467e+157", 10, MPFR_RNDN); mpfr_init(r22983); mpfr_init_set_str(r22984, "1", 10, MPFR_RNDN); mpfr_init(r22985); mpfr_init(r22986); mpfr_init(r22987); mpfr_init(r22988); mpfr_init(r22989); mpfr_init(r22990); mpfr_init(r22991); mpfr_init(r22992); mpfr_init(r22993); mpfr_init(r22994); mpfr_init(r22995); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22981, y_re, MPFR_RNDN); ; mpfr_set_si(r22983, mpfr_cmp(r22981, r22982) <= 0, MPFR_RNDN); ; mpfr_set_d(r22985, y_im, MPFR_RNDN); mpfr_hypot(r22986, r22985, r22981, MPFR_RNDN); mpfr_div(r22987, r22984, r22986, MPFR_RNDN); mpfr_set_d(r22988, x_im, MPFR_RNDN); mpfr_set_d(r22989, x_re, MPFR_RNDN); mpfr_mul(r22990, r22981, r22989, MPFR_RNDN); mpfr_fma(r22991, r22988, r22985, r22990, MPFR_RNDN); mpfr_div(r22992, r22991, r22986, MPFR_RNDN); mpfr_mul(r22993, r22987, r22992, MPFR_RNDN); mpfr_div(r22994, r22989, r22986, MPFR_RNDN); if (mpfr_get_si(r22983, MPFR_RNDN)) { mpfr_set(r22995, r22993, MPFR_RNDN); } else { mpfr_set(r22995, r22994, MPFR_RNDN); }; return mpfr_get_d(r22995, MPFR_RNDN); }