#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (13)"; double f_if(float c0, float w, float h, float D, float d, float M) { float r23938 = c0; float r23939 = 2; float r23940 = w; float r23941 = r23939 * r23940; float r23942 = r23938 / r23941; float r23943 = d; float r23944 = r23943 * r23943; float r23945 = r23938 * r23944; float r23946 = h; float r23947 = r23940 * r23946; float r23948 = D; float r23949 = r23948 * r23948; float r23950 = r23947 * r23949; float r23951 = r23945 / r23950; float r23952 = r23951 * r23951; float r23953 = M; float r23954 = r23953 * r23953; float r23955 = r23952 - r23954; float r23956 = sqrt(r23955); float r23957 = r23951 + r23956; float r23958 = r23942 * r23957; return r23958; } double f_id(double c0, double w, double h, double D, double d, double M) { double r23959 = c0; double r23960 = 2; double r23961 = w; double r23962 = r23960 * r23961; double r23963 = r23959 / r23962; double r23964 = d; double r23965 = r23964 * r23964; double r23966 = r23959 * r23965; double r23967 = h; double r23968 = r23961 * r23967; double r23969 = D; double r23970 = r23969 * r23969; double r23971 = r23968 * r23970; double r23972 = r23966 / r23971; double r23973 = r23972 * r23972; double r23974 = M; double r23975 = r23974 * r23974; double r23976 = r23973 - r23975; double r23977 = sqrt(r23976); double r23978 = r23972 + r23977; double r23979 = r23963 * r23978; return r23979; } double f_of(float c0, float w, float h, float D, float d, float M) { float r23980 = c0; float r23981 = 2; float r23982 = w; float r23983 = r23981 * r23982; float r23984 = r23980 / r23983; float r23985 = r23980 / r23982; float r23986 = h; float r23987 = r23985 / r23986; float r23988 = d; float r23989 = D; float r23990 = r23988 / r23989; float r23991 = r23990 * r23990; float r23992 = r23987 * r23991; float r23993 = M; float r23994 = r23993 * r23993; float r23995 = -r23994; float r23996 = fma(r23992, r23992, r23995); float r23997 = sqrt(r23996); float r23998 = r23997 + r23992; float r23999 = 3; float r24000 = pow(r23998, r23999); float r24001 = cbrt(r24000); float r24002 = r23984 * r24001; float r24003 = +inf.0; bool r24004 = r24002 <= r24003; float r24005 = 0; float r24006 = r24004 ? r24002 : r24005; return r24006; } double f_od(double c0, double w, double h, double D, double d, double M) { double r24007 = c0; double r24008 = 2; double r24009 = w; double r24010 = r24008 * r24009; double r24011 = r24007 / r24010; double r24012 = r24007 / r24009; double r24013 = h; double r24014 = r24012 / r24013; double r24015 = d; double r24016 = D; double r24017 = r24015 / r24016; double r24018 = r24017 * r24017; double r24019 = r24014 * r24018; double r24020 = M; double r24021 = r24020 * r24020; double r24022 = -r24021; double r24023 = fma(r24019, r24019, r24022); double r24024 = sqrt(r24023); double r24025 = r24024 + r24019; double r24026 = 3; double r24027 = pow(r24025, r24026); double r24028 = cbrt(r24027); double r24029 = r24011 * r24028; double r24030 = +inf.0; bool r24031 = r24029 <= r24030; double r24032 = 0; double r24033 = r24031 ? r24029 : r24032; return r24033; } 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 r24034, r24035, r24036, r24037, r24038, r24039, r24040, r24041, r24042, r24043, r24044, r24045, r24046, r24047, r24048, r24049, r24050, r24051, r24052, r24053, r24054; void setup_mpfr_f_im() { mpfr_set_default_prec(7312); mpfr_init(r24034); mpfr_init_set_str(r24035, "2", 10, MPFR_RNDN); 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(r24047); mpfr_init(r24048); mpfr_init(r24049); mpfr_init(r24050); mpfr_init(r24051); mpfr_init(r24052); mpfr_init(r24053); mpfr_init(r24054); } double f_im(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r24034, c0, MPFR_RNDN); ; mpfr_set_d(r24036, w, MPFR_RNDN); mpfr_mul(r24037, r24035, r24036, MPFR_RNDN); mpfr_div(r24038, r24034, r24037, MPFR_RNDN); mpfr_set_d(r24039, d, MPFR_RNDN); mpfr_mul(r24040, r24039, r24039, MPFR_RNDN); mpfr_mul(r24041, r24034, r24040, MPFR_RNDN); mpfr_set_d(r24042, h, MPFR_RNDN); mpfr_mul(r24043, r24036, r24042, MPFR_RNDN); mpfr_set_d(r24044, D, MPFR_RNDN); mpfr_mul(r24045, r24044, r24044, MPFR_RNDN); mpfr_mul(r24046, r24043, r24045, MPFR_RNDN); mpfr_div(r24047, r24041, r24046, MPFR_RNDN); mpfr_mul(r24048, r24047, r24047, MPFR_RNDN); mpfr_set_d(r24049, M, MPFR_RNDN); mpfr_mul(r24050, r24049, r24049, MPFR_RNDN); mpfr_sub(r24051, r24048, r24050, MPFR_RNDN); mpfr_sqrt(r24052, r24051, MPFR_RNDN); mpfr_add(r24053, r24047, r24052, MPFR_RNDN); mpfr_mul(r24054, r24038, r24053, MPFR_RNDN); return mpfr_get_d(r24054, MPFR_RNDN); } static mpfr_t r24055, r24056, r24057, r24058, r24059, r24060, r24061, r24062, r24063, r24064, r24065, r24066, r24067, r24068, r24069, r24070, r24071, r24072, r24073, r24074, r24075, r24076, r24077, r24078, r24079, r24080, r24081; void setup_mpfr_f_fm() { mpfr_set_default_prec(7312); mpfr_init(r24055); mpfr_init_set_str(r24056, "2", 10, MPFR_RNDN); mpfr_init(r24057); mpfr_init(r24058); mpfr_init(r24059); mpfr_init(r24060); mpfr_init(r24061); mpfr_init(r24062); mpfr_init(r24063); mpfr_init(r24064); mpfr_init(r24065); mpfr_init(r24066); mpfr_init(r24067); mpfr_init(r24068); mpfr_init(r24069); mpfr_init(r24070); mpfr_init(r24071); mpfr_init(r24072); mpfr_init(r24073); mpfr_init_set_str(r24074, "3", 10, MPFR_RNDN); mpfr_init(r24075); mpfr_init(r24076); mpfr_init(r24077); mpfr_init_set_str(r24078, "+inf.0", 10, MPFR_RNDN); mpfr_init(r24079); mpfr_init_set_str(r24080, "0", 10, MPFR_RNDN); mpfr_init(r24081); } double f_fm(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r24055, c0, MPFR_RNDN); ; mpfr_set_d(r24057, w, MPFR_RNDN); mpfr_mul(r24058, r24056, r24057, MPFR_RNDN); mpfr_div(r24059, r24055, r24058, MPFR_RNDN); mpfr_div(r24060, r24055, r24057, MPFR_RNDN); mpfr_set_d(r24061, h, MPFR_RNDN); mpfr_div(r24062, r24060, r24061, MPFR_RNDN); mpfr_set_d(r24063, d, MPFR_RNDN); mpfr_set_d(r24064, D, MPFR_RNDN); mpfr_div(r24065, r24063, r24064, MPFR_RNDN); mpfr_mul(r24066, r24065, r24065, MPFR_RNDN); mpfr_mul(r24067, r24062, r24066, MPFR_RNDN); mpfr_set_d(r24068, M, MPFR_RNDN); mpfr_mul(r24069, r24068, r24068, MPFR_RNDN); mpfr_neg(r24070, r24069, MPFR_RNDN); mpfr_fma(r24071, r24067, r24067, r24070, MPFR_RNDN); mpfr_sqrt(r24072, r24071, MPFR_RNDN); mpfr_add(r24073, r24072, r24067, MPFR_RNDN); ; mpfr_pow(r24075, r24073, r24074, MPFR_RNDN); mpfr_cbrt(r24076, r24075, MPFR_RNDN); mpfr_mul(r24077, r24059, r24076, MPFR_RNDN); ; mpfr_set_si(r24079, mpfr_cmp(r24077, r24078) <= 0, MPFR_RNDN); ; if (mpfr_get_si(r24079, MPFR_RNDN)) { mpfr_set(r24081, r24077, MPFR_RNDN); } else { mpfr_set(r24081, r24080, MPFR_RNDN); }; return mpfr_get_d(r24081, MPFR_RNDN); } static mpfr_t r24082, r24083, r24084, r24085, r24086, r24087, r24088, r24089, r24090, r24091, r24092, r24093, r24094, r24095, r24096, r24097, r24098, r24099, r24100, r24101, r24102, r24103, r24104, r24105, r24106, r24107, r24108; void setup_mpfr_f_dm() { mpfr_set_default_prec(7312); mpfr_init(r24082); mpfr_init_set_str(r24083, "2", 10, MPFR_RNDN); mpfr_init(r24084); mpfr_init(r24085); mpfr_init(r24086); mpfr_init(r24087); mpfr_init(r24088); mpfr_init(r24089); mpfr_init(r24090); mpfr_init(r24091); mpfr_init(r24092); mpfr_init(r24093); mpfr_init(r24094); mpfr_init(r24095); mpfr_init(r24096); mpfr_init(r24097); mpfr_init(r24098); mpfr_init(r24099); mpfr_init(r24100); mpfr_init_set_str(r24101, "3", 10, MPFR_RNDN); mpfr_init(r24102); mpfr_init(r24103); mpfr_init(r24104); mpfr_init_set_str(r24105, "+inf.0", 10, MPFR_RNDN); mpfr_init(r24106); mpfr_init_set_str(r24107, "0", 10, MPFR_RNDN); mpfr_init(r24108); } double f_dm(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r24082, c0, MPFR_RNDN); ; mpfr_set_d(r24084, w, MPFR_RNDN); mpfr_mul(r24085, r24083, r24084, MPFR_RNDN); mpfr_div(r24086, r24082, r24085, MPFR_RNDN); mpfr_div(r24087, r24082, r24084, MPFR_RNDN); mpfr_set_d(r24088, h, MPFR_RNDN); mpfr_div(r24089, r24087, r24088, MPFR_RNDN); mpfr_set_d(r24090, d, MPFR_RNDN); mpfr_set_d(r24091, D, MPFR_RNDN); mpfr_div(r24092, r24090, r24091, MPFR_RNDN); mpfr_mul(r24093, r24092, r24092, MPFR_RNDN); mpfr_mul(r24094, r24089, r24093, MPFR_RNDN); mpfr_set_d(r24095, M, MPFR_RNDN); mpfr_mul(r24096, r24095, r24095, MPFR_RNDN); mpfr_neg(r24097, r24096, MPFR_RNDN); mpfr_fma(r24098, r24094, r24094, r24097, MPFR_RNDN); mpfr_sqrt(r24099, r24098, MPFR_RNDN); mpfr_add(r24100, r24099, r24094, MPFR_RNDN); ; mpfr_pow(r24102, r24100, r24101, MPFR_RNDN); mpfr_cbrt(r24103, r24102, MPFR_RNDN); mpfr_mul(r24104, r24086, r24103, MPFR_RNDN); ; mpfr_set_si(r24106, mpfr_cmp(r24104, r24105) <= 0, MPFR_RNDN); ; if (mpfr_get_si(r24106, MPFR_RNDN)) { mpfr_set(r24108, r24104, MPFR_RNDN); } else { mpfr_set(r24108, r24107, MPFR_RNDN); }; return mpfr_get_d(r24108, MPFR_RNDN); }