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