#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18872 = 1.0f; float r18873 = 0.10499349236488342f; float r18874 = x; float r18875 = r18874 * r18874; float r18876 = r18873 * r18875; float r18877 = r18872 + r18876; float r18878 = 0.04240605980157852f; float r18879 = r18875 * r18875; float r18880 = r18878 * r18879; float r18881 = r18877 + r18880; float r18882 = 0.007264418061822653f; float r18883 = r18879 * r18875; float r18884 = r18882 * r18883; float r18885 = r18881 + r18884; float r18886 = 0.0005064033903181553f; float r18887 = r18883 * r18875; float r18888 = r18886 * r18887; float r18889 = r18885 + r18888; float r18890 = 0.00017899709928315133f; float r18891 = r18887 * r18875; float r18892 = r18890 * r18891; float r18893 = r18889 + r18892; float r18894 = 0.7715470790863037f; float r18895 = r18894 * r18875; float r18896 = r18872 + r18895; float r18897 = 0.29097387194633484f; float r18898 = r18897 * r18879; float r18899 = r18896 + r18898; float r18900 = 0.0694555789232254f; float r18901 = r18900 * r18883; float r18902 = r18899 + r18901; float r18903 = 0.01400054432451725f; float r18904 = r18903 * r18887; float r18905 = r18902 + r18904; float r18906 = 0.0008327945251949131f; float r18907 = r18906 * r18891; float r18908 = r18905 + r18907; float r18909 = 2.0f; float r18910 = r18909 * r18890; float r18911 = r18891 * r18875; float r18912 = r18910 * r18911; float r18913 = r18908 + r18912; float r18914 = r18893 / r18913; float r18915 = r18914 * r18874; return r18915; } double f_id(double x) { double r18916 = 1.0; double r18917 = 0.10499349236488342; double r18918 = x; double r18919 = r18918 * r18918; double r18920 = r18917 * r18919; double r18921 = r18916 + r18920; double r18922 = 0.04240605980157852; double r18923 = r18919 * r18919; double r18924 = r18922 * r18923; double r18925 = r18921 + r18924; double r18926 = 0.007264418061822653; double r18927 = r18923 * r18919; double r18928 = r18926 * r18927; double r18929 = r18925 + r18928; double r18930 = 0.0005064033903181553; double r18931 = r18927 * r18919; double r18932 = r18930 * r18931; double r18933 = r18929 + r18932; double r18934 = 0.00017899709928315133; double r18935 = r18931 * r18919; double r18936 = r18934 * r18935; double r18937 = r18933 + r18936; double r18938 = 0.7715470790863037; double r18939 = r18938 * r18919; double r18940 = r18916 + r18939; double r18941 = 0.29097387194633484; double r18942 = r18941 * r18923; double r18943 = r18940 + r18942; double r18944 = 0.0694555789232254; double r18945 = r18944 * r18927; double r18946 = r18943 + r18945; double r18947 = 0.01400054432451725; double r18948 = r18947 * r18931; double r18949 = r18946 + r18948; double r18950 = 0.0008327945251949131; double r18951 = r18950 * r18935; double r18952 = r18949 + r18951; double r18953 = 2.0; double r18954 = r18953 * r18934; double r18955 = r18935 * r18919; double r18956 = r18954 * r18955; double r18957 = r18952 + r18956; double r18958 = r18937 / r18957; double r18959 = r18958 * r18918; return r18959; } double f_of(float x) { float r18960 = x; float r18961 = r18960 * r18960; float r18962 = r18961 * (r18961 * r18961); float r18963 = 0.0694555789232254f; float r18964 = 0.01400054432451725f; float r18965 = r18961 * r18964; float r18966 = r18963 + r18965; float r18967 = r18962 * r18966; float r18968 = 0.7715470790863037f; float r18969 = r18968 * r18960; float r18970 = r18969 * r18960; float r18971 = 1.0f; float r18972 = r18970 + r18971; float r18973 = 0.29097387194633484f; float r18974 = r18960 * r18973; float r18975 = r18960 * (r18960 * r18960); float r18976 = r18974 * r18975; float r18977 = r18972 + r18976; float r18978 = r18967 + r18977; float r18979 = 2.0f; float r18980 = r18960 * r18979; float r18981 = 0.00017899709928315133f; float r18982 = r18960 * r18981; float r18983 = r18980 * r18982; float r18984 = 0.0008327945251949131f; float r18985 = r18983 + r18984; float r18986 = r18961 * r18961; float r18987 = r18986 * r18962; float r18988 = r18985 * r18987; float r18989 = r18978 + r18988; float r18990 = 0.0005064033903181553f; float r18991 = r18961 * r18990; float r18992 = 0.007264418061822653f; float r18993 = r18991 + r18992; float r18994 = r18962 * r18993; float r18995 = r18986 * r18986; float r18996 = r18982 * r18960; float r18997 = r18995 * r18996; float r18998 = 0.10499349236488342f; float r18999 = r18960 * r18998; float r19000 = r18960 * r18999; float r19001 = 0.04240605980157852f; float r19002 = r18960 * r19001; float r19003 = r18975 * r19002; float r19004 = r19000 + r19003; float r19005 = r19004 + r18971; float r19006 = r18997 + r19005; float r19007 = r18994 + r19006; float r19008 = r18989 / r19007; float r19009 = r18960 / r19008; return r19009; } double f_od(double x) { double r19010 = x; double r19011 = r19010 * r19010; double r19012 = r19011 * (r19011 * r19011); double r19013 = 0.0694555789232254; double r19014 = 0.01400054432451725; double r19015 = r19011 * r19014; double r19016 = r19013 + r19015; double r19017 = r19012 * r19016; double r19018 = 0.7715470790863037; double r19019 = r19018 * r19010; double r19020 = r19019 * r19010; double r19021 = 1.0; double r19022 = r19020 + r19021; double r19023 = 0.29097387194633484; double r19024 = r19010 * r19023; double r19025 = r19010 * (r19010 * r19010); double r19026 = r19024 * r19025; double r19027 = r19022 + r19026; double r19028 = r19017 + r19027; double r19029 = 2.0; double r19030 = r19010 * r19029; double r19031 = 0.00017899709928315133; double r19032 = r19010 * r19031; double r19033 = r19030 * r19032; double r19034 = 0.0008327945251949131; double r19035 = r19033 + r19034; double r19036 = r19011 * r19011; double r19037 = r19036 * r19012; double r19038 = r19035 * r19037; double r19039 = r19028 + r19038; double r19040 = 0.0005064033903181553; double r19041 = r19011 * r19040; double r19042 = 0.007264418061822653; double r19043 = r19041 + r19042; double r19044 = r19012 * r19043; double r19045 = r19036 * r19036; double r19046 = r19032 * r19010; double r19047 = r19045 * r19046; double r19048 = 0.10499349236488342; double r19049 = r19010 * r19048; double r19050 = r19010 * r19049; double r19051 = 0.04240605980157852; double r19052 = r19010 * r19051; double r19053 = r19025 * r19052; double r19054 = r19050 + r19053; double r19055 = r19054 + r19021; double r19056 = r19047 + r19055; double r19057 = r19044 + r19056; double r19058 = r19039 / r19057; double r19059 = r19010 / r19058; return r19059; } 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 r19060, r19061, r19062, r19063, r19064, r19065, r19066, r19067, r19068, r19069, r19070, 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, r19097, r19098, r19099, r19100, r19101, r19102, r19103; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19060, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19061, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19062); mpfr_init(r19063); mpfr_init(r19064); mpfr_init(r19065); mpfr_init_set_str(r19066, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19067); mpfr_init(r19068); mpfr_init(r19069); mpfr_init_set_str(r19070, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19071); mpfr_init(r19072); mpfr_init(r19073); mpfr_init_set_str(r19074, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19075); mpfr_init(r19076); mpfr_init(r19077); mpfr_init_set_str(r19078, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19079); mpfr_init(r19080); mpfr_init(r19081); mpfr_init_set_str(r19082, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19083); mpfr_init(r19084); mpfr_init_set_str(r19085, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19086); mpfr_init(r19087); mpfr_init_set_str(r19088, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19089); mpfr_init(r19090); mpfr_init_set_str(r19091, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19092); mpfr_init(r19093); mpfr_init_set_str(r19094, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19095); mpfr_init(r19096); mpfr_init_set_str(r19097, "2", 10, MPFR_RNDN); mpfr_init(r19098); mpfr_init(r19099); mpfr_init(r19100); mpfr_init(r19101); mpfr_init(r19102); mpfr_init(r19103); } double f_im(double x) { ; ; mpfr_set_d(r19062, x, MPFR_RNDN); mpfr_mul(r19063, r19062, r19062, MPFR_RNDN); mpfr_mul(r19064, r19061, r19063, MPFR_RNDN); mpfr_add(r19065, r19060, r19064, MPFR_RNDN); ; mpfr_mul(r19067, r19063, r19063, MPFR_RNDN); mpfr_mul(r19068, r19066, r19067, MPFR_RNDN); mpfr_add(r19069, r19065, r19068, MPFR_RNDN); ; mpfr_mul(r19071, r19067, r19063, MPFR_RNDN); mpfr_mul(r19072, r19070, r19071, MPFR_RNDN); mpfr_add(r19073, r19069, r19072, MPFR_RNDN); ; mpfr_mul(r19075, r19071, r19063, MPFR_RNDN); mpfr_mul(r19076, r19074, r19075, MPFR_RNDN); mpfr_add(r19077, r19073, r19076, MPFR_RNDN); ; mpfr_mul(r19079, r19075, r19063, MPFR_RNDN); mpfr_mul(r19080, r19078, r19079, MPFR_RNDN); mpfr_add(r19081, r19077, r19080, MPFR_RNDN); ; mpfr_mul(r19083, r19082, r19063, MPFR_RNDN); mpfr_add(r19084, r19060, r19083, MPFR_RNDN); ; mpfr_mul(r19086, r19085, r19067, MPFR_RNDN); mpfr_add(r19087, r19084, r19086, MPFR_RNDN); ; mpfr_mul(r19089, r19088, r19071, MPFR_RNDN); mpfr_add(r19090, r19087, r19089, MPFR_RNDN); ; mpfr_mul(r19092, r19091, r19075, MPFR_RNDN); mpfr_add(r19093, r19090, r19092, MPFR_RNDN); ; mpfr_mul(r19095, r19094, r19079, MPFR_RNDN); mpfr_add(r19096, r19093, r19095, MPFR_RNDN); ; mpfr_mul(r19098, r19097, r19078, MPFR_RNDN); mpfr_mul(r19099, r19079, r19063, MPFR_RNDN); mpfr_mul(r19100, r19098, r19099, MPFR_RNDN); mpfr_add(r19101, r19096, r19100, MPFR_RNDN); mpfr_div(r19102, r19081, r19101, MPFR_RNDN); mpfr_mul(r19103, r19102, r19062, MPFR_RNDN); return mpfr_get_d(r19103, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19104); mpfr_init(r19105); mpfr_init(r19106); mpfr_init_set_str(r19107, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19108, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19109); mpfr_init(r19110); mpfr_init(r19111); mpfr_init_set_str(r19112, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19113); mpfr_init(r19114); mpfr_init_set_str(r19115, "1", 10, MPFR_RNDN); mpfr_init(r19116); mpfr_init_set_str(r19117, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19118); mpfr_init(r19119); mpfr_init(r19120); mpfr_init(r19121); mpfr_init(r19122); mpfr_init_set_str(r19123, "2", 10, MPFR_RNDN); mpfr_init(r19124); mpfr_init_set_str(r19125, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19126); mpfr_init(r19127); mpfr_init_set_str(r19128, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19129); mpfr_init(r19130); mpfr_init(r19131); mpfr_init(r19132); mpfr_init(r19133); mpfr_init_set_str(r19134, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19135); mpfr_init_set_str(r19136, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init_set_str(r19142, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19143); mpfr_init(r19144); mpfr_init_set_str(r19145, "0.0424060604", 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); } double f_fm(double x) { mpfr_set_d(r19104, x, MPFR_RNDN); mpfr_sqr(r19105, r19104, MPFR_RNDN); mpfr_mul(r19106, r19105, r19105, MPFR_RNDN); mpfr_mul(r19106, r19106, r19105, MPFR_RNDN); ; ; mpfr_mul(r19109, r19105, r19108, MPFR_RNDN); mpfr_add(r19110, r19107, r19109, MPFR_RNDN); mpfr_mul(r19111, r19106, r19110, MPFR_RNDN); ; mpfr_mul(r19113, r19112, r19104, MPFR_RNDN); mpfr_mul(r19114, r19113, r19104, MPFR_RNDN); ; mpfr_add(r19116, r19114, r19115, MPFR_RNDN); ; mpfr_mul(r19118, r19104, r19117, MPFR_RNDN); mpfr_mul(r19119, r19104, r19104, MPFR_RNDN); mpfr_mul(r19119, r19119, r19104, MPFR_RNDN); mpfr_mul(r19120, r19118, r19119, MPFR_RNDN); mpfr_add(r19121, r19116, r19120, MPFR_RNDN); mpfr_add(r19122, r19111, r19121, MPFR_RNDN); ; mpfr_mul(r19124, r19104, r19123, MPFR_RNDN); ; mpfr_mul(r19126, r19104, r19125, MPFR_RNDN); mpfr_mul(r19127, r19124, r19126, MPFR_RNDN); ; mpfr_add(r19129, r19127, r19128, MPFR_RNDN); mpfr_mul(r19130, r19105, r19105, MPFR_RNDN); mpfr_mul(r19131, r19130, r19106, MPFR_RNDN); mpfr_mul(r19132, r19129, r19131, MPFR_RNDN); mpfr_add(r19133, r19122, r19132, MPFR_RNDN); ; mpfr_mul(r19135, r19105, r19134, MPFR_RNDN); ; mpfr_add(r19137, r19135, r19136, MPFR_RNDN); mpfr_mul(r19138, r19106, r19137, MPFR_RNDN); mpfr_mul(r19139, r19130, r19130, MPFR_RNDN); mpfr_mul(r19140, r19126, r19104, MPFR_RNDN); mpfr_mul(r19141, r19139, r19140, MPFR_RNDN); ; mpfr_mul(r19143, r19104, r19142, MPFR_RNDN); mpfr_mul(r19144, r19104, r19143, MPFR_RNDN); ; mpfr_mul(r19146, r19104, r19145, MPFR_RNDN); mpfr_mul(r19147, r19119, r19146, MPFR_RNDN); mpfr_add(r19148, r19144, r19147, MPFR_RNDN); mpfr_add(r19149, r19148, r19115, MPFR_RNDN); mpfr_add(r19150, r19141, r19149, MPFR_RNDN); mpfr_add(r19151, r19138, r19150, MPFR_RNDN); mpfr_div(r19152, r19133, r19151, MPFR_RNDN); mpfr_div(r19153, r19104, r19152, MPFR_RNDN); return mpfr_get_d(r19153, MPFR_RNDN); } static mpfr_t 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, r19184, r19185, r19186, r19187, r19188, r19189, r19190, r19191, r19192, r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201, r19202, r19203; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19154); mpfr_init(r19155); mpfr_init(r19156); mpfr_init_set_str(r19157, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19158, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19159); mpfr_init(r19160); mpfr_init(r19161); mpfr_init_set_str(r19162, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19163); mpfr_init(r19164); mpfr_init_set_str(r19165, "1", 10, MPFR_RNDN); mpfr_init(r19166); mpfr_init_set_str(r19167, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19168); mpfr_init(r19169); mpfr_init(r19170); mpfr_init(r19171); mpfr_init(r19172); mpfr_init_set_str(r19173, "2", 10, MPFR_RNDN); mpfr_init(r19174); mpfr_init_set_str(r19175, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19176); mpfr_init(r19177); mpfr_init_set_str(r19178, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19179); mpfr_init(r19180); mpfr_init(r19181); mpfr_init(r19182); mpfr_init(r19183); mpfr_init_set_str(r19184, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19185); mpfr_init_set_str(r19186, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19187); mpfr_init(r19188); mpfr_init(r19189); mpfr_init(r19190); mpfr_init(r19191); mpfr_init_set_str(r19192, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19193); mpfr_init(r19194); mpfr_init_set_str(r19195, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19196); mpfr_init(r19197); mpfr_init(r19198); mpfr_init(r19199); mpfr_init(r19200); mpfr_init(r19201); mpfr_init(r19202); mpfr_init(r19203); } double f_dm(double x) { mpfr_set_d(r19154, x, MPFR_RNDN); mpfr_sqr(r19155, r19154, MPFR_RNDN); mpfr_mul(r19156, r19155, r19155, MPFR_RNDN); mpfr_mul(r19156, r19156, r19155, MPFR_RNDN); ; ; mpfr_mul(r19159, r19155, r19158, MPFR_RNDN); mpfr_add(r19160, r19157, r19159, MPFR_RNDN); mpfr_mul(r19161, r19156, r19160, MPFR_RNDN); ; mpfr_mul(r19163, r19162, r19154, MPFR_RNDN); mpfr_mul(r19164, r19163, r19154, MPFR_RNDN); ; mpfr_add(r19166, r19164, r19165, MPFR_RNDN); ; mpfr_mul(r19168, r19154, r19167, MPFR_RNDN); mpfr_mul(r19169, r19154, r19154, MPFR_RNDN); mpfr_mul(r19169, r19169, r19154, MPFR_RNDN); mpfr_mul(r19170, r19168, r19169, MPFR_RNDN); mpfr_add(r19171, r19166, r19170, MPFR_RNDN); mpfr_add(r19172, r19161, r19171, MPFR_RNDN); ; mpfr_mul(r19174, r19154, r19173, MPFR_RNDN); ; mpfr_mul(r19176, r19154, r19175, MPFR_RNDN); mpfr_mul(r19177, r19174, r19176, MPFR_RNDN); ; mpfr_add(r19179, r19177, r19178, MPFR_RNDN); mpfr_mul(r19180, r19155, r19155, MPFR_RNDN); mpfr_mul(r19181, r19180, r19156, MPFR_RNDN); mpfr_mul(r19182, r19179, r19181, MPFR_RNDN); mpfr_add(r19183, r19172, r19182, MPFR_RNDN); ; mpfr_mul(r19185, r19155, r19184, MPFR_RNDN); ; mpfr_add(r19187, r19185, r19186, MPFR_RNDN); mpfr_mul(r19188, r19156, r19187, MPFR_RNDN); mpfr_mul(r19189, r19180, r19180, MPFR_RNDN); mpfr_mul(r19190, r19176, r19154, MPFR_RNDN); mpfr_mul(r19191, r19189, r19190, MPFR_RNDN); ; mpfr_mul(r19193, r19154, r19192, MPFR_RNDN); mpfr_mul(r19194, r19154, r19193, MPFR_RNDN); ; mpfr_mul(r19196, r19154, r19195, MPFR_RNDN); mpfr_mul(r19197, r19169, r19196, MPFR_RNDN); mpfr_add(r19198, r19194, r19197, MPFR_RNDN); mpfr_add(r19199, r19198, r19165, MPFR_RNDN); mpfr_add(r19200, r19191, r19199, MPFR_RNDN); mpfr_add(r19201, r19188, r19200, MPFR_RNDN); mpfr_div(r19202, r19183, r19201, MPFR_RNDN); mpfr_div(r19203, r19154, r19202, MPFR_RNDN); return mpfr_get_d(r19203, MPFR_RNDN); }