#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 r10779 = d; float r10780 = h; float r10781 = r10779 / r10780; float r10782 = 1; float r10783 = 2; float r10784 = r10782 / r10783; float r10785 = pow(r10781, r10784); float r10786 = l; float r10787 = r10779 / r10786; float r10788 = pow(r10787, r10784); float r10789 = r10785 * r10788; float r10790 = M; float r10791 = D; float r10792 = r10790 * r10791; float r10793 = r10783 * r10779; float r10794 = r10792 / r10793; float r10795 = pow(r10794, r10783); float r10796 = r10784 * r10795; float r10797 = r10780 / r10786; float r10798 = r10796 * r10797; float r10799 = r10782 - r10798; float r10800 = r10789 * r10799; return r10800; } double f_id(double d, double h, double l, double M, double D) { double r10801 = d; double r10802 = h; double r10803 = r10801 / r10802; double r10804 = 1; double r10805 = 2; double r10806 = r10804 / r10805; double r10807 = pow(r10803, r10806); double r10808 = l; double r10809 = r10801 / r10808; double r10810 = pow(r10809, r10806); double r10811 = r10807 * r10810; double r10812 = M; double r10813 = D; double r10814 = r10812 * r10813; double r10815 = r10805 * r10801; double r10816 = r10814 / r10815; double r10817 = pow(r10816, r10805); double r10818 = r10806 * r10817; double r10819 = r10802 / r10808; double r10820 = r10818 * r10819; double r10821 = r10804 - r10820; double r10822 = r10811 * r10821; return r10822; } double f_of(float d, float h, float l, float M, float D) { float r10823 = h; float r10824 = -7.243473966147417e-108; bool r10825 = r10823 <= r10824; float r10826 = 1; float r10827 = M; float r10828 = D; float r10829 = r10827 * r10828; float r10830 = d; float r10831 = 2; float r10832 = r10830 * r10831; float r10833 = r10829 / r10832; float r10834 = pow(r10833, r10831); float r10835 = r10823 * r10834; float r10836 = l; float r10837 = r10836 * r10831; float r10838 = r10835 / r10837; float r10839 = r10826 - r10838; float r10840 = 1/2; float r10841 = -1; float r10842 = r10841 / r10836; float r10843 = log(r10842); float r10844 = r10841 / r10830; float r10845 = log(r10844); float r10846 = r10843 - r10845; float r10847 = r10840 * r10846; float r10848 = exp(r10847); float r10849 = r10830 / r10823; float r10850 = r10826 / r10831; float r10851 = pow(r10849, r10850); float r10852 = r10848 * r10851; float r10853 = r10839 * r10852; float r10854 = 5.02327006811133e-310; bool r10855 = r10823 <= r10854; float r10856 = r10828 / r10831; float r10857 = r10827 / r10830; float r10858 = r10856 * r10857; float r10859 = r10823 * r10858; float r10860 = r10859 * r10858; float r10861 = r10860 / r10837; float r10862 = r10826 - r10861; float r10863 = exp(r10840); float r10864 = r10841 / r10823; float r10865 = log(r10864); float r10866 = r10865 - r10845; float r10867 = pow(r10863, r10866); float r10868 = r10830 / r10836; float r10869 = pow(r10868, r10850); float r10870 = r10867 * r10869; float r10871 = r10862 * r10870; float r10872 = 5.704880922779612e-67; bool r10873 = r10823 <= r10872; float r10874 = r10826 / r10823; float r10875 = pow(r10874, r10850); float r10876 = pow(r10830, r10850); float r10877 = r10875 * r10876; float r10878 = r10877 * r10869; float r10879 = r10878 * r10839; float r10880 = r10826 / r10836; float r10881 = pow(r10880, r10850); float r10882 = r10876 * r10881; float r10883 = r10882 * r10851; float r10884 = r10883 * r10839; float r10885 = r10873 ? r10879 : r10884; float r10886 = r10855 ? r10871 : r10885; float r10887 = r10825 ? r10853 : r10886; return r10887; } double f_od(double d, double h, double l, double M, double D) { double r10888 = h; double r10889 = -7.243473966147417e-108; bool r10890 = r10888 <= r10889; double r10891 = 1; double r10892 = M; double r10893 = D; double r10894 = r10892 * r10893; double r10895 = d; double r10896 = 2; double r10897 = r10895 * r10896; double r10898 = r10894 / r10897; double r10899 = pow(r10898, r10896); double r10900 = r10888 * r10899; double r10901 = l; double r10902 = r10901 * r10896; double r10903 = r10900 / r10902; double r10904 = r10891 - r10903; double r10905 = 1/2; double r10906 = -1; double r10907 = r10906 / r10901; double r10908 = log(r10907); double r10909 = r10906 / r10895; double r10910 = log(r10909); double r10911 = r10908 - r10910; double r10912 = r10905 * r10911; double r10913 = exp(r10912); double r10914 = r10895 / r10888; double r10915 = r10891 / r10896; double r10916 = pow(r10914, r10915); double r10917 = r10913 * r10916; double r10918 = r10904 * r10917; double r10919 = 5.02327006811133e-310; bool r10920 = r10888 <= r10919; double r10921 = r10893 / r10896; double r10922 = r10892 / r10895; double r10923 = r10921 * r10922; double r10924 = r10888 * r10923; double r10925 = r10924 * r10923; double r10926 = r10925 / r10902; double r10927 = r10891 - r10926; double r10928 = exp(r10905); double r10929 = r10906 / r10888; double r10930 = log(r10929); double r10931 = r10930 - r10910; double r10932 = pow(r10928, r10931); double r10933 = r10895 / r10901; double r10934 = pow(r10933, r10915); double r10935 = r10932 * r10934; double r10936 = r10927 * r10935; double r10937 = 5.704880922779612e-67; bool r10938 = r10888 <= r10937; double r10939 = r10891 / r10888; double r10940 = pow(r10939, r10915); double r10941 = pow(r10895, r10915); double r10942 = r10940 * r10941; double r10943 = r10942 * r10934; double r10944 = r10943 * r10904; double r10945 = r10891 / r10901; double r10946 = pow(r10945, r10915); double r10947 = r10941 * r10946; double r10948 = r10947 * r10916; double r10949 = r10948 * r10904; double r10950 = r10938 ? r10944 : r10949; double r10951 = r10920 ? r10936 : r10950; double r10952 = r10890 ? r10918 : r10951; return r10952; } 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 r10953, r10954, r10955, r10956, r10957, r10958, r10959, r10960, r10961, r10962, r10963, r10964, r10965, r10966, r10967, r10968, r10969, r10970, r10971, r10972, r10973, r10974; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10953); mpfr_init(r10954); mpfr_init(r10955); mpfr_init_set_str(r10956, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10957, "2", 10, MPFR_RNDN); mpfr_init(r10958); mpfr_init(r10959); mpfr_init(r10960); mpfr_init(r10961); mpfr_init(r10962); mpfr_init(r10963); mpfr_init(r10964); mpfr_init(r10965); mpfr_init(r10966); mpfr_init(r10967); mpfr_init(r10968); mpfr_init(r10969); mpfr_init(r10970); mpfr_init(r10971); mpfr_init(r10972); mpfr_init(r10973); mpfr_init(r10974); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r10953, d, MPFR_RNDN); mpfr_set_d(r10954, h, MPFR_RNDN); mpfr_div(r10955, r10953, r10954, MPFR_RNDN); ; ; mpfr_div(r10958, r10956, r10957, MPFR_RNDN); mpfr_pow(r10959, r10955, r10958, MPFR_RNDN); mpfr_set_d(r10960, l, MPFR_RNDN); mpfr_div(r10961, r10953, r10960, MPFR_RNDN); mpfr_pow(r10962, r10961, r10958, MPFR_RNDN); mpfr_mul(r10963, r10959, r10962, MPFR_RNDN); mpfr_set_d(r10964, M, MPFR_RNDN); mpfr_set_d(r10965, D, MPFR_RNDN); mpfr_mul(r10966, r10964, r10965, MPFR_RNDN); mpfr_mul(r10967, r10957, r10953, MPFR_RNDN); mpfr_div(r10968, r10966, r10967, MPFR_RNDN); mpfr_pow(r10969, r10968, r10957, MPFR_RNDN); mpfr_mul(r10970, r10958, r10969, MPFR_RNDN); mpfr_div(r10971, r10954, r10960, MPFR_RNDN); mpfr_mul(r10972, r10970, r10971, MPFR_RNDN); mpfr_sub(r10973, r10956, r10972, MPFR_RNDN); mpfr_mul(r10974, r10963, r10973, MPFR_RNDN); return mpfr_get_d(r10974, MPFR_RNDN); } static mpfr_t r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993, r10994, r10995, r10996, r10997, r10998, r10999, r11000, r11001, r11002, r11003, r11004, r11005, r11006, r11007, r11008, r11009, r11010, r11011, r11012, r11013, r11014, r11015, r11016, r11017, r11018, r11019, r11020, r11021, r11022, r11023, r11024, r11025, r11026, r11027, r11028, r11029, r11030, r11031, r11032, r11033, r11034, r11035, r11036, r11037, r11038, r11039; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10975); mpfr_init_set_str(r10976, "-7.243473966147417e-108", 10, MPFR_RNDN); mpfr_init(r10977); mpfr_init_set_str(r10978, "1", 10, MPFR_RNDN); mpfr_init(r10979); mpfr_init(r10980); mpfr_init(r10981); mpfr_init(r10982); mpfr_init_set_str(r10983, "2", 10, MPFR_RNDN); mpfr_init(r10984); mpfr_init(r10985); mpfr_init(r10986); mpfr_init(r10987); mpfr_init(r10988); mpfr_init(r10989); mpfr_init(r10990); mpfr_init(r10991); mpfr_init_set_str(r10992, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r10993, "-1", 10, MPFR_RNDN); mpfr_init(r10994); mpfr_init(r10995); mpfr_init(r10996); mpfr_init(r10997); mpfr_init(r10998); mpfr_init(r10999); mpfr_init(r11000); mpfr_init(r11001); mpfr_init(r11002); mpfr_init(r11003); mpfr_init(r11004); mpfr_init(r11005); mpfr_init_set_str(r11006, "5.02327006811133e-310", 10, MPFR_RNDN); mpfr_init(r11007); mpfr_init(r11008); mpfr_init(r11009); mpfr_init(r11010); mpfr_init(r11011); mpfr_init(r11012); mpfr_init(r11013); mpfr_init(r11014); mpfr_init(r11015); mpfr_init(r11016); mpfr_init(r11017); mpfr_init(r11018); mpfr_init(r11019); mpfr_init(r11020); mpfr_init(r11021); mpfr_init(r11022); mpfr_init(r11023); mpfr_init_set_str(r11024, "5.704880922779612e-67", 10, MPFR_RNDN); mpfr_init(r11025); mpfr_init(r11026); mpfr_init(r11027); mpfr_init(r11028); mpfr_init(r11029); mpfr_init(r11030); mpfr_init(r11031); mpfr_init(r11032); mpfr_init(r11033); mpfr_init(r11034); mpfr_init(r11035); mpfr_init(r11036); mpfr_init(r11037); mpfr_init(r11038); mpfr_init(r11039); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r10975, h, MPFR_RNDN); ; mpfr_set_si(r10977, mpfr_cmp(r10975, r10976) <= 0, MPFR_RNDN); ; mpfr_set_d(r10979, M, MPFR_RNDN); mpfr_set_d(r10980, D, MPFR_RNDN); mpfr_mul(r10981, r10979, r10980, MPFR_RNDN); mpfr_set_d(r10982, d, MPFR_RNDN); ; mpfr_mul(r10984, r10982, r10983, MPFR_RNDN); mpfr_div(r10985, r10981, r10984, MPFR_RNDN); mpfr_pow(r10986, r10985, r10983, MPFR_RNDN); mpfr_mul(r10987, r10975, r10986, MPFR_RNDN); mpfr_set_d(r10988, l, MPFR_RNDN); mpfr_mul(r10989, r10988, r10983, MPFR_RNDN); mpfr_div(r10990, r10987, r10989, MPFR_RNDN); mpfr_sub(r10991, r10978, r10990, MPFR_RNDN); ; ; mpfr_div(r10994, r10993, r10988, MPFR_RNDN); mpfr_log(r10995, r10994, MPFR_RNDN); mpfr_div(r10996, r10993, r10982, MPFR_RNDN); mpfr_log(r10997, r10996, MPFR_RNDN); mpfr_sub(r10998, r10995, r10997, MPFR_RNDN); mpfr_mul(r10999, r10992, r10998, MPFR_RNDN); mpfr_exp(r11000, r10999, MPFR_RNDN); mpfr_div(r11001, r10982, r10975, MPFR_RNDN); mpfr_div(r11002, r10978, r10983, MPFR_RNDN); mpfr_pow(r11003, r11001, r11002, MPFR_RNDN); mpfr_mul(r11004, r11000, r11003, MPFR_RNDN); mpfr_mul(r11005, r10991, r11004, MPFR_RNDN); ; mpfr_set_si(r11007, mpfr_cmp(r10975, r11006) <= 0, MPFR_RNDN); mpfr_div(r11008, r10980, r10983, MPFR_RNDN); mpfr_div(r11009, r10979, r10982, MPFR_RNDN); mpfr_mul(r11010, r11008, r11009, MPFR_RNDN); mpfr_mul(r11011, r10975, r11010, MPFR_RNDN); mpfr_mul(r11012, r11011, r11010, MPFR_RNDN); mpfr_div(r11013, r11012, r10989, MPFR_RNDN); mpfr_sub(r11014, r10978, r11013, MPFR_RNDN); mpfr_exp(r11015, r10992, MPFR_RNDN); mpfr_div(r11016, r10993, r10975, MPFR_RNDN); mpfr_log(r11017, r11016, MPFR_RNDN); mpfr_sub(r11018, r11017, r10997, MPFR_RNDN); mpfr_pow(r11019, r11015, r11018, MPFR_RNDN); mpfr_div(r11020, r10982, r10988, MPFR_RNDN); mpfr_pow(r11021, r11020, r11002, MPFR_RNDN); mpfr_mul(r11022, r11019, r11021, MPFR_RNDN); mpfr_mul(r11023, r11014, r11022, MPFR_RNDN); ; mpfr_set_si(r11025, mpfr_cmp(r10975, r11024) <= 0, MPFR_RNDN); mpfr_div(r11026, r10978, r10975, MPFR_RNDN); mpfr_pow(r11027, r11026, r11002, MPFR_RNDN); mpfr_pow(r11028, r10982, r11002, MPFR_RNDN); mpfr_mul(r11029, r11027, r11028, MPFR_RNDN); mpfr_mul(r11030, r11029, r11021, MPFR_RNDN); mpfr_mul(r11031, r11030, r10991, MPFR_RNDN); mpfr_div(r11032, r10978, r10988, MPFR_RNDN); mpfr_pow(r11033, r11032, r11002, MPFR_RNDN); mpfr_mul(r11034, r11028, r11033, MPFR_RNDN); mpfr_mul(r11035, r11034, r11003, MPFR_RNDN); mpfr_mul(r11036, r11035, r10991, MPFR_RNDN); if (mpfr_get_si(r11025, MPFR_RNDN)) { mpfr_set(r11037, r11031, MPFR_RNDN); } else { mpfr_set(r11037, r11036, MPFR_RNDN); }; if (mpfr_get_si(r11007, MPFR_RNDN)) { mpfr_set(r11038, r11023, MPFR_RNDN); } else { mpfr_set(r11038, r11037, MPFR_RNDN); }; if (mpfr_get_si(r10977, MPFR_RNDN)) { mpfr_set(r11039, r11005, MPFR_RNDN); } else { mpfr_set(r11039, r11038, MPFR_RNDN); }; return mpfr_get_d(r11039, MPFR_RNDN); } static mpfr_t r11040, r11041, r11042, r11043, r11044, r11045, r11046, r11047, r11048, r11049, r11050, r11051, r11052, r11053, r11054, r11055, r11056, r11057, r11058, r11059, r11060, r11061, r11062, r11063, r11064, r11065, r11066, r11067, r11068, r11069, r11070, r11071, r11072, r11073, r11074, r11075, r11076, r11077, r11078, r11079, r11080, r11081, r11082, r11083, r11084, r11085, r11086, r11087, r11088, r11089, r11090, r11091, r11092, r11093, r11094, r11095, r11096, r11097, r11098, r11099, r11100, r11101, r11102, r11103, r11104; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11040); mpfr_init_set_str(r11041, "-7.243473966147417e-108", 10, MPFR_RNDN); mpfr_init(r11042); mpfr_init_set_str(r11043, "1", 10, MPFR_RNDN); mpfr_init(r11044); mpfr_init(r11045); mpfr_init(r11046); mpfr_init(r11047); mpfr_init_set_str(r11048, "2", 10, MPFR_RNDN); mpfr_init(r11049); mpfr_init(r11050); mpfr_init(r11051); mpfr_init(r11052); mpfr_init(r11053); mpfr_init(r11054); mpfr_init(r11055); mpfr_init(r11056); mpfr_init_set_str(r11057, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r11058, "-1", 10, MPFR_RNDN); mpfr_init(r11059); mpfr_init(r11060); mpfr_init(r11061); mpfr_init(r11062); mpfr_init(r11063); mpfr_init(r11064); mpfr_init(r11065); mpfr_init(r11066); mpfr_init(r11067); mpfr_init(r11068); mpfr_init(r11069); mpfr_init(r11070); mpfr_init_set_str(r11071, "5.02327006811133e-310", 10, MPFR_RNDN); mpfr_init(r11072); mpfr_init(r11073); mpfr_init(r11074); mpfr_init(r11075); mpfr_init(r11076); mpfr_init(r11077); mpfr_init(r11078); mpfr_init(r11079); mpfr_init(r11080); mpfr_init(r11081); mpfr_init(r11082); mpfr_init(r11083); mpfr_init(r11084); mpfr_init(r11085); mpfr_init(r11086); mpfr_init(r11087); mpfr_init(r11088); mpfr_init_set_str(r11089, "5.704880922779612e-67", 10, MPFR_RNDN); mpfr_init(r11090); mpfr_init(r11091); mpfr_init(r11092); mpfr_init(r11093); mpfr_init(r11094); mpfr_init(r11095); mpfr_init(r11096); mpfr_init(r11097); mpfr_init(r11098); mpfr_init(r11099); mpfr_init(r11100); mpfr_init(r11101); mpfr_init(r11102); mpfr_init(r11103); mpfr_init(r11104); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r11040, h, MPFR_RNDN); ; mpfr_set_si(r11042, mpfr_cmp(r11040, r11041) <= 0, MPFR_RNDN); ; mpfr_set_d(r11044, M, MPFR_RNDN); mpfr_set_d(r11045, D, MPFR_RNDN); mpfr_mul(r11046, r11044, r11045, MPFR_RNDN); mpfr_set_d(r11047, d, MPFR_RNDN); ; mpfr_mul(r11049, r11047, r11048, MPFR_RNDN); mpfr_div(r11050, r11046, r11049, MPFR_RNDN); mpfr_pow(r11051, r11050, r11048, MPFR_RNDN); mpfr_mul(r11052, r11040, r11051, MPFR_RNDN); mpfr_set_d(r11053, l, MPFR_RNDN); mpfr_mul(r11054, r11053, r11048, MPFR_RNDN); mpfr_div(r11055, r11052, r11054, MPFR_RNDN); mpfr_sub(r11056, r11043, r11055, MPFR_RNDN); ; ; mpfr_div(r11059, r11058, r11053, MPFR_RNDN); mpfr_log(r11060, r11059, MPFR_RNDN); mpfr_div(r11061, r11058, r11047, MPFR_RNDN); mpfr_log(r11062, r11061, MPFR_RNDN); mpfr_sub(r11063, r11060, r11062, MPFR_RNDN); mpfr_mul(r11064, r11057, r11063, MPFR_RNDN); mpfr_exp(r11065, r11064, MPFR_RNDN); mpfr_div(r11066, r11047, r11040, MPFR_RNDN); mpfr_div(r11067, r11043, r11048, MPFR_RNDN); mpfr_pow(r11068, r11066, r11067, MPFR_RNDN); mpfr_mul(r11069, r11065, r11068, MPFR_RNDN); mpfr_mul(r11070, r11056, r11069, MPFR_RNDN); ; mpfr_set_si(r11072, mpfr_cmp(r11040, r11071) <= 0, MPFR_RNDN); mpfr_div(r11073, r11045, r11048, MPFR_RNDN); mpfr_div(r11074, r11044, r11047, MPFR_RNDN); mpfr_mul(r11075, r11073, r11074, MPFR_RNDN); mpfr_mul(r11076, r11040, r11075, MPFR_RNDN); mpfr_mul(r11077, r11076, r11075, MPFR_RNDN); mpfr_div(r11078, r11077, r11054, MPFR_RNDN); mpfr_sub(r11079, r11043, r11078, MPFR_RNDN); mpfr_exp(r11080, r11057, MPFR_RNDN); mpfr_div(r11081, r11058, r11040, MPFR_RNDN); mpfr_log(r11082, r11081, MPFR_RNDN); mpfr_sub(r11083, r11082, r11062, MPFR_RNDN); mpfr_pow(r11084, r11080, r11083, MPFR_RNDN); mpfr_div(r11085, r11047, r11053, MPFR_RNDN); mpfr_pow(r11086, r11085, r11067, MPFR_RNDN); mpfr_mul(r11087, r11084, r11086, MPFR_RNDN); mpfr_mul(r11088, r11079, r11087, MPFR_RNDN); ; mpfr_set_si(r11090, mpfr_cmp(r11040, r11089) <= 0, MPFR_RNDN); mpfr_div(r11091, r11043, r11040, MPFR_RNDN); mpfr_pow(r11092, r11091, r11067, MPFR_RNDN); mpfr_pow(r11093, r11047, r11067, MPFR_RNDN); mpfr_mul(r11094, r11092, r11093, MPFR_RNDN); mpfr_mul(r11095, r11094, r11086, MPFR_RNDN); mpfr_mul(r11096, r11095, r11056, MPFR_RNDN); mpfr_div(r11097, r11043, r11053, MPFR_RNDN); mpfr_pow(r11098, r11097, r11067, MPFR_RNDN); mpfr_mul(r11099, r11093, r11098, MPFR_RNDN); mpfr_mul(r11100, r11099, r11068, MPFR_RNDN); mpfr_mul(r11101, r11100, r11056, MPFR_RNDN); if (mpfr_get_si(r11090, MPFR_RNDN)) { mpfr_set(r11102, r11096, MPFR_RNDN); } else { mpfr_set(r11102, r11101, MPFR_RNDN); }; if (mpfr_get_si(r11072, MPFR_RNDN)) { mpfr_set(r11103, r11088, MPFR_RNDN); } else { mpfr_set(r11103, r11102, MPFR_RNDN); }; if (mpfr_get_si(r11042, MPFR_RNDN)) { mpfr_set(r11104, r11070, MPFR_RNDN); } else { mpfr_set(r11104, r11103, MPFR_RNDN); }; return mpfr_get_d(r11104, MPFR_RNDN); }