#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r10803 = 1.0f; float r10804 = 0.1049934947f; float r10805 = x; float r10806 = r10805 * r10805; float r10807 = r10804 * r10806; float r10808 = r10803 + r10807; float r10809 = 0.0424060604f; float r10810 = r10806 * r10806; float r10811 = r10809 * r10810; float r10812 = r10808 + r10811; float r10813 = 0.0072644182f; float r10814 = r10810 * r10806; float r10815 = r10813 * r10814; float r10816 = r10812 + r10815; float r10817 = 0.0005064034f; float r10818 = r10814 * r10806; float r10819 = r10817 * r10818; float r10820 = r10816 + r10819; float r10821 = 0.0001789971f; float r10822 = r10818 * r10806; float r10823 = r10821 * r10822; float r10824 = r10820 + r10823; float r10825 = 0.7715471019f; float r10826 = r10825 * r10806; float r10827 = r10803 + r10826; float r10828 = 0.2909738639f; float r10829 = r10828 * r10810; float r10830 = r10827 + r10829; float r10831 = 0.0694555761f; float r10832 = r10831 * r10814; float r10833 = r10830 + r10832; float r10834 = 0.0140005442f; float r10835 = r10834 * r10818; float r10836 = r10833 + r10835; float r10837 = 0.0008327945f; float r10838 = r10837 * r10822; float r10839 = r10836 + r10838; float r10840 = 2.0f; float r10841 = r10840 * r10821; float r10842 = r10822 * r10806; float r10843 = r10841 * r10842; float r10844 = r10839 + r10843; float r10845 = r10824 / r10844; float r10846 = r10845 * r10805; return r10846; } double f_id(double x) { double r10847 = 1.0; double r10848 = 0.1049934947; double r10849 = x; double r10850 = r10849 * r10849; double r10851 = r10848 * r10850; double r10852 = r10847 + r10851; double r10853 = 0.0424060604; double r10854 = r10850 * r10850; double r10855 = r10853 * r10854; double r10856 = r10852 + r10855; double r10857 = 0.0072644182; double r10858 = r10854 * r10850; double r10859 = r10857 * r10858; double r10860 = r10856 + r10859; double r10861 = 0.0005064034; double r10862 = r10858 * r10850; double r10863 = r10861 * r10862; double r10864 = r10860 + r10863; double r10865 = 0.0001789971; double r10866 = r10862 * r10850; double r10867 = r10865 * r10866; double r10868 = r10864 + r10867; double r10869 = 0.7715471019; double r10870 = r10869 * r10850; double r10871 = r10847 + r10870; double r10872 = 0.2909738639; double r10873 = r10872 * r10854; double r10874 = r10871 + r10873; double r10875 = 0.0694555761; double r10876 = r10875 * r10858; double r10877 = r10874 + r10876; double r10878 = 0.0140005442; double r10879 = r10878 * r10862; double r10880 = r10877 + r10879; double r10881 = 0.0008327945; double r10882 = r10881 * r10866; double r10883 = r10880 + r10882; double r10884 = 2.0; double r10885 = r10884 * r10865; double r10886 = r10866 * r10850; double r10887 = r10885 * r10886; double r10888 = r10883 + r10887; double r10889 = r10868 / r10888; double r10890 = r10889 * r10849; return r10890; } double f_of(float x) { float r10891 = 0.15298196345929327f; float r10892 = x; float r10893 = 5.0f; float r10894 = pow(r10892, r10893); float r10895 = r10891 / r10894; float r10896 = 1.0f; float r10897 = r10896 / r10892; float r10898 = 0.2514179000665375f; float r10899 = r10898 / r10892; float r10900 = 0.5f; float r10901 = fma(r10897, r10899, r10900); float r10902 = r10901 * r10897; float r10903 = r10895 + r10902; float r10904 = -0.0004036445015576861f; bool r10905 = r10903 <= r10904; float r10906 = 0.265709700396151f; float r10907 = 4.0f; float r10908 = pow(r10892, r10907); float r10909 = fma(r10906, r10908, r10896); float r10910 = 0.6665536072f; float r10911 = r10910 * r10892; float r10912 = r10892 * r10911; float r10913 = r10909 - r10912; float r10914 = r10892 * r10913; float r10915 = 8.03349330015686e-09f; bool r10916 = r10903 <= r10915; float r10917 = r10892 * r10892; float r10918 = 0.0001789971f; float r10919 = r10917 * r10918; float r10920 = r10917 * r10917; float r10921 = r10920 * r10920; float r10922 = r10919 * r10921; float r10923 = 0.0005064034f; float r10924 = 0.0072644182f; float r10925 = r10924 * r10917; float r10926 = r10925 * r10920; float r10927 = 0.0424060604f; float r10928 = 0.1049934947f; float r10929 = fma(r10928, r10917, r10896); float r10930 = fma(r10920, r10927, r10929); float r10931 = r10926 + r10930; float r10932 = fma(r10923, r10921, r10931); float r10933 = r10922 + r10932; float r10934 = 0.0140005442f; float r10935 = 0.2909738639f; float r10936 = 0.7715471019f; float r10937 = fma(r10917, r10936, r10896); float r10938 = fma(r10935, r10920, r10937); float r10939 = 0.0694555761f; float r10940 = r10939 * r10917; float r10941 = r10920 * r10940; float r10942 = r10938 + r10941; float r10943 = fma(r10921, r10934, r10942); float r10944 = 0.0008327945f; float r10945 = r10917 * r10944; float r10946 = r10945 * r10921; float r10947 = 2.0f; float r10948 = r10918 * r10947; float r10949 = r10948 * r10920; float r10950 = r10949 * r10921; float r10951 = r10946 + r10950; float r10952 = r10943 + r10951; float r10953 = r10933 / r10952; float r10954 = log(r10953); float r10955 = exp(r10954); float r10956 = r10892 * r10955; float r10957 = r10916 ? r10903 : r10956; float r10958 = r10905 ? r10914 : r10957; return r10958; } double f_od(double x) { double r10959 = 0.15298196345929327; double r10960 = x; double r10961 = 5.0; double r10962 = pow(r10960, r10961); double r10963 = r10959 / r10962; double r10964 = 1.0; double r10965 = r10964 / r10960; double r10966 = 0.2514179000665375; double r10967 = r10966 / r10960; double r10968 = 0.5; double r10969 = fma(r10965, r10967, r10968); double r10970 = r10969 * r10965; double r10971 = r10963 + r10970; double r10972 = -0.0004036445015576861; bool r10973 = r10971 <= r10972; double r10974 = 0.265709700396151; double r10975 = 4.0; double r10976 = pow(r10960, r10975); double r10977 = fma(r10974, r10976, r10964); double r10978 = 0.6665536072; double r10979 = r10978 * r10960; double r10980 = r10960 * r10979; double r10981 = r10977 - r10980; double r10982 = r10960 * r10981; double r10983 = 8.03349330015686e-09; bool r10984 = r10971 <= r10983; double r10985 = r10960 * r10960; double r10986 = 0.0001789971; double r10987 = r10985 * r10986; double r10988 = r10985 * r10985; double r10989 = r10988 * r10988; double r10990 = r10987 * r10989; double r10991 = 0.0005064034; double r10992 = 0.0072644182; double r10993 = r10992 * r10985; double r10994 = r10993 * r10988; double r10995 = 0.0424060604; double r10996 = 0.1049934947; double r10997 = fma(r10996, r10985, r10964); double r10998 = fma(r10988, r10995, r10997); double r10999 = r10994 + r10998; double r11000 = fma(r10991, r10989, r10999); double r11001 = r10990 + r11000; double r11002 = 0.0140005442; double r11003 = 0.2909738639; double r11004 = 0.7715471019; double r11005 = fma(r10985, r11004, r10964); double r11006 = fma(r11003, r10988, r11005); double r11007 = 0.0694555761; double r11008 = r11007 * r10985; double r11009 = r10988 * r11008; double r11010 = r11006 + r11009; double r11011 = fma(r10989, r11002, r11010); double r11012 = 0.0008327945; double r11013 = r10985 * r11012; double r11014 = r11013 * r10989; double r11015 = 2.0; double r11016 = r10986 * r11015; double r11017 = r11016 * r10988; double r11018 = r11017 * r10989; double r11019 = r11014 + r11018; double r11020 = r11011 + r11019; double r11021 = r11001 / r11020; double r11022 = log(r11021); double r11023 = exp(r11022); double r11024 = r10960 * r11023; double r11025 = r10984 ? r10971 : r11024; double r11026 = r10973 ? r10982 : r11025; return r11026; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r11027, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11028, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11029); mpfr_init(r11030); mpfr_init(r11031); mpfr_init(r11032); mpfr_init_set_str(r11033, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11034); mpfr_init(r11035); mpfr_init(r11036); mpfr_init_set_str(r11037, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11038); mpfr_init(r11039); mpfr_init(r11040); mpfr_init_set_str(r11041, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11042); mpfr_init(r11043); mpfr_init(r11044); mpfr_init_set_str(r11045, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11046); mpfr_init(r11047); mpfr_init(r11048); mpfr_init_set_str(r11049, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11050); mpfr_init(r11051); mpfr_init_set_str(r11052, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11053); mpfr_init(r11054); mpfr_init_set_str(r11055, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11056); mpfr_init(r11057); mpfr_init_set_str(r11058, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11059); mpfr_init(r11060); mpfr_init_set_str(r11061, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11062); mpfr_init(r11063); mpfr_init_set_str(r11064, "2", 10, MPFR_RNDN); mpfr_init(r11065); mpfr_init(r11066); mpfr_init(r11067); mpfr_init(r11068); mpfr_init(r11069); mpfr_init(r11070); } double f_im(double x) { ; ; mpfr_set_d(r11029, x, MPFR_RNDN); mpfr_mul(r11030, r11029, r11029, MPFR_RNDN); mpfr_mul(r11031, r11028, r11030, MPFR_RNDN); mpfr_add(r11032, r11027, r11031, MPFR_RNDN); ; mpfr_mul(r11034, r11030, r11030, MPFR_RNDN); mpfr_mul(r11035, r11033, r11034, MPFR_RNDN); mpfr_add(r11036, r11032, r11035, MPFR_RNDN); ; mpfr_mul(r11038, r11034, r11030, MPFR_RNDN); mpfr_mul(r11039, r11037, r11038, MPFR_RNDN); mpfr_add(r11040, r11036, r11039, MPFR_RNDN); ; mpfr_mul(r11042, r11038, r11030, MPFR_RNDN); mpfr_mul(r11043, r11041, r11042, MPFR_RNDN); mpfr_add(r11044, r11040, r11043, MPFR_RNDN); ; mpfr_mul(r11046, r11042, r11030, MPFR_RNDN); mpfr_mul(r11047, r11045, r11046, MPFR_RNDN); mpfr_add(r11048, r11044, r11047, MPFR_RNDN); ; mpfr_mul(r11050, r11049, r11030, MPFR_RNDN); mpfr_add(r11051, r11027, r11050, MPFR_RNDN); ; mpfr_mul(r11053, r11052, r11034, MPFR_RNDN); mpfr_add(r11054, r11051, r11053, MPFR_RNDN); ; mpfr_mul(r11056, r11055, r11038, MPFR_RNDN); mpfr_add(r11057, r11054, r11056, MPFR_RNDN); ; mpfr_mul(r11059, r11058, r11042, MPFR_RNDN); mpfr_add(r11060, r11057, r11059, MPFR_RNDN); ; mpfr_mul(r11062, r11061, r11046, MPFR_RNDN); mpfr_add(r11063, r11060, r11062, MPFR_RNDN); ; mpfr_mul(r11065, r11064, r11045, MPFR_RNDN); mpfr_mul(r11066, r11046, r11030, MPFR_RNDN); mpfr_mul(r11067, r11065, r11066, MPFR_RNDN); mpfr_add(r11068, r11063, r11067, MPFR_RNDN); mpfr_div(r11069, r11048, r11068, MPFR_RNDN); mpfr_mul(r11070, r11069, r11029, MPFR_RNDN); return mpfr_get_d(r11070, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r11071, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r11072); mpfr_init_set_str(r11073, "5", 10, MPFR_RNDN); mpfr_init(r11074); mpfr_init(r11075); mpfr_init_set_str(r11076, "1", 10, MPFR_RNDN); mpfr_init(r11077); mpfr_init_set_str(r11078, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r11079); mpfr_init_set_str(r11080, "0.5", 10, MPFR_RNDN); mpfr_init(r11081); mpfr_init(r11082); mpfr_init(r11083); mpfr_init_set_str(r11084, "-0.0004036445015576861", 10, MPFR_RNDN); mpfr_init(r11085); mpfr_init_set_str(r11086, "0.265709700396151", 10, MPFR_RNDN); mpfr_init_set_str(r11087, "4", 10, MPFR_RNDN); mpfr_init(r11088); mpfr_init(r11089); mpfr_init_set_str(r11090, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r11091); mpfr_init(r11092); mpfr_init(r11093); mpfr_init(r11094); mpfr_init_set_str(r11095, "8.03349330015686e-09", 10, MPFR_RNDN); mpfr_init(r11096); mpfr_init(r11097); mpfr_init_set_str(r11098, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11099); mpfr_init(r11100); mpfr_init(r11101); mpfr_init(r11102); mpfr_init_set_str(r11103, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r11104, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11105); mpfr_init(r11106); mpfr_init_set_str(r11107, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r11108, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11109); mpfr_init(r11110); mpfr_init(r11111); mpfr_init(r11112); mpfr_init(r11113); mpfr_init_set_str(r11114, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r11115, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r11116, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11117); mpfr_init(r11118); mpfr_init_set_str(r11119, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11120); mpfr_init(r11121); mpfr_init(r11122); mpfr_init(r11123); mpfr_init_set_str(r11124, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11125); mpfr_init(r11126); mpfr_init_set_str(r11127, "2", 10, MPFR_RNDN); mpfr_init(r11128); mpfr_init(r11129); mpfr_init(r11130); mpfr_init(r11131); mpfr_init(r11132); mpfr_init(r11133); mpfr_init(r11134); mpfr_init(r11135); mpfr_init(r11136); mpfr_init(r11137); mpfr_init(r11138); } double f_fm(double x) { ; mpfr_set_d(r11072, x, MPFR_RNDN); ; mpfr_pow(r11074, r11072, r11073, MPFR_RNDN); mpfr_div(r11075, r11071, r11074, MPFR_RNDN); ; mpfr_div(r11077, r11076, r11072, MPFR_RNDN); ; mpfr_div(r11079, r11078, r11072, MPFR_RNDN); ; mpfr_fma(r11081, r11077, r11079, r11080, MPFR_RNDN); mpfr_mul(r11082, r11081, r11077, MPFR_RNDN); mpfr_add(r11083, r11075, r11082, MPFR_RNDN); ; mpfr_set_si(r11085, mpfr_cmp(r11083, r11084) <= 0, MPFR_RNDN); ; ; mpfr_pow(r11088, r11072, r11087, MPFR_RNDN); mpfr_fma(r11089, r11086, r11088, r11076, MPFR_RNDN); ; mpfr_mul(r11091, r11090, r11072, MPFR_RNDN); mpfr_mul(r11092, r11072, r11091, MPFR_RNDN); mpfr_sub(r11093, r11089, r11092, MPFR_RNDN); mpfr_mul(r11094, r11072, r11093, MPFR_RNDN); ; mpfr_set_si(r11096, mpfr_cmp(r11083, r11095) <= 0, MPFR_RNDN); mpfr_mul(r11097, r11072, r11072, MPFR_RNDN); ; mpfr_mul(r11099, r11097, r11098, MPFR_RNDN); mpfr_mul(r11100, r11097, r11097, MPFR_RNDN); mpfr_mul(r11101, r11100, r11100, MPFR_RNDN); mpfr_mul(r11102, r11099, r11101, MPFR_RNDN); ; ; mpfr_mul(r11105, r11104, r11097, MPFR_RNDN); mpfr_mul(r11106, r11105, r11100, MPFR_RNDN); ; ; mpfr_fma(r11109, r11108, r11097, r11076, MPFR_RNDN); mpfr_fma(r11110, r11100, r11107, r11109, MPFR_RNDN); mpfr_add(r11111, r11106, r11110, MPFR_RNDN); mpfr_fma(r11112, r11103, r11101, r11111, MPFR_RNDN); mpfr_add(r11113, r11102, r11112, MPFR_RNDN); ; ; ; mpfr_fma(r11117, r11097, r11116, r11076, MPFR_RNDN); mpfr_fma(r11118, r11115, r11100, r11117, MPFR_RNDN); ; mpfr_mul(r11120, r11119, r11097, MPFR_RNDN); mpfr_mul(r11121, r11100, r11120, MPFR_RNDN); mpfr_add(r11122, r11118, r11121, MPFR_RNDN); mpfr_fma(r11123, r11101, r11114, r11122, MPFR_RNDN); ; mpfr_mul(r11125, r11097, r11124, MPFR_RNDN); mpfr_mul(r11126, r11125, r11101, MPFR_RNDN); ; mpfr_mul(r11128, r11098, r11127, MPFR_RNDN); mpfr_mul(r11129, r11128, r11100, MPFR_RNDN); mpfr_mul(r11130, r11129, r11101, MPFR_RNDN); mpfr_add(r11131, r11126, r11130, MPFR_RNDN); mpfr_add(r11132, r11123, r11131, MPFR_RNDN); mpfr_div(r11133, r11113, r11132, MPFR_RNDN); mpfr_log(r11134, r11133, MPFR_RNDN); mpfr_exp(r11135, r11134, MPFR_RNDN); mpfr_mul(r11136, r11072, r11135, MPFR_RNDN); if (mpfr_get_si(r11096, MPFR_RNDN)) { mpfr_set(r11137, r11083, MPFR_RNDN); } else { mpfr_set(r11137, r11136, MPFR_RNDN); }; if (mpfr_get_si(r11085, MPFR_RNDN)) { mpfr_set(r11138, r11094, MPFR_RNDN); } else { mpfr_set(r11138, r11137, MPFR_RNDN); }; return mpfr_get_d(r11138, MPFR_RNDN); } static mpfr_t r11139, r11140, r11141, r11142, r11143, r11144, r11145, r11146, r11147, r11148, r11149, r11150, r11151, r11152, r11153, r11154, r11155, r11156, r11157, r11158, r11159, r11160, r11161, r11162, r11163, r11164, r11165, r11166, r11167, r11168, r11169, r11170, r11171, r11172, r11173, r11174, r11175, r11176, r11177, r11178, r11179, r11180, r11181, r11182, r11183, r11184, r11185, r11186, r11187, r11188, r11189, r11190, r11191, r11192, r11193, r11194, r11195, r11196, r11197, r11198, r11199, r11200, r11201, r11202, r11203, r11204, r11205, r11206; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r11139, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r11140); mpfr_init_set_str(r11141, "5", 10, MPFR_RNDN); mpfr_init(r11142); mpfr_init(r11143); mpfr_init_set_str(r11144, "1", 10, MPFR_RNDN); mpfr_init(r11145); mpfr_init_set_str(r11146, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r11147); mpfr_init_set_str(r11148, "0.5", 10, MPFR_RNDN); mpfr_init(r11149); mpfr_init(r11150); mpfr_init(r11151); mpfr_init_set_str(r11152, "-0.0004036445015576861", 10, MPFR_RNDN); mpfr_init(r11153); mpfr_init_set_str(r11154, "0.265709700396151", 10, MPFR_RNDN); mpfr_init_set_str(r11155, "4", 10, MPFR_RNDN); mpfr_init(r11156); mpfr_init(r11157); mpfr_init_set_str(r11158, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r11159); mpfr_init(r11160); mpfr_init(r11161); mpfr_init(r11162); mpfr_init_set_str(r11163, "8.03349330015686e-09", 10, MPFR_RNDN); mpfr_init(r11164); mpfr_init(r11165); mpfr_init_set_str(r11166, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11167); mpfr_init(r11168); mpfr_init(r11169); mpfr_init(r11170); mpfr_init_set_str(r11171, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r11172, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11173); mpfr_init(r11174); mpfr_init_set_str(r11175, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r11176, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11177); mpfr_init(r11178); mpfr_init(r11179); mpfr_init(r11180); mpfr_init(r11181); mpfr_init_set_str(r11182, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r11183, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r11184, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11185); mpfr_init(r11186); mpfr_init_set_str(r11187, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11188); mpfr_init(r11189); mpfr_init(r11190); mpfr_init(r11191); mpfr_init_set_str(r11192, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11193); mpfr_init(r11194); mpfr_init_set_str(r11195, "2", 10, MPFR_RNDN); mpfr_init(r11196); mpfr_init(r11197); mpfr_init(r11198); mpfr_init(r11199); mpfr_init(r11200); mpfr_init(r11201); mpfr_init(r11202); mpfr_init(r11203); mpfr_init(r11204); mpfr_init(r11205); mpfr_init(r11206); } double f_dm(double x) { ; mpfr_set_d(r11140, x, MPFR_RNDN); ; mpfr_pow(r11142, r11140, r11141, MPFR_RNDN); mpfr_div(r11143, r11139, r11142, MPFR_RNDN); ; mpfr_div(r11145, r11144, r11140, MPFR_RNDN); ; mpfr_div(r11147, r11146, r11140, MPFR_RNDN); ; mpfr_fma(r11149, r11145, r11147, r11148, MPFR_RNDN); mpfr_mul(r11150, r11149, r11145, MPFR_RNDN); mpfr_add(r11151, r11143, r11150, MPFR_RNDN); ; mpfr_set_si(r11153, mpfr_cmp(r11151, r11152) <= 0, MPFR_RNDN); ; ; mpfr_pow(r11156, r11140, r11155, MPFR_RNDN); mpfr_fma(r11157, r11154, r11156, r11144, MPFR_RNDN); ; mpfr_mul(r11159, r11158, r11140, MPFR_RNDN); mpfr_mul(r11160, r11140, r11159, MPFR_RNDN); mpfr_sub(r11161, r11157, r11160, MPFR_RNDN); mpfr_mul(r11162, r11140, r11161, MPFR_RNDN); ; mpfr_set_si(r11164, mpfr_cmp(r11151, r11163) <= 0, MPFR_RNDN); mpfr_mul(r11165, r11140, r11140, MPFR_RNDN); ; mpfr_mul(r11167, r11165, r11166, MPFR_RNDN); mpfr_mul(r11168, r11165, r11165, MPFR_RNDN); mpfr_mul(r11169, r11168, r11168, MPFR_RNDN); mpfr_mul(r11170, r11167, r11169, MPFR_RNDN); ; ; mpfr_mul(r11173, r11172, r11165, MPFR_RNDN); mpfr_mul(r11174, r11173, r11168, MPFR_RNDN); ; ; mpfr_fma(r11177, r11176, r11165, r11144, MPFR_RNDN); mpfr_fma(r11178, r11168, r11175, r11177, MPFR_RNDN); mpfr_add(r11179, r11174, r11178, MPFR_RNDN); mpfr_fma(r11180, r11171, r11169, r11179, MPFR_RNDN); mpfr_add(r11181, r11170, r11180, MPFR_RNDN); ; ; ; mpfr_fma(r11185, r11165, r11184, r11144, MPFR_RNDN); mpfr_fma(r11186, r11183, r11168, r11185, MPFR_RNDN); ; mpfr_mul(r11188, r11187, r11165, MPFR_RNDN); mpfr_mul(r11189, r11168, r11188, MPFR_RNDN); mpfr_add(r11190, r11186, r11189, MPFR_RNDN); mpfr_fma(r11191, r11169, r11182, r11190, MPFR_RNDN); ; mpfr_mul(r11193, r11165, r11192, MPFR_RNDN); mpfr_mul(r11194, r11193, r11169, MPFR_RNDN); ; mpfr_mul(r11196, r11166, r11195, MPFR_RNDN); mpfr_mul(r11197, r11196, r11168, MPFR_RNDN); mpfr_mul(r11198, r11197, r11169, MPFR_RNDN); mpfr_add(r11199, r11194, r11198, MPFR_RNDN); mpfr_add(r11200, r11191, r11199, MPFR_RNDN); mpfr_div(r11201, r11181, r11200, MPFR_RNDN); mpfr_log(r11202, r11201, MPFR_RNDN); mpfr_exp(r11203, r11202, MPFR_RNDN); mpfr_mul(r11204, r11140, r11203, MPFR_RNDN); if (mpfr_get_si(r11164, MPFR_RNDN)) { mpfr_set(r11205, r11151, MPFR_RNDN); } else { mpfr_set(r11205, r11204, MPFR_RNDN); }; if (mpfr_get_si(r11153, MPFR_RNDN)) { mpfr_set(r11206, r11162, MPFR_RNDN); } else { mpfr_set(r11206, r11205, MPFR_RNDN); }; return mpfr_get_d(r11206, MPFR_RNDN); }