#include #include #include #include #include char *name = "Complex division, real part"; double f_if(float a, float b, float c, float d) { float r23868 = a; float r23869 = c; float r23870 = r23868 * r23869; float r23871 = b; float r23872 = d; float r23873 = r23871 * r23872; float r23874 = r23870 + r23873; float r23875 = r23869 * r23869; float r23876 = r23872 * r23872; float r23877 = r23875 + r23876; float r23878 = r23874 / r23877; return r23878; } double f_id(double a, double b, double c, double d) { double r23879 = a; double r23880 = c; double r23881 = r23879 * r23880; double r23882 = b; double r23883 = d; double r23884 = r23882 * r23883; double r23885 = r23881 + r23884; double r23886 = r23880 * r23880; double r23887 = r23883 * r23883; double r23888 = r23886 + r23887; double r23889 = r23885 / r23888; return r23889; } double f_of(float a, float b, float c, float d) { float r23890 = c; float r23891 = -3.508226146289393e+102; bool r23892 = r23890 <= r23891; float r23893 = a; float r23894 = -r23893; float r23895 = r23890 * r23890; float r23896 = d; float r23897 = r23896 * r23896; float r23898 = r23895 + r23897; float r23899 = sqrt(r23898); float r23900 = r23894 / r23899; float r23901 = r23893 * r23890; float r23902 = b; float r23903 = r23902 * r23896; float r23904 = r23901 + r23903; float r23905 = r23904 / r23899; float r23906 = r23905 / r23899; float r23907 = r23892 ? r23900 : r23906; return r23907; } double f_od(double a, double b, double c, double d) { double r23908 = c; double r23909 = -3.508226146289393e+102; bool r23910 = r23908 <= r23909; double r23911 = a; double r23912 = -r23911; double r23913 = r23908 * r23908; double r23914 = d; double r23915 = r23914 * r23914; double r23916 = r23913 + r23915; double r23917 = sqrt(r23916); double r23918 = r23912 / r23917; double r23919 = r23911 * r23908; double r23920 = b; double r23921 = r23920 * r23914; double r23922 = r23919 + r23921; double r23923 = r23922 / r23917; double r23924 = r23923 / r23917; double r23925 = r23910 ? r23918 : r23924; return r23925; } 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 r23926, r23927, r23928, r23929, r23930, r23931, r23932, r23933, r23934, r23935, r23936; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23926); mpfr_init(r23927); mpfr_init(r23928); mpfr_init(r23929); mpfr_init(r23930); mpfr_init(r23931); mpfr_init(r23932); mpfr_init(r23933); mpfr_init(r23934); mpfr_init(r23935); mpfr_init(r23936); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r23926, a, MPFR_RNDN); mpfr_set_d(r23927, c, MPFR_RNDN); mpfr_mul(r23928, r23926, r23927, MPFR_RNDN); mpfr_set_d(r23929, b, MPFR_RNDN); mpfr_set_d(r23930, d, MPFR_RNDN); mpfr_mul(r23931, r23929, r23930, MPFR_RNDN); mpfr_add(r23932, r23928, r23931, MPFR_RNDN); mpfr_mul(r23933, r23927, r23927, MPFR_RNDN); mpfr_mul(r23934, r23930, r23930, MPFR_RNDN); mpfr_add(r23935, r23933, r23934, MPFR_RNDN); mpfr_div(r23936, r23932, r23935, MPFR_RNDN); return mpfr_get_d(r23936, MPFR_RNDN); } static mpfr_t r23937, r23938, r23939, r23940, r23941, r23942, r23943, r23944, r23945, r23946, r23947, r23948, r23949, r23950, r23951, r23952, r23953, r23954; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23937); mpfr_init_set_str(r23938, "-3.508226146289393e+102", 10, MPFR_RNDN); mpfr_init(r23939); mpfr_init(r23940); mpfr_init(r23941); mpfr_init(r23942); mpfr_init(r23943); mpfr_init(r23944); mpfr_init(r23945); mpfr_init(r23946); mpfr_init(r23947); mpfr_init(r23948); mpfr_init(r23949); mpfr_init(r23950); mpfr_init(r23951); mpfr_init(r23952); mpfr_init(r23953); mpfr_init(r23954); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r23937, c, MPFR_RNDN); ; mpfr_set_si(r23939, mpfr_cmp(r23937, r23938) <= 0, MPFR_RNDN); mpfr_set_d(r23940, a, MPFR_RNDN); mpfr_neg(r23941, r23940, MPFR_RNDN); mpfr_mul(r23942, r23937, r23937, MPFR_RNDN); mpfr_set_d(r23943, d, MPFR_RNDN); mpfr_mul(r23944, r23943, r23943, MPFR_RNDN); mpfr_add(r23945, r23942, r23944, MPFR_RNDN); mpfr_sqrt(r23946, r23945, MPFR_RNDN); mpfr_div(r23947, r23941, r23946, MPFR_RNDN); mpfr_mul(r23948, r23940, r23937, MPFR_RNDN); mpfr_set_d(r23949, b, MPFR_RNDN); mpfr_mul(r23950, r23949, r23943, MPFR_RNDN); mpfr_add(r23951, r23948, r23950, MPFR_RNDN); mpfr_div(r23952, r23951, r23946, MPFR_RNDN); mpfr_div(r23953, r23952, r23946, MPFR_RNDN); if (mpfr_get_si(r23939, MPFR_RNDN)) { mpfr_set(r23954, r23947, MPFR_RNDN); } else { mpfr_set(r23954, r23953, MPFR_RNDN); }; return mpfr_get_d(r23954, MPFR_RNDN); } static mpfr_t r23955, r23956, r23957, r23958, r23959, r23960, r23961, r23962, r23963, r23964, r23965, r23966, r23967, r23968, r23969, r23970, r23971, r23972; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23955); mpfr_init_set_str(r23956, "-3.508226146289393e+102", 10, MPFR_RNDN); mpfr_init(r23957); mpfr_init(r23958); mpfr_init(r23959); mpfr_init(r23960); mpfr_init(r23961); mpfr_init(r23962); mpfr_init(r23963); mpfr_init(r23964); mpfr_init(r23965); mpfr_init(r23966); mpfr_init(r23967); mpfr_init(r23968); mpfr_init(r23969); mpfr_init(r23970); mpfr_init(r23971); mpfr_init(r23972); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r23955, c, MPFR_RNDN); ; mpfr_set_si(r23957, mpfr_cmp(r23955, r23956) <= 0, MPFR_RNDN); mpfr_set_d(r23958, a, MPFR_RNDN); mpfr_neg(r23959, r23958, MPFR_RNDN); mpfr_mul(r23960, r23955, r23955, MPFR_RNDN); mpfr_set_d(r23961, d, MPFR_RNDN); mpfr_mul(r23962, r23961, r23961, MPFR_RNDN); mpfr_add(r23963, r23960, r23962, MPFR_RNDN); mpfr_sqrt(r23964, r23963, MPFR_RNDN); mpfr_div(r23965, r23959, r23964, MPFR_RNDN); mpfr_mul(r23966, r23958, r23955, MPFR_RNDN); mpfr_set_d(r23967, b, MPFR_RNDN); mpfr_mul(r23968, r23967, r23961, MPFR_RNDN); mpfr_add(r23969, r23966, r23968, MPFR_RNDN); mpfr_div(r23970, r23969, r23964, MPFR_RNDN); mpfr_div(r23971, r23970, r23964, MPFR_RNDN); if (mpfr_get_si(r23957, MPFR_RNDN)) { mpfr_set(r23972, r23965, MPFR_RNDN); } else { mpfr_set(r23972, r23971, MPFR_RNDN); }; return mpfr_get_d(r23972, MPFR_RNDN); }