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