#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r27778 = 1; float r27779 = 0.1049934947; float r27780 = x; float r27781 = r27780 * r27780; float r27782 = r27779 * r27781; float r27783 = r27778 + r27782; float r27784 = 0.0424060604; float r27785 = r27781 * r27781; float r27786 = r27784 * r27785; float r27787 = r27783 + r27786; float r27788 = 0.0072644182; float r27789 = r27785 * r27781; float r27790 = r27788 * r27789; float r27791 = r27787 + r27790; float r27792 = 0.0005064034; float r27793 = r27789 * r27781; float r27794 = r27792 * r27793; float r27795 = r27791 + r27794; float r27796 = 0.0001789971; float r27797 = r27793 * r27781; float r27798 = r27796 * r27797; float r27799 = r27795 + r27798; float r27800 = 0.7715471019; float r27801 = r27800 * r27781; float r27802 = r27778 + r27801; float r27803 = 0.2909738639; float r27804 = r27803 * r27785; float r27805 = r27802 + r27804; float r27806 = 0.0694555761; float r27807 = r27806 * r27789; float r27808 = r27805 + r27807; float r27809 = 0.0140005442; float r27810 = r27809 * r27793; float r27811 = r27808 + r27810; float r27812 = 0.0008327945; float r27813 = r27812 * r27797; float r27814 = r27811 + r27813; float r27815 = 2; float r27816 = r27815 * r27796; float r27817 = r27797 * r27781; float r27818 = r27816 * r27817; float r27819 = r27814 + r27818; float r27820 = r27799 / r27819; float r27821 = r27820 * r27780; return r27821; } double f_id(double x) { double r27822 = 1; double r27823 = 0.1049934947; double r27824 = x; double r27825 = r27824 * r27824; double r27826 = r27823 * r27825; double r27827 = r27822 + r27826; double r27828 = 0.0424060604; double r27829 = r27825 * r27825; double r27830 = r27828 * r27829; double r27831 = r27827 + r27830; double r27832 = 0.0072644182; double r27833 = r27829 * r27825; double r27834 = r27832 * r27833; double r27835 = r27831 + r27834; double r27836 = 0.0005064034; double r27837 = r27833 * r27825; double r27838 = r27836 * r27837; double r27839 = r27835 + r27838; double r27840 = 0.0001789971; double r27841 = r27837 * r27825; double r27842 = r27840 * r27841; double r27843 = r27839 + r27842; double r27844 = 0.7715471019; double r27845 = r27844 * r27825; double r27846 = r27822 + r27845; double r27847 = 0.2909738639; double r27848 = r27847 * r27829; double r27849 = r27846 + r27848; double r27850 = 0.0694555761; double r27851 = r27850 * r27833; double r27852 = r27849 + r27851; double r27853 = 0.0140005442; double r27854 = r27853 * r27837; double r27855 = r27852 + r27854; double r27856 = 0.0008327945; double r27857 = r27856 * r27841; double r27858 = r27855 + r27857; double r27859 = 2; double r27860 = r27859 * r27840; double r27861 = r27841 * r27825; double r27862 = r27860 * r27861; double r27863 = r27858 + r27862; double r27864 = r27843 / r27863; double r27865 = r27864 * r27824; return r27865; } double f_of(float x) { float r27866 = 0.5; float r27867 = x; float r27868 = r27866 / r27867; float r27869 = 0.2514179000665375; float r27870 = r27869 * r27867; float r27871 = 4; float r27872 = pow(r27867, r27871); float r27873 = r27870 / r27872; float r27874 = 0.15298196345929327; float r27875 = r27874 * r27867; float r27876 = 6; float r27877 = pow(r27867, r27876); float r27878 = r27875 / r27877; float r27879 = r27873 + r27878; float r27880 = r27868 + r27879; float r27881 = -2.0383110714971574e-09; bool r27882 = r27880 <= r27881; float r27883 = 1; float r27884 = 0.1049934947; float r27885 = r27867 * r27867; float r27886 = r27884 * r27885; float r27887 = r27883 + r27886; float r27888 = 0.0424060604; float r27889 = r27885 * r27885; float r27890 = r27888 * r27889; float r27891 = r27887 + r27890; float r27892 = 0.0072644182; float r27893 = r27889 * r27885; float r27894 = r27892 * r27893; float r27895 = r27891 + r27894; float r27896 = 0.0005064034; float r27897 = r27893 * r27885; float r27898 = r27896 * r27897; float r27899 = r27895 + r27898; float r27900 = 0.0001789971; float r27901 = r27897 * r27885; float r27902 = r27900 * r27901; float r27903 = r27899 + r27902; float r27904 = 0.7715471019; float r27905 = r27904 * r27885; float r27906 = r27883 + r27905; float r27907 = 0.2909738639; float r27908 = r27907 * r27889; float r27909 = r27906 + r27908; float r27910 = 0.0694555761; float r27911 = r27910 * r27893; float r27912 = r27909 + r27911; float r27913 = 0.0140005442; float r27914 = r27913 * r27897; float r27915 = r27912 + r27914; float r27916 = 0.0008327945; float r27917 = r27916 * r27901; float r27918 = r27915 + r27917; float r27919 = 2; float r27920 = r27919 * r27900; float r27921 = r27901 * r27885; float r27922 = r27920 * r27921; float r27923 = r27918 + r27922; float r27924 = r27903 / r27923; float r27925 = r27924 * r27867; float r27926 = 1.903280832029775e-24; bool r27927 = r27880 <= r27926; float r27928 = r27885 * r27884; float r27929 = r27928 + r27883; float r27930 = r27867 * r27892; float r27931 = r27930 * r27867; float r27932 = r27888 + r27931; float r27933 = r27889 * r27932; float r27934 = 3; float r27935 = r27934 + r27883; float r27936 = pow(r27885, r27935); float r27937 = r27885 * r27900; float r27938 = r27937 + r27896; float r27939 = r27936 * r27938; float r27940 = r27933 + r27939; float r27941 = r27929 + r27940; float r27942 = r27907 * r27867; float r27943 = pow(r27867, r27934); float r27944 = r27942 * r27943; float r27945 = r27883 + r27944; float r27946 = r27904 * r27867; float r27947 = r27946 * r27867; float r27948 = r27945 + r27947; float r27949 = r27943 * r27885; float r27950 = r27949 * r27949; float r27951 = r27900 * r27919; float r27952 = r27951 * r27885; float r27953 = r27916 + r27952; float r27954 = r27950 * r27953; float r27955 = r27943 * r27943; float r27956 = r27867 * r27913; float r27957 = r27956 * r27867; float r27958 = r27957 + r27910; float r27959 = r27955 * r27958; float r27960 = r27954 + r27959; float r27961 = r27948 + r27960; float r27962 = r27941 / r27961; float r27963 = exp(r27962); float r27964 = log(r27963); float r27965 = r27964 * r27867; float r27966 = r27927 ? r27880 : r27965; float r27967 = r27882 ? r27925 : r27966; return r27967; } double f_od(double x) { double r27968 = 0.5; double r27969 = x; double r27970 = r27968 / r27969; double r27971 = 0.2514179000665375; double r27972 = r27971 * r27969; double r27973 = 4; double r27974 = pow(r27969, r27973); double r27975 = r27972 / r27974; double r27976 = 0.15298196345929327; double r27977 = r27976 * r27969; double r27978 = 6; double r27979 = pow(r27969, r27978); double r27980 = r27977 / r27979; double r27981 = r27975 + r27980; double r27982 = r27970 + r27981; double r27983 = -2.0383110714971574e-09; bool r27984 = r27982 <= r27983; double r27985 = 1; double r27986 = 0.1049934947; double r27987 = r27969 * r27969; double r27988 = r27986 * r27987; double r27989 = r27985 + r27988; double r27990 = 0.0424060604; double r27991 = r27987 * r27987; double r27992 = r27990 * r27991; double r27993 = r27989 + r27992; double r27994 = 0.0072644182; double r27995 = r27991 * r27987; double r27996 = r27994 * r27995; double r27997 = r27993 + r27996; double r27998 = 0.0005064034; double r27999 = r27995 * r27987; double r28000 = r27998 * r27999; double r28001 = r27997 + r28000; double r28002 = 0.0001789971; double r28003 = r27999 * r27987; double r28004 = r28002 * r28003; double r28005 = r28001 + r28004; double r28006 = 0.7715471019; double r28007 = r28006 * r27987; double r28008 = r27985 + r28007; double r28009 = 0.2909738639; double r28010 = r28009 * r27991; double r28011 = r28008 + r28010; double r28012 = 0.0694555761; double r28013 = r28012 * r27995; double r28014 = r28011 + r28013; double r28015 = 0.0140005442; double r28016 = r28015 * r27999; double r28017 = r28014 + r28016; double r28018 = 0.0008327945; double r28019 = r28018 * r28003; double r28020 = r28017 + r28019; double r28021 = 2; double r28022 = r28021 * r28002; double r28023 = r28003 * r27987; double r28024 = r28022 * r28023; double r28025 = r28020 + r28024; double r28026 = r28005 / r28025; double r28027 = r28026 * r27969; double r28028 = 1.903280832029775e-24; bool r28029 = r27982 <= r28028; double r28030 = r27987 * r27986; double r28031 = r28030 + r27985; double r28032 = r27969 * r27994; double r28033 = r28032 * r27969; double r28034 = r27990 + r28033; double r28035 = r27991 * r28034; double r28036 = 3; double r28037 = r28036 + r27985; double r28038 = pow(r27987, r28037); double r28039 = r27987 * r28002; double r28040 = r28039 + r27998; double r28041 = r28038 * r28040; double r28042 = r28035 + r28041; double r28043 = r28031 + r28042; double r28044 = r28009 * r27969; double r28045 = pow(r27969, r28036); double r28046 = r28044 * r28045; double r28047 = r27985 + r28046; double r28048 = r28006 * r27969; double r28049 = r28048 * r27969; double r28050 = r28047 + r28049; double r28051 = r28045 * r27987; double r28052 = r28051 * r28051; double r28053 = r28002 * r28021; double r28054 = r28053 * r27987; double r28055 = r28018 + r28054; double r28056 = r28052 * r28055; double r28057 = r28045 * r28045; double r28058 = r27969 * r28015; double r28059 = r28058 * r27969; double r28060 = r28059 + r28012; double r28061 = r28057 * r28060; double r28062 = r28056 + r28061; double r28063 = r28050 + r28062; double r28064 = r28043 / r28063; double r28065 = exp(r28064); double r28066 = log(r28065); double r28067 = r28066 * r27969; double r28068 = r28029 ? r27982 : r28067; double r28069 = r27984 ? r28027 : r28068; return r28069; } 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 r28070, r28071, r28072, r28073, r28074, r28075, r28076, r28077, r28078, r28079, r28080, r28081, r28082, r28083, r28084, r28085, r28086, r28087, r28088, r28089, r28090, r28091, r28092, r28093, r28094, r28095, r28096, r28097, r28098, r28099, r28100, r28101, r28102, r28103, r28104, r28105, r28106, r28107, r28108, r28109, r28110, r28111, r28112, r28113; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r28070, "1", 10, MPFR_RNDN); mpfr_init_set_str(r28071, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r28072); mpfr_init(r28073); mpfr_init(r28074); mpfr_init(r28075); mpfr_init_set_str(r28076, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r28077); mpfr_init(r28078); mpfr_init(r28079); mpfr_init_set_str(r28080, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r28081); mpfr_init(r28082); mpfr_init(r28083); mpfr_init_set_str(r28084, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r28085); mpfr_init(r28086); mpfr_init(r28087); mpfr_init_set_str(r28088, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28089); mpfr_init(r28090); mpfr_init(r28091); mpfr_init_set_str(r28092, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28093); mpfr_init(r28094); mpfr_init_set_str(r28095, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r28096); mpfr_init(r28097); mpfr_init_set_str(r28098, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28099); mpfr_init(r28100); mpfr_init_set_str(r28101, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r28102); mpfr_init(r28103); mpfr_init_set_str(r28104, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28105); mpfr_init(r28106); mpfr_init_set_str(r28107, "2", 10, MPFR_RNDN); mpfr_init(r28108); mpfr_init(r28109); mpfr_init(r28110); mpfr_init(r28111); mpfr_init(r28112); mpfr_init(r28113); } double f_im(double x) { ; ; mpfr_set_d(r28072, x, MPFR_RNDN); mpfr_mul(r28073, r28072, r28072, MPFR_RNDN); mpfr_mul(r28074, r28071, r28073, MPFR_RNDN); mpfr_add(r28075, r28070, r28074, MPFR_RNDN); ; mpfr_mul(r28077, r28073, r28073, MPFR_RNDN); mpfr_mul(r28078, r28076, r28077, MPFR_RNDN); mpfr_add(r28079, r28075, r28078, MPFR_RNDN); ; mpfr_mul(r28081, r28077, r28073, MPFR_RNDN); mpfr_mul(r28082, r28080, r28081, MPFR_RNDN); mpfr_add(r28083, r28079, r28082, MPFR_RNDN); ; mpfr_mul(r28085, r28081, r28073, MPFR_RNDN); mpfr_mul(r28086, r28084, r28085, MPFR_RNDN); mpfr_add(r28087, r28083, r28086, MPFR_RNDN); ; mpfr_mul(r28089, r28085, r28073, MPFR_RNDN); mpfr_mul(r28090, r28088, r28089, MPFR_RNDN); mpfr_add(r28091, r28087, r28090, MPFR_RNDN); ; mpfr_mul(r28093, r28092, r28073, MPFR_RNDN); mpfr_add(r28094, r28070, r28093, MPFR_RNDN); ; mpfr_mul(r28096, r28095, r28077, MPFR_RNDN); mpfr_add(r28097, r28094, r28096, MPFR_RNDN); ; mpfr_mul(r28099, r28098, r28081, MPFR_RNDN); mpfr_add(r28100, r28097, r28099, MPFR_RNDN); ; mpfr_mul(r28102, r28101, r28085, MPFR_RNDN); mpfr_add(r28103, r28100, r28102, MPFR_RNDN); ; mpfr_mul(r28105, r28104, r28089, MPFR_RNDN); mpfr_add(r28106, r28103, r28105, MPFR_RNDN); ; mpfr_mul(r28108, r28107, r28088, MPFR_RNDN); mpfr_mul(r28109, r28089, r28073, MPFR_RNDN); mpfr_mul(r28110, r28108, r28109, MPFR_RNDN); mpfr_add(r28111, r28106, r28110, MPFR_RNDN); mpfr_div(r28112, r28091, r28111, MPFR_RNDN); mpfr_mul(r28113, r28112, r28072, MPFR_RNDN); return mpfr_get_d(r28113, MPFR_RNDN); } static mpfr_t r28114, r28115, r28116, r28117, r28118, r28119, r28120, r28121, r28122, r28123, r28124, r28125, r28126, r28127, r28128, r28129, r28130, r28131, r28132, r28133, r28134, r28135, r28136, r28137, r28138, r28139, r28140, r28141, r28142, r28143, r28144, r28145, r28146, r28147, r28148, r28149, r28150, r28151, r28152, r28153, r28154, r28155, r28156, r28157, r28158, r28159, r28160, r28161, r28162, r28163, r28164, r28165, r28166, r28167, r28168, r28169, r28170, r28171, r28172, r28173, r28174, r28175, r28176, r28177, r28178, r28179, r28180, r28181, r28182, r28183, r28184, r28185, r28186, r28187, r28188, r28189, r28190, r28191, r28192, r28193, r28194, r28195, r28196, r28197, r28198, r28199, r28200, r28201, r28202, r28203, r28204, r28205, r28206, r28207, r28208, r28209, r28210, r28211, r28212, r28213, r28214, r28215; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r28114, "0.5", 10, MPFR_RNDN); mpfr_init(r28115); mpfr_init(r28116); mpfr_init_set_str(r28117, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r28118); mpfr_init_set_str(r28119, "4", 10, MPFR_RNDN); mpfr_init(r28120); mpfr_init(r28121); mpfr_init_set_str(r28122, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r28123); mpfr_init_set_str(r28124, "6", 10, MPFR_RNDN); mpfr_init(r28125); mpfr_init(r28126); mpfr_init(r28127); mpfr_init(r28128); mpfr_init_set_str(r28129, "-2.0383110714971574e-09", 10, MPFR_RNDN); mpfr_init(r28130); mpfr_init_set_str(r28131, "1", 10, MPFR_RNDN); mpfr_init_set_str(r28132, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r28133); mpfr_init(r28134); mpfr_init(r28135); mpfr_init_set_str(r28136, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r28137); mpfr_init(r28138); mpfr_init(r28139); mpfr_init_set_str(r28140, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r28141); mpfr_init(r28142); mpfr_init(r28143); mpfr_init_set_str(r28144, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r28145); mpfr_init(r28146); mpfr_init(r28147); mpfr_init_set_str(r28148, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28149); mpfr_init(r28150); mpfr_init(r28151); mpfr_init_set_str(r28152, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28153); mpfr_init(r28154); mpfr_init_set_str(r28155, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r28156); mpfr_init(r28157); mpfr_init_set_str(r28158, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28159); mpfr_init(r28160); mpfr_init_set_str(r28161, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r28162); mpfr_init(r28163); mpfr_init_set_str(r28164, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28165); mpfr_init(r28166); mpfr_init_set_str(r28167, "2", 10, MPFR_RNDN); mpfr_init(r28168); mpfr_init(r28169); mpfr_init(r28170); mpfr_init(r28171); mpfr_init(r28172); mpfr_init(r28173); mpfr_init_set_str(r28174, "1.903280832029775e-24", 10, MPFR_RNDN); mpfr_init(r28175); mpfr_init(r28176); mpfr_init(r28177); mpfr_init(r28178); mpfr_init(r28179); mpfr_init(r28180); mpfr_init(r28181); mpfr_init_set_str(r28182, "3", 10, MPFR_RNDN); mpfr_init(r28183); mpfr_init(r28184); mpfr_init(r28185); mpfr_init(r28186); mpfr_init(r28187); mpfr_init(r28188); mpfr_init(r28189); mpfr_init(r28190); mpfr_init(r28191); mpfr_init(r28192); mpfr_init(r28193); mpfr_init(r28194); mpfr_init(r28195); mpfr_init(r28196); mpfr_init(r28197); mpfr_init(r28198); mpfr_init(r28199); mpfr_init(r28200); mpfr_init(r28201); mpfr_init(r28202); mpfr_init(r28203); mpfr_init(r28204); mpfr_init(r28205); mpfr_init(r28206); mpfr_init(r28207); mpfr_init(r28208); mpfr_init(r28209); mpfr_init(r28210); mpfr_init(r28211); mpfr_init(r28212); mpfr_init(r28213); mpfr_init(r28214); mpfr_init(r28215); } double f_fm(double x) { ; mpfr_set_d(r28115, x, MPFR_RNDN); mpfr_div(r28116, r28114, r28115, MPFR_RNDN); ; mpfr_mul(r28118, r28117, r28115, MPFR_RNDN); ; mpfr_pow(r28120, r28115, r28119, MPFR_RNDN); mpfr_div(r28121, r28118, r28120, MPFR_RNDN); ; mpfr_mul(r28123, r28122, r28115, MPFR_RNDN); ; mpfr_pow(r28125, r28115, r28124, MPFR_RNDN); mpfr_div(r28126, r28123, r28125, MPFR_RNDN); mpfr_add(r28127, r28121, r28126, MPFR_RNDN); mpfr_add(r28128, r28116, r28127, MPFR_RNDN); ; mpfr_set_si(r28130, mpfr_cmp(r28128, r28129) <= 0, MPFR_RNDN); ; ; mpfr_mul(r28133, r28115, r28115, MPFR_RNDN); mpfr_mul(r28134, r28132, r28133, MPFR_RNDN); mpfr_add(r28135, r28131, r28134, MPFR_RNDN); ; mpfr_mul(r28137, r28133, r28133, MPFR_RNDN); mpfr_mul(r28138, r28136, r28137, MPFR_RNDN); mpfr_add(r28139, r28135, r28138, MPFR_RNDN); ; mpfr_mul(r28141, r28137, r28133, MPFR_RNDN); mpfr_mul(r28142, r28140, r28141, MPFR_RNDN); mpfr_add(r28143, r28139, r28142, MPFR_RNDN); ; mpfr_mul(r28145, r28141, r28133, MPFR_RNDN); mpfr_mul(r28146, r28144, r28145, MPFR_RNDN); mpfr_add(r28147, r28143, r28146, MPFR_RNDN); ; mpfr_mul(r28149, r28145, r28133, MPFR_RNDN); mpfr_mul(r28150, r28148, r28149, MPFR_RNDN); mpfr_add(r28151, r28147, r28150, MPFR_RNDN); ; mpfr_mul(r28153, r28152, r28133, MPFR_RNDN); mpfr_add(r28154, r28131, r28153, MPFR_RNDN); ; mpfr_mul(r28156, r28155, r28137, MPFR_RNDN); mpfr_add(r28157, r28154, r28156, MPFR_RNDN); ; mpfr_mul(r28159, r28158, r28141, MPFR_RNDN); mpfr_add(r28160, r28157, r28159, MPFR_RNDN); ; mpfr_mul(r28162, r28161, r28145, MPFR_RNDN); mpfr_add(r28163, r28160, r28162, MPFR_RNDN); ; mpfr_mul(r28165, r28164, r28149, MPFR_RNDN); mpfr_add(r28166, r28163, r28165, MPFR_RNDN); ; mpfr_mul(r28168, r28167, r28148, MPFR_RNDN); mpfr_mul(r28169, r28149, r28133, MPFR_RNDN); mpfr_mul(r28170, r28168, r28169, MPFR_RNDN); mpfr_add(r28171, r28166, r28170, MPFR_RNDN); mpfr_div(r28172, r28151, r28171, MPFR_RNDN); mpfr_mul(r28173, r28172, r28115, MPFR_RNDN); ; mpfr_set_si(r28175, mpfr_cmp(r28128, r28174) <= 0, MPFR_RNDN); mpfr_mul(r28176, r28133, r28132, MPFR_RNDN); mpfr_add(r28177, r28176, r28131, MPFR_RNDN); mpfr_mul(r28178, r28115, r28140, MPFR_RNDN); mpfr_mul(r28179, r28178, r28115, MPFR_RNDN); mpfr_add(r28180, r28136, r28179, MPFR_RNDN); mpfr_mul(r28181, r28137, r28180, MPFR_RNDN); ; mpfr_add(r28183, r28182, r28131, MPFR_RNDN); mpfr_pow(r28184, r28133, r28183, MPFR_RNDN); mpfr_mul(r28185, r28133, r28148, MPFR_RNDN); mpfr_add(r28186, r28185, r28144, MPFR_RNDN); mpfr_mul(r28187, r28184, r28186, MPFR_RNDN); mpfr_add(r28188, r28181, r28187, MPFR_RNDN); mpfr_add(r28189, r28177, r28188, MPFR_RNDN); mpfr_mul(r28190, r28155, r28115, MPFR_RNDN); mpfr_pow(r28191, r28115, r28182, MPFR_RNDN); mpfr_mul(r28192, r28190, r28191, MPFR_RNDN); mpfr_add(r28193, r28131, r28192, MPFR_RNDN); mpfr_mul(r28194, r28152, r28115, MPFR_RNDN); mpfr_mul(r28195, r28194, r28115, MPFR_RNDN); mpfr_add(r28196, r28193, r28195, MPFR_RNDN); mpfr_mul(r28197, r28191, r28133, MPFR_RNDN); mpfr_mul(r28198, r28197, r28197, MPFR_RNDN); mpfr_mul(r28199, r28148, r28167, MPFR_RNDN); mpfr_mul(r28200, r28199, r28133, MPFR_RNDN); mpfr_add(r28201, r28164, r28200, MPFR_RNDN); mpfr_mul(r28202, r28198, r28201, MPFR_RNDN); mpfr_mul(r28203, r28191, r28191, MPFR_RNDN); mpfr_mul(r28204, r28115, r28161, MPFR_RNDN); mpfr_mul(r28205, r28204, r28115, MPFR_RNDN); mpfr_add(r28206, r28205, r28158, MPFR_RNDN); mpfr_mul(r28207, r28203, r28206, MPFR_RNDN); mpfr_add(r28208, r28202, r28207, MPFR_RNDN); mpfr_add(r28209, r28196, r28208, MPFR_RNDN); mpfr_div(r28210, r28189, r28209, MPFR_RNDN); mpfr_exp(r28211, r28210, MPFR_RNDN); mpfr_log(r28212, r28211, MPFR_RNDN); mpfr_mul(r28213, r28212, r28115, MPFR_RNDN); if (mpfr_get_si(r28175, MPFR_RNDN)) { mpfr_set(r28214, r28128, MPFR_RNDN); } else { mpfr_set(r28214, r28213, MPFR_RNDN); }; if (mpfr_get_si(r28130, MPFR_RNDN)) { mpfr_set(r28215, r28173, MPFR_RNDN); } else { mpfr_set(r28215, r28214, MPFR_RNDN); }; return mpfr_get_d(r28215, MPFR_RNDN); } static mpfr_t r28216, r28217, r28218, r28219, r28220, r28221, r28222, r28223, r28224, r28225, r28226, r28227, r28228, r28229, r28230, r28231, r28232, r28233, r28234, r28235, r28236, r28237, r28238, r28239, r28240, r28241, r28242, r28243, r28244, r28245, r28246, r28247, r28248, r28249, r28250, r28251, r28252, r28253, r28254, r28255, r28256, r28257, r28258, r28259, r28260, r28261, r28262, r28263, r28264, r28265, r28266, r28267, r28268, r28269, r28270, r28271, r28272, r28273, r28274, r28275, r28276, r28277, r28278, r28279, r28280, r28281, r28282, r28283, r28284, r28285, r28286, r28287, r28288, r28289, r28290, r28291, r28292, r28293, r28294, r28295, r28296, r28297, r28298, r28299, r28300, r28301, r28302, r28303, r28304, r28305, r28306, r28307, r28308, r28309, r28310, r28311, r28312, r28313, r28314, r28315, r28316, r28317; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r28216, "0.5", 10, MPFR_RNDN); mpfr_init(r28217); mpfr_init(r28218); mpfr_init_set_str(r28219, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r28220); mpfr_init_set_str(r28221, "4", 10, MPFR_RNDN); mpfr_init(r28222); mpfr_init(r28223); mpfr_init_set_str(r28224, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r28225); mpfr_init_set_str(r28226, "6", 10, MPFR_RNDN); mpfr_init(r28227); mpfr_init(r28228); mpfr_init(r28229); mpfr_init(r28230); mpfr_init_set_str(r28231, "-2.0383110714971574e-09", 10, MPFR_RNDN); mpfr_init(r28232); mpfr_init_set_str(r28233, "1", 10, MPFR_RNDN); mpfr_init_set_str(r28234, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r28235); mpfr_init(r28236); mpfr_init(r28237); mpfr_init_set_str(r28238, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r28239); mpfr_init(r28240); mpfr_init(r28241); mpfr_init_set_str(r28242, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r28243); mpfr_init(r28244); mpfr_init(r28245); mpfr_init_set_str(r28246, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r28247); mpfr_init(r28248); mpfr_init(r28249); mpfr_init_set_str(r28250, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28251); mpfr_init(r28252); mpfr_init(r28253); mpfr_init_set_str(r28254, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28255); mpfr_init(r28256); mpfr_init_set_str(r28257, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r28258); mpfr_init(r28259); mpfr_init_set_str(r28260, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28261); mpfr_init(r28262); mpfr_init_set_str(r28263, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r28264); mpfr_init(r28265); mpfr_init_set_str(r28266, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28267); mpfr_init(r28268); mpfr_init_set_str(r28269, "2", 10, MPFR_RNDN); mpfr_init(r28270); mpfr_init(r28271); mpfr_init(r28272); mpfr_init(r28273); mpfr_init(r28274); mpfr_init(r28275); mpfr_init_set_str(r28276, "1.903280832029775e-24", 10, MPFR_RNDN); mpfr_init(r28277); mpfr_init(r28278); mpfr_init(r28279); mpfr_init(r28280); mpfr_init(r28281); mpfr_init(r28282); mpfr_init(r28283); mpfr_init_set_str(r28284, "3", 10, MPFR_RNDN); mpfr_init(r28285); mpfr_init(r28286); mpfr_init(r28287); mpfr_init(r28288); mpfr_init(r28289); mpfr_init(r28290); mpfr_init(r28291); mpfr_init(r28292); mpfr_init(r28293); mpfr_init(r28294); mpfr_init(r28295); mpfr_init(r28296); mpfr_init(r28297); mpfr_init(r28298); mpfr_init(r28299); mpfr_init(r28300); mpfr_init(r28301); mpfr_init(r28302); mpfr_init(r28303); mpfr_init(r28304); mpfr_init(r28305); mpfr_init(r28306); mpfr_init(r28307); mpfr_init(r28308); mpfr_init(r28309); mpfr_init(r28310); mpfr_init(r28311); mpfr_init(r28312); mpfr_init(r28313); mpfr_init(r28314); mpfr_init(r28315); mpfr_init(r28316); mpfr_init(r28317); } double f_dm(double x) { ; mpfr_set_d(r28217, x, MPFR_RNDN); mpfr_div(r28218, r28216, r28217, MPFR_RNDN); ; mpfr_mul(r28220, r28219, r28217, MPFR_RNDN); ; mpfr_pow(r28222, r28217, r28221, MPFR_RNDN); mpfr_div(r28223, r28220, r28222, MPFR_RNDN); ; mpfr_mul(r28225, r28224, r28217, MPFR_RNDN); ; mpfr_pow(r28227, r28217, r28226, MPFR_RNDN); mpfr_div(r28228, r28225, r28227, MPFR_RNDN); mpfr_add(r28229, r28223, r28228, MPFR_RNDN); mpfr_add(r28230, r28218, r28229, MPFR_RNDN); ; mpfr_set_si(r28232, mpfr_cmp(r28230, r28231) <= 0, MPFR_RNDN); ; ; mpfr_mul(r28235, r28217, r28217, MPFR_RNDN); mpfr_mul(r28236, r28234, r28235, MPFR_RNDN); mpfr_add(r28237, r28233, r28236, MPFR_RNDN); ; mpfr_mul(r28239, r28235, r28235, MPFR_RNDN); mpfr_mul(r28240, r28238, r28239, MPFR_RNDN); mpfr_add(r28241, r28237, r28240, MPFR_RNDN); ; mpfr_mul(r28243, r28239, r28235, MPFR_RNDN); mpfr_mul(r28244, r28242, r28243, MPFR_RNDN); mpfr_add(r28245, r28241, r28244, MPFR_RNDN); ; mpfr_mul(r28247, r28243, r28235, MPFR_RNDN); mpfr_mul(r28248, r28246, r28247, MPFR_RNDN); mpfr_add(r28249, r28245, r28248, MPFR_RNDN); ; mpfr_mul(r28251, r28247, r28235, MPFR_RNDN); mpfr_mul(r28252, r28250, r28251, MPFR_RNDN); mpfr_add(r28253, r28249, r28252, MPFR_RNDN); ; mpfr_mul(r28255, r28254, r28235, MPFR_RNDN); mpfr_add(r28256, r28233, r28255, MPFR_RNDN); ; mpfr_mul(r28258, r28257, r28239, MPFR_RNDN); mpfr_add(r28259, r28256, r28258, MPFR_RNDN); ; mpfr_mul(r28261, r28260, r28243, MPFR_RNDN); mpfr_add(r28262, r28259, r28261, MPFR_RNDN); ; mpfr_mul(r28264, r28263, r28247, MPFR_RNDN); mpfr_add(r28265, r28262, r28264, MPFR_RNDN); ; mpfr_mul(r28267, r28266, r28251, MPFR_RNDN); mpfr_add(r28268, r28265, r28267, MPFR_RNDN); ; mpfr_mul(r28270, r28269, r28250, MPFR_RNDN); mpfr_mul(r28271, r28251, r28235, MPFR_RNDN); mpfr_mul(r28272, r28270, r28271, MPFR_RNDN); mpfr_add(r28273, r28268, r28272, MPFR_RNDN); mpfr_div(r28274, r28253, r28273, MPFR_RNDN); mpfr_mul(r28275, r28274, r28217, MPFR_RNDN); ; mpfr_set_si(r28277, mpfr_cmp(r28230, r28276) <= 0, MPFR_RNDN); mpfr_mul(r28278, r28235, r28234, MPFR_RNDN); mpfr_add(r28279, r28278, r28233, MPFR_RNDN); mpfr_mul(r28280, r28217, r28242, MPFR_RNDN); mpfr_mul(r28281, r28280, r28217, MPFR_RNDN); mpfr_add(r28282, r28238, r28281, MPFR_RNDN); mpfr_mul(r28283, r28239, r28282, MPFR_RNDN); ; mpfr_add(r28285, r28284, r28233, MPFR_RNDN); mpfr_pow(r28286, r28235, r28285, MPFR_RNDN); mpfr_mul(r28287, r28235, r28250, MPFR_RNDN); mpfr_add(r28288, r28287, r28246, MPFR_RNDN); mpfr_mul(r28289, r28286, r28288, MPFR_RNDN); mpfr_add(r28290, r28283, r28289, MPFR_RNDN); mpfr_add(r28291, r28279, r28290, MPFR_RNDN); mpfr_mul(r28292, r28257, r28217, MPFR_RNDN); mpfr_pow(r28293, r28217, r28284, MPFR_RNDN); mpfr_mul(r28294, r28292, r28293, MPFR_RNDN); mpfr_add(r28295, r28233, r28294, MPFR_RNDN); mpfr_mul(r28296, r28254, r28217, MPFR_RNDN); mpfr_mul(r28297, r28296, r28217, MPFR_RNDN); mpfr_add(r28298, r28295, r28297, MPFR_RNDN); mpfr_mul(r28299, r28293, r28235, MPFR_RNDN); mpfr_mul(r28300, r28299, r28299, MPFR_RNDN); mpfr_mul(r28301, r28250, r28269, MPFR_RNDN); mpfr_mul(r28302, r28301, r28235, MPFR_RNDN); mpfr_add(r28303, r28266, r28302, MPFR_RNDN); mpfr_mul(r28304, r28300, r28303, MPFR_RNDN); mpfr_mul(r28305, r28293, r28293, MPFR_RNDN); mpfr_mul(r28306, r28217, r28263, MPFR_RNDN); mpfr_mul(r28307, r28306, r28217, MPFR_RNDN); mpfr_add(r28308, r28307, r28260, MPFR_RNDN); mpfr_mul(r28309, r28305, r28308, MPFR_RNDN); mpfr_add(r28310, r28304, r28309, MPFR_RNDN); mpfr_add(r28311, r28298, r28310, MPFR_RNDN); mpfr_div(r28312, r28291, r28311, MPFR_RNDN); mpfr_exp(r28313, r28312, MPFR_RNDN); mpfr_log(r28314, r28313, MPFR_RNDN); mpfr_mul(r28315, r28314, r28217, MPFR_RNDN); if (mpfr_get_si(r28277, MPFR_RNDN)) { mpfr_set(r28316, r28230, MPFR_RNDN); } else { mpfr_set(r28316, r28315, MPFR_RNDN); }; if (mpfr_get_si(r28232, MPFR_RNDN)) { mpfr_set(r28317, r28275, MPFR_RNDN); } else { mpfr_set(r28317, r28316, MPFR_RNDN); }; return mpfr_get_d(r28317, MPFR_RNDN); }