#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r10713 = 1; float r10714 = 0.1049934947; float r10715 = x; float r10716 = r10715 * r10715; float r10717 = r10714 * r10716; float r10718 = r10713 + r10717; float r10719 = 0.0424060604; float r10720 = r10716 * r10716; float r10721 = r10719 * r10720; float r10722 = r10718 + r10721; float r10723 = 0.0072644182; float r10724 = r10720 * r10716; float r10725 = r10723 * r10724; float r10726 = r10722 + r10725; float r10727 = 0.0005064034; float r10728 = r10724 * r10716; float r10729 = r10727 * r10728; float r10730 = r10726 + r10729; float r10731 = 0.0001789971; float r10732 = r10728 * r10716; float r10733 = r10731 * r10732; float r10734 = r10730 + r10733; float r10735 = 0.7715471019; float r10736 = r10735 * r10716; float r10737 = r10713 + r10736; float r10738 = 0.2909738639; float r10739 = r10738 * r10720; float r10740 = r10737 + r10739; float r10741 = 0.0694555761; float r10742 = r10741 * r10724; float r10743 = r10740 + r10742; float r10744 = 0.0140005442; float r10745 = r10744 * r10728; float r10746 = r10743 + r10745; float r10747 = 0.0008327945; float r10748 = r10747 * r10732; float r10749 = r10746 + r10748; float r10750 = 2; float r10751 = r10750 * r10731; float r10752 = r10732 * r10716; float r10753 = r10751 * r10752; float r10754 = r10749 + r10753; float r10755 = r10734 / r10754; float r10756 = r10755 * r10715; return r10756; } double f_id(double x) { double r10757 = 1; double r10758 = 0.1049934947; double r10759 = x; double r10760 = r10759 * r10759; double r10761 = r10758 * r10760; double r10762 = r10757 + r10761; double r10763 = 0.0424060604; double r10764 = r10760 * r10760; double r10765 = r10763 * r10764; double r10766 = r10762 + r10765; double r10767 = 0.0072644182; double r10768 = r10764 * r10760; double r10769 = r10767 * r10768; double r10770 = r10766 + r10769; double r10771 = 0.0005064034; double r10772 = r10768 * r10760; double r10773 = r10771 * r10772; double r10774 = r10770 + r10773; double r10775 = 0.0001789971; double r10776 = r10772 * r10760; double r10777 = r10775 * r10776; double r10778 = r10774 + r10777; double r10779 = 0.7715471019; double r10780 = r10779 * r10760; double r10781 = r10757 + r10780; double r10782 = 0.2909738639; double r10783 = r10782 * r10764; double r10784 = r10781 + r10783; double r10785 = 0.0694555761; double r10786 = r10785 * r10768; double r10787 = r10784 + r10786; double r10788 = 0.0140005442; double r10789 = r10788 * r10772; double r10790 = r10787 + r10789; double r10791 = 0.0008327945; double r10792 = r10791 * r10776; double r10793 = r10790 + r10792; double r10794 = 2; double r10795 = r10794 * r10775; double r10796 = r10776 * r10760; double r10797 = r10795 * r10796; double r10798 = r10793 + r10797; double r10799 = r10778 / r10798; double r10800 = r10799 * r10759; return r10800; } double f_of(float x) { float r10801 = x; float r10802 = -42794565.356408514; bool r10803 = r10801 <= r10802; float r10804 = 52046178.931661196; bool r10805 = r10801 <= r10804; float r10806 = !r10805; bool r10807 = r10803 || r10806; float r10808 = 0.5; float r10809 = r10808 / r10801; float r10810 = 43733.204511252174; float r10811 = 15; float r10812 = pow(r10801, r10811); float r10813 = r10810 / r10812; float r10814 = r10809 - r10813; float r10815 = 666.049723509856; float r10816 = 9; float r10817 = pow(r10801, r10816); float r10818 = r10815 / r10817; float r10819 = 71.24974274308389; float r10820 = 7; float r10821 = pow(r10801, r10820); float r10822 = r10819 / r10821; float r10823 = 20.29199970278848; float r10824 = 5; float r10825 = pow(r10801, r10824); float r10826 = r10823 / r10825; float r10827 = r10822 + r10826; float r10828 = r10818 - r10827; float r10829 = r10814 - r10828; float r10830 = 1; float r10831 = 0.1049934947; float r10832 = r10801 * r10801; float r10833 = r10831 * r10832; float r10834 = r10830 + r10833; float r10835 = r10832 * r10832; float r10836 = 0.0072644182; float r10837 = r10832 * r10836; float r10838 = 0.0424060604; float r10839 = r10837 + r10838; float r10840 = r10835 * r10839; float r10841 = 3; float r10842 = r10841 + r10830; float r10843 = pow(r10832, r10842); float r10844 = 0.0001789971; float r10845 = r10832 * r10844; float r10846 = 0.0005064034; float r10847 = r10845 + r10846; float r10848 = r10843 * r10847; float r10849 = r10840 + r10848; float r10850 = r10834 + r10849; float r10851 = 0.2909738639; float r10852 = -r10851; float r10853 = r10835 * r10852; float r10854 = 0.7715471019; float r10855 = r10832 * r10854; float r10856 = r10830 + r10855; float r10857 = r10853 - r10856; float r10858 = pow(r10832, r10841); float r10859 = r10858 * r10835; float r10860 = 0.0008327945; float r10861 = 2; float r10862 = r10844 * r10861; float r10863 = r10832 * r10862; float r10864 = r10860 + r10863; float r10865 = r10859 * r10864; float r10866 = 0.0694555761; float r10867 = 0.0140005442; float r10868 = r10832 * r10867; float r10869 = r10866 + r10868; float r10870 = r10858 * r10869; float r10871 = r10865 + r10870; float r10872 = r10857 - r10871; float r10873 = r10850 / r10872; float r10874 = -r10801; float r10875 = r10873 * r10874; float r10876 = r10807 ? r10829 : r10875; return r10876; } double f_od(double x) { double r10877 = x; double r10878 = -42794565.356408514; bool r10879 = r10877 <= r10878; double r10880 = 52046178.931661196; bool r10881 = r10877 <= r10880; double r10882 = !r10881; bool r10883 = r10879 || r10882; double r10884 = 0.5; double r10885 = r10884 / r10877; double r10886 = 43733.204511252174; double r10887 = 15; double r10888 = pow(r10877, r10887); double r10889 = r10886 / r10888; double r10890 = r10885 - r10889; double r10891 = 666.049723509856; double r10892 = 9; double r10893 = pow(r10877, r10892); double r10894 = r10891 / r10893; double r10895 = 71.24974274308389; double r10896 = 7; double r10897 = pow(r10877, r10896); double r10898 = r10895 / r10897; double r10899 = 20.29199970278848; double r10900 = 5; double r10901 = pow(r10877, r10900); double r10902 = r10899 / r10901; double r10903 = r10898 + r10902; double r10904 = r10894 - r10903; double r10905 = r10890 - r10904; double r10906 = 1; double r10907 = 0.1049934947; double r10908 = r10877 * r10877; double r10909 = r10907 * r10908; double r10910 = r10906 + r10909; double r10911 = r10908 * r10908; double r10912 = 0.0072644182; double r10913 = r10908 * r10912; double r10914 = 0.0424060604; double r10915 = r10913 + r10914; double r10916 = r10911 * r10915; double r10917 = 3; double r10918 = r10917 + r10906; double r10919 = pow(r10908, r10918); double r10920 = 0.0001789971; double r10921 = r10908 * r10920; double r10922 = 0.0005064034; double r10923 = r10921 + r10922; double r10924 = r10919 * r10923; double r10925 = r10916 + r10924; double r10926 = r10910 + r10925; double r10927 = 0.2909738639; double r10928 = -r10927; double r10929 = r10911 * r10928; double r10930 = 0.7715471019; double r10931 = r10908 * r10930; double r10932 = r10906 + r10931; double r10933 = r10929 - r10932; double r10934 = pow(r10908, r10917); double r10935 = r10934 * r10911; double r10936 = 0.0008327945; double r10937 = 2; double r10938 = r10920 * r10937; double r10939 = r10908 * r10938; double r10940 = r10936 + r10939; double r10941 = r10935 * r10940; double r10942 = 0.0694555761; double r10943 = 0.0140005442; double r10944 = r10908 * r10943; double r10945 = r10942 + r10944; double r10946 = r10934 * r10945; double r10947 = r10941 + r10946; double r10948 = r10933 - r10947; double r10949 = r10926 / r10948; double r10950 = -r10877; double r10951 = r10949 * r10950; double r10952 = r10883 ? r10905 : 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, r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993, r10994, r10995, r10996; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r10953, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10954, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r10955); mpfr_init(r10956); mpfr_init(r10957); mpfr_init(r10958); mpfr_init_set_str(r10959, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r10960); mpfr_init(r10961); mpfr_init(r10962); mpfr_init_set_str(r10963, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r10964); mpfr_init(r10965); mpfr_init(r10966); mpfr_init_set_str(r10967, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r10968); mpfr_init(r10969); mpfr_init(r10970); mpfr_init_set_str(r10971, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r10972); mpfr_init(r10973); mpfr_init(r10974); mpfr_init_set_str(r10975, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r10976); mpfr_init(r10977); mpfr_init_set_str(r10978, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r10979); mpfr_init(r10980); mpfr_init_set_str(r10981, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r10982); mpfr_init(r10983); mpfr_init_set_str(r10984, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r10985); mpfr_init(r10986); mpfr_init_set_str(r10987, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r10988); mpfr_init(r10989); mpfr_init_set_str(r10990, "2", 10, MPFR_RNDN); mpfr_init(r10991); mpfr_init(r10992); mpfr_init(r10993); mpfr_init(r10994); mpfr_init(r10995); mpfr_init(r10996); } double f_im(double x) { ; ; mpfr_set_d(r10955, x, MPFR_RNDN); mpfr_mul(r10956, r10955, r10955, MPFR_RNDN); mpfr_mul(r10957, r10954, r10956, MPFR_RNDN); mpfr_add(r10958, r10953, r10957, MPFR_RNDN); ; mpfr_mul(r10960, r10956, r10956, MPFR_RNDN); mpfr_mul(r10961, r10959, r10960, MPFR_RNDN); mpfr_add(r10962, r10958, r10961, MPFR_RNDN); ; mpfr_mul(r10964, r10960, r10956, MPFR_RNDN); mpfr_mul(r10965, r10963, r10964, MPFR_RNDN); mpfr_add(r10966, r10962, r10965, MPFR_RNDN); ; mpfr_mul(r10968, r10964, r10956, MPFR_RNDN); mpfr_mul(r10969, r10967, r10968, MPFR_RNDN); mpfr_add(r10970, r10966, r10969, MPFR_RNDN); ; mpfr_mul(r10972, r10968, r10956, MPFR_RNDN); mpfr_mul(r10973, r10971, r10972, MPFR_RNDN); mpfr_add(r10974, r10970, r10973, MPFR_RNDN); ; mpfr_mul(r10976, r10975, r10956, MPFR_RNDN); mpfr_add(r10977, r10953, r10976, MPFR_RNDN); ; mpfr_mul(r10979, r10978, r10960, MPFR_RNDN); mpfr_add(r10980, r10977, r10979, MPFR_RNDN); ; mpfr_mul(r10982, r10981, r10964, MPFR_RNDN); mpfr_add(r10983, r10980, r10982, MPFR_RNDN); ; mpfr_mul(r10985, r10984, r10968, MPFR_RNDN); mpfr_add(r10986, r10983, r10985, MPFR_RNDN); ; mpfr_mul(r10988, r10987, r10972, MPFR_RNDN); mpfr_add(r10989, r10986, r10988, MPFR_RNDN); ; mpfr_mul(r10991, r10990, r10971, MPFR_RNDN); mpfr_mul(r10992, r10972, r10956, MPFR_RNDN); mpfr_mul(r10993, r10991, r10992, MPFR_RNDN); mpfr_add(r10994, r10989, r10993, MPFR_RNDN); mpfr_div(r10995, r10974, r10994, MPFR_RNDN); mpfr_mul(r10996, r10995, r10955, MPFR_RNDN); return mpfr_get_d(r10996, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10997); mpfr_init_set_str(r10998, "-42794565.356408514", 10, MPFR_RNDN); mpfr_init(r10999); mpfr_init_set_str(r11000, "52046178.931661196", 10, MPFR_RNDN); mpfr_init(r11001); mpfr_init(r11002); mpfr_init(r11003); mpfr_init_set_str(r11004, "0.5", 10, MPFR_RNDN); mpfr_init(r11005); mpfr_init_set_str(r11006, "43733.204511252174", 10, MPFR_RNDN); mpfr_init_set_str(r11007, "15", 10, MPFR_RNDN); mpfr_init(r11008); mpfr_init(r11009); mpfr_init(r11010); mpfr_init_set_str(r11011, "666.049723509856", 10, MPFR_RNDN); mpfr_init_set_str(r11012, "9", 10, MPFR_RNDN); mpfr_init(r11013); mpfr_init(r11014); mpfr_init_set_str(r11015, "71.24974274308389", 10, MPFR_RNDN); mpfr_init_set_str(r11016, "7", 10, MPFR_RNDN); mpfr_init(r11017); mpfr_init(r11018); mpfr_init_set_str(r11019, "20.29199970278848", 10, MPFR_RNDN); mpfr_init_set_str(r11020, "5", 10, MPFR_RNDN); mpfr_init(r11021); mpfr_init(r11022); mpfr_init(r11023); mpfr_init(r11024); mpfr_init(r11025); mpfr_init_set_str(r11026, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11027, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11028); mpfr_init(r11029); mpfr_init(r11030); mpfr_init(r11031); mpfr_init_set_str(r11032, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11033); mpfr_init_set_str(r11034, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11035); mpfr_init(r11036); mpfr_init_set_str(r11037, "3", 10, MPFR_RNDN); mpfr_init(r11038); mpfr_init(r11039); mpfr_init_set_str(r11040, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11041); mpfr_init_set_str(r11042, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11043); mpfr_init(r11044); mpfr_init(r11045); mpfr_init(r11046); mpfr_init_set_str(r11047, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11048); mpfr_init(r11049); mpfr_init_set_str(r11050, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11051); mpfr_init(r11052); mpfr_init(r11053); mpfr_init(r11054); mpfr_init(r11055); mpfr_init_set_str(r11056, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r11057, "2", 10, MPFR_RNDN); mpfr_init(r11058); mpfr_init(r11059); mpfr_init(r11060); mpfr_init(r11061); mpfr_init_set_str(r11062, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r11063, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11064); mpfr_init(r11065); mpfr_init(r11066); mpfr_init(r11067); mpfr_init(r11068); mpfr_init(r11069); mpfr_init(r11070); mpfr_init(r11071); mpfr_init(r11072); } double f_fm(double x) { mpfr_set_d(r10997, x, MPFR_RNDN); ; mpfr_set_si(r10999, mpfr_cmp(r10997, r10998) <= 0, MPFR_RNDN); ; mpfr_set_si(r11001, mpfr_cmp(r10997, r11000) <= 0, MPFR_RNDN); mpfr_set_si(r11002, !mpfr_get_si(r11001, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11003, mpfr_get_si(r10999, MPFR_RNDN) || mpfr_get_si(r11002, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r11005, r11004, r10997, MPFR_RNDN); ; ; mpfr_pow(r11008, r10997, r11007, MPFR_RNDN); mpfr_div(r11009, r11006, r11008, MPFR_RNDN); mpfr_sub(r11010, r11005, r11009, MPFR_RNDN); ; ; mpfr_pow(r11013, r10997, r11012, MPFR_RNDN); mpfr_div(r11014, r11011, r11013, MPFR_RNDN); ; ; mpfr_pow(r11017, r10997, r11016, MPFR_RNDN); mpfr_div(r11018, r11015, r11017, MPFR_RNDN); ; ; mpfr_pow(r11021, r10997, r11020, MPFR_RNDN); mpfr_div(r11022, r11019, r11021, MPFR_RNDN); mpfr_add(r11023, r11018, r11022, MPFR_RNDN); mpfr_sub(r11024, r11014, r11023, MPFR_RNDN); mpfr_sub(r11025, r11010, r11024, MPFR_RNDN); ; ; mpfr_mul(r11028, r10997, r10997, MPFR_RNDN); mpfr_mul(r11029, r11027, r11028, MPFR_RNDN); mpfr_add(r11030, r11026, r11029, MPFR_RNDN); mpfr_mul(r11031, r11028, r11028, MPFR_RNDN); ; mpfr_mul(r11033, r11028, r11032, MPFR_RNDN); ; mpfr_add(r11035, r11033, r11034, MPFR_RNDN); mpfr_mul(r11036, r11031, r11035, MPFR_RNDN); ; mpfr_add(r11038, r11037, r11026, MPFR_RNDN); mpfr_pow(r11039, r11028, r11038, MPFR_RNDN); ; mpfr_mul(r11041, r11028, r11040, MPFR_RNDN); ; mpfr_add(r11043, r11041, r11042, MPFR_RNDN); mpfr_mul(r11044, r11039, r11043, MPFR_RNDN); mpfr_add(r11045, r11036, r11044, MPFR_RNDN); mpfr_add(r11046, r11030, r11045, MPFR_RNDN); ; mpfr_neg(r11048, r11047, MPFR_RNDN); mpfr_mul(r11049, r11031, r11048, MPFR_RNDN); ; mpfr_mul(r11051, r11028, r11050, MPFR_RNDN); mpfr_add(r11052, r11026, r11051, MPFR_RNDN); mpfr_sub(r11053, r11049, r11052, MPFR_RNDN); mpfr_pow(r11054, r11028, r11037, MPFR_RNDN); mpfr_mul(r11055, r11054, r11031, MPFR_RNDN); ; ; mpfr_mul(r11058, r11040, r11057, MPFR_RNDN); mpfr_mul(r11059, r11028, r11058, MPFR_RNDN); mpfr_add(r11060, r11056, r11059, MPFR_RNDN); mpfr_mul(r11061, r11055, r11060, MPFR_RNDN); ; ; mpfr_mul(r11064, r11028, r11063, MPFR_RNDN); mpfr_add(r11065, r11062, r11064, MPFR_RNDN); mpfr_mul(r11066, r11054, r11065, MPFR_RNDN); mpfr_add(r11067, r11061, r11066, MPFR_RNDN); mpfr_sub(r11068, r11053, r11067, MPFR_RNDN); mpfr_div(r11069, r11046, r11068, MPFR_RNDN); mpfr_neg(r11070, r10997, MPFR_RNDN); mpfr_mul(r11071, r11069, r11070, MPFR_RNDN); if (mpfr_get_si(r11003, MPFR_RNDN)) { mpfr_set(r11072, r11025, MPFR_RNDN); } else { mpfr_set(r11072, r11071, MPFR_RNDN); }; return mpfr_get_d(r11072, MPFR_RNDN); } static mpfr_t 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, r11105, r11106, r11107, r11108, r11109, r11110, r11111, r11112, r11113, r11114, r11115, r11116, r11117, r11118, r11119, r11120, r11121, r11122, r11123, r11124, r11125, r11126, r11127, r11128, r11129, r11130, r11131, r11132, r11133, r11134, r11135, r11136, r11137, r11138, r11139, r11140, r11141, r11142, r11143, r11144, r11145, r11146, r11147, r11148; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11073); mpfr_init_set_str(r11074, "-42794565.356408514", 10, MPFR_RNDN); mpfr_init(r11075); mpfr_init_set_str(r11076, "52046178.931661196", 10, MPFR_RNDN); mpfr_init(r11077); mpfr_init(r11078); mpfr_init(r11079); mpfr_init_set_str(r11080, "0.5", 10, MPFR_RNDN); mpfr_init(r11081); mpfr_init_set_str(r11082, "43733.204511252174", 10, MPFR_RNDN); mpfr_init_set_str(r11083, "15", 10, MPFR_RNDN); mpfr_init(r11084); mpfr_init(r11085); mpfr_init(r11086); mpfr_init_set_str(r11087, "666.049723509856", 10, MPFR_RNDN); mpfr_init_set_str(r11088, "9", 10, MPFR_RNDN); mpfr_init(r11089); mpfr_init(r11090); mpfr_init_set_str(r11091, "71.24974274308389", 10, MPFR_RNDN); mpfr_init_set_str(r11092, "7", 10, MPFR_RNDN); mpfr_init(r11093); mpfr_init(r11094); mpfr_init_set_str(r11095, "20.29199970278848", 10, MPFR_RNDN); mpfr_init_set_str(r11096, "5", 10, MPFR_RNDN); mpfr_init(r11097); mpfr_init(r11098); mpfr_init(r11099); mpfr_init(r11100); mpfr_init(r11101); mpfr_init_set_str(r11102, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11103, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11104); mpfr_init(r11105); mpfr_init(r11106); mpfr_init(r11107); mpfr_init_set_str(r11108, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11109); mpfr_init_set_str(r11110, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11111); mpfr_init(r11112); mpfr_init_set_str(r11113, "3", 10, MPFR_RNDN); mpfr_init(r11114); mpfr_init(r11115); mpfr_init_set_str(r11116, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11117); mpfr_init_set_str(r11118, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11119); mpfr_init(r11120); mpfr_init(r11121); mpfr_init(r11122); mpfr_init_set_str(r11123, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11124); mpfr_init(r11125); mpfr_init_set_str(r11126, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11127); mpfr_init(r11128); mpfr_init(r11129); mpfr_init(r11130); mpfr_init(r11131); mpfr_init_set_str(r11132, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r11133, "2", 10, MPFR_RNDN); mpfr_init(r11134); mpfr_init(r11135); mpfr_init(r11136); mpfr_init(r11137); mpfr_init_set_str(r11138, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r11139, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11140); mpfr_init(r11141); mpfr_init(r11142); mpfr_init(r11143); mpfr_init(r11144); mpfr_init(r11145); mpfr_init(r11146); mpfr_init(r11147); mpfr_init(r11148); } double f_dm(double x) { mpfr_set_d(r11073, x, MPFR_RNDN); ; mpfr_set_si(r11075, mpfr_cmp(r11073, r11074) <= 0, MPFR_RNDN); ; mpfr_set_si(r11077, mpfr_cmp(r11073, r11076) <= 0, MPFR_RNDN); mpfr_set_si(r11078, !mpfr_get_si(r11077, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11079, mpfr_get_si(r11075, MPFR_RNDN) || mpfr_get_si(r11078, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r11081, r11080, r11073, MPFR_RNDN); ; ; mpfr_pow(r11084, r11073, r11083, MPFR_RNDN); mpfr_div(r11085, r11082, r11084, MPFR_RNDN); mpfr_sub(r11086, r11081, r11085, MPFR_RNDN); ; ; mpfr_pow(r11089, r11073, r11088, MPFR_RNDN); mpfr_div(r11090, r11087, r11089, MPFR_RNDN); ; ; mpfr_pow(r11093, r11073, r11092, MPFR_RNDN); mpfr_div(r11094, r11091, r11093, MPFR_RNDN); ; ; mpfr_pow(r11097, r11073, r11096, MPFR_RNDN); mpfr_div(r11098, r11095, r11097, MPFR_RNDN); mpfr_add(r11099, r11094, r11098, MPFR_RNDN); mpfr_sub(r11100, r11090, r11099, MPFR_RNDN); mpfr_sub(r11101, r11086, r11100, MPFR_RNDN); ; ; mpfr_mul(r11104, r11073, r11073, MPFR_RNDN); mpfr_mul(r11105, r11103, r11104, MPFR_RNDN); mpfr_add(r11106, r11102, r11105, MPFR_RNDN); mpfr_mul(r11107, r11104, r11104, MPFR_RNDN); ; mpfr_mul(r11109, r11104, r11108, MPFR_RNDN); ; mpfr_add(r11111, r11109, r11110, MPFR_RNDN); mpfr_mul(r11112, r11107, r11111, MPFR_RNDN); ; mpfr_add(r11114, r11113, r11102, MPFR_RNDN); mpfr_pow(r11115, r11104, r11114, MPFR_RNDN); ; mpfr_mul(r11117, r11104, r11116, MPFR_RNDN); ; mpfr_add(r11119, r11117, r11118, MPFR_RNDN); mpfr_mul(r11120, r11115, r11119, MPFR_RNDN); mpfr_add(r11121, r11112, r11120, MPFR_RNDN); mpfr_add(r11122, r11106, r11121, MPFR_RNDN); ; mpfr_neg(r11124, r11123, MPFR_RNDN); mpfr_mul(r11125, r11107, r11124, MPFR_RNDN); ; mpfr_mul(r11127, r11104, r11126, MPFR_RNDN); mpfr_add(r11128, r11102, r11127, MPFR_RNDN); mpfr_sub(r11129, r11125, r11128, MPFR_RNDN); mpfr_pow(r11130, r11104, r11113, MPFR_RNDN); mpfr_mul(r11131, r11130, r11107, MPFR_RNDN); ; ; mpfr_mul(r11134, r11116, r11133, MPFR_RNDN); mpfr_mul(r11135, r11104, r11134, MPFR_RNDN); mpfr_add(r11136, r11132, r11135, MPFR_RNDN); mpfr_mul(r11137, r11131, r11136, MPFR_RNDN); ; ; mpfr_mul(r11140, r11104, r11139, MPFR_RNDN); mpfr_add(r11141, r11138, r11140, MPFR_RNDN); mpfr_mul(r11142, r11130, r11141, MPFR_RNDN); mpfr_add(r11143, r11137, r11142, MPFR_RNDN); mpfr_sub(r11144, r11129, r11143, MPFR_RNDN); mpfr_div(r11145, r11122, r11144, MPFR_RNDN); mpfr_neg(r11146, r11073, MPFR_RNDN); mpfr_mul(r11147, r11145, r11146, MPFR_RNDN); if (mpfr_get_si(r11079, MPFR_RNDN)) { mpfr_set(r11148, r11101, MPFR_RNDN); } else { mpfr_set(r11148, r11147, MPFR_RNDN); }; return mpfr_get_d(r11148, MPFR_RNDN); }