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