#include #include #include #include #include char *name = "Jmat.Real.erf"; double f_if(float x) { float r18845 = 1.0f; float r18846 = 0.3275911f; float r18847 = x; float r18848 = fabs(r18847); float r18849 = r18846 * r18848; float r18850 = r18845 + r18849; float r18851 = r18845 / r18850; float r18852 = 0.254829592f; float r18853 = -0.284496736f; float r18854 = 1.421413741f; float r18855 = -1.453152027f; float r18856 = 1.061405429f; float r18857 = r18851 * r18856; float r18858 = r18855 + r18857; float r18859 = r18851 * r18858; float r18860 = r18854 + r18859; float r18861 = r18851 * r18860; float r18862 = r18853 + r18861; float r18863 = r18851 * r18862; float r18864 = r18852 + r18863; float r18865 = r18851 * r18864; float r18866 = r18848 * r18848; float r18867 = -r18866; float r18868 = exp(r18867); float r18869 = r18865 * r18868; float r18870 = r18845 - r18869; return r18870; } double f_id(double x) { double r18871 = 1.0; double r18872 = 0.3275911; double r18873 = x; double r18874 = fabs(r18873); double r18875 = r18872 * r18874; double r18876 = r18871 + r18875; double r18877 = r18871 / r18876; double r18878 = 0.254829592; double r18879 = -0.284496736; double r18880 = 1.421413741; double r18881 = -1.453152027; double r18882 = 1.061405429; double r18883 = r18877 * r18882; double r18884 = r18881 + r18883; double r18885 = r18877 * r18884; double r18886 = r18880 + r18885; double r18887 = r18877 * r18886; double r18888 = r18879 + r18887; double r18889 = r18877 * r18888; double r18890 = r18878 + r18889; double r18891 = r18877 * r18890; double r18892 = r18874 * r18874; double r18893 = -r18892; double r18894 = exp(r18893); double r18895 = r18891 * r18894; double r18896 = r18871 - r18895; return r18896; } double f_of(float x) { float r18897 = x; float r18898 = -3.878602931916025e-24f; bool r18899 = r18897 <= r18898; float r18900 = 1.0f; float r18901 = 1.453152027f; float r18902 = fabs(r18897); float r18903 = r18902 * r18902; float r18904 = -r18903; float r18905 = exp(r18904); float r18906 = 0.3275911f; float r18907 = r18906 * r18902; float r18908 = r18900 + r18907; float r18909 = 4.0f; float r18910 = pow(r18908, r18909); float r18911 = r18905 / r18910; float r18912 = r18901 * r18911; float r18913 = 0.284496736f; float r18914 = r18908 * r18908; float r18915 = r18905 / r18914; float r18916 = r18913 * r18915; float r18917 = r18912 + r18916; float r18918 = r18900 + r18917; float r18919 = 3.0f; float r18920 = pow(r18918, r18919); float r18921 = 0.254829592f; float r18922 = r18905 / r18908; float r18923 = r18921 * r18922; float r18924 = 1.061405429f; float r18925 = 5.0f; float r18926 = pow(r18908, r18925); float r18927 = cbrt(r18926); float r18928 = r18927 * (r18927 * r18927); float r18929 = r18905 / r18928; float r18930 = r18924 * r18929; float r18931 = 1.421413741f; float r18932 = pow(r18908, r18919); float r18933 = r18905 / r18932; float r18934 = r18931 * r18933; float r18935 = r18930 + r18934; float r18936 = r18923 + r18935; float r18937 = pow(r18936, r18919); float r18938 = r18920 - r18937; float r18939 = r18918 * r18918; float r18940 = r18936 * r18936; float r18941 = r18918 * r18936; float r18942 = r18940 + r18941; float r18943 = r18939 + r18942; float r18944 = r18938 / r18943; float r18945 = 1.1454841341683684e-10f; bool r18946 = r18897 <= r18945; float r18947 = r18900 * r18900; float r18948 = r18917 * r18917; float r18949 = r18947 - r18948; float r18950 = r18905 / r18926; float r18951 = r18924 * r18950; float r18952 = r18951 + r18934; float r18953 = r18923 - r18952; float r18954 = r18949 * r18953; float r18955 = r18900 - r18917; float r18956 = r18923 * r18923; float r18957 = r18952 * r18952; float r18958 = r18956 - r18957; float r18959 = r18955 * r18958; float r18960 = r18954 - r18959; float r18961 = r18955 * r18953; float r18962 = r18960 / r18961; float r18963 = r18900 / r18908; float r18964 = -0.284496736f; float r18965 = -1.453152027f; float r18966 = sqrt(r18963); float r18967 = r18966 * r18966; float r18968 = r18967 * r18924; float r18969 = r18965 + r18968; float r18970 = r18963 * r18969; float r18971 = r18931 + r18970; float r18972 = r18963 * r18971; float r18973 = r18964 + r18972; float r18974 = r18963 * r18973; float r18975 = r18921 + r18974; float r18976 = r18963 * r18975; float r18977 = r18902 * r18902; float r18978 = -r18977; float r18979 = exp(r18978); float r18980 = r18976 * r18979; float r18981 = r18900 - r18980; float r18982 = r18946 ? r18962 : r18981; float r18983 = r18899 ? r18944 : r18982; return r18983; } double f_od(double x) { double r18984 = x; double r18985 = -3.878602931916025e-24; bool r18986 = r18984 <= r18985; double r18987 = 1.0; double r18988 = 1.453152027; double r18989 = fabs(r18984); double r18990 = r18989 * r18989; double r18991 = -r18990; double r18992 = exp(r18991); double r18993 = 0.3275911; double r18994 = r18993 * r18989; double r18995 = r18987 + r18994; double r18996 = 4.0; double r18997 = pow(r18995, r18996); double r18998 = r18992 / r18997; double r18999 = r18988 * r18998; double r19000 = 0.284496736; double r19001 = r18995 * r18995; double r19002 = r18992 / r19001; double r19003 = r19000 * r19002; double r19004 = r18999 + r19003; double r19005 = r18987 + r19004; double r19006 = 3.0; double r19007 = pow(r19005, r19006); double r19008 = 0.254829592; double r19009 = r18992 / r18995; double r19010 = r19008 * r19009; double r19011 = 1.061405429; double r19012 = 5.0; double r19013 = pow(r18995, r19012); double r19014 = cbrt(r19013); double r19015 = r19014 * (r19014 * r19014); double r19016 = r18992 / r19015; double r19017 = r19011 * r19016; double r19018 = 1.421413741; double r19019 = pow(r18995, r19006); double r19020 = r18992 / r19019; double r19021 = r19018 * r19020; double r19022 = r19017 + r19021; double r19023 = r19010 + r19022; double r19024 = pow(r19023, r19006); double r19025 = r19007 - r19024; double r19026 = r19005 * r19005; double r19027 = r19023 * r19023; double r19028 = r19005 * r19023; double r19029 = r19027 + r19028; double r19030 = r19026 + r19029; double r19031 = r19025 / r19030; double r19032 = 1.1454841341683684e-10; bool r19033 = r18984 <= r19032; double r19034 = r18987 * r18987; double r19035 = r19004 * r19004; double r19036 = r19034 - r19035; double r19037 = r18992 / r19013; double r19038 = r19011 * r19037; double r19039 = r19038 + r19021; double r19040 = r19010 - r19039; double r19041 = r19036 * r19040; double r19042 = r18987 - r19004; double r19043 = r19010 * r19010; double r19044 = r19039 * r19039; double r19045 = r19043 - r19044; double r19046 = r19042 * r19045; double r19047 = r19041 - r19046; double r19048 = r19042 * r19040; double r19049 = r19047 / r19048; double r19050 = r18987 / r18995; double r19051 = -0.284496736; double r19052 = -1.453152027; double r19053 = sqrt(r19050); double r19054 = r19053 * r19053; double r19055 = r19054 * r19011; double r19056 = r19052 + r19055; double r19057 = r19050 * r19056; double r19058 = r19018 + r19057; double r19059 = r19050 * r19058; double r19060 = r19051 + r19059; double r19061 = r19050 * r19060; double r19062 = r19008 + r19061; double r19063 = r19050 * r19062; double r19064 = r18989 * r18989; double r19065 = -r19064; double r19066 = exp(r19065); double r19067 = r19063 * r19066; double r19068 = r18987 - r19067; double r19069 = r19033 ? r19049 : r19068; double r19070 = r18986 ? r19031 : r19069; return r19070; } 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 r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086, r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19071, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19072, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19073); mpfr_init(r19074); mpfr_init(r19075); mpfr_init(r19076); mpfr_init(r19077); mpfr_init_set_str(r19078, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r19079, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19080, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r19081, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r19082, "1.061405429", 10, MPFR_RNDN); mpfr_init(r19083); mpfr_init(r19084); mpfr_init(r19085); mpfr_init(r19086); mpfr_init(r19087); mpfr_init(r19088); mpfr_init(r19089); mpfr_init(r19090); mpfr_init(r19091); mpfr_init(r19092); mpfr_init(r19093); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); } double f_im(double x) { ; ; mpfr_set_d(r19073, x, MPFR_RNDN); mpfr_abs(r19074, r19073, MPFR_RNDN); mpfr_mul(r19075, r19072, r19074, MPFR_RNDN); mpfr_add(r19076, r19071, r19075, MPFR_RNDN); mpfr_div(r19077, r19071, r19076, MPFR_RNDN); ; ; ; ; ; mpfr_mul(r19083, r19077, r19082, MPFR_RNDN); mpfr_add(r19084, r19081, r19083, MPFR_RNDN); mpfr_mul(r19085, r19077, r19084, MPFR_RNDN); mpfr_add(r19086, r19080, r19085, MPFR_RNDN); mpfr_mul(r19087, r19077, r19086, MPFR_RNDN); mpfr_add(r19088, r19079, r19087, MPFR_RNDN); mpfr_mul(r19089, r19077, r19088, MPFR_RNDN); mpfr_add(r19090, r19078, r19089, MPFR_RNDN); mpfr_mul(r19091, r19077, r19090, MPFR_RNDN); mpfr_mul(r19092, r19074, r19074, MPFR_RNDN); mpfr_neg(r19093, r19092, MPFR_RNDN); mpfr_exp(r19094, r19093, MPFR_RNDN); mpfr_mul(r19095, r19091, r19094, MPFR_RNDN); mpfr_sub(r19096, r19071, r19095, MPFR_RNDN); return mpfr_get_d(r19096, MPFR_RNDN); } static mpfr_t r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117, r19118, r19119, r19120, r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138, r19139, r19140, r19141, r19142, r19143, r19144, r19145, r19146, r19147, r19148, r19149, r19150, r19151, r19152, r19153, r19154, r19155, r19156, r19157, r19158, r19159, r19160, r19161, r19162, r19163, r19164, r19165, r19166, r19167, r19168, r19169, r19170, r19171, r19172, r19173, r19174, r19175, r19176, r19177, r19178, r19179, r19180, r19181, r19182, r19183; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19097); mpfr_init_set_str(r19098, "-3.878602931916025e-24", 10, MPFR_RNDN); mpfr_init(r19099); mpfr_init_set_str(r19100, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19101, "1.453152027", 10, MPFR_RNDN); mpfr_init(r19102); mpfr_init(r19103); mpfr_init(r19104); mpfr_init(r19105); mpfr_init_set_str(r19106, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19107); mpfr_init(r19108); mpfr_init_set_str(r19109, "4", 10, MPFR_RNDN); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init_set_str(r19113, "0.284496736", 10, MPFR_RNDN); mpfr_init(r19114); mpfr_init(r19115); mpfr_init(r19116); mpfr_init(r19117); mpfr_init(r19118); mpfr_init_set_str(r19119, "3", 10, MPFR_RNDN); mpfr_init(r19120); mpfr_init_set_str(r19121, "0.254829592", 10, MPFR_RNDN); mpfr_init(r19122); mpfr_init(r19123); mpfr_init_set_str(r19124, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r19125, "5", 10, MPFR_RNDN); mpfr_init(r19126); mpfr_init(r19127); mpfr_init(r19128); mpfr_init(r19129); mpfr_init(r19130); mpfr_init_set_str(r19131, "1.421413741", 10, MPFR_RNDN); mpfr_init(r19132); mpfr_init(r19133); mpfr_init(r19134); mpfr_init(r19135); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init(r19142); mpfr_init(r19143); mpfr_init(r19144); mpfr_init_set_str(r19145, "1.1454841341683684e-10", 10, MPFR_RNDN); mpfr_init(r19146); mpfr_init(r19147); mpfr_init(r19148); mpfr_init(r19149); mpfr_init(r19150); mpfr_init(r19151); mpfr_init(r19152); mpfr_init(r19153); mpfr_init(r19154); mpfr_init(r19155); mpfr_init(r19156); mpfr_init(r19157); mpfr_init(r19158); mpfr_init(r19159); mpfr_init(r19160); mpfr_init(r19161); mpfr_init(r19162); mpfr_init(r19163); mpfr_init_set_str(r19164, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19165, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19166); mpfr_init(r19167); mpfr_init(r19168); mpfr_init(r19169); mpfr_init(r19170); mpfr_init(r19171); mpfr_init(r19172); mpfr_init(r19173); mpfr_init(r19174); mpfr_init(r19175); mpfr_init(r19176); mpfr_init(r19177); mpfr_init(r19178); mpfr_init(r19179); mpfr_init(r19180); mpfr_init(r19181); mpfr_init(r19182); mpfr_init(r19183); } double f_fm(double x) { mpfr_set_d(r19097, x, MPFR_RNDN); ; mpfr_set_si(r19099, mpfr_cmp(r19097, r19098) <= 0, MPFR_RNDN); ; ; mpfr_abs(r19102, r19097, MPFR_RNDN); mpfr_sqr(r19103, r19102, MPFR_RNDN); mpfr_neg(r19104, r19103, MPFR_RNDN); mpfr_exp(r19105, r19104, MPFR_RNDN); ; mpfr_mul(r19107, r19106, r19102, MPFR_RNDN); mpfr_add(r19108, r19100, r19107, MPFR_RNDN); ; mpfr_pow(r19110, r19108, r19109, MPFR_RNDN); mpfr_div(r19111, r19105, r19110, MPFR_RNDN); mpfr_mul(r19112, r19101, r19111, MPFR_RNDN); ; mpfr_sqr(r19114, r19108, MPFR_RNDN); mpfr_div(r19115, r19105, r19114, MPFR_RNDN); mpfr_mul(r19116, r19113, r19115, MPFR_RNDN); mpfr_add(r19117, r19112, r19116, MPFR_RNDN); mpfr_add(r19118, r19100, r19117, MPFR_RNDN); ; mpfr_pow(r19120, r19118, r19119, MPFR_RNDN); ; mpfr_div(r19122, r19105, r19108, MPFR_RNDN); mpfr_mul(r19123, r19121, r19122, MPFR_RNDN); ; ; mpfr_pow(r19126, r19108, r19125, MPFR_RNDN); mpfr_cbrt(r19127, r19126, MPFR_RNDN); mpfr_mul(r19128, r19127, r19127, MPFR_RNDN); mpfr_mul(r19128, r19128, r19127, MPFR_RNDN); mpfr_div(r19129, r19105, r19128, MPFR_RNDN); mpfr_mul(r19130, r19124, r19129, MPFR_RNDN); ; mpfr_pow(r19132, r19108, r19119, MPFR_RNDN); mpfr_div(r19133, r19105, r19132, MPFR_RNDN); mpfr_mul(r19134, r19131, r19133, MPFR_RNDN); mpfr_add(r19135, r19130, r19134, MPFR_RNDN); mpfr_add(r19136, r19123, r19135, MPFR_RNDN); mpfr_pow(r19137, r19136, r19119, MPFR_RNDN); mpfr_sub(r19138, r19120, r19137, MPFR_RNDN); mpfr_sqr(r19139, r19118, MPFR_RNDN); mpfr_sqr(r19140, r19136, MPFR_RNDN); mpfr_mul(r19141, r19118, r19136, MPFR_RNDN); mpfr_add(r19142, r19140, r19141, MPFR_RNDN); mpfr_add(r19143, r19139, r19142, MPFR_RNDN); mpfr_div(r19144, r19138, r19143, MPFR_RNDN); ; mpfr_set_si(r19146, mpfr_cmp(r19097, r19145) <= 0, MPFR_RNDN); mpfr_sqr(r19147, r19100, MPFR_RNDN); mpfr_sqr(r19148, r19117, MPFR_RNDN); mpfr_sub(r19149, r19147, r19148, MPFR_RNDN); mpfr_div(r19150, r19105, r19126, MPFR_RNDN); mpfr_mul(r19151, r19124, r19150, MPFR_RNDN); mpfr_add(r19152, r19151, r19134, MPFR_RNDN); mpfr_sub(r19153, r19123, r19152, MPFR_RNDN); mpfr_mul(r19154, r19149, r19153, MPFR_RNDN); mpfr_sub(r19155, r19100, r19117, MPFR_RNDN); mpfr_sqr(r19156, r19123, MPFR_RNDN); mpfr_sqr(r19157, r19152, MPFR_RNDN); mpfr_sub(r19158, r19156, r19157, MPFR_RNDN); mpfr_mul(r19159, r19155, r19158, MPFR_RNDN); mpfr_sub(r19160, r19154, r19159, MPFR_RNDN); mpfr_mul(r19161, r19155, r19153, MPFR_RNDN); mpfr_div(r19162, r19160, r19161, MPFR_RNDN); mpfr_div(r19163, r19100, r19108, MPFR_RNDN); ; ; mpfr_sqrt(r19166, r19163, MPFR_RNDN); mpfr_sqr(r19167, r19166, MPFR_RNDN); mpfr_mul(r19168, r19167, r19124, MPFR_RNDN); mpfr_add(r19169, r19165, r19168, MPFR_RNDN); mpfr_mul(r19170, r19163, r19169, MPFR_RNDN); mpfr_add(r19171, r19131, r19170, MPFR_RNDN); mpfr_mul(r19172, r19163, r19171, MPFR_RNDN); mpfr_add(r19173, r19164, r19172, MPFR_RNDN); mpfr_mul(r19174, r19163, r19173, MPFR_RNDN); mpfr_add(r19175, r19121, r19174, MPFR_RNDN); mpfr_mul(r19176, r19163, r19175, MPFR_RNDN); mpfr_mul(r19177, r19102, r19102, MPFR_RNDN); mpfr_neg(r19178, r19177, MPFR_RNDN); mpfr_exp(r19179, r19178, MPFR_RNDN); mpfr_mul(r19180, r19176, r19179, MPFR_RNDN); mpfr_sub(r19181, r19100, r19180, MPFR_RNDN); if (mpfr_get_si(r19146, MPFR_RNDN)) { mpfr_set(r19182, r19162, MPFR_RNDN); } else { mpfr_set(r19182, r19181, MPFR_RNDN); }; if (mpfr_get_si(r19099, MPFR_RNDN)) { mpfr_set(r19183, r19144, MPFR_RNDN); } else { mpfr_set(r19183, r19182, MPFR_RNDN); }; return mpfr_get_d(r19183, MPFR_RNDN); } static mpfr_t r19184, r19185, r19186, r19187, r19188, r19189, r19190, r19191, r19192, r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201, r19202, r19203, r19204, r19205, r19206, r19207, r19208, r19209, r19210, r19211, r19212, r19213, r19214, r19215, r19216, r19217, r19218, r19219, r19220, r19221, r19222, r19223, r19224, r19225, r19226, r19227, r19228, r19229, r19230, r19231, r19232, r19233, r19234, r19235, r19236, r19237, r19238, r19239, r19240, r19241, r19242, r19243, r19244, r19245, r19246, r19247, r19248, r19249, r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265, r19266, r19267, r19268, r19269, r19270; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19184); mpfr_init_set_str(r19185, "-3.878602931916025e-24", 10, MPFR_RNDN); mpfr_init(r19186); mpfr_init_set_str(r19187, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19188, "1.453152027", 10, MPFR_RNDN); mpfr_init(r19189); mpfr_init(r19190); mpfr_init(r19191); mpfr_init(r19192); mpfr_init_set_str(r19193, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19194); mpfr_init(r19195); mpfr_init_set_str(r19196, "4", 10, MPFR_RNDN); mpfr_init(r19197); mpfr_init(r19198); mpfr_init(r19199); mpfr_init_set_str(r19200, "0.284496736", 10, MPFR_RNDN); mpfr_init(r19201); mpfr_init(r19202); mpfr_init(r19203); mpfr_init(r19204); mpfr_init(r19205); mpfr_init_set_str(r19206, "3", 10, MPFR_RNDN); mpfr_init(r19207); mpfr_init_set_str(r19208, "0.254829592", 10, MPFR_RNDN); mpfr_init(r19209); mpfr_init(r19210); mpfr_init_set_str(r19211, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r19212, "5", 10, MPFR_RNDN); mpfr_init(r19213); mpfr_init(r19214); mpfr_init(r19215); mpfr_init(r19216); mpfr_init(r19217); mpfr_init_set_str(r19218, "1.421413741", 10, MPFR_RNDN); mpfr_init(r19219); mpfr_init(r19220); mpfr_init(r19221); mpfr_init(r19222); mpfr_init(r19223); mpfr_init(r19224); mpfr_init(r19225); mpfr_init(r19226); mpfr_init(r19227); mpfr_init(r19228); mpfr_init(r19229); mpfr_init(r19230); mpfr_init(r19231); mpfr_init_set_str(r19232, "1.1454841341683684e-10", 10, MPFR_RNDN); mpfr_init(r19233); mpfr_init(r19234); mpfr_init(r19235); mpfr_init(r19236); mpfr_init(r19237); mpfr_init(r19238); mpfr_init(r19239); mpfr_init(r19240); mpfr_init(r19241); mpfr_init(r19242); mpfr_init(r19243); mpfr_init(r19244); mpfr_init(r19245); mpfr_init(r19246); mpfr_init(r19247); mpfr_init(r19248); mpfr_init(r19249); mpfr_init(r19250); mpfr_init_set_str(r19251, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19252, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19253); mpfr_init(r19254); mpfr_init(r19255); mpfr_init(r19256); mpfr_init(r19257); mpfr_init(r19258); mpfr_init(r19259); mpfr_init(r19260); mpfr_init(r19261); mpfr_init(r19262); mpfr_init(r19263); mpfr_init(r19264); mpfr_init(r19265); mpfr_init(r19266); mpfr_init(r19267); mpfr_init(r19268); mpfr_init(r19269); mpfr_init(r19270); } double f_dm(double x) { mpfr_set_d(r19184, x, MPFR_RNDN); ; mpfr_set_si(r19186, mpfr_cmp(r19184, r19185) <= 0, MPFR_RNDN); ; ; mpfr_abs(r19189, r19184, MPFR_RNDN); mpfr_sqr(r19190, r19189, MPFR_RNDN); mpfr_neg(r19191, r19190, MPFR_RNDN); mpfr_exp(r19192, r19191, MPFR_RNDN); ; mpfr_mul(r19194, r19193, r19189, MPFR_RNDN); mpfr_add(r19195, r19187, r19194, MPFR_RNDN); ; mpfr_pow(r19197, r19195, r19196, MPFR_RNDN); mpfr_div(r19198, r19192, r19197, MPFR_RNDN); mpfr_mul(r19199, r19188, r19198, MPFR_RNDN); ; mpfr_sqr(r19201, r19195, MPFR_RNDN); mpfr_div(r19202, r19192, r19201, MPFR_RNDN); mpfr_mul(r19203, r19200, r19202, MPFR_RNDN); mpfr_add(r19204, r19199, r19203, MPFR_RNDN); mpfr_add(r19205, r19187, r19204, MPFR_RNDN); ; mpfr_pow(r19207, r19205, r19206, MPFR_RNDN); ; mpfr_div(r19209, r19192, r19195, MPFR_RNDN); mpfr_mul(r19210, r19208, r19209, MPFR_RNDN); ; ; mpfr_pow(r19213, r19195, r19212, MPFR_RNDN); mpfr_cbrt(r19214, r19213, MPFR_RNDN); mpfr_mul(r19215, r19214, r19214, MPFR_RNDN); mpfr_mul(r19215, r19215, r19214, MPFR_RNDN); mpfr_div(r19216, r19192, r19215, MPFR_RNDN); mpfr_mul(r19217, r19211, r19216, MPFR_RNDN); ; mpfr_pow(r19219, r19195, r19206, MPFR_RNDN); mpfr_div(r19220, r19192, r19219, MPFR_RNDN); mpfr_mul(r19221, r19218, r19220, MPFR_RNDN); mpfr_add(r19222, r19217, r19221, MPFR_RNDN); mpfr_add(r19223, r19210, r19222, MPFR_RNDN); mpfr_pow(r19224, r19223, r19206, MPFR_RNDN); mpfr_sub(r19225, r19207, r19224, MPFR_RNDN); mpfr_sqr(r19226, r19205, MPFR_RNDN); mpfr_sqr(r19227, r19223, MPFR_RNDN); mpfr_mul(r19228, r19205, r19223, MPFR_RNDN); mpfr_add(r19229, r19227, r19228, MPFR_RNDN); mpfr_add(r19230, r19226, r19229, MPFR_RNDN); mpfr_div(r19231, r19225, r19230, MPFR_RNDN); ; mpfr_set_si(r19233, mpfr_cmp(r19184, r19232) <= 0, MPFR_RNDN); mpfr_sqr(r19234, r19187, MPFR_RNDN); mpfr_sqr(r19235, r19204, MPFR_RNDN); mpfr_sub(r19236, r19234, r19235, MPFR_RNDN); mpfr_div(r19237, r19192, r19213, MPFR_RNDN); mpfr_mul(r19238, r19211, r19237, MPFR_RNDN); mpfr_add(r19239, r19238, r19221, MPFR_RNDN); mpfr_sub(r19240, r19210, r19239, MPFR_RNDN); mpfr_mul(r19241, r19236, r19240, MPFR_RNDN); mpfr_sub(r19242, r19187, r19204, MPFR_RNDN); mpfr_sqr(r19243, r19210, MPFR_RNDN); mpfr_sqr(r19244, r19239, MPFR_RNDN); mpfr_sub(r19245, r19243, r19244, MPFR_RNDN); mpfr_mul(r19246, r19242, r19245, MPFR_RNDN); mpfr_sub(r19247, r19241, r19246, MPFR_RNDN); mpfr_mul(r19248, r19242, r19240, MPFR_RNDN); mpfr_div(r19249, r19247, r19248, MPFR_RNDN); mpfr_div(r19250, r19187, r19195, MPFR_RNDN); ; ; mpfr_sqrt(r19253, r19250, MPFR_RNDN); mpfr_sqr(r19254, r19253, MPFR_RNDN); mpfr_mul(r19255, r19254, r19211, MPFR_RNDN); mpfr_add(r19256, r19252, r19255, MPFR_RNDN); mpfr_mul(r19257, r19250, r19256, MPFR_RNDN); mpfr_add(r19258, r19218, r19257, MPFR_RNDN); mpfr_mul(r19259, r19250, r19258, MPFR_RNDN); mpfr_add(r19260, r19251, r19259, MPFR_RNDN); mpfr_mul(r19261, r19250, r19260, MPFR_RNDN); mpfr_add(r19262, r19208, r19261, MPFR_RNDN); mpfr_mul(r19263, r19250, r19262, MPFR_RNDN); mpfr_mul(r19264, r19189, r19189, MPFR_RNDN); mpfr_neg(r19265, r19264, MPFR_RNDN); mpfr_exp(r19266, r19265, MPFR_RNDN); mpfr_mul(r19267, r19263, r19266, MPFR_RNDN); mpfr_sub(r19268, r19187, r19267, MPFR_RNDN); if (mpfr_get_si(r19233, MPFR_RNDN)) { mpfr_set(r19269, r19249, MPFR_RNDN); } else { mpfr_set(r19269, r19268, MPFR_RNDN); }; if (mpfr_get_si(r19186, MPFR_RNDN)) { mpfr_set(r19270, r19231, MPFR_RNDN); } else { mpfr_set(r19270, r19269, MPFR_RNDN); }; return mpfr_get_d(r19270, MPFR_RNDN); }