#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (12)"; double f_if(float d, float h, float l, float M, float D) { float r26843 = d; float r26844 = h; float r26845 = r26843 / r26844; float r26846 = 1; float r26847 = 2; float r26848 = r26846 / r26847; float r26849 = pow(r26845, r26848); float r26850 = l; float r26851 = r26843 / r26850; float r26852 = pow(r26851, r26848); float r26853 = r26849 * r26852; float r26854 = M; float r26855 = D; float r26856 = r26854 * r26855; float r26857 = r26847 * r26843; float r26858 = r26856 / r26857; float r26859 = pow(r26858, r26847); float r26860 = r26848 * r26859; float r26861 = r26844 / r26850; float r26862 = r26860 * r26861; float r26863 = r26846 - r26862; float r26864 = r26853 * r26863; return r26864; } double f_id(double d, double h, double l, double M, double D) { double r26865 = d; double r26866 = h; double r26867 = r26865 / r26866; double r26868 = 1; double r26869 = 2; double r26870 = r26868 / r26869; double r26871 = pow(r26867, r26870); double r26872 = l; double r26873 = r26865 / r26872; double r26874 = pow(r26873, r26870); double r26875 = r26871 * r26874; double r26876 = M; double r26877 = D; double r26878 = r26876 * r26877; double r26879 = r26869 * r26865; double r26880 = r26878 / r26879; double r26881 = pow(r26880, r26869); double r26882 = r26870 * r26881; double r26883 = r26866 / r26872; double r26884 = r26882 * r26883; double r26885 = r26868 - r26884; double r26886 = r26875 * r26885; return r26886; } double f_of(float d, float h, float l, float M, float D) { float r26887 = h; float r26888 = -6.21156791415509e-93; bool r26889 = r26887 <= r26888; float r26890 = 1; float r26891 = l; float r26892 = r26890 / r26891; float r26893 = 2; float r26894 = r26890 / r26893; float r26895 = M; float r26896 = D; float r26897 = r26895 * r26896; float r26898 = d; float r26899 = r26898 * r26893; float r26900 = r26897 / r26899; float r26901 = pow(r26900, r26893); float r26902 = r26894 * r26901; float r26903 = r26887 * r26902; float r26904 = r26892 * r26903; float r26905 = r26890 - r26904; float r26906 = r26898 / r26887; float r26907 = pow(r26906, r26894); float r26908 = 1/2; float r26909 = -1; float r26910 = r26909 / r26891; float r26911 = log(r26910); float r26912 = r26909 / r26898; float r26913 = log(r26912); float r26914 = r26911 - r26913; float r26915 = r26908 * r26914; float r26916 = exp(r26915); float r26917 = r26907 * r26916; float r26918 = r26905 * r26917; float r26919 = -5.7840309184133e-311; bool r26920 = r26887 <= r26919; float r26921 = r26909 / r26887; float r26922 = log(r26921); float r26923 = r26922 - r26913; float r26924 = r26923 * r26908; float r26925 = exp(r26924); float r26926 = r26898 / r26891; float r26927 = pow(r26926, r26894); float r26928 = r26925 * r26927; float r26929 = r26928 * r26905; float r26930 = 1.2379471543476946e-36; bool r26931 = r26887 <= r26930; float r26932 = 9.612895863884516e+116; bool r26933 = r26887 <= r26932; float r26934 = !r26933; bool r26935 = r26931 || r26934; float r26936 = r26890 / r26887; float r26937 = pow(r26936, r26894); float r26938 = pow(r26898, r26894); float r26939 = r26937 * r26938; float r26940 = r26927 * r26939; float r26941 = r26940 * r26905; float r26942 = pow(r26892, r26894); float r26943 = r26938 * r26942; float r26944 = r26943 * r26907; float r26945 = r26905 * r26944; float r26946 = r26935 ? r26941 : r26945; float r26947 = r26920 ? r26929 : r26946; float r26948 = r26889 ? r26918 : r26947; return r26948; } double f_od(double d, double h, double l, double M, double D) { double r26949 = h; double r26950 = -6.21156791415509e-93; bool r26951 = r26949 <= r26950; double r26952 = 1; double r26953 = l; double r26954 = r26952 / r26953; double r26955 = 2; double r26956 = r26952 / r26955; double r26957 = M; double r26958 = D; double r26959 = r26957 * r26958; double r26960 = d; double r26961 = r26960 * r26955; double r26962 = r26959 / r26961; double r26963 = pow(r26962, r26955); double r26964 = r26956 * r26963; double r26965 = r26949 * r26964; double r26966 = r26954 * r26965; double r26967 = r26952 - r26966; double r26968 = r26960 / r26949; double r26969 = pow(r26968, r26956); double r26970 = 1/2; double r26971 = -1; double r26972 = r26971 / r26953; double r26973 = log(r26972); double r26974 = r26971 / r26960; double r26975 = log(r26974); double r26976 = r26973 - r26975; double r26977 = r26970 * r26976; double r26978 = exp(r26977); double r26979 = r26969 * r26978; double r26980 = r26967 * r26979; double r26981 = -5.7840309184133e-311; bool r26982 = r26949 <= r26981; double r26983 = r26971 / r26949; double r26984 = log(r26983); double r26985 = r26984 - r26975; double r26986 = r26985 * r26970; double r26987 = exp(r26986); double r26988 = r26960 / r26953; double r26989 = pow(r26988, r26956); double r26990 = r26987 * r26989; double r26991 = r26990 * r26967; double r26992 = 1.2379471543476946e-36; bool r26993 = r26949 <= r26992; double r26994 = 9.612895863884516e+116; bool r26995 = r26949 <= r26994; double r26996 = !r26995; bool r26997 = r26993 || r26996; double r26998 = r26952 / r26949; double r26999 = pow(r26998, r26956); double r27000 = pow(r26960, r26956); double r27001 = r26999 * r27000; double r27002 = r26989 * r27001; double r27003 = r27002 * r26967; double r27004 = pow(r26954, r26956); double r27005 = r27000 * r27004; double r27006 = r27005 * r26969; double r27007 = r26967 * r27006; double r27008 = r26997 ? r27003 : r27007; double r27009 = r26982 ? r26991 : r27008; double r27010 = r26951 ? r26980 : r27009; return r27010; } 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 r27011, r27012, r27013, r27014, r27015, r27016, r27017, r27018, r27019, r27020, r27021, r27022, r27023, r27024, r27025, r27026, r27027, r27028, r27029, r27030, r27031, r27032; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r27011); mpfr_init(r27012); mpfr_init(r27013); mpfr_init_set_str(r27014, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27015, "2", 10, MPFR_RNDN); mpfr_init(r27016); mpfr_init(r27017); mpfr_init(r27018); mpfr_init(r27019); mpfr_init(r27020); mpfr_init(r27021); mpfr_init(r27022); mpfr_init(r27023); mpfr_init(r27024); mpfr_init(r27025); mpfr_init(r27026); mpfr_init(r27027); mpfr_init(r27028); mpfr_init(r27029); mpfr_init(r27030); mpfr_init(r27031); mpfr_init(r27032); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r27011, d, MPFR_RNDN); mpfr_set_d(r27012, h, MPFR_RNDN); mpfr_div(r27013, r27011, r27012, MPFR_RNDN); ; ; mpfr_div(r27016, r27014, r27015, MPFR_RNDN); mpfr_pow(r27017, r27013, r27016, MPFR_RNDN); mpfr_set_d(r27018, l, MPFR_RNDN); mpfr_div(r27019, r27011, r27018, MPFR_RNDN); mpfr_pow(r27020, r27019, r27016, MPFR_RNDN); mpfr_mul(r27021, r27017, r27020, MPFR_RNDN); mpfr_set_d(r27022, M, MPFR_RNDN); mpfr_set_d(r27023, D, MPFR_RNDN); mpfr_mul(r27024, r27022, r27023, MPFR_RNDN); mpfr_mul(r27025, r27015, r27011, MPFR_RNDN); mpfr_div(r27026, r27024, r27025, MPFR_RNDN); mpfr_pow(r27027, r27026, r27015, MPFR_RNDN); mpfr_mul(r27028, r27016, r27027, MPFR_RNDN); mpfr_div(r27029, r27012, r27018, MPFR_RNDN); mpfr_mul(r27030, r27028, r27029, MPFR_RNDN); mpfr_sub(r27031, r27014, r27030, MPFR_RNDN); mpfr_mul(r27032, r27021, r27031, MPFR_RNDN); return mpfr_get_d(r27032, MPFR_RNDN); } static mpfr_t r27033, r27034, r27035, r27036, r27037, r27038, r27039, r27040, r27041, r27042, r27043, r27044, r27045, r27046, r27047, r27048, r27049, r27050, r27051, r27052, r27053, r27054, r27055, r27056, r27057, r27058, r27059, r27060, r27061, r27062, r27063, r27064, r27065, r27066, r27067, r27068, r27069, r27070, r27071, r27072, r27073, r27074, r27075, r27076, r27077, r27078, r27079, r27080, r27081, r27082, r27083, r27084, r27085, r27086, r27087, r27088, r27089, r27090, r27091, r27092, r27093, r27094; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27033); mpfr_init_set_str(r27034, "-6.21156791415509e-93", 10, MPFR_RNDN); mpfr_init(r27035); mpfr_init_set_str(r27036, "1", 10, MPFR_RNDN); mpfr_init(r27037); mpfr_init(r27038); mpfr_init_set_str(r27039, "2", 10, MPFR_RNDN); mpfr_init(r27040); mpfr_init(r27041); mpfr_init(r27042); mpfr_init(r27043); mpfr_init(r27044); mpfr_init(r27045); mpfr_init(r27046); mpfr_init(r27047); mpfr_init(r27048); mpfr_init(r27049); mpfr_init(r27050); mpfr_init(r27051); mpfr_init(r27052); mpfr_init(r27053); mpfr_init_set_str(r27054, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r27055, "-1", 10, MPFR_RNDN); mpfr_init(r27056); mpfr_init(r27057); mpfr_init(r27058); mpfr_init(r27059); mpfr_init(r27060); mpfr_init(r27061); mpfr_init(r27062); mpfr_init(r27063); mpfr_init(r27064); mpfr_init_set_str(r27065, "-5.7840309184133e-311", 10, MPFR_RNDN); mpfr_init(r27066); mpfr_init(r27067); mpfr_init(r27068); mpfr_init(r27069); mpfr_init(r27070); mpfr_init(r27071); mpfr_init(r27072); mpfr_init(r27073); mpfr_init(r27074); mpfr_init(r27075); mpfr_init_set_str(r27076, "1.2379471543476946e-36", 10, MPFR_RNDN); mpfr_init(r27077); mpfr_init_set_str(r27078, "9.612895863884516e+116", 10, MPFR_RNDN); mpfr_init(r27079); mpfr_init(r27080); mpfr_init(r27081); mpfr_init(r27082); mpfr_init(r27083); mpfr_init(r27084); mpfr_init(r27085); mpfr_init(r27086); mpfr_init(r27087); mpfr_init(r27088); mpfr_init(r27089); mpfr_init(r27090); mpfr_init(r27091); mpfr_init(r27092); mpfr_init(r27093); mpfr_init(r27094); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r27033, h, MPFR_RNDN); ; mpfr_set_si(r27035, mpfr_cmp(r27033, r27034) <= 0, MPFR_RNDN); ; mpfr_set_d(r27037, l, MPFR_RNDN); mpfr_div(r27038, r27036, r27037, MPFR_RNDN); ; mpfr_div(r27040, r27036, r27039, MPFR_RNDN); mpfr_set_d(r27041, M, MPFR_RNDN); mpfr_set_d(r27042, D, MPFR_RNDN); mpfr_mul(r27043, r27041, r27042, MPFR_RNDN); mpfr_set_d(r27044, d, MPFR_RNDN); mpfr_mul(r27045, r27044, r27039, MPFR_RNDN); mpfr_div(r27046, r27043, r27045, MPFR_RNDN); mpfr_pow(r27047, r27046, r27039, MPFR_RNDN); mpfr_mul(r27048, r27040, r27047, MPFR_RNDN); mpfr_mul(r27049, r27033, r27048, MPFR_RNDN); mpfr_mul(r27050, r27038, r27049, MPFR_RNDN); mpfr_sub(r27051, r27036, r27050, MPFR_RNDN); mpfr_div(r27052, r27044, r27033, MPFR_RNDN); mpfr_pow(r27053, r27052, r27040, MPFR_RNDN); ; ; mpfr_div(r27056, r27055, r27037, MPFR_RNDN); mpfr_log(r27057, r27056, MPFR_RNDN); mpfr_div(r27058, r27055, r27044, MPFR_RNDN); mpfr_log(r27059, r27058, MPFR_RNDN); mpfr_sub(r27060, r27057, r27059, MPFR_RNDN); mpfr_mul(r27061, r27054, r27060, MPFR_RNDN); mpfr_exp(r27062, r27061, MPFR_RNDN); mpfr_mul(r27063, r27053, r27062, MPFR_RNDN); mpfr_mul(r27064, r27051, r27063, MPFR_RNDN); ; mpfr_set_si(r27066, mpfr_cmp(r27033, r27065) <= 0, MPFR_RNDN); mpfr_div(r27067, r27055, r27033, MPFR_RNDN); mpfr_log(r27068, r27067, MPFR_RNDN); mpfr_sub(r27069, r27068, r27059, MPFR_RNDN); mpfr_mul(r27070, r27069, r27054, MPFR_RNDN); mpfr_exp(r27071, r27070, MPFR_RNDN); mpfr_div(r27072, r27044, r27037, MPFR_RNDN); mpfr_pow(r27073, r27072, r27040, MPFR_RNDN); mpfr_mul(r27074, r27071, r27073, MPFR_RNDN); mpfr_mul(r27075, r27074, r27051, MPFR_RNDN); ; mpfr_set_si(r27077, mpfr_cmp(r27033, r27076) <= 0, MPFR_RNDN); ; mpfr_set_si(r27079, mpfr_cmp(r27033, r27078) <= 0, MPFR_RNDN); mpfr_set_si(r27080, !mpfr_get_si(r27079, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27081, mpfr_get_si(r27077, MPFR_RNDN) || mpfr_get_si(r27080, MPFR_RNDN), MPFR_RNDN); mpfr_div(r27082, r27036, r27033, MPFR_RNDN); mpfr_pow(r27083, r27082, r27040, MPFR_RNDN); mpfr_pow(r27084, r27044, r27040, MPFR_RNDN); mpfr_mul(r27085, r27083, r27084, MPFR_RNDN); mpfr_mul(r27086, r27073, r27085, MPFR_RNDN); mpfr_mul(r27087, r27086, r27051, MPFR_RNDN); mpfr_pow(r27088, r27038, r27040, MPFR_RNDN); mpfr_mul(r27089, r27084, r27088, MPFR_RNDN); mpfr_mul(r27090, r27089, r27053, MPFR_RNDN); mpfr_mul(r27091, r27051, r27090, MPFR_RNDN); if (mpfr_get_si(r27081, MPFR_RNDN)) { mpfr_set(r27092, r27087, MPFR_RNDN); } else { mpfr_set(r27092, r27091, MPFR_RNDN); }; if (mpfr_get_si(r27066, MPFR_RNDN)) { mpfr_set(r27093, r27075, MPFR_RNDN); } else { mpfr_set(r27093, r27092, MPFR_RNDN); }; if (mpfr_get_si(r27035, MPFR_RNDN)) { mpfr_set(r27094, r27064, MPFR_RNDN); } else { mpfr_set(r27094, r27093, MPFR_RNDN); }; return mpfr_get_d(r27094, MPFR_RNDN); } static mpfr_t r27095, r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120, r27121, r27122, r27123, r27124, r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27095); mpfr_init_set_str(r27096, "-6.21156791415509e-93", 10, MPFR_RNDN); mpfr_init(r27097); mpfr_init_set_str(r27098, "1", 10, MPFR_RNDN); mpfr_init(r27099); mpfr_init(r27100); mpfr_init_set_str(r27101, "2", 10, MPFR_RNDN); mpfr_init(r27102); mpfr_init(r27103); mpfr_init(r27104); mpfr_init(r27105); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init(r27109); mpfr_init(r27110); mpfr_init(r27111); mpfr_init(r27112); mpfr_init(r27113); mpfr_init(r27114); mpfr_init(r27115); mpfr_init_set_str(r27116, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r27117, "-1", 10, MPFR_RNDN); mpfr_init(r27118); mpfr_init(r27119); mpfr_init(r27120); mpfr_init(r27121); mpfr_init(r27122); mpfr_init(r27123); mpfr_init(r27124); mpfr_init(r27125); mpfr_init(r27126); mpfr_init_set_str(r27127, "-5.7840309184133e-311", 10, MPFR_RNDN); mpfr_init(r27128); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init(r27135); mpfr_init(r27136); mpfr_init(r27137); mpfr_init_set_str(r27138, "1.2379471543476946e-36", 10, MPFR_RNDN); mpfr_init(r27139); mpfr_init_set_str(r27140, "9.612895863884516e+116", 10, MPFR_RNDN); mpfr_init(r27141); mpfr_init(r27142); mpfr_init(r27143); mpfr_init(r27144); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init(r27152); mpfr_init(r27153); mpfr_init(r27154); mpfr_init(r27155); mpfr_init(r27156); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r27095, h, MPFR_RNDN); ; mpfr_set_si(r27097, mpfr_cmp(r27095, r27096) <= 0, MPFR_RNDN); ; mpfr_set_d(r27099, l, MPFR_RNDN); mpfr_div(r27100, r27098, r27099, MPFR_RNDN); ; mpfr_div(r27102, r27098, r27101, MPFR_RNDN); mpfr_set_d(r27103, M, MPFR_RNDN); mpfr_set_d(r27104, D, MPFR_RNDN); mpfr_mul(r27105, r27103, r27104, MPFR_RNDN); mpfr_set_d(r27106, d, MPFR_RNDN); mpfr_mul(r27107, r27106, r27101, MPFR_RNDN); mpfr_div(r27108, r27105, r27107, MPFR_RNDN); mpfr_pow(r27109, r27108, r27101, MPFR_RNDN); mpfr_mul(r27110, r27102, r27109, MPFR_RNDN); mpfr_mul(r27111, r27095, r27110, MPFR_RNDN); mpfr_mul(r27112, r27100, r27111, MPFR_RNDN); mpfr_sub(r27113, r27098, r27112, MPFR_RNDN); mpfr_div(r27114, r27106, r27095, MPFR_RNDN); mpfr_pow(r27115, r27114, r27102, MPFR_RNDN); ; ; mpfr_div(r27118, r27117, r27099, MPFR_RNDN); mpfr_log(r27119, r27118, MPFR_RNDN); mpfr_div(r27120, r27117, r27106, MPFR_RNDN); mpfr_log(r27121, r27120, MPFR_RNDN); mpfr_sub(r27122, r27119, r27121, MPFR_RNDN); mpfr_mul(r27123, r27116, r27122, MPFR_RNDN); mpfr_exp(r27124, r27123, MPFR_RNDN); mpfr_mul(r27125, r27115, r27124, MPFR_RNDN); mpfr_mul(r27126, r27113, r27125, MPFR_RNDN); ; mpfr_set_si(r27128, mpfr_cmp(r27095, r27127) <= 0, MPFR_RNDN); mpfr_div(r27129, r27117, r27095, MPFR_RNDN); mpfr_log(r27130, r27129, MPFR_RNDN); mpfr_sub(r27131, r27130, r27121, MPFR_RNDN); mpfr_mul(r27132, r27131, r27116, MPFR_RNDN); mpfr_exp(r27133, r27132, MPFR_RNDN); mpfr_div(r27134, r27106, r27099, MPFR_RNDN); mpfr_pow(r27135, r27134, r27102, MPFR_RNDN); mpfr_mul(r27136, r27133, r27135, MPFR_RNDN); mpfr_mul(r27137, r27136, r27113, MPFR_RNDN); ; mpfr_set_si(r27139, mpfr_cmp(r27095, r27138) <= 0, MPFR_RNDN); ; mpfr_set_si(r27141, mpfr_cmp(r27095, r27140) <= 0, MPFR_RNDN); mpfr_set_si(r27142, !mpfr_get_si(r27141, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27143, mpfr_get_si(r27139, MPFR_RNDN) || mpfr_get_si(r27142, MPFR_RNDN), MPFR_RNDN); mpfr_div(r27144, r27098, r27095, MPFR_RNDN); mpfr_pow(r27145, r27144, r27102, MPFR_RNDN); mpfr_pow(r27146, r27106, r27102, MPFR_RNDN); mpfr_mul(r27147, r27145, r27146, MPFR_RNDN); mpfr_mul(r27148, r27135, r27147, MPFR_RNDN); mpfr_mul(r27149, r27148, r27113, MPFR_RNDN); mpfr_pow(r27150, r27100, r27102, MPFR_RNDN); mpfr_mul(r27151, r27146, r27150, MPFR_RNDN); mpfr_mul(r27152, r27151, r27115, MPFR_RNDN); mpfr_mul(r27153, r27113, r27152, MPFR_RNDN); if (mpfr_get_si(r27143, MPFR_RNDN)) { mpfr_set(r27154, r27149, MPFR_RNDN); } else { mpfr_set(r27154, r27153, MPFR_RNDN); }; if (mpfr_get_si(r27128, MPFR_RNDN)) { mpfr_set(r27155, r27137, MPFR_RNDN); } else { mpfr_set(r27155, r27154, MPFR_RNDN); }; if (mpfr_get_si(r27097, MPFR_RNDN)) { mpfr_set(r27156, r27126, MPFR_RNDN); } else { mpfr_set(r27156, r27155, MPFR_RNDN); }; return mpfr_get_d(r27156, MPFR_RNDN); }