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