#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 r10679 = d; float r10680 = h; float r10681 = r10679 / r10680; float r10682 = 1; float r10683 = 2; float r10684 = r10682 / r10683; float r10685 = pow(r10681, r10684); float r10686 = l; float r10687 = r10679 / r10686; float r10688 = pow(r10687, r10684); float r10689 = r10685 * r10688; float r10690 = M; float r10691 = D; float r10692 = r10690 * r10691; float r10693 = r10683 * r10679; float r10694 = r10692 / r10693; float r10695 = pow(r10694, r10683); float r10696 = r10684 * r10695; float r10697 = r10680 / r10686; float r10698 = r10696 * r10697; float r10699 = r10682 - r10698; float r10700 = r10689 * r10699; return r10700; } double f_id(double d, double h, double l, double M, double D) { double r10701 = d; double r10702 = h; double r10703 = r10701 / r10702; double r10704 = 1; double r10705 = 2; double r10706 = r10704 / r10705; double r10707 = pow(r10703, r10706); double r10708 = l; double r10709 = r10701 / r10708; double r10710 = pow(r10709, r10706); double r10711 = r10707 * r10710; double r10712 = M; double r10713 = D; double r10714 = r10712 * r10713; double r10715 = r10705 * r10701; double r10716 = r10714 / r10715; double r10717 = pow(r10716, r10705); double r10718 = r10706 * r10717; double r10719 = r10702 / r10708; double r10720 = r10718 * r10719; double r10721 = r10704 - r10720; double r10722 = r10711 * r10721; return r10722; } double f_of(float d, float h, float l, float M, float D) { float r10723 = h; float r10724 = -3.1210736997942694e-116; bool r10725 = r10723 <= r10724; float r10726 = 1; float r10727 = M; float r10728 = D; float r10729 = r10727 * r10728; float r10730 = d; float r10731 = 2; float r10732 = r10730 * r10731; float r10733 = r10729 / r10732; float r10734 = pow(r10733, r10731); float r10735 = r10723 * r10734; float r10736 = l; float r10737 = r10736 * r10731; float r10738 = r10735 / r10737; float r10739 = r10726 - r10738; float r10740 = 1/2; float r10741 = -1; float r10742 = r10741 / r10736; float r10743 = log(r10742); float r10744 = r10741 / r10730; float r10745 = log(r10744); float r10746 = r10743 - r10745; float r10747 = r10740 * r10746; float r10748 = exp(r10747); float r10749 = r10730 / r10723; float r10750 = r10726 / r10731; float r10751 = pow(r10749, r10750); float r10752 = r10748 * r10751; float r10753 = r10739 * r10752; float r10754 = 5.02327006811133e-310; bool r10755 = r10723 <= r10754; float r10756 = r10741 / r10723; float r10757 = log(r10756); float r10758 = r10757 - r10745; float r10759 = r10758 * r10740; float r10760 = exp(r10759); float r10761 = r10730 / r10736; float r10762 = pow(r10761, r10750); float r10763 = r10760 * r10762; float r10764 = r10734 * r10750; float r10765 = r10723 / r10736; float r10766 = r10764 * r10765; float r10767 = r10726 - r10766; float r10768 = r10763 * r10767; float r10769 = 5.704880922779612e-67; bool r10770 = r10723 <= r10769; float r10771 = pow(r10730, r10750); float r10772 = r10726 / r10723; float r10773 = pow(r10772, r10750); float r10774 = r10771 * r10773; float r10775 = r10774 * r10762; float r10776 = r10775 * r10739; float r10777 = r10726 / r10736; float r10778 = pow(r10777, r10750); float r10779 = r10778 * r10771; float r10780 = r10751 * r10779; float r10781 = r10739 * r10780; float r10782 = r10770 ? r10776 : r10781; float r10783 = r10755 ? r10768 : r10782; float r10784 = r10725 ? r10753 : r10783; return r10784; } double f_od(double d, double h, double l, double M, double D) { double r10785 = h; double r10786 = -3.1210736997942694e-116; bool r10787 = r10785 <= r10786; double r10788 = 1; double r10789 = M; double r10790 = D; double r10791 = r10789 * r10790; double r10792 = d; double r10793 = 2; double r10794 = r10792 * r10793; double r10795 = r10791 / r10794; double r10796 = pow(r10795, r10793); double r10797 = r10785 * r10796; double r10798 = l; double r10799 = r10798 * r10793; double r10800 = r10797 / r10799; double r10801 = r10788 - r10800; double r10802 = 1/2; double r10803 = -1; double r10804 = r10803 / r10798; double r10805 = log(r10804); double r10806 = r10803 / r10792; double r10807 = log(r10806); double r10808 = r10805 - r10807; double r10809 = r10802 * r10808; double r10810 = exp(r10809); double r10811 = r10792 / r10785; double r10812 = r10788 / r10793; double r10813 = pow(r10811, r10812); double r10814 = r10810 * r10813; double r10815 = r10801 * r10814; double r10816 = 5.02327006811133e-310; bool r10817 = r10785 <= r10816; double r10818 = r10803 / r10785; double r10819 = log(r10818); double r10820 = r10819 - r10807; double r10821 = r10820 * r10802; double r10822 = exp(r10821); double r10823 = r10792 / r10798; double r10824 = pow(r10823, r10812); double r10825 = r10822 * r10824; double r10826 = r10796 * r10812; double r10827 = r10785 / r10798; double r10828 = r10826 * r10827; double r10829 = r10788 - r10828; double r10830 = r10825 * r10829; double r10831 = 5.704880922779612e-67; bool r10832 = r10785 <= r10831; double r10833 = pow(r10792, r10812); double r10834 = r10788 / r10785; double r10835 = pow(r10834, r10812); double r10836 = r10833 * r10835; double r10837 = r10836 * r10824; double r10838 = r10837 * r10801; double r10839 = r10788 / r10798; double r10840 = pow(r10839, r10812); double r10841 = r10840 * r10833; double r10842 = r10813 * r10841; double r10843 = r10801 * r10842; double r10844 = r10832 ? r10838 : r10843; double r10845 = r10817 ? r10830 : r10844; double r10846 = r10787 ? r10815 : r10845; return r10846; } 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 r10847, r10848, r10849, r10850, r10851, r10852, r10853, r10854, r10855, r10856, r10857, r10858, r10859, r10860, r10861, r10862, r10863, r10864, r10865, r10866, r10867, r10868; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10847); mpfr_init(r10848); mpfr_init(r10849); mpfr_init_set_str(r10850, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10851, "2", 10, MPFR_RNDN); mpfr_init(r10852); mpfr_init(r10853); mpfr_init(r10854); mpfr_init(r10855); mpfr_init(r10856); mpfr_init(r10857); mpfr_init(r10858); mpfr_init(r10859); mpfr_init(r10860); mpfr_init(r10861); mpfr_init(r10862); mpfr_init(r10863); mpfr_init(r10864); mpfr_init(r10865); mpfr_init(r10866); mpfr_init(r10867); mpfr_init(r10868); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r10847, d, MPFR_RNDN); mpfr_set_d(r10848, h, MPFR_RNDN); mpfr_div(r10849, r10847, r10848, MPFR_RNDN); ; ; mpfr_div(r10852, r10850, r10851, MPFR_RNDN); mpfr_pow(r10853, r10849, r10852, MPFR_RNDN); mpfr_set_d(r10854, l, MPFR_RNDN); mpfr_div(r10855, r10847, r10854, MPFR_RNDN); mpfr_pow(r10856, r10855, r10852, MPFR_RNDN); mpfr_mul(r10857, r10853, r10856, MPFR_RNDN); mpfr_set_d(r10858, M, MPFR_RNDN); mpfr_set_d(r10859, D, MPFR_RNDN); mpfr_mul(r10860, r10858, r10859, MPFR_RNDN); mpfr_mul(r10861, r10851, r10847, MPFR_RNDN); mpfr_div(r10862, r10860, r10861, MPFR_RNDN); mpfr_pow(r10863, r10862, r10851, MPFR_RNDN); mpfr_mul(r10864, r10852, r10863, MPFR_RNDN); mpfr_div(r10865, r10848, r10854, MPFR_RNDN); mpfr_mul(r10866, r10864, r10865, MPFR_RNDN); mpfr_sub(r10867, r10850, r10866, MPFR_RNDN); mpfr_mul(r10868, r10857, r10867, MPFR_RNDN); return mpfr_get_d(r10868, MPFR_RNDN); } static mpfr_t r10869, r10870, r10871, r10872, r10873, r10874, r10875, r10876, r10877, r10878, r10879, r10880, r10881, r10882, r10883, r10884, r10885, r10886, r10887, r10888, r10889, r10890, r10891, r10892, r10893, r10894, r10895, r10896, r10897, r10898, r10899, r10900, r10901, r10902, r10903, r10904, r10905, r10906, r10907, r10908, r10909, r10910, r10911, r10912, r10913, r10914, r10915, r10916, r10917, r10918, r10919, r10920, r10921, r10922, r10923, r10924, r10925, r10926, r10927, r10928, r10929, r10930; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10869); mpfr_init_set_str(r10870, "-3.1210736997942694e-116", 10, MPFR_RNDN); mpfr_init(r10871); mpfr_init_set_str(r10872, "1", 10, MPFR_RNDN); mpfr_init(r10873); mpfr_init(r10874); mpfr_init(r10875); mpfr_init(r10876); mpfr_init_set_str(r10877, "2", 10, MPFR_RNDN); mpfr_init(r10878); mpfr_init(r10879); mpfr_init(r10880); mpfr_init(r10881); mpfr_init(r10882); mpfr_init(r10883); mpfr_init(r10884); mpfr_init(r10885); mpfr_init_set_str(r10886, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r10887, "-1", 10, MPFR_RNDN); mpfr_init(r10888); mpfr_init(r10889); mpfr_init(r10890); mpfr_init(r10891); mpfr_init(r10892); mpfr_init(r10893); mpfr_init(r10894); mpfr_init(r10895); mpfr_init(r10896); mpfr_init(r10897); mpfr_init(r10898); mpfr_init(r10899); mpfr_init_set_str(r10900, "5.02327006811133e-310", 10, MPFR_RNDN); mpfr_init(r10901); mpfr_init(r10902); mpfr_init(r10903); mpfr_init(r10904); mpfr_init(r10905); mpfr_init(r10906); mpfr_init(r10907); mpfr_init(r10908); mpfr_init(r10909); mpfr_init(r10910); mpfr_init(r10911); mpfr_init(r10912); mpfr_init(r10913); mpfr_init(r10914); mpfr_init_set_str(r10915, "5.704880922779612e-67", 10, MPFR_RNDN); mpfr_init(r10916); mpfr_init(r10917); mpfr_init(r10918); mpfr_init(r10919); mpfr_init(r10920); mpfr_init(r10921); mpfr_init(r10922); mpfr_init(r10923); mpfr_init(r10924); mpfr_init(r10925); mpfr_init(r10926); mpfr_init(r10927); mpfr_init(r10928); mpfr_init(r10929); mpfr_init(r10930); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r10869, h, MPFR_RNDN); ; mpfr_set_si(r10871, mpfr_cmp(r10869, r10870) <= 0, MPFR_RNDN); ; mpfr_set_d(r10873, M, MPFR_RNDN); mpfr_set_d(r10874, D, MPFR_RNDN); mpfr_mul(r10875, r10873, r10874, MPFR_RNDN); mpfr_set_d(r10876, d, MPFR_RNDN); ; mpfr_mul(r10878, r10876, r10877, MPFR_RNDN); mpfr_div(r10879, r10875, r10878, MPFR_RNDN); mpfr_pow(r10880, r10879, r10877, MPFR_RNDN); mpfr_mul(r10881, r10869, r10880, MPFR_RNDN); mpfr_set_d(r10882, l, MPFR_RNDN); mpfr_mul(r10883, r10882, r10877, MPFR_RNDN); mpfr_div(r10884, r10881, r10883, MPFR_RNDN); mpfr_sub(r10885, r10872, r10884, MPFR_RNDN); ; ; mpfr_div(r10888, r10887, r10882, MPFR_RNDN); mpfr_log(r10889, r10888, MPFR_RNDN); mpfr_div(r10890, r10887, r10876, MPFR_RNDN); mpfr_log(r10891, r10890, MPFR_RNDN); mpfr_sub(r10892, r10889, r10891, MPFR_RNDN); mpfr_mul(r10893, r10886, r10892, MPFR_RNDN); mpfr_exp(r10894, r10893, MPFR_RNDN); mpfr_div(r10895, r10876, r10869, MPFR_RNDN); mpfr_div(r10896, r10872, r10877, MPFR_RNDN); mpfr_pow(r10897, r10895, r10896, MPFR_RNDN); mpfr_mul(r10898, r10894, r10897, MPFR_RNDN); mpfr_mul(r10899, r10885, r10898, MPFR_RNDN); ; mpfr_set_si(r10901, mpfr_cmp(r10869, r10900) <= 0, MPFR_RNDN); mpfr_div(r10902, r10887, r10869, MPFR_RNDN); mpfr_log(r10903, r10902, MPFR_RNDN); mpfr_sub(r10904, r10903, r10891, MPFR_RNDN); mpfr_mul(r10905, r10904, r10886, MPFR_RNDN); mpfr_exp(r10906, r10905, MPFR_RNDN); mpfr_div(r10907, r10876, r10882, MPFR_RNDN); mpfr_pow(r10908, r10907, r10896, MPFR_RNDN); mpfr_mul(r10909, r10906, r10908, MPFR_RNDN); mpfr_mul(r10910, r10880, r10896, MPFR_RNDN); mpfr_div(r10911, r10869, r10882, MPFR_RNDN); mpfr_mul(r10912, r10910, r10911, MPFR_RNDN); mpfr_sub(r10913, r10872, r10912, MPFR_RNDN); mpfr_mul(r10914, r10909, r10913, MPFR_RNDN); ; mpfr_set_si(r10916, mpfr_cmp(r10869, r10915) <= 0, MPFR_RNDN); mpfr_pow(r10917, r10876, r10896, MPFR_RNDN); mpfr_div(r10918, r10872, r10869, MPFR_RNDN); mpfr_pow(r10919, r10918, r10896, MPFR_RNDN); mpfr_mul(r10920, r10917, r10919, MPFR_RNDN); mpfr_mul(r10921, r10920, r10908, MPFR_RNDN); mpfr_mul(r10922, r10921, r10885, MPFR_RNDN); mpfr_div(r10923, r10872, r10882, MPFR_RNDN); mpfr_pow(r10924, r10923, r10896, MPFR_RNDN); mpfr_mul(r10925, r10924, r10917, MPFR_RNDN); mpfr_mul(r10926, r10897, r10925, MPFR_RNDN); mpfr_mul(r10927, r10885, r10926, MPFR_RNDN); if (mpfr_get_si(r10916, MPFR_RNDN)) { mpfr_set(r10928, r10922, MPFR_RNDN); } else { mpfr_set(r10928, r10927, MPFR_RNDN); }; if (mpfr_get_si(r10901, MPFR_RNDN)) { mpfr_set(r10929, r10914, MPFR_RNDN); } else { mpfr_set(r10929, r10928, MPFR_RNDN); }; if (mpfr_get_si(r10871, MPFR_RNDN)) { mpfr_set(r10930, r10899, MPFR_RNDN); } else { mpfr_set(r10930, r10929, MPFR_RNDN); }; return mpfr_get_d(r10930, MPFR_RNDN); } static mpfr_t r10931, r10932, r10933, r10934, r10935, r10936, r10937, r10938, r10939, r10940, r10941, r10942, r10943, r10944, r10945, r10946, r10947, r10948, r10949, r10950, r10951, r10952, r10953, r10954, r10955, r10956, r10957, r10958, r10959, r10960, r10961, r10962, r10963, r10964, r10965, r10966, r10967, r10968, r10969, r10970, r10971, r10972, r10973, r10974, r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10931); mpfr_init_set_str(r10932, "-3.1210736997942694e-116", 10, MPFR_RNDN); mpfr_init(r10933); mpfr_init_set_str(r10934, "1", 10, MPFR_RNDN); mpfr_init(r10935); mpfr_init(r10936); mpfr_init(r10937); mpfr_init(r10938); mpfr_init_set_str(r10939, "2", 10, MPFR_RNDN); mpfr_init(r10940); mpfr_init(r10941); mpfr_init(r10942); mpfr_init(r10943); mpfr_init(r10944); mpfr_init(r10945); mpfr_init(r10946); mpfr_init(r10947); mpfr_init_set_str(r10948, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r10949, "-1", 10, MPFR_RNDN); mpfr_init(r10950); mpfr_init(r10951); mpfr_init(r10952); mpfr_init(r10953); mpfr_init(r10954); mpfr_init(r10955); mpfr_init(r10956); mpfr_init(r10957); mpfr_init(r10958); mpfr_init(r10959); mpfr_init(r10960); mpfr_init(r10961); mpfr_init_set_str(r10962, "5.02327006811133e-310", 10, MPFR_RNDN); 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); mpfr_init(r10975); mpfr_init(r10976); mpfr_init_set_str(r10977, "5.704880922779612e-67", 10, MPFR_RNDN); mpfr_init(r10978); mpfr_init(r10979); mpfr_init(r10980); mpfr_init(r10981); mpfr_init(r10982); mpfr_init(r10983); 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(r10992); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r10931, h, MPFR_RNDN); ; mpfr_set_si(r10933, mpfr_cmp(r10931, r10932) <= 0, MPFR_RNDN); ; mpfr_set_d(r10935, M, MPFR_RNDN); mpfr_set_d(r10936, D, MPFR_RNDN); mpfr_mul(r10937, r10935, r10936, MPFR_RNDN); mpfr_set_d(r10938, d, MPFR_RNDN); ; mpfr_mul(r10940, r10938, r10939, MPFR_RNDN); mpfr_div(r10941, r10937, r10940, MPFR_RNDN); mpfr_pow(r10942, r10941, r10939, MPFR_RNDN); mpfr_mul(r10943, r10931, r10942, MPFR_RNDN); mpfr_set_d(r10944, l, MPFR_RNDN); mpfr_mul(r10945, r10944, r10939, MPFR_RNDN); mpfr_div(r10946, r10943, r10945, MPFR_RNDN); mpfr_sub(r10947, r10934, r10946, MPFR_RNDN); ; ; mpfr_div(r10950, r10949, r10944, MPFR_RNDN); mpfr_log(r10951, r10950, MPFR_RNDN); mpfr_div(r10952, r10949, r10938, MPFR_RNDN); mpfr_log(r10953, r10952, MPFR_RNDN); mpfr_sub(r10954, r10951, r10953, MPFR_RNDN); mpfr_mul(r10955, r10948, r10954, MPFR_RNDN); mpfr_exp(r10956, r10955, MPFR_RNDN); mpfr_div(r10957, r10938, r10931, MPFR_RNDN); mpfr_div(r10958, r10934, r10939, MPFR_RNDN); mpfr_pow(r10959, r10957, r10958, MPFR_RNDN); mpfr_mul(r10960, r10956, r10959, MPFR_RNDN); mpfr_mul(r10961, r10947, r10960, MPFR_RNDN); ; mpfr_set_si(r10963, mpfr_cmp(r10931, r10962) <= 0, MPFR_RNDN); mpfr_div(r10964, r10949, r10931, MPFR_RNDN); mpfr_log(r10965, r10964, MPFR_RNDN); mpfr_sub(r10966, r10965, r10953, MPFR_RNDN); mpfr_mul(r10967, r10966, r10948, MPFR_RNDN); mpfr_exp(r10968, r10967, MPFR_RNDN); mpfr_div(r10969, r10938, r10944, MPFR_RNDN); mpfr_pow(r10970, r10969, r10958, MPFR_RNDN); mpfr_mul(r10971, r10968, r10970, MPFR_RNDN); mpfr_mul(r10972, r10942, r10958, MPFR_RNDN); mpfr_div(r10973, r10931, r10944, MPFR_RNDN); mpfr_mul(r10974, r10972, r10973, MPFR_RNDN); mpfr_sub(r10975, r10934, r10974, MPFR_RNDN); mpfr_mul(r10976, r10971, r10975, MPFR_RNDN); ; mpfr_set_si(r10978, mpfr_cmp(r10931, r10977) <= 0, MPFR_RNDN); mpfr_pow(r10979, r10938, r10958, MPFR_RNDN); mpfr_div(r10980, r10934, r10931, MPFR_RNDN); mpfr_pow(r10981, r10980, r10958, MPFR_RNDN); mpfr_mul(r10982, r10979, r10981, MPFR_RNDN); mpfr_mul(r10983, r10982, r10970, MPFR_RNDN); mpfr_mul(r10984, r10983, r10947, MPFR_RNDN); mpfr_div(r10985, r10934, r10944, MPFR_RNDN); mpfr_pow(r10986, r10985, r10958, MPFR_RNDN); mpfr_mul(r10987, r10986, r10979, MPFR_RNDN); mpfr_mul(r10988, r10959, r10987, MPFR_RNDN); mpfr_mul(r10989, r10947, r10988, MPFR_RNDN); if (mpfr_get_si(r10978, MPFR_RNDN)) { mpfr_set(r10990, r10984, MPFR_RNDN); } else { mpfr_set(r10990, r10989, MPFR_RNDN); }; if (mpfr_get_si(r10963, MPFR_RNDN)) { mpfr_set(r10991, r10976, MPFR_RNDN); } else { mpfr_set(r10991, r10990, MPFR_RNDN); }; if (mpfr_get_si(r10933, MPFR_RNDN)) { mpfr_set(r10992, r10961, MPFR_RNDN); } else { mpfr_set(r10992, r10991, MPFR_RNDN); }; return mpfr_get_d(r10992, MPFR_RNDN); }