#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (9a)"; double f_if(float w0, float M, float D, float h, float l, float d) { float r25925 = w0; float r25926 = 1; float r25927 = M; float r25928 = D; float r25929 = r25927 * r25928; float r25930 = 2; float r25931 = d; float r25932 = r25930 * r25931; float r25933 = r25929 / r25932; float r25934 = pow(r25933, r25930); float r25935 = h; float r25936 = l; float r25937 = r25935 / r25936; float r25938 = r25934 * r25937; float r25939 = r25926 - r25938; float r25940 = sqrt(r25939); float r25941 = r25925 * r25940; return r25941; } double f_id(double w0, double M, double D, double h, double l, double d) { double r25942 = w0; double r25943 = 1; double r25944 = M; double r25945 = D; double r25946 = r25944 * r25945; double r25947 = 2; double r25948 = d; double r25949 = r25947 * r25948; double r25950 = r25946 / r25949; double r25951 = pow(r25950, r25947); double r25952 = h; double r25953 = l; double r25954 = r25952 / r25953; double r25955 = r25951 * r25954; double r25956 = r25943 - r25955; double r25957 = sqrt(r25956); double r25958 = r25942 * r25957; return r25958; } double f_of(float w0, float M, float D, float h, float l, float d) { float r25959 = 1; float r25960 = M; float r25961 = d; float r25962 = r25960 / r25961; float r25963 = h; float r25964 = D; float r25965 = r25963 * r25964; float r25966 = 1/2; float r25967 = r25965 * r25966; float r25968 = r25962 * r25967; float r25969 = 2; float r25970 = r25964 / r25969; float r25971 = r25970 * r25962; float r25972 = l; float r25973 = r25971 / r25972; float r25974 = r25968 * r25973; float r25975 = r25959 - r25974; float r25976 = sqrt(r25975); float r25977 = 4.712989610511589e+137; bool r25978 = r25976 <= r25977; float r25979 = w0; float r25980 = r25976 * r25979; float r25981 = r25960 * r25964; float r25982 = r25969 * r25961; float r25983 = r25981 / r25982; float r25984 = cbrt(r25983); float r25985 = r25984 * r25984; float r25986 = pow(r25985, r25969); float r25987 = pow(r25984, r25969); float r25988 = r25987 * r25963; float r25989 = r25986 * r25988; float r25990 = r25959 / r25972; float r25991 = r25989 * r25990; float r25992 = r25959 - r25991; float r25993 = sqrt(r25992); float r25994 = r25979 * r25993; float r25995 = r25978 ? r25980 : r25994; return r25995; } double f_od(double w0, double M, double D, double h, double l, double d) { double r25996 = 1; double r25997 = M; double r25998 = d; double r25999 = r25997 / r25998; double r26000 = h; double r26001 = D; double r26002 = r26000 * r26001; double r26003 = 1/2; double r26004 = r26002 * r26003; double r26005 = r25999 * r26004; double r26006 = 2; double r26007 = r26001 / r26006; double r26008 = r26007 * r25999; double r26009 = l; double r26010 = r26008 / r26009; double r26011 = r26005 * r26010; double r26012 = r25996 - r26011; double r26013 = sqrt(r26012); double r26014 = 4.712989610511589e+137; bool r26015 = r26013 <= r26014; double r26016 = w0; double r26017 = r26013 * r26016; double r26018 = r25997 * r26001; double r26019 = r26006 * r25998; double r26020 = r26018 / r26019; double r26021 = cbrt(r26020); double r26022 = r26021 * r26021; double r26023 = pow(r26022, r26006); double r26024 = pow(r26021, r26006); double r26025 = r26024 * r26000; double r26026 = r26023 * r26025; double r26027 = r25996 / r26009; double r26028 = r26026 * r26027; double r26029 = r25996 - r26028; double r26030 = sqrt(r26029); double r26031 = r26016 * r26030; double r26032 = r26015 ? r26017 : r26031; return r26032; } 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 r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26033); mpfr_init_set_str(r26034, "1", 10, MPFR_RNDN); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init_set_str(r26038, "2", 10, MPFR_RNDN); mpfr_init(r26039); mpfr_init(r26040); mpfr_init(r26041); mpfr_init(r26042); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); mpfr_init(r26048); mpfr_init(r26049); } double f_im(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r26033, w0, MPFR_RNDN); ; mpfr_set_d(r26035, M, MPFR_RNDN); mpfr_set_d(r26036, D, MPFR_RNDN); mpfr_mul(r26037, r26035, r26036, MPFR_RNDN); ; mpfr_set_d(r26039, d, MPFR_RNDN); mpfr_mul(r26040, r26038, r26039, MPFR_RNDN); mpfr_div(r26041, r26037, r26040, MPFR_RNDN); mpfr_pow(r26042, r26041, r26038, MPFR_RNDN); mpfr_set_d(r26043, h, MPFR_RNDN); mpfr_set_d(r26044, l, MPFR_RNDN); mpfr_div(r26045, r26043, r26044, MPFR_RNDN); mpfr_mul(r26046, r26042, r26045, MPFR_RNDN); mpfr_sub(r26047, r26034, r26046, MPFR_RNDN); mpfr_sqrt(r26048, r26047, MPFR_RNDN); mpfr_mul(r26049, r26033, r26048, MPFR_RNDN); return mpfr_get_d(r26049, MPFR_RNDN); } static mpfr_t r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26050, "1", 10, MPFR_RNDN); mpfr_init(r26051); mpfr_init(r26052); mpfr_init(r26053); mpfr_init(r26054); mpfr_init(r26055); mpfr_init(r26056); mpfr_init_set_str(r26057, "1/2", 10, MPFR_RNDN); mpfr_init(r26058); mpfr_init(r26059); mpfr_init_set_str(r26060, "2", 10, MPFR_RNDN); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); mpfr_init(r26065); mpfr_init(r26066); mpfr_init(r26067); mpfr_init_set_str(r26068, "4.712989610511589e+137", 10, MPFR_RNDN); mpfr_init(r26069); mpfr_init(r26070); mpfr_init(r26071); mpfr_init(r26072); mpfr_init(r26073); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init(r26084); mpfr_init(r26085); mpfr_init(r26086); } double f_fm(double w0, double M, double D, double h, double l, double d) { ; mpfr_set_d(r26051, M, MPFR_RNDN); mpfr_set_d(r26052, d, MPFR_RNDN); mpfr_div(r26053, r26051, r26052, MPFR_RNDN); mpfr_set_d(r26054, h, MPFR_RNDN); mpfr_set_d(r26055, D, MPFR_RNDN); mpfr_mul(r26056, r26054, r26055, MPFR_RNDN); ; mpfr_mul(r26058, r26056, r26057, MPFR_RNDN); mpfr_mul(r26059, r26053, r26058, MPFR_RNDN); ; mpfr_div(r26061, r26055, r26060, MPFR_RNDN); mpfr_mul(r26062, r26061, r26053, MPFR_RNDN); mpfr_set_d(r26063, l, MPFR_RNDN); mpfr_div(r26064, r26062, r26063, MPFR_RNDN); mpfr_mul(r26065, r26059, r26064, MPFR_RNDN); mpfr_sub(r26066, r26050, r26065, MPFR_RNDN); mpfr_sqrt(r26067, r26066, MPFR_RNDN); ; mpfr_set_si(r26069, mpfr_cmp(r26067, r26068) <= 0, MPFR_RNDN); mpfr_set_d(r26070, w0, MPFR_RNDN); mpfr_mul(r26071, r26067, r26070, MPFR_RNDN); mpfr_mul(r26072, r26051, r26055, MPFR_RNDN); mpfr_mul(r26073, r26060, r26052, MPFR_RNDN); mpfr_div(r26074, r26072, r26073, MPFR_RNDN); mpfr_cbrt(r26075, r26074, MPFR_RNDN); mpfr_mul(r26076, r26075, r26075, MPFR_RNDN); mpfr_pow(r26077, r26076, r26060, MPFR_RNDN); mpfr_pow(r26078, r26075, r26060, MPFR_RNDN); mpfr_mul(r26079, r26078, r26054, MPFR_RNDN); mpfr_mul(r26080, r26077, r26079, MPFR_RNDN); mpfr_div(r26081, r26050, r26063, MPFR_RNDN); mpfr_mul(r26082, r26080, r26081, MPFR_RNDN); mpfr_sub(r26083, r26050, r26082, MPFR_RNDN); mpfr_sqrt(r26084, r26083, MPFR_RNDN); mpfr_mul(r26085, r26070, r26084, MPFR_RNDN); if (mpfr_get_si(r26069, MPFR_RNDN)) { mpfr_set(r26086, r26071, MPFR_RNDN); } else { mpfr_set(r26086, r26085, MPFR_RNDN); }; return mpfr_get_d(r26086, MPFR_RNDN); } static mpfr_t r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26087, "1", 10, MPFR_RNDN); mpfr_init(r26088); mpfr_init(r26089); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); mpfr_init(r26093); mpfr_init_set_str(r26094, "1/2", 10, MPFR_RNDN); mpfr_init(r26095); mpfr_init(r26096); mpfr_init_set_str(r26097, "2", 10, MPFR_RNDN); mpfr_init(r26098); mpfr_init(r26099); mpfr_init(r26100); mpfr_init(r26101); mpfr_init(r26102); mpfr_init(r26103); mpfr_init(r26104); mpfr_init_set_str(r26105, "4.712989610511589e+137", 10, MPFR_RNDN); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init(r26109); mpfr_init(r26110); mpfr_init(r26111); mpfr_init(r26112); mpfr_init(r26113); mpfr_init(r26114); mpfr_init(r26115); mpfr_init(r26116); mpfr_init(r26117); mpfr_init(r26118); mpfr_init(r26119); mpfr_init(r26120); mpfr_init(r26121); mpfr_init(r26122); mpfr_init(r26123); } double f_dm(double w0, double M, double D, double h, double l, double d) { ; mpfr_set_d(r26088, M, MPFR_RNDN); mpfr_set_d(r26089, d, MPFR_RNDN); mpfr_div(r26090, r26088, r26089, MPFR_RNDN); mpfr_set_d(r26091, h, MPFR_RNDN); mpfr_set_d(r26092, D, MPFR_RNDN); mpfr_mul(r26093, r26091, r26092, MPFR_RNDN); ; mpfr_mul(r26095, r26093, r26094, MPFR_RNDN); mpfr_mul(r26096, r26090, r26095, MPFR_RNDN); ; mpfr_div(r26098, r26092, r26097, MPFR_RNDN); mpfr_mul(r26099, r26098, r26090, MPFR_RNDN); mpfr_set_d(r26100, l, MPFR_RNDN); mpfr_div(r26101, r26099, r26100, MPFR_RNDN); mpfr_mul(r26102, r26096, r26101, MPFR_RNDN); mpfr_sub(r26103, r26087, r26102, MPFR_RNDN); mpfr_sqrt(r26104, r26103, MPFR_RNDN); ; mpfr_set_si(r26106, mpfr_cmp(r26104, r26105) <= 0, MPFR_RNDN); mpfr_set_d(r26107, w0, MPFR_RNDN); mpfr_mul(r26108, r26104, r26107, MPFR_RNDN); mpfr_mul(r26109, r26088, r26092, MPFR_RNDN); mpfr_mul(r26110, r26097, r26089, MPFR_RNDN); mpfr_div(r26111, r26109, r26110, MPFR_RNDN); mpfr_cbrt(r26112, r26111, MPFR_RNDN); mpfr_mul(r26113, r26112, r26112, MPFR_RNDN); mpfr_pow(r26114, r26113, r26097, MPFR_RNDN); mpfr_pow(r26115, r26112, r26097, MPFR_RNDN); mpfr_mul(r26116, r26115, r26091, MPFR_RNDN); mpfr_mul(r26117, r26114, r26116, MPFR_RNDN); mpfr_div(r26118, r26087, r26100, MPFR_RNDN); mpfr_mul(r26119, r26117, r26118, MPFR_RNDN); mpfr_sub(r26120, r26087, r26119, MPFR_RNDN); mpfr_sqrt(r26121, r26120, MPFR_RNDN); mpfr_mul(r26122, r26107, r26121, MPFR_RNDN); if (mpfr_get_si(r26106, MPFR_RNDN)) { mpfr_set(r26123, r26108, MPFR_RNDN); } else { mpfr_set(r26123, r26122, MPFR_RNDN); }; return mpfr_get_d(r26123, MPFR_RNDN); }