#include #include #include #include #include char *name = "expq3 (problem 3.4.2)"; double f_if(float a, float b, float eps) { float r23914 = eps; float r23915 = a; float r23916 = b; float r23917 = r23915 + r23916; float r23918 = r23917 * r23914; float r23919 = exp(r23918); float r23920 = 1; float r23921 = r23919 - r23920; float r23922 = r23914 * r23921; float r23923 = r23915 * r23914; float r23924 = exp(r23923); float r23925 = r23924 - r23920; float r23926 = r23916 * r23914; float r23927 = exp(r23926); float r23928 = r23927 - r23920; float r23929 = r23925 * r23928; float r23930 = r23922 / r23929; return r23930; } double f_id(double a, double b, double eps) { double r23931 = eps; double r23932 = a; double r23933 = b; double r23934 = r23932 + r23933; double r23935 = r23934 * r23931; double r23936 = exp(r23935); double r23937 = 1; double r23938 = r23936 - r23937; double r23939 = r23931 * r23938; double r23940 = r23932 * r23931; double r23941 = exp(r23940); double r23942 = r23941 - r23937; double r23943 = r23933 * r23931; double r23944 = exp(r23943); double r23945 = r23944 - r23937; double r23946 = r23942 * r23945; double r23947 = r23939 / r23946; return r23947; } double f_of(float a, float b, float eps) { float r23948 = eps; float r23949 = b; float r23950 = r23949 * r23948; float r23951 = expm1(r23950); float r23952 = a; float r23953 = r23948 * r23952; float r23954 = expm1(r23953); float r23955 = r23951 * r23954; float r23956 = r23949 + r23952; float r23957 = r23956 * r23948; float r23958 = expm1(r23957); float r23959 = r23955 / r23958; float r23960 = r23948 / r23959; float r23961 = -5.270434816906832e-10; bool r23962 = r23960 <= r23961; float r23963 = 1; float r23964 = r23963 / r23949; float r23965 = r23963 / r23952; float r23966 = r23964 + r23965; float r23967 = 2.088519104255703e-20; bool r23968 = r23960 <= r23967; float r23969 = r23968 ? r23960 : r23966; float r23970 = r23962 ? r23966 : r23969; return r23970; } double f_od(double a, double b, double eps) { double r23971 = eps; double r23972 = b; double r23973 = r23972 * r23971; double r23974 = expm1(r23973); double r23975 = a; double r23976 = r23971 * r23975; double r23977 = expm1(r23976); double r23978 = r23974 * r23977; double r23979 = r23972 + r23975; double r23980 = r23979 * r23971; double r23981 = expm1(r23980); double r23982 = r23978 / r23981; double r23983 = r23971 / r23982; double r23984 = -5.270434816906832e-10; bool r23985 = r23983 <= r23984; double r23986 = 1; double r23987 = r23986 / r23972; double r23988 = r23986 / r23975; double r23989 = r23987 + r23988; double r23990 = 2.088519104255703e-20; bool r23991 = r23983 <= r23990; double r23992 = r23991 ? r23983 : r23989; double r23993 = r23985 ? r23989 : r23992; return r23993; } 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 r23994, r23995, r23996, r23997, r23998, r23999, r24000, r24001, r24002, r24003, r24004, r24005, r24006, r24007, r24008, r24009, r24010; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r23994); mpfr_init(r23995); mpfr_init(r23996); mpfr_init(r23997); mpfr_init(r23998); mpfr_init(r23999); mpfr_init_set_str(r24000, "1", 10, MPFR_RNDN); mpfr_init(r24001); mpfr_init(r24002); mpfr_init(r24003); mpfr_init(r24004); mpfr_init(r24005); mpfr_init(r24006); mpfr_init(r24007); mpfr_init(r24008); mpfr_init(r24009); mpfr_init(r24010); } double f_im(double a, double b, double eps) { mpfr_set_d(r23994, eps, MPFR_RNDN); mpfr_set_d(r23995, a, MPFR_RNDN); mpfr_set_d(r23996, b, MPFR_RNDN); mpfr_add(r23997, r23995, r23996, MPFR_RNDN); mpfr_mul(r23998, r23997, r23994, MPFR_RNDN); mpfr_exp(r23999, r23998, MPFR_RNDN); ; mpfr_sub(r24001, r23999, r24000, MPFR_RNDN); mpfr_mul(r24002, r23994, r24001, MPFR_RNDN); mpfr_mul(r24003, r23995, r23994, MPFR_RNDN); mpfr_exp(r24004, r24003, MPFR_RNDN); mpfr_sub(r24005, r24004, r24000, MPFR_RNDN); mpfr_mul(r24006, r23996, r23994, MPFR_RNDN); mpfr_exp(r24007, r24006, MPFR_RNDN); mpfr_sub(r24008, r24007, r24000, MPFR_RNDN); mpfr_mul(r24009, r24005, r24008, MPFR_RNDN); mpfr_div(r24010, r24002, r24009, MPFR_RNDN); return mpfr_get_d(r24010, MPFR_RNDN); } static mpfr_t r24011, r24012, r24013, r24014, r24015, r24016, r24017, r24018, r24019, r24020, r24021, r24022, r24023, r24024, r24025, r24026, r24027, r24028, r24029, r24030, r24031, r24032, r24033; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r24011); mpfr_init(r24012); mpfr_init(r24013); mpfr_init(r24014); mpfr_init(r24015); mpfr_init(r24016); mpfr_init(r24017); mpfr_init(r24018); mpfr_init(r24019); mpfr_init(r24020); mpfr_init(r24021); mpfr_init(r24022); mpfr_init(r24023); mpfr_init_set_str(r24024, "-5.270434816906832e-10", 10, MPFR_RNDN); mpfr_init(r24025); mpfr_init_set_str(r24026, "1", 10, MPFR_RNDN); mpfr_init(r24027); mpfr_init(r24028); mpfr_init(r24029); mpfr_init_set_str(r24030, "2.088519104255703e-20", 10, MPFR_RNDN); mpfr_init(r24031); mpfr_init(r24032); mpfr_init(r24033); } double f_fm(double a, double b, double eps) { mpfr_set_d(r24011, eps, MPFR_RNDN); mpfr_set_d(r24012, b, MPFR_RNDN); mpfr_mul(r24013, r24012, r24011, MPFR_RNDN); mpfr_expm1(r24014, r24013, MPFR_RNDN); mpfr_set_d(r24015, a, MPFR_RNDN); mpfr_mul(r24016, r24011, r24015, MPFR_RNDN); mpfr_expm1(r24017, r24016, MPFR_RNDN); mpfr_mul(r24018, r24014, r24017, MPFR_RNDN); mpfr_add(r24019, r24012, r24015, MPFR_RNDN); mpfr_mul(r24020, r24019, r24011, MPFR_RNDN); mpfr_expm1(r24021, r24020, MPFR_RNDN); mpfr_div(r24022, r24018, r24021, MPFR_RNDN); mpfr_div(r24023, r24011, r24022, MPFR_RNDN); ; mpfr_set_si(r24025, mpfr_cmp(r24023, r24024) <= 0, MPFR_RNDN); ; mpfr_div(r24027, r24026, r24012, MPFR_RNDN); mpfr_div(r24028, r24026, r24015, MPFR_RNDN); mpfr_add(r24029, r24027, r24028, MPFR_RNDN); ; mpfr_set_si(r24031, mpfr_cmp(r24023, r24030) <= 0, MPFR_RNDN); if (mpfr_get_si(r24031, MPFR_RNDN)) { mpfr_set(r24032, r24023, MPFR_RNDN); } else { mpfr_set(r24032, r24029, MPFR_RNDN); }; if (mpfr_get_si(r24025, MPFR_RNDN)) { mpfr_set(r24033, r24029, MPFR_RNDN); } else { mpfr_set(r24033, r24032, MPFR_RNDN); }; return mpfr_get_d(r24033, MPFR_RNDN); } static mpfr_t r24034, r24035, r24036, r24037, r24038, r24039, r24040, r24041, r24042, r24043, r24044, r24045, r24046, r24047, r24048, r24049, r24050, r24051, r24052, r24053, r24054, r24055, r24056; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r24034); mpfr_init(r24035); mpfr_init(r24036); mpfr_init(r24037); mpfr_init(r24038); mpfr_init(r24039); mpfr_init(r24040); mpfr_init(r24041); mpfr_init(r24042); mpfr_init(r24043); mpfr_init(r24044); mpfr_init(r24045); mpfr_init(r24046); mpfr_init_set_str(r24047, "-5.270434816906832e-10", 10, MPFR_RNDN); mpfr_init(r24048); mpfr_init_set_str(r24049, "1", 10, MPFR_RNDN); mpfr_init(r24050); mpfr_init(r24051); mpfr_init(r24052); mpfr_init_set_str(r24053, "2.088519104255703e-20", 10, MPFR_RNDN); mpfr_init(r24054); mpfr_init(r24055); mpfr_init(r24056); } double f_dm(double a, double b, double eps) { mpfr_set_d(r24034, eps, MPFR_RNDN); mpfr_set_d(r24035, b, MPFR_RNDN); mpfr_mul(r24036, r24035, r24034, MPFR_RNDN); mpfr_expm1(r24037, r24036, MPFR_RNDN); mpfr_set_d(r24038, a, MPFR_RNDN); mpfr_mul(r24039, r24034, r24038, MPFR_RNDN); mpfr_expm1(r24040, r24039, MPFR_RNDN); mpfr_mul(r24041, r24037, r24040, MPFR_RNDN); mpfr_add(r24042, r24035, r24038, MPFR_RNDN); mpfr_mul(r24043, r24042, r24034, MPFR_RNDN); mpfr_expm1(r24044, r24043, MPFR_RNDN); mpfr_div(r24045, r24041, r24044, MPFR_RNDN); mpfr_div(r24046, r24034, r24045, MPFR_RNDN); ; mpfr_set_si(r24048, mpfr_cmp(r24046, r24047) <= 0, MPFR_RNDN); ; mpfr_div(r24050, r24049, r24035, MPFR_RNDN); mpfr_div(r24051, r24049, r24038, MPFR_RNDN); mpfr_add(r24052, r24050, r24051, MPFR_RNDN); ; mpfr_set_si(r24054, mpfr_cmp(r24046, r24053) <= 0, MPFR_RNDN); if (mpfr_get_si(r24054, MPFR_RNDN)) { mpfr_set(r24055, r24046, MPFR_RNDN); } else { mpfr_set(r24055, r24052, MPFR_RNDN); }; if (mpfr_get_si(r24048, MPFR_RNDN)) { mpfr_set(r24056, r24052, MPFR_RNDN); } else { mpfr_set(r24056, r24055, MPFR_RNDN); }; return mpfr_get_d(r24056, MPFR_RNDN); }