#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 r22856 = x_re; float r22857 = y_re; float r22858 = r22856 * r22857; float r22859 = x_im; float r22860 = y_im; float r22861 = r22859 * r22860; float r22862 = r22858 + r22861; float r22863 = r22857 * r22857; float r22864 = r22860 * r22860; float r22865 = r22863 + r22864; float r22866 = r22862 / r22865; return r22866; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r22867 = x_re; double r22868 = y_re; double r22869 = r22867 * r22868; double r22870 = x_im; double r22871 = y_im; double r22872 = r22870 * r22871; double r22873 = r22869 + r22872; double r22874 = r22868 * r22868; double r22875 = r22871 * r22871; double r22876 = r22874 + r22875; double r22877 = r22873 / r22876; return r22877; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r22878 = y_re; float r22879 = -1.609443697480765e+108; bool r22880 = r22878 <= r22879; float r22881 = y_im; float r22882 = r22881 / r22878; float r22883 = x_im; float r22884 = -r22883; float r22885 = x_re; float r22886 = -r22885; float r22887 = fma(r22882, r22884, r22886); float r22888 = hypot(r22878, r22881); float r22889 = r22887 / r22888; float r22890 = 714.2898411289591; bool r22891 = r22878 <= r22890; float r22892 = 1; float r22893 = r22892 / r22888; float r22894 = r22878 * r22885; float r22895 = fma(r22883, r22881, r22894); float r22896 = r22895 / r22888; float r22897 = r22893 * r22896; float r22898 = fma(r22882, r22883, r22885); float r22899 = r22898 / r22888; float r22900 = r22891 ? r22897 : r22899; float r22901 = r22880 ? r22889 : r22900; return r22901; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r22902 = y_re; double r22903 = -1.609443697480765e+108; bool r22904 = r22902 <= r22903; double r22905 = y_im; double r22906 = r22905 / r22902; double r22907 = x_im; double r22908 = -r22907; double r22909 = x_re; double r22910 = -r22909; double r22911 = fma(r22906, r22908, r22910); double r22912 = hypot(r22902, r22905); double r22913 = r22911 / r22912; double r22914 = 714.2898411289591; bool r22915 = r22902 <= r22914; double r22916 = 1; double r22917 = r22916 / r22912; double r22918 = r22902 * r22909; double r22919 = fma(r22907, r22905, r22918); double r22920 = r22919 / r22912; double r22921 = r22917 * r22920; double r22922 = fma(r22906, r22907, r22909); double r22923 = r22922 / r22912; double r22924 = r22915 ? r22921 : r22923; double r22925 = r22904 ? r22913 : r22924; return r22925; } 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 r22926, r22927, r22928, r22929, r22930, r22931, r22932, r22933, r22934, r22935, r22936; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22926); mpfr_init(r22927); mpfr_init(r22928); mpfr_init(r22929); mpfr_init(r22930); mpfr_init(r22931); mpfr_init(r22932); mpfr_init(r22933); mpfr_init(r22934); mpfr_init(r22935); mpfr_init(r22936); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22926, x_re, MPFR_RNDN); mpfr_set_d(r22927, y_re, MPFR_RNDN); mpfr_mul(r22928, r22926, r22927, MPFR_RNDN); mpfr_set_d(r22929, x_im, MPFR_RNDN); mpfr_set_d(r22930, y_im, MPFR_RNDN); mpfr_mul(r22931, r22929, r22930, MPFR_RNDN); mpfr_add(r22932, r22928, r22931, MPFR_RNDN); mpfr_mul(r22933, r22927, r22927, MPFR_RNDN); mpfr_mul(r22934, r22930, r22930, MPFR_RNDN); mpfr_add(r22935, r22933, r22934, MPFR_RNDN); mpfr_div(r22936, r22932, r22935, MPFR_RNDN); return mpfr_get_d(r22936, MPFR_RNDN); } static mpfr_t r22937, r22938, r22939, r22940, r22941, r22942, r22943, r22944, r22945, r22946, r22947, r22948, r22949, r22950, r22951, r22952, r22953, r22954, r22955, r22956, r22957, r22958, r22959, r22960; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22937); mpfr_init_set_str(r22938, "-1.609443697480765e+108", 10, MPFR_RNDN); mpfr_init(r22939); mpfr_init(r22940); mpfr_init(r22941); mpfr_init(r22942); mpfr_init(r22943); mpfr_init(r22944); mpfr_init(r22945); mpfr_init(r22946); mpfr_init(r22947); mpfr_init(r22948); mpfr_init_set_str(r22949, "714.2898411289591", 10, MPFR_RNDN); mpfr_init(r22950); mpfr_init_set_str(r22951, "1", 10, MPFR_RNDN); mpfr_init(r22952); mpfr_init(r22953); mpfr_init(r22954); mpfr_init(r22955); mpfr_init(r22956); mpfr_init(r22957); mpfr_init(r22958); mpfr_init(r22959); mpfr_init(r22960); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22937, y_re, MPFR_RNDN); ; mpfr_set_si(r22939, mpfr_cmp(r22937, r22938) <= 0, MPFR_RNDN); mpfr_set_d(r22940, y_im, MPFR_RNDN); mpfr_div(r22941, r22940, r22937, MPFR_RNDN); mpfr_set_d(r22942, x_im, MPFR_RNDN); mpfr_neg(r22943, r22942, MPFR_RNDN); mpfr_set_d(r22944, x_re, MPFR_RNDN); mpfr_neg(r22945, r22944, MPFR_RNDN); mpfr_fma(r22946, r22941, r22943, r22945, MPFR_RNDN); mpfr_hypot(r22947, r22937, r22940, MPFR_RNDN); mpfr_div(r22948, r22946, r22947, MPFR_RNDN); ; mpfr_set_si(r22950, mpfr_cmp(r22937, r22949) <= 0, MPFR_RNDN); ; mpfr_div(r22952, r22951, r22947, MPFR_RNDN); mpfr_mul(r22953, r22937, r22944, MPFR_RNDN); mpfr_fma(r22954, r22942, r22940, r22953, MPFR_RNDN); mpfr_div(r22955, r22954, r22947, MPFR_RNDN); mpfr_mul(r22956, r22952, r22955, MPFR_RNDN); mpfr_fma(r22957, r22941, r22942, r22944, MPFR_RNDN); mpfr_div(r22958, r22957, r22947, MPFR_RNDN); if (mpfr_get_si(r22950, MPFR_RNDN)) { mpfr_set(r22959, r22956, MPFR_RNDN); } else { mpfr_set(r22959, r22958, MPFR_RNDN); }; if (mpfr_get_si(r22939, MPFR_RNDN)) { mpfr_set(r22960, r22948, MPFR_RNDN); } else { mpfr_set(r22960, r22959, MPFR_RNDN); }; return mpfr_get_d(r22960, MPFR_RNDN); } static mpfr_t r22961, r22962, r22963, r22964, r22965, r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980, r22981, r22982, r22983, r22984; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22961); mpfr_init_set_str(r22962, "-1.609443697480765e+108", 10, MPFR_RNDN); mpfr_init(r22963); mpfr_init(r22964); mpfr_init(r22965); mpfr_init(r22966); mpfr_init(r22967); mpfr_init(r22968); mpfr_init(r22969); mpfr_init(r22970); mpfr_init(r22971); mpfr_init(r22972); mpfr_init_set_str(r22973, "714.2898411289591", 10, MPFR_RNDN); mpfr_init(r22974); mpfr_init_set_str(r22975, "1", 10, MPFR_RNDN); mpfr_init(r22976); mpfr_init(r22977); mpfr_init(r22978); mpfr_init(r22979); mpfr_init(r22980); mpfr_init(r22981); mpfr_init(r22982); mpfr_init(r22983); mpfr_init(r22984); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22961, y_re, MPFR_RNDN); ; mpfr_set_si(r22963, mpfr_cmp(r22961, r22962) <= 0, MPFR_RNDN); mpfr_set_d(r22964, y_im, MPFR_RNDN); mpfr_div(r22965, r22964, r22961, MPFR_RNDN); mpfr_set_d(r22966, x_im, MPFR_RNDN); mpfr_neg(r22967, r22966, MPFR_RNDN); mpfr_set_d(r22968, x_re, MPFR_RNDN); mpfr_neg(r22969, r22968, MPFR_RNDN); mpfr_fma(r22970, r22965, r22967, r22969, MPFR_RNDN); mpfr_hypot(r22971, r22961, r22964, MPFR_RNDN); mpfr_div(r22972, r22970, r22971, MPFR_RNDN); ; mpfr_set_si(r22974, mpfr_cmp(r22961, r22973) <= 0, MPFR_RNDN); ; mpfr_div(r22976, r22975, r22971, MPFR_RNDN); mpfr_mul(r22977, r22961, r22968, MPFR_RNDN); mpfr_fma(r22978, r22966, r22964, r22977, MPFR_RNDN); mpfr_div(r22979, r22978, r22971, MPFR_RNDN); mpfr_mul(r22980, r22976, r22979, MPFR_RNDN); mpfr_fma(r22981, r22965, r22966, r22968, MPFR_RNDN); mpfr_div(r22982, r22981, r22971, MPFR_RNDN); if (mpfr_get_si(r22974, MPFR_RNDN)) { mpfr_set(r22983, r22980, MPFR_RNDN); } else { mpfr_set(r22983, r22982, MPFR_RNDN); }; if (mpfr_get_si(r22963, MPFR_RNDN)) { mpfr_set(r22984, r22972, MPFR_RNDN); } else { mpfr_set(r22984, r22983, MPFR_RNDN); }; return mpfr_get_d(r22984, MPFR_RNDN); }