#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18955 = 1.0f; float r18956 = 0.10499349236488342f; float r18957 = x; float r18958 = r18957 * r18957; float r18959 = r18956 * r18958; float r18960 = r18955 + r18959; float r18961 = 0.04240605980157852f; float r18962 = r18958 * r18958; float r18963 = r18961 * r18962; float r18964 = r18960 + r18963; float r18965 = 0.007264418061822653f; float r18966 = r18962 * r18958; float r18967 = r18965 * r18966; float r18968 = r18964 + r18967; float r18969 = 0.0005064033903181553f; float r18970 = r18966 * r18958; float r18971 = r18969 * r18970; float r18972 = r18968 + r18971; float r18973 = 0.00017899709928315133f; float r18974 = r18970 * r18958; float r18975 = r18973 * r18974; float r18976 = r18972 + r18975; float r18977 = 0.7715470790863037f; float r18978 = r18977 * r18958; float r18979 = r18955 + r18978; float r18980 = 0.29097387194633484f; float r18981 = r18980 * r18962; float r18982 = r18979 + r18981; float r18983 = 0.0694555789232254f; float r18984 = r18983 * r18966; float r18985 = r18982 + r18984; float r18986 = 0.01400054432451725f; float r18987 = r18986 * r18970; float r18988 = r18985 + r18987; float r18989 = 0.0008327945251949131f; float r18990 = r18989 * r18974; float r18991 = r18988 + r18990; float r18992 = 2.0f; float r18993 = r18992 * r18973; float r18994 = r18974 * r18958; float r18995 = r18993 * r18994; float r18996 = r18991 + r18995; float r18997 = r18976 / r18996; float r18998 = r18997 * r18957; return r18998; } double f_id(double x) { double r18999 = 1.0; double r19000 = 0.10499349236488342; double r19001 = x; double r19002 = r19001 * r19001; double r19003 = r19000 * r19002; double r19004 = r18999 + r19003; double r19005 = 0.04240605980157852; double r19006 = r19002 * r19002; double r19007 = r19005 * r19006; double r19008 = r19004 + r19007; double r19009 = 0.007264418061822653; double r19010 = r19006 * r19002; double r19011 = r19009 * r19010; double r19012 = r19008 + r19011; double r19013 = 0.0005064033903181553; double r19014 = r19010 * r19002; double r19015 = r19013 * r19014; double r19016 = r19012 + r19015; double r19017 = 0.00017899709928315133; double r19018 = r19014 * r19002; double r19019 = r19017 * r19018; double r19020 = r19016 + r19019; double r19021 = 0.7715470790863037; double r19022 = r19021 * r19002; double r19023 = r18999 + r19022; double r19024 = 0.29097387194633484; double r19025 = r19024 * r19006; double r19026 = r19023 + r19025; double r19027 = 0.0694555789232254; double r19028 = r19027 * r19010; double r19029 = r19026 + r19028; double r19030 = 0.01400054432451725; double r19031 = r19030 * r19014; double r19032 = r19029 + r19031; double r19033 = 0.0008327945251949131; double r19034 = r19033 * r19018; double r19035 = r19032 + r19034; double r19036 = 2.0; double r19037 = r19036 * r19017; double r19038 = r19018 * r19002; double r19039 = r19037 * r19038; double r19040 = r19035 + r19039; double r19041 = r19020 / r19040; double r19042 = r19041 * r19001; return r19042; } double f_of(float x) { float r19043 = x; float r19044 = r19043 * r19043; float r19045 = r19044 * (r19044 * r19044); float r19046 = r19043 * (r19043 * r19043); float r19047 = 0.00017899709928315133f; float r19048 = r19043 * r19047; float r19049 = r19046 * r19048; float r19050 = 0.007264418061822653f; float r19051 = r19049 + r19050; float r19052 = 0.0005064033903181553f; float r19053 = r19043 * r19052; float r19054 = r19043 * r19053; float r19055 = r19051 + r19054; float r19056 = r19045 * r19055; float r19057 = 0.10499349236488342f; float r19058 = r19057 * r19043; float r19059 = r19043 * r19058; float r19060 = 1.0f; float r19061 = 0.04240605980157852f; float r19062 = r19061 * r19043; float r19063 = r19062 * r19046; float r19064 = r19060 + r19063; float r19065 = r19059 + r19064; float r19066 = r19056 + r19065; float r19067 = 0.01400054432451725f; float r19068 = r19067 * r19043; float r19069 = r19043 * r19068; float r19070 = 0.0694555789232254f; float r19071 = r19069 + r19070; float r19072 = r19045 * r19071; float r19073 = 0.29097387194633484f; float r19074 = r19043 * r19073; float r19075 = r19074 * r19046; float r19076 = r19072 + r19075; float r19077 = r19044 * r19044; float r19078 = r19077 * r19045; float r19079 = 2.0f; float r19080 = r19079 * r19047; float r19081 = r19044 * r19080; float r19082 = 0.0008327945251949131f; float r19083 = r19081 + r19082; float r19084 = r19078 * r19083; float r19085 = 0.7715470790863037f; float r19086 = r19044 * r19085; float r19087 = r19060 + r19086; float r19088 = r19084 + r19087; float r19089 = r19076 + r19088; float r19090 = r19043 / r19089; float r19091 = r19066 * r19090; return r19091; } double f_od(double x) { double r19092 = x; double r19093 = r19092 * r19092; double r19094 = r19093 * (r19093 * r19093); double r19095 = r19092 * (r19092 * r19092); double r19096 = 0.00017899709928315133; double r19097 = r19092 * r19096; double r19098 = r19095 * r19097; double r19099 = 0.007264418061822653; double r19100 = r19098 + r19099; double r19101 = 0.0005064033903181553; double r19102 = r19092 * r19101; double r19103 = r19092 * r19102; double r19104 = r19100 + r19103; double r19105 = r19094 * r19104; double r19106 = 0.10499349236488342; double r19107 = r19106 * r19092; double r19108 = r19092 * r19107; double r19109 = 1.0; double r19110 = 0.04240605980157852; double r19111 = r19110 * r19092; double r19112 = r19111 * r19095; double r19113 = r19109 + r19112; double r19114 = r19108 + r19113; double r19115 = r19105 + r19114; double r19116 = 0.01400054432451725; double r19117 = r19116 * r19092; double r19118 = r19092 * r19117; double r19119 = 0.0694555789232254; double r19120 = r19118 + r19119; double r19121 = r19094 * r19120; double r19122 = 0.29097387194633484; double r19123 = r19092 * r19122; double r19124 = r19123 * r19095; double r19125 = r19121 + r19124; double r19126 = r19093 * r19093; double r19127 = r19126 * r19094; double r19128 = 2.0; double r19129 = r19128 * r19096; double r19130 = r19093 * r19129; double r19131 = 0.0008327945251949131; double r19132 = r19130 + r19131; double r19133 = r19127 * r19132; double r19134 = 0.7715470790863037; double r19135 = r19093 * r19134; double r19136 = r19109 + r19135; double r19137 = r19133 + r19136; double r19138 = r19125 + r19137; double r19139 = r19092 / r19138; double r19140 = r19115 * r19139; return r19140; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19141, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19142, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19143); mpfr_init(r19144); mpfr_init(r19145); mpfr_init(r19146); mpfr_init_set_str(r19147, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19148); mpfr_init(r19149); mpfr_init(r19150); mpfr_init_set_str(r19151, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19152); mpfr_init(r19153); mpfr_init(r19154); mpfr_init_set_str(r19155, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19156); mpfr_init(r19157); mpfr_init(r19158); mpfr_init_set_str(r19159, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19160); mpfr_init(r19161); mpfr_init(r19162); mpfr_init_set_str(r19163, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19164); mpfr_init(r19165); mpfr_init_set_str(r19166, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19167); mpfr_init(r19168); mpfr_init_set_str(r19169, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19170); mpfr_init(r19171); mpfr_init_set_str(r19172, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19173); mpfr_init(r19174); mpfr_init_set_str(r19175, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19176); mpfr_init(r19177); mpfr_init_set_str(r19178, "2", 10, MPFR_RNDN); mpfr_init(r19179); mpfr_init(r19180); mpfr_init(r19181); mpfr_init(r19182); mpfr_init(r19183); mpfr_init(r19184); } double f_im(double x) { ; ; mpfr_set_d(r19143, x, MPFR_RNDN); mpfr_mul(r19144, r19143, r19143, MPFR_RNDN); mpfr_mul(r19145, r19142, r19144, MPFR_RNDN); mpfr_add(r19146, r19141, r19145, MPFR_RNDN); ; mpfr_mul(r19148, r19144, r19144, MPFR_RNDN); mpfr_mul(r19149, r19147, r19148, MPFR_RNDN); mpfr_add(r19150, r19146, r19149, MPFR_RNDN); ; mpfr_mul(r19152, r19148, r19144, MPFR_RNDN); mpfr_mul(r19153, r19151, r19152, MPFR_RNDN); mpfr_add(r19154, r19150, r19153, MPFR_RNDN); ; mpfr_mul(r19156, r19152, r19144, MPFR_RNDN); mpfr_mul(r19157, r19155, r19156, MPFR_RNDN); mpfr_add(r19158, r19154, r19157, MPFR_RNDN); ; mpfr_mul(r19160, r19156, r19144, MPFR_RNDN); mpfr_mul(r19161, r19159, r19160, MPFR_RNDN); mpfr_add(r19162, r19158, r19161, MPFR_RNDN); ; mpfr_mul(r19164, r19163, r19144, MPFR_RNDN); mpfr_add(r19165, r19141, r19164, MPFR_RNDN); ; mpfr_mul(r19167, r19166, r19148, MPFR_RNDN); mpfr_add(r19168, r19165, r19167, MPFR_RNDN); ; mpfr_mul(r19170, r19169, r19152, MPFR_RNDN); mpfr_add(r19171, r19168, r19170, MPFR_RNDN); ; mpfr_mul(r19173, r19172, r19156, MPFR_RNDN); mpfr_add(r19174, r19171, r19173, MPFR_RNDN); ; mpfr_mul(r19176, r19175, r19160, MPFR_RNDN); mpfr_add(r19177, r19174, r19176, MPFR_RNDN); ; mpfr_mul(r19179, r19178, r19159, MPFR_RNDN); mpfr_mul(r19180, r19160, r19144, MPFR_RNDN); mpfr_mul(r19181, r19179, r19180, MPFR_RNDN); mpfr_add(r19182, r19177, r19181, MPFR_RNDN); mpfr_div(r19183, r19162, r19182, MPFR_RNDN); mpfr_mul(r19184, r19183, r19143, MPFR_RNDN); return mpfr_get_d(r19184, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19185); mpfr_init(r19186); mpfr_init(r19187); mpfr_init(r19188); mpfr_init_set_str(r19189, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19190); mpfr_init(r19191); mpfr_init_set_str(r19192, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19193); mpfr_init_set_str(r19194, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19195); mpfr_init(r19196); mpfr_init(r19197); mpfr_init(r19198); mpfr_init_set_str(r19199, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19200); mpfr_init(r19201); mpfr_init_set_str(r19202, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19203, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19204); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); mpfr_init(r19208); mpfr_init_set_str(r19209, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19210); mpfr_init(r19211); mpfr_init_set_str(r19212, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19213); mpfr_init(r19214); mpfr_init_set_str(r19215, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19216); mpfr_init(r19217); mpfr_init(r19218); mpfr_init(r19219); mpfr_init(r19220); mpfr_init_set_str(r19221, "2", 10, MPFR_RNDN); mpfr_init(r19222); mpfr_init(r19223); mpfr_init_set_str(r19224, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19225); mpfr_init(r19226); mpfr_init_set_str(r19227, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19228); mpfr_init(r19229); mpfr_init(r19230); mpfr_init(r19231); mpfr_init(r19232); mpfr_init(r19233); } double f_fm(double x) { mpfr_set_d(r19185, x, MPFR_RNDN); mpfr_sqr(r19186, r19185, MPFR_RNDN); mpfr_mul(r19187, r19186, r19186, MPFR_RNDN); mpfr_mul(r19187, r19187, r19186, MPFR_RNDN); mpfr_mul(r19188, r19185, r19185, MPFR_RNDN); mpfr_mul(r19188, r19188, r19185, MPFR_RNDN); ; mpfr_mul(r19190, r19185, r19189, MPFR_RNDN); mpfr_mul(r19191, r19188, r19190, MPFR_RNDN); ; mpfr_add(r19193, r19191, r19192, MPFR_RNDN); ; mpfr_mul(r19195, r19185, r19194, MPFR_RNDN); mpfr_mul(r19196, r19185, r19195, MPFR_RNDN); mpfr_add(r19197, r19193, r19196, MPFR_RNDN); mpfr_mul(r19198, r19187, r19197, MPFR_RNDN); ; mpfr_mul(r19200, r19199, r19185, MPFR_RNDN); mpfr_mul(r19201, r19185, r19200, MPFR_RNDN); ; ; mpfr_mul(r19204, r19203, r19185, MPFR_RNDN); mpfr_mul(r19205, r19204, r19188, MPFR_RNDN); mpfr_add(r19206, r19202, r19205, MPFR_RNDN); mpfr_add(r19207, r19201, r19206, MPFR_RNDN); mpfr_add(r19208, r19198, r19207, MPFR_RNDN); ; mpfr_mul(r19210, r19209, r19185, MPFR_RNDN); mpfr_mul(r19211, r19185, r19210, MPFR_RNDN); ; mpfr_add(r19213, r19211, r19212, MPFR_RNDN); mpfr_mul(r19214, r19187, r19213, MPFR_RNDN); ; mpfr_mul(r19216, r19185, r19215, MPFR_RNDN); mpfr_mul(r19217, r19216, r19188, MPFR_RNDN); mpfr_add(r19218, r19214, r19217, MPFR_RNDN); mpfr_mul(r19219, r19186, r19186, MPFR_RNDN); mpfr_mul(r19220, r19219, r19187, MPFR_RNDN); ; mpfr_mul(r19222, r19221, r19189, MPFR_RNDN); mpfr_mul(r19223, r19186, r19222, MPFR_RNDN); ; mpfr_add(r19225, r19223, r19224, MPFR_RNDN); mpfr_mul(r19226, r19220, r19225, MPFR_RNDN); ; mpfr_mul(r19228, r19186, r19227, MPFR_RNDN); mpfr_add(r19229, r19202, r19228, MPFR_RNDN); mpfr_add(r19230, r19226, r19229, MPFR_RNDN); mpfr_add(r19231, r19218, r19230, MPFR_RNDN); mpfr_div(r19232, r19185, r19231, MPFR_RNDN); mpfr_mul(r19233, r19208, r19232, MPFR_RNDN); return mpfr_get_d(r19233, MPFR_RNDN); } static mpfr_t 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, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19234); mpfr_init(r19235); mpfr_init(r19236); mpfr_init(r19237); mpfr_init_set_str(r19238, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19239); mpfr_init(r19240); mpfr_init_set_str(r19241, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19242); mpfr_init_set_str(r19243, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19244); mpfr_init(r19245); mpfr_init(r19246); mpfr_init(r19247); mpfr_init_set_str(r19248, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19249); mpfr_init(r19250); mpfr_init_set_str(r19251, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19252, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19253); mpfr_init(r19254); mpfr_init(r19255); mpfr_init(r19256); mpfr_init(r19257); mpfr_init_set_str(r19258, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19259); mpfr_init(r19260); mpfr_init_set_str(r19261, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19262); mpfr_init(r19263); mpfr_init_set_str(r19264, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19265); mpfr_init(r19266); mpfr_init(r19267); mpfr_init(r19268); mpfr_init(r19269); mpfr_init_set_str(r19270, "2", 10, MPFR_RNDN); mpfr_init(r19271); mpfr_init(r19272); mpfr_init_set_str(r19273, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19274); mpfr_init(r19275); mpfr_init_set_str(r19276, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19277); mpfr_init(r19278); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init(r19282); } double f_dm(double x) { mpfr_set_d(r19234, x, MPFR_RNDN); mpfr_sqr(r19235, r19234, MPFR_RNDN); mpfr_mul(r19236, r19235, r19235, MPFR_RNDN); mpfr_mul(r19236, r19236, r19235, MPFR_RNDN); mpfr_mul(r19237, r19234, r19234, MPFR_RNDN); mpfr_mul(r19237, r19237, r19234, MPFR_RNDN); ; mpfr_mul(r19239, r19234, r19238, MPFR_RNDN); mpfr_mul(r19240, r19237, r19239, MPFR_RNDN); ; mpfr_add(r19242, r19240, r19241, MPFR_RNDN); ; mpfr_mul(r19244, r19234, r19243, MPFR_RNDN); mpfr_mul(r19245, r19234, r19244, MPFR_RNDN); mpfr_add(r19246, r19242, r19245, MPFR_RNDN); mpfr_mul(r19247, r19236, r19246, MPFR_RNDN); ; mpfr_mul(r19249, r19248, r19234, MPFR_RNDN); mpfr_mul(r19250, r19234, r19249, MPFR_RNDN); ; ; mpfr_mul(r19253, r19252, r19234, MPFR_RNDN); mpfr_mul(r19254, r19253, r19237, MPFR_RNDN); mpfr_add(r19255, r19251, r19254, MPFR_RNDN); mpfr_add(r19256, r19250, r19255, MPFR_RNDN); mpfr_add(r19257, r19247, r19256, MPFR_RNDN); ; mpfr_mul(r19259, r19258, r19234, MPFR_RNDN); mpfr_mul(r19260, r19234, r19259, MPFR_RNDN); ; mpfr_add(r19262, r19260, r19261, MPFR_RNDN); mpfr_mul(r19263, r19236, r19262, MPFR_RNDN); ; mpfr_mul(r19265, r19234, r19264, MPFR_RNDN); mpfr_mul(r19266, r19265, r19237, MPFR_RNDN); mpfr_add(r19267, r19263, r19266, MPFR_RNDN); mpfr_mul(r19268, r19235, r19235, MPFR_RNDN); mpfr_mul(r19269, r19268, r19236, MPFR_RNDN); ; mpfr_mul(r19271, r19270, r19238, MPFR_RNDN); mpfr_mul(r19272, r19235, r19271, MPFR_RNDN); ; mpfr_add(r19274, r19272, r19273, MPFR_RNDN); mpfr_mul(r19275, r19269, r19274, MPFR_RNDN); ; mpfr_mul(r19277, r19235, r19276, MPFR_RNDN); mpfr_add(r19278, r19251, r19277, MPFR_RNDN); mpfr_add(r19279, r19275, r19278, MPFR_RNDN); mpfr_add(r19280, r19267, r19279, MPFR_RNDN); mpfr_div(r19281, r19234, r19280, MPFR_RNDN); mpfr_mul(r19282, r19257, r19281, MPFR_RNDN); return mpfr_get_d(r19282, MPFR_RNDN); }