#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 r26712 = d; float r26713 = h; float r26714 = r26712 / r26713; float r26715 = 1; float r26716 = 2; float r26717 = r26715 / r26716; float r26718 = pow(r26714, r26717); float r26719 = l; float r26720 = r26712 / r26719; float r26721 = pow(r26720, r26717); float r26722 = r26718 * r26721; float r26723 = M; float r26724 = D; float r26725 = r26723 * r26724; float r26726 = r26716 * r26712; float r26727 = r26725 / r26726; float r26728 = pow(r26727, r26716); float r26729 = r26717 * r26728; float r26730 = r26713 / r26719; float r26731 = r26729 * r26730; float r26732 = r26715 - r26731; float r26733 = r26722 * r26732; return r26733; } double f_id(double d, double h, double l, double M, double D) { double r26734 = d; double r26735 = h; double r26736 = r26734 / r26735; double r26737 = 1; double r26738 = 2; double r26739 = r26737 / r26738; double r26740 = pow(r26736, r26739); double r26741 = l; double r26742 = r26734 / r26741; double r26743 = pow(r26742, r26739); double r26744 = r26740 * r26743; double r26745 = M; double r26746 = D; double r26747 = r26745 * r26746; double r26748 = r26738 * r26734; double r26749 = r26747 / r26748; double r26750 = pow(r26749, r26738); double r26751 = r26739 * r26750; double r26752 = r26735 / r26741; double r26753 = r26751 * r26752; double r26754 = r26737 - r26753; double r26755 = r26744 * r26754; return r26755; } double f_of(float d, float h, float l, float M, float D) { float r26756 = d; float r26757 = -2.1958774546394325e+108; bool r26758 = r26756 <= r26757; float r26759 = 1; float r26760 = h; float r26761 = M; float r26762 = D; float r26763 = r26761 * r26762; float r26764 = 2; float r26765 = r26764 * r26756; float r26766 = r26763 / r26765; float r26767 = pow(r26766, r26764); float r26768 = r26760 * r26767; float r26769 = l; float r26770 = r26769 * r26764; float r26771 = r26768 / r26770; float r26772 = r26759 - r26771; float r26773 = r26756 / r26769; float r26774 = r26759 / r26764; float r26775 = pow(r26773, r26774); float r26776 = -1; float r26777 = r26776 / r26760; float r26778 = log(r26777); float r26779 = r26776 / r26756; float r26780 = log(r26779); float r26781 = r26778 - r26780; float r26782 = 1/2; float r26783 = r26781 * r26782; float r26784 = exp(r26783); float r26785 = r26775 * r26784; float r26786 = r26772 * r26785; float r26787 = -3.3486659915059e-310; bool r26788 = r26756 <= r26787; float r26789 = r26762 / r26756; float r26790 = r26761 / r26764; float r26791 = r26789 * r26790; float r26792 = r26760 * r26791; float r26793 = r26792 * r26791; float r26794 = r26793 / r26770; float r26795 = r26759 - r26794; float r26796 = exp(r26782); float r26797 = r26776 / r26769; float r26798 = log(r26797); float r26799 = r26798 - r26780; float r26800 = pow(r26796, r26799); float r26801 = r26756 / r26760; float r26802 = pow(r26801, r26774); float r26803 = r26800 * r26802; float r26804 = r26795 * r26803; float r26805 = 3.0702186674933227e+77; bool r26806 = r26756 <= r26805; float r26807 = log(r26756); float r26808 = log(r26769); float r26809 = r26807 - r26808; float r26810 = pow(r26796, r26809); float r26811 = r26810 * r26802; float r26812 = r26760 * r26766; float r26813 = r26812 * r26766; float r26814 = r26813 / r26770; float r26815 = r26759 - r26814; float r26816 = r26811 * r26815; float r26817 = r26759 / r26760; float r26818 = pow(r26817, r26774); float r26819 = pow(r26756, r26774); float r26820 = r26818 * r26819; float r26821 = r26775 * r26820; float r26822 = r26774 * r26767; float r26823 = r26760 / r26769; float r26824 = r26822 * r26823; float r26825 = r26759 - r26824; float r26826 = r26821 * r26825; float r26827 = r26806 ? r26816 : r26826; float r26828 = r26788 ? r26804 : r26827; float r26829 = r26758 ? r26786 : r26828; return r26829; } double f_od(double d, double h, double l, double M, double D) { double r26830 = d; double r26831 = -2.1958774546394325e+108; bool r26832 = r26830 <= r26831; double r26833 = 1; double r26834 = h; double r26835 = M; double r26836 = D; double r26837 = r26835 * r26836; double r26838 = 2; double r26839 = r26838 * r26830; double r26840 = r26837 / r26839; double r26841 = pow(r26840, r26838); double r26842 = r26834 * r26841; double r26843 = l; double r26844 = r26843 * r26838; double r26845 = r26842 / r26844; double r26846 = r26833 - r26845; double r26847 = r26830 / r26843; double r26848 = r26833 / r26838; double r26849 = pow(r26847, r26848); double r26850 = -1; double r26851 = r26850 / r26834; double r26852 = log(r26851); double r26853 = r26850 / r26830; double r26854 = log(r26853); double r26855 = r26852 - r26854; double r26856 = 1/2; double r26857 = r26855 * r26856; double r26858 = exp(r26857); double r26859 = r26849 * r26858; double r26860 = r26846 * r26859; double r26861 = -3.3486659915059e-310; bool r26862 = r26830 <= r26861; double r26863 = r26836 / r26830; double r26864 = r26835 / r26838; double r26865 = r26863 * r26864; double r26866 = r26834 * r26865; double r26867 = r26866 * r26865; double r26868 = r26867 / r26844; double r26869 = r26833 - r26868; double r26870 = exp(r26856); double r26871 = r26850 / r26843; double r26872 = log(r26871); double r26873 = r26872 - r26854; double r26874 = pow(r26870, r26873); double r26875 = r26830 / r26834; double r26876 = pow(r26875, r26848); double r26877 = r26874 * r26876; double r26878 = r26869 * r26877; double r26879 = 3.0702186674933227e+77; bool r26880 = r26830 <= r26879; double r26881 = log(r26830); double r26882 = log(r26843); double r26883 = r26881 - r26882; double r26884 = pow(r26870, r26883); double r26885 = r26884 * r26876; double r26886 = r26834 * r26840; double r26887 = r26886 * r26840; double r26888 = r26887 / r26844; double r26889 = r26833 - r26888; double r26890 = r26885 * r26889; double r26891 = r26833 / r26834; double r26892 = pow(r26891, r26848); double r26893 = pow(r26830, r26848); double r26894 = r26892 * r26893; double r26895 = r26849 * r26894; double r26896 = r26848 * r26841; double r26897 = r26834 / r26843; double r26898 = r26896 * r26897; double r26899 = r26833 - r26898; double r26900 = r26895 * r26899; double r26901 = r26880 ? r26890 : r26900; double r26902 = r26862 ? r26878 : r26901; double r26903 = r26832 ? r26860 : r26902; return r26903; } 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 r26904, r26905, r26906, r26907, r26908, r26909, r26910, r26911, r26912, r26913, r26914, r26915, r26916, r26917, r26918, r26919, r26920, r26921, r26922, r26923, r26924, r26925; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26904); mpfr_init(r26905); mpfr_init(r26906); mpfr_init_set_str(r26907, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26908, "2", 10, MPFR_RNDN); mpfr_init(r26909); mpfr_init(r26910); mpfr_init(r26911); mpfr_init(r26912); mpfr_init(r26913); mpfr_init(r26914); mpfr_init(r26915); mpfr_init(r26916); mpfr_init(r26917); mpfr_init(r26918); mpfr_init(r26919); mpfr_init(r26920); mpfr_init(r26921); mpfr_init(r26922); mpfr_init(r26923); mpfr_init(r26924); mpfr_init(r26925); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r26904, d, MPFR_RNDN); mpfr_set_d(r26905, h, MPFR_RNDN); mpfr_div(r26906, r26904, r26905, MPFR_RNDN); ; ; mpfr_div(r26909, r26907, r26908, MPFR_RNDN); mpfr_pow(r26910, r26906, r26909, MPFR_RNDN); mpfr_set_d(r26911, l, MPFR_RNDN); mpfr_div(r26912, r26904, r26911, MPFR_RNDN); mpfr_pow(r26913, r26912, r26909, MPFR_RNDN); mpfr_mul(r26914, r26910, r26913, MPFR_RNDN); mpfr_set_d(r26915, M, MPFR_RNDN); mpfr_set_d(r26916, D, MPFR_RNDN); mpfr_mul(r26917, r26915, r26916, MPFR_RNDN); mpfr_mul(r26918, r26908, r26904, MPFR_RNDN); mpfr_div(r26919, r26917, r26918, MPFR_RNDN); mpfr_pow(r26920, r26919, r26908, MPFR_RNDN); mpfr_mul(r26921, r26909, r26920, MPFR_RNDN); mpfr_div(r26922, r26905, r26911, MPFR_RNDN); mpfr_mul(r26923, r26921, r26922, MPFR_RNDN); mpfr_sub(r26924, r26907, r26923, MPFR_RNDN); mpfr_mul(r26925, r26914, r26924, MPFR_RNDN); return mpfr_get_d(r26925, MPFR_RNDN); } static mpfr_t r26926, r26927, r26928, r26929, r26930, r26931, r26932, r26933, r26934, r26935, r26936, r26937, r26938, r26939, r26940, r26941, r26942, r26943, r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982, r26983, r26984, r26985, r26986, r26987, r26988, r26989, r26990, r26991, r26992, r26993, r26994, r26995, r26996, r26997, r26998, r26999; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26926); mpfr_init_set_str(r26927, "-2.1958774546394325e+108", 10, MPFR_RNDN); mpfr_init(r26928); mpfr_init_set_str(r26929, "1", 10, MPFR_RNDN); mpfr_init(r26930); mpfr_init(r26931); mpfr_init(r26932); mpfr_init(r26933); mpfr_init_set_str(r26934, "2", 10, MPFR_RNDN); mpfr_init(r26935); mpfr_init(r26936); mpfr_init(r26937); mpfr_init(r26938); mpfr_init(r26939); mpfr_init(r26940); mpfr_init(r26941); mpfr_init(r26942); mpfr_init(r26943); mpfr_init(r26944); mpfr_init(r26945); mpfr_init_set_str(r26946, "-1", 10, MPFR_RNDN); mpfr_init(r26947); mpfr_init(r26948); mpfr_init(r26949); mpfr_init(r26950); mpfr_init(r26951); mpfr_init_set_str(r26952, "1/2", 10, MPFR_RNDN); mpfr_init(r26953); mpfr_init(r26954); mpfr_init(r26955); mpfr_init(r26956); mpfr_init_set_str(r26957, "-3.3486659915059e-310", 10, MPFR_RNDN); mpfr_init(r26958); mpfr_init(r26959); mpfr_init(r26960); mpfr_init(r26961); mpfr_init(r26962); mpfr_init(r26963); mpfr_init(r26964); mpfr_init(r26965); mpfr_init(r26966); mpfr_init(r26967); mpfr_init(r26968); mpfr_init(r26969); mpfr_init(r26970); mpfr_init(r26971); mpfr_init(r26972); mpfr_init(r26973); mpfr_init(r26974); mpfr_init_set_str(r26975, "3.0702186674933227e+77", 10, MPFR_RNDN); mpfr_init(r26976); mpfr_init(r26977); mpfr_init(r26978); mpfr_init(r26979); mpfr_init(r26980); mpfr_init(r26981); mpfr_init(r26982); mpfr_init(r26983); mpfr_init(r26984); mpfr_init(r26985); mpfr_init(r26986); mpfr_init(r26987); mpfr_init(r26988); mpfr_init(r26989); mpfr_init(r26990); mpfr_init(r26991); mpfr_init(r26992); mpfr_init(r26993); mpfr_init(r26994); mpfr_init(r26995); mpfr_init(r26996); mpfr_init(r26997); mpfr_init(r26998); mpfr_init(r26999); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r26926, d, MPFR_RNDN); ; mpfr_set_si(r26928, mpfr_cmp(r26926, r26927) <= 0, MPFR_RNDN); ; mpfr_set_d(r26930, h, MPFR_RNDN); mpfr_set_d(r26931, M, MPFR_RNDN); mpfr_set_d(r26932, D, MPFR_RNDN); mpfr_mul(r26933, r26931, r26932, MPFR_RNDN); ; mpfr_mul(r26935, r26934, r26926, MPFR_RNDN); mpfr_div(r26936, r26933, r26935, MPFR_RNDN); mpfr_pow(r26937, r26936, r26934, MPFR_RNDN); mpfr_mul(r26938, r26930, r26937, MPFR_RNDN); mpfr_set_d(r26939, l, MPFR_RNDN); mpfr_mul(r26940, r26939, r26934, MPFR_RNDN); mpfr_div(r26941, r26938, r26940, MPFR_RNDN); mpfr_sub(r26942, r26929, r26941, MPFR_RNDN); mpfr_div(r26943, r26926, r26939, MPFR_RNDN); mpfr_div(r26944, r26929, r26934, MPFR_RNDN); mpfr_pow(r26945, r26943, r26944, MPFR_RNDN); ; mpfr_div(r26947, r26946, r26930, MPFR_RNDN); mpfr_log(r26948, r26947, MPFR_RNDN); mpfr_div(r26949, r26946, r26926, MPFR_RNDN); mpfr_log(r26950, r26949, MPFR_RNDN); mpfr_sub(r26951, r26948, r26950, MPFR_RNDN); ; mpfr_mul(r26953, r26951, r26952, MPFR_RNDN); mpfr_exp(r26954, r26953, MPFR_RNDN); mpfr_mul(r26955, r26945, r26954, MPFR_RNDN); mpfr_mul(r26956, r26942, r26955, MPFR_RNDN); ; mpfr_set_si(r26958, mpfr_cmp(r26926, r26957) <= 0, MPFR_RNDN); mpfr_div(r26959, r26932, r26926, MPFR_RNDN); mpfr_div(r26960, r26931, r26934, MPFR_RNDN); mpfr_mul(r26961, r26959, r26960, MPFR_RNDN); mpfr_mul(r26962, r26930, r26961, MPFR_RNDN); mpfr_mul(r26963, r26962, r26961, MPFR_RNDN); mpfr_div(r26964, r26963, r26940, MPFR_RNDN); mpfr_sub(r26965, r26929, r26964, MPFR_RNDN); mpfr_exp(r26966, r26952, MPFR_RNDN); mpfr_div(r26967, r26946, r26939, MPFR_RNDN); mpfr_log(r26968, r26967, MPFR_RNDN); mpfr_sub(r26969, r26968, r26950, MPFR_RNDN); mpfr_pow(r26970, r26966, r26969, MPFR_RNDN); mpfr_div(r26971, r26926, r26930, MPFR_RNDN); mpfr_pow(r26972, r26971, r26944, MPFR_RNDN); mpfr_mul(r26973, r26970, r26972, MPFR_RNDN); mpfr_mul(r26974, r26965, r26973, MPFR_RNDN); ; mpfr_set_si(r26976, mpfr_cmp(r26926, r26975) <= 0, MPFR_RNDN); mpfr_log(r26977, r26926, MPFR_RNDN); mpfr_log(r26978, r26939, MPFR_RNDN); mpfr_sub(r26979, r26977, r26978, MPFR_RNDN); mpfr_pow(r26980, r26966, r26979, MPFR_RNDN); mpfr_mul(r26981, r26980, r26972, MPFR_RNDN); mpfr_mul(r26982, r26930, r26936, MPFR_RNDN); mpfr_mul(r26983, r26982, r26936, MPFR_RNDN); mpfr_div(r26984, r26983, r26940, MPFR_RNDN); mpfr_sub(r26985, r26929, r26984, MPFR_RNDN); mpfr_mul(r26986, r26981, r26985, MPFR_RNDN); mpfr_div(r26987, r26929, r26930, MPFR_RNDN); mpfr_pow(r26988, r26987, r26944, MPFR_RNDN); mpfr_pow(r26989, r26926, r26944, MPFR_RNDN); mpfr_mul(r26990, r26988, r26989, MPFR_RNDN); mpfr_mul(r26991, r26945, r26990, MPFR_RNDN); mpfr_mul(r26992, r26944, r26937, MPFR_RNDN); mpfr_div(r26993, r26930, r26939, MPFR_RNDN); mpfr_mul(r26994, r26992, r26993, MPFR_RNDN); mpfr_sub(r26995, r26929, r26994, MPFR_RNDN); mpfr_mul(r26996, r26991, r26995, MPFR_RNDN); if (mpfr_get_si(r26976, MPFR_RNDN)) { mpfr_set(r26997, r26986, MPFR_RNDN); } else { mpfr_set(r26997, r26996, MPFR_RNDN); }; if (mpfr_get_si(r26958, MPFR_RNDN)) { mpfr_set(r26998, r26974, MPFR_RNDN); } else { mpfr_set(r26998, r26997, MPFR_RNDN); }; if (mpfr_get_si(r26928, MPFR_RNDN)) { mpfr_set(r26999, r26956, MPFR_RNDN); } else { mpfr_set(r26999, r26998, MPFR_RNDN); }; return mpfr_get_d(r26999, MPFR_RNDN); } static mpfr_t r27000, r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010, r27011, r27012, r27013, r27014, r27015, r27016, r27017, r27018, r27019, r27020, r27021, r27022, r27023, r27024, r27025, r27026, r27027, r27028, r27029, r27030, r27031, r27032, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27000); mpfr_init_set_str(r27001, "-2.1958774546394325e+108", 10, MPFR_RNDN); mpfr_init(r27002); mpfr_init_set_str(r27003, "1", 10, MPFR_RNDN); mpfr_init(r27004); mpfr_init(r27005); mpfr_init(r27006); mpfr_init(r27007); mpfr_init_set_str(r27008, "2", 10, MPFR_RNDN); mpfr_init(r27009); mpfr_init(r27010); mpfr_init(r27011); mpfr_init(r27012); mpfr_init(r27013); mpfr_init(r27014); mpfr_init(r27015); mpfr_init(r27016); mpfr_init(r27017); mpfr_init(r27018); mpfr_init(r27019); mpfr_init_set_str(r27020, "-1", 10, MPFR_RNDN); mpfr_init(r27021); mpfr_init(r27022); mpfr_init(r27023); mpfr_init(r27024); mpfr_init(r27025); mpfr_init_set_str(r27026, "1/2", 10, MPFR_RNDN); mpfr_init(r27027); mpfr_init(r27028); mpfr_init(r27029); mpfr_init(r27030); mpfr_init_set_str(r27031, "-3.3486659915059e-310", 10, MPFR_RNDN); mpfr_init(r27032); mpfr_init(r27033); mpfr_init(r27034); mpfr_init(r27035); mpfr_init(r27036); mpfr_init(r27037); mpfr_init(r27038); mpfr_init(r27039); 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_set_str(r27049, "3.0702186674933227e+77", 10, MPFR_RNDN); mpfr_init(r27050); mpfr_init(r27051); mpfr_init(r27052); mpfr_init(r27053); mpfr_init(r27054); mpfr_init(r27055); 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(r27065); mpfr_init(r27066); mpfr_init(r27067); mpfr_init(r27068); mpfr_init(r27069); mpfr_init(r27070); mpfr_init(r27071); mpfr_init(r27072); mpfr_init(r27073); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r27000, d, MPFR_RNDN); ; mpfr_set_si(r27002, mpfr_cmp(r27000, r27001) <= 0, MPFR_RNDN); ; mpfr_set_d(r27004, h, MPFR_RNDN); mpfr_set_d(r27005, M, MPFR_RNDN); mpfr_set_d(r27006, D, MPFR_RNDN); mpfr_mul(r27007, r27005, r27006, MPFR_RNDN); ; mpfr_mul(r27009, r27008, r27000, MPFR_RNDN); mpfr_div(r27010, r27007, r27009, MPFR_RNDN); mpfr_pow(r27011, r27010, r27008, MPFR_RNDN); mpfr_mul(r27012, r27004, r27011, MPFR_RNDN); mpfr_set_d(r27013, l, MPFR_RNDN); mpfr_mul(r27014, r27013, r27008, MPFR_RNDN); mpfr_div(r27015, r27012, r27014, MPFR_RNDN); mpfr_sub(r27016, r27003, r27015, MPFR_RNDN); mpfr_div(r27017, r27000, r27013, MPFR_RNDN); mpfr_div(r27018, r27003, r27008, MPFR_RNDN); mpfr_pow(r27019, r27017, r27018, MPFR_RNDN); ; mpfr_div(r27021, r27020, r27004, MPFR_RNDN); mpfr_log(r27022, r27021, MPFR_RNDN); mpfr_div(r27023, r27020, r27000, MPFR_RNDN); mpfr_log(r27024, r27023, MPFR_RNDN); mpfr_sub(r27025, r27022, r27024, MPFR_RNDN); ; mpfr_mul(r27027, r27025, r27026, MPFR_RNDN); mpfr_exp(r27028, r27027, MPFR_RNDN); mpfr_mul(r27029, r27019, r27028, MPFR_RNDN); mpfr_mul(r27030, r27016, r27029, MPFR_RNDN); ; mpfr_set_si(r27032, mpfr_cmp(r27000, r27031) <= 0, MPFR_RNDN); mpfr_div(r27033, r27006, r27000, MPFR_RNDN); mpfr_div(r27034, r27005, r27008, MPFR_RNDN); mpfr_mul(r27035, r27033, r27034, MPFR_RNDN); mpfr_mul(r27036, r27004, r27035, MPFR_RNDN); mpfr_mul(r27037, r27036, r27035, MPFR_RNDN); mpfr_div(r27038, r27037, r27014, MPFR_RNDN); mpfr_sub(r27039, r27003, r27038, MPFR_RNDN); mpfr_exp(r27040, r27026, MPFR_RNDN); mpfr_div(r27041, r27020, r27013, MPFR_RNDN); mpfr_log(r27042, r27041, MPFR_RNDN); mpfr_sub(r27043, r27042, r27024, MPFR_RNDN); mpfr_pow(r27044, r27040, r27043, MPFR_RNDN); mpfr_div(r27045, r27000, r27004, MPFR_RNDN); mpfr_pow(r27046, r27045, r27018, MPFR_RNDN); mpfr_mul(r27047, r27044, r27046, MPFR_RNDN); mpfr_mul(r27048, r27039, r27047, MPFR_RNDN); ; mpfr_set_si(r27050, mpfr_cmp(r27000, r27049) <= 0, MPFR_RNDN); mpfr_log(r27051, r27000, MPFR_RNDN); mpfr_log(r27052, r27013, MPFR_RNDN); mpfr_sub(r27053, r27051, r27052, MPFR_RNDN); mpfr_pow(r27054, r27040, r27053, MPFR_RNDN); mpfr_mul(r27055, r27054, r27046, MPFR_RNDN); mpfr_mul(r27056, r27004, r27010, MPFR_RNDN); mpfr_mul(r27057, r27056, r27010, MPFR_RNDN); mpfr_div(r27058, r27057, r27014, MPFR_RNDN); mpfr_sub(r27059, r27003, r27058, MPFR_RNDN); mpfr_mul(r27060, r27055, r27059, MPFR_RNDN); mpfr_div(r27061, r27003, r27004, MPFR_RNDN); mpfr_pow(r27062, r27061, r27018, MPFR_RNDN); mpfr_pow(r27063, r27000, r27018, MPFR_RNDN); mpfr_mul(r27064, r27062, r27063, MPFR_RNDN); mpfr_mul(r27065, r27019, r27064, MPFR_RNDN); mpfr_mul(r27066, r27018, r27011, MPFR_RNDN); mpfr_div(r27067, r27004, r27013, MPFR_RNDN); mpfr_mul(r27068, r27066, r27067, MPFR_RNDN); mpfr_sub(r27069, r27003, r27068, MPFR_RNDN); mpfr_mul(r27070, r27065, r27069, MPFR_RNDN); if (mpfr_get_si(r27050, MPFR_RNDN)) { mpfr_set(r27071, r27060, MPFR_RNDN); } else { mpfr_set(r27071, r27070, MPFR_RNDN); }; if (mpfr_get_si(r27032, MPFR_RNDN)) { mpfr_set(r27072, r27048, MPFR_RNDN); } else { mpfr_set(r27072, r27071, MPFR_RNDN); }; if (mpfr_get_si(r27002, MPFR_RNDN)) { mpfr_set(r27073, r27030, MPFR_RNDN); } else { mpfr_set(r27073, r27072, MPFR_RNDN); }; return mpfr_get_d(r27073, MPFR_RNDN); }