#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 r21961 = c0; float r21962 = 2; float r21963 = w; float r21964 = r21962 * r21963; float r21965 = r21961 / r21964; float r21966 = d; float r21967 = r21966 * r21966; float r21968 = r21961 * r21967; float r21969 = h; float r21970 = r21963 * r21969; float r21971 = D; float r21972 = r21971 * r21971; float r21973 = r21970 * r21972; float r21974 = r21968 / r21973; float r21975 = r21974 * r21974; float r21976 = M; float r21977 = r21976 * r21976; float r21978 = r21975 - r21977; float r21979 = sqrt(r21978); float r21980 = r21974 + r21979; float r21981 = r21965 * r21980; return r21981; } double f_id(double c0, double w, double h, double D, double d, double M) { double r21982 = c0; double r21983 = 2; double r21984 = w; double r21985 = r21983 * r21984; double r21986 = r21982 / r21985; double r21987 = d; double r21988 = r21987 * r21987; double r21989 = r21982 * r21988; double r21990 = h; double r21991 = r21984 * r21990; double r21992 = D; double r21993 = r21992 * r21992; double r21994 = r21991 * r21993; double r21995 = r21989 / r21994; double r21996 = r21995 * r21995; double r21997 = M; double r21998 = r21997 * r21997; double r21999 = r21996 - r21998; double r22000 = sqrt(r21999); double r22001 = r21995 + r22000; double r22002 = r21986 * r22001; return r22002; } double f_of(float c0, float w, float h, float D, float d, float __attribute__((unused)) M) { float r22003 = D; float r22004 = 3.973420093857094e-14; bool r22005 = r22003 <= r22004; float r22006 = 0; float r22007 = 3.634866208990275e+21; bool r22008 = r22003 <= r22007; float r22009 = c0; float r22010 = w; float r22011 = r22010 + r22010; float r22012 = r22009 / r22011; float r22013 = d; float r22014 = r22013 * r22009; float r22015 = r22013 + r22013; float r22016 = r22014 * r22015; float r22017 = r22003 * r22003; float r22018 = h; float r22019 = r22010 * r22018; float r22020 = r22017 * r22019; float r22021 = r22016 / r22020; float r22022 = r22012 * r22021; float r22023 = r22008 ? r22022 : r22006; float r22024 = r22005 ? r22006 : r22023; return r22024; } double f_od(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { double r22025 = D; double r22026 = 3.973420093857094e-14; bool r22027 = r22025 <= r22026; double r22028 = 0; double r22029 = 3.634866208990275e+21; bool r22030 = r22025 <= r22029; double r22031 = c0; double r22032 = w; double r22033 = r22032 + r22032; double r22034 = r22031 / r22033; double r22035 = d; double r22036 = r22035 * r22031; double r22037 = r22035 + r22035; double r22038 = r22036 * r22037; double r22039 = r22025 * r22025; double r22040 = h; double r22041 = r22032 * r22040; double r22042 = r22039 * r22041; double r22043 = r22038 / r22042; double r22044 = r22034 * r22043; double r22045 = r22030 ? r22044 : r22028; double r22046 = r22027 ? r22028 : r22045; return r22046; } 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 r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067; void setup_mpfr_f_im() { mpfr_set_default_prec(7056); mpfr_init(r22047); mpfr_init_set_str(r22048, "2", 10, MPFR_RNDN); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); } double f_im(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r22047, c0, MPFR_RNDN); ; mpfr_set_d(r22049, w, MPFR_RNDN); mpfr_mul(r22050, r22048, r22049, MPFR_RNDN); mpfr_div(r22051, r22047, r22050, MPFR_RNDN); mpfr_set_d(r22052, d, MPFR_RNDN); mpfr_mul(r22053, r22052, r22052, MPFR_RNDN); mpfr_mul(r22054, r22047, r22053, MPFR_RNDN); mpfr_set_d(r22055, h, MPFR_RNDN); mpfr_mul(r22056, r22049, r22055, MPFR_RNDN); mpfr_set_d(r22057, D, MPFR_RNDN); mpfr_mul(r22058, r22057, r22057, MPFR_RNDN); mpfr_mul(r22059, r22056, r22058, MPFR_RNDN); mpfr_div(r22060, r22054, r22059, MPFR_RNDN); mpfr_mul(r22061, r22060, r22060, MPFR_RNDN); mpfr_set_d(r22062, M, MPFR_RNDN); mpfr_mul(r22063, r22062, r22062, MPFR_RNDN); mpfr_sub(r22064, r22061, r22063, MPFR_RNDN); mpfr_sqrt(r22065, r22064, MPFR_RNDN); mpfr_add(r22066, r22060, r22065, MPFR_RNDN); mpfr_mul(r22067, r22051, r22066, MPFR_RNDN); return mpfr_get_d(r22067, MPFR_RNDN); } static mpfr_t r22068, r22069, r22070, r22071, r22072, r22073, r22074, r22075, r22076, r22077, r22078, r22079, r22080, r22081, r22082, r22083, r22084, r22085, r22086, r22087, r22088, r22089; void setup_mpfr_f_fm() { mpfr_set_default_prec(7056); mpfr_init(r22068); mpfr_init_set_str(r22069, "3.973420093857094e-14", 10, MPFR_RNDN); mpfr_init(r22070); mpfr_init_set_str(r22071, "0", 10, MPFR_RNDN); mpfr_init_set_str(r22072, "3.634866208990275e+21", 10, MPFR_RNDN); mpfr_init(r22073); mpfr_init(r22074); mpfr_init(r22075); mpfr_init(r22076); mpfr_init(r22077); mpfr_init(r22078); mpfr_init(r22079); mpfr_init(r22080); mpfr_init(r22081); mpfr_init(r22082); mpfr_init(r22083); mpfr_init(r22084); mpfr_init(r22085); mpfr_init(r22086); mpfr_init(r22087); mpfr_init(r22088); mpfr_init(r22089); } double f_fm(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { mpfr_set_d(r22068, D, MPFR_RNDN); ; mpfr_set_si(r22070, mpfr_cmp(r22068, r22069) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r22073, mpfr_cmp(r22068, r22072) <= 0, MPFR_RNDN); mpfr_set_d(r22074, c0, MPFR_RNDN); mpfr_set_d(r22075, w, MPFR_RNDN); mpfr_add(r22076, r22075, r22075, MPFR_RNDN); mpfr_div(r22077, r22074, r22076, MPFR_RNDN); mpfr_set_d(r22078, d, MPFR_RNDN); mpfr_mul(r22079, r22078, r22074, MPFR_RNDN); mpfr_add(r22080, r22078, r22078, MPFR_RNDN); mpfr_mul(r22081, r22079, r22080, MPFR_RNDN); mpfr_mul(r22082, r22068, r22068, MPFR_RNDN); mpfr_set_d(r22083, h, MPFR_RNDN); mpfr_mul(r22084, r22075, r22083, MPFR_RNDN); mpfr_mul(r22085, r22082, r22084, MPFR_RNDN); mpfr_div(r22086, r22081, r22085, MPFR_RNDN); mpfr_mul(r22087, r22077, r22086, MPFR_RNDN); if (mpfr_get_si(r22073, MPFR_RNDN)) { mpfr_set(r22088, r22087, MPFR_RNDN); } else { mpfr_set(r22088, r22071, MPFR_RNDN); }; if (mpfr_get_si(r22070, MPFR_RNDN)) { mpfr_set(r22089, r22071, MPFR_RNDN); } else { mpfr_set(r22089, r22088, MPFR_RNDN); }; return mpfr_get_d(r22089, MPFR_RNDN); } static mpfr_t r22090, r22091, r22092, r22093, r22094, r22095, r22096, r22097, r22098, r22099, r22100, r22101, r22102, r22103, r22104, r22105, r22106, r22107, r22108, r22109, r22110, r22111; void setup_mpfr_f_dm() { mpfr_set_default_prec(7056); mpfr_init(r22090); mpfr_init_set_str(r22091, "3.973420093857094e-14", 10, MPFR_RNDN); mpfr_init(r22092); mpfr_init_set_str(r22093, "0", 10, MPFR_RNDN); mpfr_init_set_str(r22094, "3.634866208990275e+21", 10, MPFR_RNDN); mpfr_init(r22095); mpfr_init(r22096); mpfr_init(r22097); mpfr_init(r22098); mpfr_init(r22099); mpfr_init(r22100); mpfr_init(r22101); mpfr_init(r22102); mpfr_init(r22103); mpfr_init(r22104); mpfr_init(r22105); mpfr_init(r22106); mpfr_init(r22107); mpfr_init(r22108); mpfr_init(r22109); mpfr_init(r22110); mpfr_init(r22111); } double f_dm(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { mpfr_set_d(r22090, D, MPFR_RNDN); ; mpfr_set_si(r22092, mpfr_cmp(r22090, r22091) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r22095, mpfr_cmp(r22090, r22094) <= 0, MPFR_RNDN); mpfr_set_d(r22096, c0, MPFR_RNDN); mpfr_set_d(r22097, w, MPFR_RNDN); mpfr_add(r22098, r22097, r22097, MPFR_RNDN); mpfr_div(r22099, r22096, r22098, MPFR_RNDN); mpfr_set_d(r22100, d, MPFR_RNDN); mpfr_mul(r22101, r22100, r22096, MPFR_RNDN); mpfr_add(r22102, r22100, r22100, MPFR_RNDN); mpfr_mul(r22103, r22101, r22102, MPFR_RNDN); mpfr_mul(r22104, r22090, r22090, MPFR_RNDN); mpfr_set_d(r22105, h, MPFR_RNDN); mpfr_mul(r22106, r22097, r22105, MPFR_RNDN); mpfr_mul(r22107, r22104, r22106, MPFR_RNDN); mpfr_div(r22108, r22103, r22107, MPFR_RNDN); mpfr_mul(r22109, r22099, r22108, MPFR_RNDN); if (mpfr_get_si(r22095, MPFR_RNDN)) { mpfr_set(r22110, r22109, MPFR_RNDN); } else { mpfr_set(r22110, r22093, MPFR_RNDN); }; if (mpfr_get_si(r22092, MPFR_RNDN)) { mpfr_set(r22111, r22093, MPFR_RNDN); } else { mpfr_set(r22111, r22110, MPFR_RNDN); }; return mpfr_get_d(r22111, MPFR_RNDN); }