#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 r29811 = d; float r29812 = h; float r29813 = r29811 / r29812; float r29814 = 1.0f; float r29815 = 2.0f; float r29816 = r29814 / r29815; float r29817 = pow(r29813, r29816); float r29818 = l; float r29819 = r29811 / r29818; float r29820 = pow(r29819, r29816); float r29821 = r29817 * r29820; float r29822 = M; float r29823 = D; float r29824 = r29822 * r29823; float r29825 = r29815 * r29811; float r29826 = r29824 / r29825; float r29827 = r29826 * r29826; float r29828 = r29816 * r29827; float r29829 = r29812 / r29818; float r29830 = r29828 * r29829; float r29831 = r29814 - r29830; float r29832 = r29821 * r29831; return r29832; } double f_id(double d, double h, double l, double M, double D) { double r29833 = d; double r29834 = h; double r29835 = r29833 / r29834; double r29836 = 1.0; double r29837 = 2.0; double r29838 = r29836 / r29837; double r29839 = pow(r29835, r29838); double r29840 = l; double r29841 = r29833 / r29840; double r29842 = pow(r29841, r29838); double r29843 = r29839 * r29842; double r29844 = M; double r29845 = D; double r29846 = r29844 * r29845; double r29847 = r29837 * r29833; double r29848 = r29846 / r29847; double r29849 = r29848 * r29848; double r29850 = r29838 * r29849; double r29851 = r29834 / r29840; double r29852 = r29850 * r29851; double r29853 = r29836 - r29852; double r29854 = r29843 * r29853; return r29854; } double f_of(float d, float h, float l, float M, float D) { float r29855 = h; float r29856 = -1.759225207600335e+58f; bool r29857 = r29855 <= r29856; float r29858 = d; float r29859 = r29858 / r29855; float r29860 = 1.0f; float r29861 = 2.0f; float r29862 = r29860 / r29861; float r29863 = pow(r29859, r29862); float r29864 = l; float r29865 = r29858 / r29864; float r29866 = pow(r29865, r29862); float r29867 = r29863 * r29866; float r29868 = r29860 / r29864; float r29869 = D; float r29870 = M; float r29871 = r29869 * r29870; float r29872 = r29858 + r29858; float r29873 = r29871 / r29872; float r29874 = r29873 * r29873; float r29875 = r29861 / r29855; float r29876 = r29874 / r29875; float r29877 = r29868 * r29876; float r29878 = r29860 - r29877; float r29879 = r29867 * r29878; float r29880 = -1.3193602092139957e-197f; bool r29881 = r29855 <= r29880; float r29882 = sqrt(r29866); float r29883 = r29882 * r29882; float r29884 = r29863 * r29883; float r29885 = r29870 / r29860; float r29886 = r29869 / r29872; float r29887 = r29885 * r29886; float r29888 = r29855 / r29861; float r29889 = r29864 / r29888; float r29890 = sqrt(r29889); float r29891 = r29887 / r29890; float r29892 = r29891 * r29891; float r29893 = r29860 - r29892; float r29894 = r29884 * r29893; float r29895 = 7.060529603865913e-283f; bool r29896 = r29855 <= r29895; float r29897 = r29870 * r29869; float r29898 = r29897 / r29872; float r29899 = r29898 / r29864; float r29900 = r29870 * r29855; float r29901 = r29900 / r29872; float r29902 = r29869 / r29861; float r29903 = r29901 * r29902; float r29904 = r29899 * r29903; float r29905 = r29860 - r29904; float r29906 = r29867 * r29905; float r29907 = 8.953263541046556e-147f; bool r29908 = r29855 <= r29907; float r29909 = pow(r29858, r29862); float r29910 = r29860 / r29855; float r29911 = pow(r29910, r29862); float r29912 = r29909 * r29911; float r29913 = r29912 * r29866; float r29914 = r29913 * r29893; float r29915 = pow(r29868, r29862); float r29916 = r29909 * r29915; float r29917 = r29863 * r29916; float r29918 = r29898 / r29890; float r29919 = r29918 * r29918; float r29920 = r29860 - r29919; float r29921 = r29917 * r29920; float r29922 = r29908 ? r29914 : r29921; float r29923 = r29896 ? r29906 : r29922; float r29924 = r29881 ? r29894 : r29923; float r29925 = r29857 ? r29879 : r29924; return r29925; } double f_od(double d, double h, double l, double M, double D) { double r29926 = h; double r29927 = -1.759225207600335e+58; bool r29928 = r29926 <= r29927; double r29929 = d; double r29930 = r29929 / r29926; double r29931 = 1.0; double r29932 = 2.0; double r29933 = r29931 / r29932; double r29934 = pow(r29930, r29933); double r29935 = l; double r29936 = r29929 / r29935; double r29937 = pow(r29936, r29933); double r29938 = r29934 * r29937; double r29939 = r29931 / r29935; double r29940 = D; double r29941 = M; double r29942 = r29940 * r29941; double r29943 = r29929 + r29929; double r29944 = r29942 / r29943; double r29945 = r29944 * r29944; double r29946 = r29932 / r29926; double r29947 = r29945 / r29946; double r29948 = r29939 * r29947; double r29949 = r29931 - r29948; double r29950 = r29938 * r29949; double r29951 = -1.3193602092139957e-197; bool r29952 = r29926 <= r29951; double r29953 = sqrt(r29937); double r29954 = r29953 * r29953; double r29955 = r29934 * r29954; double r29956 = r29941 / r29931; double r29957 = r29940 / r29943; double r29958 = r29956 * r29957; double r29959 = r29926 / r29932; double r29960 = r29935 / r29959; double r29961 = sqrt(r29960); double r29962 = r29958 / r29961; double r29963 = r29962 * r29962; double r29964 = r29931 - r29963; double r29965 = r29955 * r29964; double r29966 = 7.060529603865913e-283; bool r29967 = r29926 <= r29966; double r29968 = r29941 * r29940; double r29969 = r29968 / r29943; double r29970 = r29969 / r29935; double r29971 = r29941 * r29926; double r29972 = r29971 / r29943; double r29973 = r29940 / r29932; double r29974 = r29972 * r29973; double r29975 = r29970 * r29974; double r29976 = r29931 - r29975; double r29977 = r29938 * r29976; double r29978 = 8.953263541046556e-147; bool r29979 = r29926 <= r29978; double r29980 = pow(r29929, r29933); double r29981 = r29931 / r29926; double r29982 = pow(r29981, r29933); double r29983 = r29980 * r29982; double r29984 = r29983 * r29937; double r29985 = r29984 * r29964; double r29986 = pow(r29939, r29933); double r29987 = r29980 * r29986; double r29988 = r29934 * r29987; double r29989 = r29969 / r29961; double r29990 = r29989 * r29989; double r29991 = r29931 - r29990; double r29992 = r29988 * r29991; double r29993 = r29979 ? r29985 : r29992; double r29994 = r29967 ? r29977 : r29993; double r29995 = r29952 ? r29965 : r29994; double r29996 = r29928 ? r29950 : r29995; return r29996; } 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 r29997, r29998, r29999, r30000, r30001, r30002, r30003, r30004, r30005, r30006, r30007, r30008, r30009, r30010, r30011, r30012, r30013, r30014, r30015, r30016, r30017, r30018; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r29997); mpfr_init(r29998); mpfr_init(r29999); mpfr_init_set_str(r30000, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30001, "2", 10, MPFR_RNDN); mpfr_init(r30002); mpfr_init(r30003); mpfr_init(r30004); mpfr_init(r30005); mpfr_init(r30006); mpfr_init(r30007); mpfr_init(r30008); mpfr_init(r30009); mpfr_init(r30010); mpfr_init(r30011); mpfr_init(r30012); mpfr_init(r30013); mpfr_init(r30014); mpfr_init(r30015); mpfr_init(r30016); mpfr_init(r30017); mpfr_init(r30018); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r29997, d, MPFR_RNDN); mpfr_set_d(r29998, h, MPFR_RNDN); mpfr_div(r29999, r29997, r29998, MPFR_RNDN); ; ; mpfr_div(r30002, r30000, r30001, MPFR_RNDN); mpfr_pow(r30003, r29999, r30002, MPFR_RNDN); mpfr_set_d(r30004, l, MPFR_RNDN); mpfr_div(r30005, r29997, r30004, MPFR_RNDN); mpfr_pow(r30006, r30005, r30002, MPFR_RNDN); mpfr_mul(r30007, r30003, r30006, MPFR_RNDN); mpfr_set_d(r30008, M, MPFR_RNDN); mpfr_set_d(r30009, D, MPFR_RNDN); mpfr_mul(r30010, r30008, r30009, MPFR_RNDN); mpfr_mul(r30011, r30001, r29997, MPFR_RNDN); mpfr_div(r30012, r30010, r30011, MPFR_RNDN); mpfr_sqr(r30013, r30012, MPFR_RNDN); mpfr_mul(r30014, r30002, r30013, MPFR_RNDN); mpfr_div(r30015, r29998, r30004, MPFR_RNDN); mpfr_mul(r30016, r30014, r30015, MPFR_RNDN); mpfr_sub(r30017, r30000, r30016, MPFR_RNDN); mpfr_mul(r30018, r30007, r30017, MPFR_RNDN); return mpfr_get_d(r30018, MPFR_RNDN); } static mpfr_t r30019, r30020, r30021, r30022, r30023, r30024, r30025, r30026, r30027, r30028, r30029, r30030, r30031, r30032, r30033, r30034, r30035, r30036, r30037, r30038, r30039, r30040, r30041, r30042, r30043, r30044, r30045, r30046, r30047, r30048, r30049, r30050, r30051, r30052, r30053, r30054, r30055, r30056, r30057, r30058, r30059, r30060, r30061, r30062, r30063, r30064, r30065, r30066, r30067, r30068, r30069, r30070, r30071, r30072, r30073, r30074, r30075, r30076, r30077, r30078, r30079, r30080, r30081, r30082, r30083, r30084, r30085, r30086, r30087, r30088, r30089; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r30019); mpfr_init_set_str(r30020, "-1.759225207600335e+58", 10, MPFR_RNDN); mpfr_init(r30021); mpfr_init(r30022); mpfr_init(r30023); mpfr_init_set_str(r30024, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30025, "2", 10, MPFR_RNDN); mpfr_init(r30026); mpfr_init(r30027); mpfr_init(r30028); mpfr_init(r30029); mpfr_init(r30030); mpfr_init(r30031); mpfr_init(r30032); mpfr_init(r30033); mpfr_init(r30034); mpfr_init(r30035); mpfr_init(r30036); mpfr_init(r30037); mpfr_init(r30038); mpfr_init(r30039); mpfr_init(r30040); mpfr_init(r30041); mpfr_init(r30042); mpfr_init(r30043); mpfr_init_set_str(r30044, "-1.3193602092139957e-197", 10, MPFR_RNDN); mpfr_init(r30045); mpfr_init(r30046); mpfr_init(r30047); mpfr_init(r30048); mpfr_init(r30049); mpfr_init(r30050); mpfr_init(r30051); mpfr_init(r30052); mpfr_init(r30053); mpfr_init(r30054); mpfr_init(r30055); mpfr_init(r30056); mpfr_init(r30057); mpfr_init(r30058); mpfr_init_set_str(r30059, "7.060529603865913e-283", 10, MPFR_RNDN); mpfr_init(r30060); mpfr_init(r30061); mpfr_init(r30062); mpfr_init(r30063); mpfr_init(r30064); mpfr_init(r30065); mpfr_init(r30066); mpfr_init(r30067); mpfr_init(r30068); mpfr_init(r30069); mpfr_init(r30070); mpfr_init_set_str(r30071, "8.953263541046556e-147", 10, MPFR_RNDN); mpfr_init(r30072); mpfr_init(r30073); mpfr_init(r30074); mpfr_init(r30075); mpfr_init(r30076); mpfr_init(r30077); mpfr_init(r30078); mpfr_init(r30079); mpfr_init(r30080); mpfr_init(r30081); mpfr_init(r30082); mpfr_init(r30083); mpfr_init(r30084); mpfr_init(r30085); mpfr_init(r30086); mpfr_init(r30087); mpfr_init(r30088); mpfr_init(r30089); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r30019, h, MPFR_RNDN); ; mpfr_set_si(r30021, mpfr_cmp(r30019, r30020) <= 0, MPFR_RNDN); mpfr_set_d(r30022, d, MPFR_RNDN); mpfr_div(r30023, r30022, r30019, MPFR_RNDN); ; ; mpfr_div(r30026, r30024, r30025, MPFR_RNDN); mpfr_pow(r30027, r30023, r30026, MPFR_RNDN); mpfr_set_d(r30028, l, MPFR_RNDN); mpfr_div(r30029, r30022, r30028, MPFR_RNDN); mpfr_pow(r30030, r30029, r30026, MPFR_RNDN); mpfr_mul(r30031, r30027, r30030, MPFR_RNDN); mpfr_div(r30032, r30024, r30028, MPFR_RNDN); mpfr_set_d(r30033, D, MPFR_RNDN); mpfr_set_d(r30034, M, MPFR_RNDN); mpfr_mul(r30035, r30033, r30034, MPFR_RNDN); mpfr_add(r30036, r30022, r30022, MPFR_RNDN); mpfr_div(r30037, r30035, r30036, MPFR_RNDN); mpfr_sqr(r30038, r30037, MPFR_RNDN); mpfr_div(r30039, r30025, r30019, MPFR_RNDN); mpfr_div(r30040, r30038, r30039, MPFR_RNDN); mpfr_mul(r30041, r30032, r30040, MPFR_RNDN); mpfr_sub(r30042, r30024, r30041, MPFR_RNDN); mpfr_mul(r30043, r30031, r30042, MPFR_RNDN); ; mpfr_set_si(r30045, mpfr_cmp(r30019, r30044) <= 0, MPFR_RNDN); mpfr_sqrt(r30046, r30030, MPFR_RNDN); mpfr_sqr(r30047, r30046, MPFR_RNDN); mpfr_mul(r30048, r30027, r30047, MPFR_RNDN); mpfr_div(r30049, r30034, r30024, MPFR_RNDN); mpfr_div(r30050, r30033, r30036, MPFR_RNDN); mpfr_mul(r30051, r30049, r30050, MPFR_RNDN); mpfr_div(r30052, r30019, r30025, MPFR_RNDN); mpfr_div(r30053, r30028, r30052, MPFR_RNDN); mpfr_sqrt(r30054, r30053, MPFR_RNDN); mpfr_div(r30055, r30051, r30054, MPFR_RNDN); mpfr_sqr(r30056, r30055, MPFR_RNDN); mpfr_sub(r30057, r30024, r30056, MPFR_RNDN); mpfr_mul(r30058, r30048, r30057, MPFR_RNDN); ; mpfr_set_si(r30060, mpfr_cmp(r30019, r30059) <= 0, MPFR_RNDN); mpfr_mul(r30061, r30034, r30033, MPFR_RNDN); mpfr_div(r30062, r30061, r30036, MPFR_RNDN); mpfr_div(r30063, r30062, r30028, MPFR_RNDN); mpfr_mul(r30064, r30034, r30019, MPFR_RNDN); mpfr_div(r30065, r30064, r30036, MPFR_RNDN); mpfr_div(r30066, r30033, r30025, MPFR_RNDN); mpfr_mul(r30067, r30065, r30066, MPFR_RNDN); mpfr_mul(r30068, r30063, r30067, MPFR_RNDN); mpfr_sub(r30069, r30024, r30068, MPFR_RNDN); mpfr_mul(r30070, r30031, r30069, MPFR_RNDN); ; mpfr_set_si(r30072, mpfr_cmp(r30019, r30071) <= 0, MPFR_RNDN); mpfr_pow(r30073, r30022, r30026, MPFR_RNDN); mpfr_div(r30074, r30024, r30019, MPFR_RNDN); mpfr_pow(r30075, r30074, r30026, MPFR_RNDN); mpfr_mul(r30076, r30073, r30075, MPFR_RNDN); mpfr_mul(r30077, r30076, r30030, MPFR_RNDN); mpfr_mul(r30078, r30077, r30057, MPFR_RNDN); mpfr_pow(r30079, r30032, r30026, MPFR_RNDN); mpfr_mul(r30080, r30073, r30079, MPFR_RNDN); mpfr_mul(r30081, r30027, r30080, MPFR_RNDN); mpfr_div(r30082, r30062, r30054, MPFR_RNDN); mpfr_sqr(r30083, r30082, MPFR_RNDN); mpfr_sub(r30084, r30024, r30083, MPFR_RNDN); mpfr_mul(r30085, r30081, r30084, MPFR_RNDN); if (mpfr_get_si(r30072, MPFR_RNDN)) { mpfr_set(r30086, r30078, MPFR_RNDN); } else { mpfr_set(r30086, r30085, MPFR_RNDN); }; if (mpfr_get_si(r30060, MPFR_RNDN)) { mpfr_set(r30087, r30070, MPFR_RNDN); } else { mpfr_set(r30087, r30086, MPFR_RNDN); }; if (mpfr_get_si(r30045, MPFR_RNDN)) { mpfr_set(r30088, r30058, MPFR_RNDN); } else { mpfr_set(r30088, r30087, MPFR_RNDN); }; if (mpfr_get_si(r30021, MPFR_RNDN)) { mpfr_set(r30089, r30043, MPFR_RNDN); } else { mpfr_set(r30089, r30088, MPFR_RNDN); }; return mpfr_get_d(r30089, MPFR_RNDN); } static mpfr_t r30090, r30091, r30092, r30093, r30094, r30095, r30096, r30097, r30098, r30099, r30100, r30101, r30102, r30103, r30104, r30105, r30106, r30107, r30108, r30109, r30110, r30111, r30112, r30113, r30114, r30115, r30116, r30117, r30118, r30119, r30120, r30121, r30122, r30123, r30124, r30125, r30126, r30127, r30128, r30129, r30130, r30131, r30132, r30133, r30134, r30135, r30136, r30137, r30138, r30139, r30140, r30141, r30142, r30143, r30144, r30145, r30146, r30147, r30148, r30149, r30150, r30151, r30152, r30153, r30154, r30155, r30156, r30157, r30158, r30159, r30160; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r30090); mpfr_init_set_str(r30091, "-1.759225207600335e+58", 10, MPFR_RNDN); mpfr_init(r30092); mpfr_init(r30093); mpfr_init(r30094); mpfr_init_set_str(r30095, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30096, "2", 10, MPFR_RNDN); mpfr_init(r30097); mpfr_init(r30098); mpfr_init(r30099); mpfr_init(r30100); mpfr_init(r30101); mpfr_init(r30102); mpfr_init(r30103); mpfr_init(r30104); mpfr_init(r30105); mpfr_init(r30106); mpfr_init(r30107); mpfr_init(r30108); mpfr_init(r30109); mpfr_init(r30110); mpfr_init(r30111); mpfr_init(r30112); mpfr_init(r30113); mpfr_init(r30114); mpfr_init_set_str(r30115, "-1.3193602092139957e-197", 10, MPFR_RNDN); mpfr_init(r30116); mpfr_init(r30117); mpfr_init(r30118); mpfr_init(r30119); mpfr_init(r30120); mpfr_init(r30121); mpfr_init(r30122); mpfr_init(r30123); mpfr_init(r30124); mpfr_init(r30125); mpfr_init(r30126); mpfr_init(r30127); mpfr_init(r30128); mpfr_init(r30129); mpfr_init_set_str(r30130, "7.060529603865913e-283", 10, MPFR_RNDN); mpfr_init(r30131); mpfr_init(r30132); mpfr_init(r30133); mpfr_init(r30134); mpfr_init(r30135); mpfr_init(r30136); mpfr_init(r30137); mpfr_init(r30138); mpfr_init(r30139); mpfr_init(r30140); mpfr_init(r30141); mpfr_init_set_str(r30142, "8.953263541046556e-147", 10, MPFR_RNDN); mpfr_init(r30143); mpfr_init(r30144); mpfr_init(r30145); mpfr_init(r30146); mpfr_init(r30147); mpfr_init(r30148); mpfr_init(r30149); mpfr_init(r30150); mpfr_init(r30151); mpfr_init(r30152); mpfr_init(r30153); mpfr_init(r30154); mpfr_init(r30155); mpfr_init(r30156); mpfr_init(r30157); mpfr_init(r30158); mpfr_init(r30159); mpfr_init(r30160); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r30090, h, MPFR_RNDN); ; mpfr_set_si(r30092, mpfr_cmp(r30090, r30091) <= 0, MPFR_RNDN); mpfr_set_d(r30093, d, MPFR_RNDN); mpfr_div(r30094, r30093, r30090, MPFR_RNDN); ; ; mpfr_div(r30097, r30095, r30096, MPFR_RNDN); mpfr_pow(r30098, r30094, r30097, MPFR_RNDN); mpfr_set_d(r30099, l, MPFR_RNDN); mpfr_div(r30100, r30093, r30099, MPFR_RNDN); mpfr_pow(r30101, r30100, r30097, MPFR_RNDN); mpfr_mul(r30102, r30098, r30101, MPFR_RNDN); mpfr_div(r30103, r30095, r30099, MPFR_RNDN); mpfr_set_d(r30104, D, MPFR_RNDN); mpfr_set_d(r30105, M, MPFR_RNDN); mpfr_mul(r30106, r30104, r30105, MPFR_RNDN); mpfr_add(r30107, r30093, r30093, MPFR_RNDN); mpfr_div(r30108, r30106, r30107, MPFR_RNDN); mpfr_sqr(r30109, r30108, MPFR_RNDN); mpfr_div(r30110, r30096, r30090, MPFR_RNDN); mpfr_div(r30111, r30109, r30110, MPFR_RNDN); mpfr_mul(r30112, r30103, r30111, MPFR_RNDN); mpfr_sub(r30113, r30095, r30112, MPFR_RNDN); mpfr_mul(r30114, r30102, r30113, MPFR_RNDN); ; mpfr_set_si(r30116, mpfr_cmp(r30090, r30115) <= 0, MPFR_RNDN); mpfr_sqrt(r30117, r30101, MPFR_RNDN); mpfr_sqr(r30118, r30117, MPFR_RNDN); mpfr_mul(r30119, r30098, r30118, MPFR_RNDN); mpfr_div(r30120, r30105, r30095, MPFR_RNDN); mpfr_div(r30121, r30104, r30107, MPFR_RNDN); mpfr_mul(r30122, r30120, r30121, MPFR_RNDN); mpfr_div(r30123, r30090, r30096, MPFR_RNDN); mpfr_div(r30124, r30099, r30123, MPFR_RNDN); mpfr_sqrt(r30125, r30124, MPFR_RNDN); mpfr_div(r30126, r30122, r30125, MPFR_RNDN); mpfr_sqr(r30127, r30126, MPFR_RNDN); mpfr_sub(r30128, r30095, r30127, MPFR_RNDN); mpfr_mul(r30129, r30119, r30128, MPFR_RNDN); ; mpfr_set_si(r30131, mpfr_cmp(r30090, r30130) <= 0, MPFR_RNDN); mpfr_mul(r30132, r30105, r30104, MPFR_RNDN); mpfr_div(r30133, r30132, r30107, MPFR_RNDN); mpfr_div(r30134, r30133, r30099, MPFR_RNDN); mpfr_mul(r30135, r30105, r30090, MPFR_RNDN); mpfr_div(r30136, r30135, r30107, MPFR_RNDN); mpfr_div(r30137, r30104, r30096, MPFR_RNDN); mpfr_mul(r30138, r30136, r30137, MPFR_RNDN); mpfr_mul(r30139, r30134, r30138, MPFR_RNDN); mpfr_sub(r30140, r30095, r30139, MPFR_RNDN); mpfr_mul(r30141, r30102, r30140, MPFR_RNDN); ; mpfr_set_si(r30143, mpfr_cmp(r30090, r30142) <= 0, MPFR_RNDN); mpfr_pow(r30144, r30093, r30097, MPFR_RNDN); mpfr_div(r30145, r30095, r30090, MPFR_RNDN); mpfr_pow(r30146, r30145, r30097, MPFR_RNDN); mpfr_mul(r30147, r30144, r30146, MPFR_RNDN); mpfr_mul(r30148, r30147, r30101, MPFR_RNDN); mpfr_mul(r30149, r30148, r30128, MPFR_RNDN); mpfr_pow(r30150, r30103, r30097, MPFR_RNDN); mpfr_mul(r30151, r30144, r30150, MPFR_RNDN); mpfr_mul(r30152, r30098, r30151, MPFR_RNDN); mpfr_div(r30153, r30133, r30125, MPFR_RNDN); mpfr_sqr(r30154, r30153, MPFR_RNDN); mpfr_sub(r30155, r30095, r30154, MPFR_RNDN); mpfr_mul(r30156, r30152, r30155, MPFR_RNDN); if (mpfr_get_si(r30143, MPFR_RNDN)) { mpfr_set(r30157, r30149, MPFR_RNDN); } else { mpfr_set(r30157, r30156, MPFR_RNDN); }; if (mpfr_get_si(r30131, MPFR_RNDN)) { mpfr_set(r30158, r30141, MPFR_RNDN); } else { mpfr_set(r30158, r30157, MPFR_RNDN); }; if (mpfr_get_si(r30116, MPFR_RNDN)) { mpfr_set(r30159, r30129, MPFR_RNDN); } else { mpfr_set(r30159, r30158, MPFR_RNDN); }; if (mpfr_get_si(r30092, MPFR_RNDN)) { mpfr_set(r30160, r30114, MPFR_RNDN); } else { mpfr_set(r30160, r30159, MPFR_RNDN); }; return mpfr_get_d(r30160, MPFR_RNDN); }