#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18912 = 1.0f; float r18913 = 0.10499349236488342f; float r18914 = x; float r18915 = r18914 * r18914; float r18916 = r18913 * r18915; float r18917 = r18912 + r18916; float r18918 = 0.04240605980157852f; float r18919 = r18915 * r18915; float r18920 = r18918 * r18919; float r18921 = r18917 + r18920; float r18922 = 0.007264418061822653f; float r18923 = r18919 * r18915; float r18924 = r18922 * r18923; float r18925 = r18921 + r18924; float r18926 = 0.0005064033903181553f; float r18927 = r18923 * r18915; float r18928 = r18926 * r18927; float r18929 = r18925 + r18928; float r18930 = 0.00017899709928315133f; float r18931 = r18927 * r18915; float r18932 = r18930 * r18931; float r18933 = r18929 + r18932; float r18934 = 0.7715470790863037f; float r18935 = r18934 * r18915; float r18936 = r18912 + r18935; float r18937 = 0.29097387194633484f; float r18938 = r18937 * r18919; float r18939 = r18936 + r18938; float r18940 = 0.0694555789232254f; float r18941 = r18940 * r18923; float r18942 = r18939 + r18941; float r18943 = 0.01400054432451725f; float r18944 = r18943 * r18927; float r18945 = r18942 + r18944; float r18946 = 0.0008327945251949131f; float r18947 = r18946 * r18931; float r18948 = r18945 + r18947; float r18949 = 2.0f; float r18950 = r18949 * r18930; float r18951 = r18931 * r18915; float r18952 = r18950 * r18951; float r18953 = r18948 + r18952; float r18954 = r18933 / r18953; float r18955 = r18954 * r18914; return r18955; } double f_id(double x) { double r18956 = 1.0; double r18957 = 0.10499349236488342; double r18958 = x; double r18959 = r18958 * r18958; double r18960 = r18957 * r18959; double r18961 = r18956 + r18960; double r18962 = 0.04240605980157852; double r18963 = r18959 * r18959; double r18964 = r18962 * r18963; double r18965 = r18961 + r18964; double r18966 = 0.007264418061822653; double r18967 = r18963 * r18959; double r18968 = r18966 * r18967; double r18969 = r18965 + r18968; double r18970 = 0.0005064033903181553; double r18971 = r18967 * r18959; double r18972 = r18970 * r18971; double r18973 = r18969 + r18972; double r18974 = 0.00017899709928315133; double r18975 = r18971 * r18959; double r18976 = r18974 * r18975; double r18977 = r18973 + r18976; double r18978 = 0.7715470790863037; double r18979 = r18978 * r18959; double r18980 = r18956 + r18979; double r18981 = 0.29097387194633484; double r18982 = r18981 * r18963; double r18983 = r18980 + r18982; double r18984 = 0.0694555789232254; double r18985 = r18984 * r18967; double r18986 = r18983 + r18985; double r18987 = 0.01400054432451725; double r18988 = r18987 * r18971; double r18989 = r18986 + r18988; double r18990 = 0.0008327945251949131; double r18991 = r18990 * r18975; double r18992 = r18989 + r18991; double r18993 = 2.0; double r18994 = r18993 * r18974; double r18995 = r18975 * r18959; double r18996 = r18994 * r18995; double r18997 = r18992 + r18996; double r18998 = r18977 / r18997; double r18999 = r18998 * r18958; return r18999; } double f_of(float x) { float r19000 = x; float r19001 = -11433.3388671875f; bool r19002 = r19000 <= r19001; float r19003 = 1.0f; float r19004 = 0.10499349236488342f; float r19005 = r19000 * r19000; float r19006 = r19004 / r19005; float r19007 = r19003 + r19006; float r19008 = 0.007264418061822653f; float r19009 = r19008 * r19003; float r19010 = r19000 * (r19000 * r19000); float r19011 = r19010 * r19010; float r19012 = r19009 / r19011; float r19013 = r19007 + r19012; float r19014 = 0.04240605980157852f; float r19015 = r19014 * r19003; float r19016 = r19015 / r19010; float r19017 = r19016 / r19000; float r19018 = 0.00017899709928315133f; float r19019 = 10.0f; float r19020 = pow(r19000, r19019); float r19021 = r19018 / r19020; float r19022 = r19017 + r19021; float r19023 = 0.0005064033903181553f; float r19024 = r19003 * r19023; float r19025 = r19024 / r19011; float r19026 = r19025 / r19005; float r19027 = r19022 + r19026; float r19028 = r19013 + r19027; float r19029 = r19003 / r19000; float r19030 = r19029 / r19010; float r19031 = 0.29097387194633484f; float r19032 = r19030 * r19031; float r19033 = r19003 + r19032; float r19034 = 0.0694555789232254f; float r19035 = r19034 * r19003; float r19036 = r19035 / r19011; float r19037 = r19033 + r19036; float r19038 = 0.01400054432451725f; float r19039 = r19003 * r19038; float r19040 = r19039 / r19011; float r19041 = r19040 / r19005; float r19042 = r19003 / r19010; float r19043 = r19042 * r19042; float r19044 = 6.0f; float r19045 = pow(r19000, r19044); float r19046 = 0.00035799419856630266f; float r19047 = r19045 / r19046; float r19048 = r19043 / r19047; float r19049 = r19041 + r19048; float r19050 = r19037 + r19049; float r19051 = 0.7715470790863037f; float r19052 = r19051 / r19005; float r19053 = 4.0f; float r19054 = pow(r19000, r19053); float r19055 = 0.0008327945251949131f; float r19056 = r19054 / r19055; float r19057 = r19043 / r19056; float r19058 = r19052 + r19057; float r19059 = r19050 + r19058; float r19060 = r19059 * r19000; float r19061 = r19028 / r19060; float r19062 = 33073.84375f; bool r19063 = r19000 <= r19062; float r19064 = r19010 * r19010; float r19065 = r19000 * r19000; float r19066 = r19023 * r19065; float r19067 = r19008 + r19066; float r19068 = r19064 * r19067; float r19069 = r19000 * r19014; float r19070 = r19069 * r19010; float r19071 = r19003 + r19070; float r19072 = r19004 * r19065; float r19073 = r19071 + r19072; float r19074 = r19068 + r19073; float r19075 = r19065 * r19065; float r19076 = r19064 * r19075; float r19077 = r19076 * r19018; float r19078 = r19074 + r19077; float r19079 = r19054 * r19055; float r19080 = r19064 * r19079; float r19081 = r19003 + r19080; float r19082 = r19038 * r19000; float r19083 = r19000 * r19082; float r19084 = r19083 + r19034; float r19085 = r19064 * r19084; float r19086 = r19081 + r19085; float r19087 = r19045 * r19046; float r19088 = r19087 * r19064; float r19089 = r19031 * r19075; float r19090 = r19000 * r19051; float r19091 = r19090 * r19000; float r19092 = r19089 + r19091; float r19093 = r19088 + r19092; float r19094 = r19086 + r19093; float r19095 = r19000 / r19094; float r19096 = r19078 * r19095; float r19097 = r19063 ? r19096 : r19061; float r19098 = r19002 ? r19061 : r19097; return r19098; } double f_od(double x) { double r19099 = x; double r19100 = -11433.3388671875; bool r19101 = r19099 <= r19100; double r19102 = 1.0; double r19103 = 0.10499349236488342; double r19104 = r19099 * r19099; double r19105 = r19103 / r19104; double r19106 = r19102 + r19105; double r19107 = 0.007264418061822653; double r19108 = r19107 * r19102; double r19109 = r19099 * (r19099 * r19099); double r19110 = r19109 * r19109; double r19111 = r19108 / r19110; double r19112 = r19106 + r19111; double r19113 = 0.04240605980157852; double r19114 = r19113 * r19102; double r19115 = r19114 / r19109; double r19116 = r19115 / r19099; double r19117 = 0.00017899709928315133; double r19118 = 10.0; double r19119 = pow(r19099, r19118); double r19120 = r19117 / r19119; double r19121 = r19116 + r19120; double r19122 = 0.0005064033903181553; double r19123 = r19102 * r19122; double r19124 = r19123 / r19110; double r19125 = r19124 / r19104; double r19126 = r19121 + r19125; double r19127 = r19112 + r19126; double r19128 = r19102 / r19099; double r19129 = r19128 / r19109; double r19130 = 0.29097387194633484; double r19131 = r19129 * r19130; double r19132 = r19102 + r19131; double r19133 = 0.0694555789232254; double r19134 = r19133 * r19102; double r19135 = r19134 / r19110; double r19136 = r19132 + r19135; double r19137 = 0.01400054432451725; double r19138 = r19102 * r19137; double r19139 = r19138 / r19110; double r19140 = r19139 / r19104; double r19141 = r19102 / r19109; double r19142 = r19141 * r19141; double r19143 = 6.0; double r19144 = pow(r19099, r19143); double r19145 = 0.00035799419856630266; double r19146 = r19144 / r19145; double r19147 = r19142 / r19146; double r19148 = r19140 + r19147; double r19149 = r19136 + r19148; double r19150 = 0.7715470790863037; double r19151 = r19150 / r19104; double r19152 = 4.0; double r19153 = pow(r19099, r19152); double r19154 = 0.0008327945251949131; double r19155 = r19153 / r19154; double r19156 = r19142 / r19155; double r19157 = r19151 + r19156; double r19158 = r19149 + r19157; double r19159 = r19158 * r19099; double r19160 = r19127 / r19159; double r19161 = 33073.84375; bool r19162 = r19099 <= r19161; double r19163 = r19109 * r19109; double r19164 = r19099 * r19099; double r19165 = r19122 * r19164; double r19166 = r19107 + r19165; double r19167 = r19163 * r19166; double r19168 = r19099 * r19113; double r19169 = r19168 * r19109; double r19170 = r19102 + r19169; double r19171 = r19103 * r19164; double r19172 = r19170 + r19171; double r19173 = r19167 + r19172; double r19174 = r19164 * r19164; double r19175 = r19163 * r19174; double r19176 = r19175 * r19117; double r19177 = r19173 + r19176; double r19178 = r19153 * r19154; double r19179 = r19163 * r19178; double r19180 = r19102 + r19179; double r19181 = r19137 * r19099; double r19182 = r19099 * r19181; double r19183 = r19182 + r19133; double r19184 = r19163 * r19183; double r19185 = r19180 + r19184; double r19186 = r19144 * r19145; double r19187 = r19186 * r19163; double r19188 = r19130 * r19174; double r19189 = r19099 * r19150; double r19190 = r19189 * r19099; double r19191 = r19188 + r19190; double r19192 = r19187 + r19191; double r19193 = r19185 + r19192; double r19194 = r19099 / r19193; double r19195 = r19177 * r19194; double r19196 = r19162 ? r19195 : r19160; double r19197 = r19101 ? r19160 : r19196; return r19197; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19198, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19199, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19200); mpfr_init(r19201); mpfr_init(r19202); mpfr_init(r19203); mpfr_init_set_str(r19204, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); mpfr_init_set_str(r19208, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19209); mpfr_init(r19210); mpfr_init(r19211); mpfr_init_set_str(r19212, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19213); mpfr_init(r19214); mpfr_init(r19215); mpfr_init_set_str(r19216, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19217); mpfr_init(r19218); mpfr_init(r19219); mpfr_init_set_str(r19220, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19221); mpfr_init(r19222); mpfr_init_set_str(r19223, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19224); mpfr_init(r19225); mpfr_init_set_str(r19226, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19227); mpfr_init(r19228); mpfr_init_set_str(r19229, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19230); mpfr_init(r19231); mpfr_init_set_str(r19232, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19233); mpfr_init(r19234); mpfr_init_set_str(r19235, "2", 10, MPFR_RNDN); mpfr_init(r19236); mpfr_init(r19237); mpfr_init(r19238); mpfr_init(r19239); mpfr_init(r19240); mpfr_init(r19241); } double f_im(double x) { ; ; mpfr_set_d(r19200, x, MPFR_RNDN); mpfr_mul(r19201, r19200, r19200, MPFR_RNDN); mpfr_mul(r19202, r19199, r19201, MPFR_RNDN); mpfr_add(r19203, r19198, r19202, MPFR_RNDN); ; mpfr_mul(r19205, r19201, r19201, MPFR_RNDN); mpfr_mul(r19206, r19204, r19205, MPFR_RNDN); mpfr_add(r19207, r19203, r19206, MPFR_RNDN); ; mpfr_mul(r19209, r19205, r19201, MPFR_RNDN); mpfr_mul(r19210, r19208, r19209, MPFR_RNDN); mpfr_add(r19211, r19207, r19210, MPFR_RNDN); ; mpfr_mul(r19213, r19209, r19201, MPFR_RNDN); mpfr_mul(r19214, r19212, r19213, MPFR_RNDN); mpfr_add(r19215, r19211, r19214, MPFR_RNDN); ; mpfr_mul(r19217, r19213, r19201, MPFR_RNDN); mpfr_mul(r19218, r19216, r19217, MPFR_RNDN); mpfr_add(r19219, r19215, r19218, MPFR_RNDN); ; mpfr_mul(r19221, r19220, r19201, MPFR_RNDN); mpfr_add(r19222, r19198, r19221, MPFR_RNDN); ; mpfr_mul(r19224, r19223, r19205, MPFR_RNDN); mpfr_add(r19225, r19222, r19224, MPFR_RNDN); ; mpfr_mul(r19227, r19226, r19209, MPFR_RNDN); mpfr_add(r19228, r19225, r19227, MPFR_RNDN); ; mpfr_mul(r19230, r19229, r19213, MPFR_RNDN); mpfr_add(r19231, r19228, r19230, MPFR_RNDN); ; mpfr_mul(r19233, r19232, r19217, MPFR_RNDN); mpfr_add(r19234, r19231, r19233, MPFR_RNDN); ; mpfr_mul(r19236, r19235, r19216, MPFR_RNDN); mpfr_mul(r19237, r19217, r19201, MPFR_RNDN); mpfr_mul(r19238, r19236, r19237, MPFR_RNDN); mpfr_add(r19239, r19234, r19238, MPFR_RNDN); mpfr_div(r19240, r19219, r19239, MPFR_RNDN); mpfr_mul(r19241, r19240, r19200, MPFR_RNDN); return mpfr_get_d(r19241, MPFR_RNDN); } static mpfr_t 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, r19283, r19284, r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328, r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336, r19337, r19338, r19339, r19340; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19242); mpfr_init_set_str(r19243, "-11433.339f0", 10, MPFR_RNDN); mpfr_init(r19244); mpfr_init_set_str(r19245, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19246, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19247); mpfr_init(r19248); mpfr_init(r19249); mpfr_init_set_str(r19250, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19251); mpfr_init(r19252); mpfr_init(r19253); mpfr_init(r19254); mpfr_init(r19255); mpfr_init_set_str(r19256, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19257); mpfr_init(r19258); mpfr_init(r19259); mpfr_init_set_str(r19260, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19261, "10", 10, MPFR_RNDN); mpfr_init(r19262); mpfr_init(r19263); mpfr_init(r19264); mpfr_init_set_str(r19265, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19266); mpfr_init(r19267); mpfr_init(r19268); mpfr_init(r19269); mpfr_init(r19270); mpfr_init(r19271); mpfr_init(r19272); mpfr_init_set_str(r19273, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19274); mpfr_init(r19275); mpfr_init_set_str(r19276, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19277); mpfr_init(r19278); mpfr_init(r19279); mpfr_init_set_str(r19280, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19281); mpfr_init(r19282); mpfr_init(r19283); mpfr_init(r19284); mpfr_init(r19285); mpfr_init_set_str(r19286, "6", 10, MPFR_RNDN); mpfr_init(r19287); mpfr_init_set_str(r19288, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); mpfr_init(r19292); mpfr_init_set_str(r19293, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19294); mpfr_init_set_str(r19295, "4", 10, MPFR_RNDN); mpfr_init(r19296); mpfr_init_set_str(r19297, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19298); mpfr_init(r19299); mpfr_init(r19300); mpfr_init(r19301); mpfr_init(r19302); mpfr_init(r19303); mpfr_init_set_str(r19304, "33073.844f0", 10, MPFR_RNDN); mpfr_init(r19305); mpfr_init(r19306); mpfr_init(r19307); mpfr_init(r19308); mpfr_init(r19309); mpfr_init(r19310); mpfr_init(r19311); mpfr_init(r19312); mpfr_init(r19313); mpfr_init(r19314); mpfr_init(r19315); mpfr_init(r19316); mpfr_init(r19317); mpfr_init(r19318); mpfr_init(r19319); mpfr_init(r19320); mpfr_init(r19321); mpfr_init(r19322); mpfr_init(r19323); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init(r19327); mpfr_init(r19328); mpfr_init(r19329); mpfr_init(r19330); mpfr_init(r19331); mpfr_init(r19332); mpfr_init(r19333); mpfr_init(r19334); mpfr_init(r19335); mpfr_init(r19336); mpfr_init(r19337); mpfr_init(r19338); mpfr_init(r19339); mpfr_init(r19340); } double f_fm(double x) { mpfr_set_d(r19242, x, MPFR_RNDN); ; mpfr_set_si(r19244, mpfr_cmp(r19242, r19243) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19247, r19242, r19242, MPFR_RNDN); mpfr_div(r19248, r19246, r19247, MPFR_RNDN); mpfr_add(r19249, r19245, r19248, MPFR_RNDN); ; mpfr_mul(r19251, r19250, r19245, MPFR_RNDN); mpfr_mul(r19252, r19242, r19242, MPFR_RNDN); mpfr_mul(r19252, r19252, r19242, MPFR_RNDN); mpfr_sqr(r19253, r19252, MPFR_RNDN); mpfr_div(r19254, r19251, r19253, MPFR_RNDN); mpfr_add(r19255, r19249, r19254, MPFR_RNDN); ; mpfr_mul(r19257, r19256, r19245, MPFR_RNDN); mpfr_div(r19258, r19257, r19252, MPFR_RNDN); mpfr_div(r19259, r19258, r19242, MPFR_RNDN); ; ; mpfr_pow(r19262, r19242, r19261, MPFR_RNDN); mpfr_div(r19263, r19260, r19262, MPFR_RNDN); mpfr_add(r19264, r19259, r19263, MPFR_RNDN); ; mpfr_mul(r19266, r19245, r19265, MPFR_RNDN); mpfr_div(r19267, r19266, r19253, MPFR_RNDN); mpfr_div(r19268, r19267, r19247, MPFR_RNDN); mpfr_add(r19269, r19264, r19268, MPFR_RNDN); mpfr_add(r19270, r19255, r19269, MPFR_RNDN); mpfr_div(r19271, r19245, r19242, MPFR_RNDN); mpfr_div(r19272, r19271, r19252, MPFR_RNDN); ; mpfr_mul(r19274, r19272, r19273, MPFR_RNDN); mpfr_add(r19275, r19245, r19274, MPFR_RNDN); ; mpfr_mul(r19277, r19276, r19245, MPFR_RNDN); mpfr_div(r19278, r19277, r19253, MPFR_RNDN); mpfr_add(r19279, r19275, r19278, MPFR_RNDN); ; mpfr_mul(r19281, r19245, r19280, MPFR_RNDN); mpfr_div(r19282, r19281, r19253, MPFR_RNDN); mpfr_div(r19283, r19282, r19247, MPFR_RNDN); mpfr_div(r19284, r19245, r19252, MPFR_RNDN); mpfr_mul(r19285, r19284, r19284, MPFR_RNDN); ; mpfr_pow(r19287, r19242, r19286, MPFR_RNDN); ; mpfr_div(r19289, r19287, r19288, MPFR_RNDN); mpfr_div(r19290, r19285, r19289, MPFR_RNDN); mpfr_add(r19291, r19283, r19290, MPFR_RNDN); mpfr_add(r19292, r19279, r19291, MPFR_RNDN); ; mpfr_div(r19294, r19293, r19247, MPFR_RNDN); ; mpfr_pow(r19296, r19242, r19295, MPFR_RNDN); ; mpfr_div(r19298, r19296, r19297, MPFR_RNDN); mpfr_div(r19299, r19285, r19298, MPFR_RNDN); mpfr_add(r19300, r19294, r19299, MPFR_RNDN); mpfr_add(r19301, r19292, r19300, MPFR_RNDN); mpfr_mul(r19302, r19301, r19242, MPFR_RNDN); mpfr_div(r19303, r19270, r19302, MPFR_RNDN); ; mpfr_set_si(r19305, mpfr_cmp(r19242, r19304) <= 0, MPFR_RNDN); mpfr_mul(r19306, r19252, r19252, MPFR_RNDN); mpfr_sqr(r19307, r19242, MPFR_RNDN); mpfr_mul(r19308, r19265, r19307, MPFR_RNDN); mpfr_add(r19309, r19250, r19308, MPFR_RNDN); mpfr_mul(r19310, r19306, r19309, MPFR_RNDN); mpfr_mul(r19311, r19242, r19256, MPFR_RNDN); mpfr_mul(r19312, r19311, r19252, MPFR_RNDN); mpfr_add(r19313, r19245, r19312, MPFR_RNDN); mpfr_mul(r19314, r19246, r19307, MPFR_RNDN); mpfr_add(r19315, r19313, r19314, MPFR_RNDN); mpfr_add(r19316, r19310, r19315, MPFR_RNDN); mpfr_sqr(r19317, r19307, MPFR_RNDN); mpfr_mul(r19318, r19306, r19317, MPFR_RNDN); mpfr_mul(r19319, r19318, r19260, MPFR_RNDN); mpfr_add(r19320, r19316, r19319, MPFR_RNDN); mpfr_mul(r19321, r19296, r19297, MPFR_RNDN); mpfr_mul(r19322, r19306, r19321, MPFR_RNDN); mpfr_add(r19323, r19245, r19322, MPFR_RNDN); mpfr_mul(r19324, r19280, r19242, MPFR_RNDN); mpfr_mul(r19325, r19242, r19324, MPFR_RNDN); mpfr_add(r19326, r19325, r19276, MPFR_RNDN); mpfr_mul(r19327, r19306, r19326, MPFR_RNDN); mpfr_add(r19328, r19323, r19327, MPFR_RNDN); mpfr_mul(r19329, r19287, r19288, MPFR_RNDN); mpfr_mul(r19330, r19329, r19306, MPFR_RNDN); mpfr_mul(r19331, r19273, r19317, MPFR_RNDN); mpfr_mul(r19332, r19242, r19293, MPFR_RNDN); mpfr_mul(r19333, r19332, r19242, MPFR_RNDN); mpfr_add(r19334, r19331, r19333, MPFR_RNDN); mpfr_add(r19335, r19330, r19334, MPFR_RNDN); mpfr_add(r19336, r19328, r19335, MPFR_RNDN); mpfr_div(r19337, r19242, r19336, MPFR_RNDN); mpfr_mul(r19338, r19320, r19337, MPFR_RNDN); if (mpfr_get_si(r19305, MPFR_RNDN)) { mpfr_set(r19339, r19338, MPFR_RNDN); } else { mpfr_set(r19339, r19303, MPFR_RNDN); }; if (mpfr_get_si(r19244, MPFR_RNDN)) { mpfr_set(r19340, r19303, MPFR_RNDN); } else { mpfr_set(r19340, r19339, MPFR_RNDN); }; return mpfr_get_d(r19340, MPFR_RNDN); } static mpfr_t r19341, r19342, r19343, r19344, r19345, r19346, r19347, r19348, r19349, r19350, r19351, r19352, r19353, r19354, r19355, r19356, r19357, r19358, r19359, r19360, r19361, r19362, r19363, r19364, r19365, r19366, r19367, r19368, r19369, r19370, r19371, r19372, r19373, r19374, r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414, r19415, r19416, r19417, r19418, r19419, r19420, r19421, r19422, r19423, r19424, r19425, r19426, r19427, r19428, r19429, r19430, r19431, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19341); mpfr_init_set_str(r19342, "-11433.339f0", 10, MPFR_RNDN); mpfr_init(r19343); mpfr_init_set_str(r19344, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19345, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19346); mpfr_init(r19347); mpfr_init(r19348); mpfr_init_set_str(r19349, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19350); mpfr_init(r19351); mpfr_init(r19352); mpfr_init(r19353); mpfr_init(r19354); mpfr_init_set_str(r19355, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19356); mpfr_init(r19357); mpfr_init(r19358); mpfr_init_set_str(r19359, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19360, "10", 10, MPFR_RNDN); mpfr_init(r19361); mpfr_init(r19362); mpfr_init(r19363); mpfr_init_set_str(r19364, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19365); mpfr_init(r19366); mpfr_init(r19367); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init_set_str(r19372, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19373); mpfr_init(r19374); mpfr_init_set_str(r19375, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19376); mpfr_init(r19377); mpfr_init(r19378); mpfr_init_set_str(r19379, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init(r19384); mpfr_init_set_str(r19385, "6", 10, MPFR_RNDN); mpfr_init(r19386); mpfr_init_set_str(r19387, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init_set_str(r19392, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19393); mpfr_init_set_str(r19394, "4", 10, MPFR_RNDN); mpfr_init(r19395); mpfr_init_set_str(r19396, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init_set_str(r19403, "33073.844f0", 10, MPFR_RNDN); mpfr_init(r19404); mpfr_init(r19405); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init(r19415); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); mpfr_init(r19422); mpfr_init(r19423); mpfr_init(r19424); mpfr_init(r19425); mpfr_init(r19426); mpfr_init(r19427); mpfr_init(r19428); mpfr_init(r19429); mpfr_init(r19430); mpfr_init(r19431); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); mpfr_init(r19435); mpfr_init(r19436); mpfr_init(r19437); mpfr_init(r19438); mpfr_init(r19439); } double f_dm(double x) { mpfr_set_d(r19341, x, MPFR_RNDN); ; mpfr_set_si(r19343, mpfr_cmp(r19341, r19342) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19346, r19341, r19341, MPFR_RNDN); mpfr_div(r19347, r19345, r19346, MPFR_RNDN); mpfr_add(r19348, r19344, r19347, MPFR_RNDN); ; mpfr_mul(r19350, r19349, r19344, MPFR_RNDN); mpfr_mul(r19351, r19341, r19341, MPFR_RNDN); mpfr_mul(r19351, r19351, r19341, MPFR_RNDN); mpfr_sqr(r19352, r19351, MPFR_RNDN); mpfr_div(r19353, r19350, r19352, MPFR_RNDN); mpfr_add(r19354, r19348, r19353, MPFR_RNDN); ; mpfr_mul(r19356, r19355, r19344, MPFR_RNDN); mpfr_div(r19357, r19356, r19351, MPFR_RNDN); mpfr_div(r19358, r19357, r19341, MPFR_RNDN); ; ; mpfr_pow(r19361, r19341, r19360, MPFR_RNDN); mpfr_div(r19362, r19359, r19361, MPFR_RNDN); mpfr_add(r19363, r19358, r19362, MPFR_RNDN); ; mpfr_mul(r19365, r19344, r19364, MPFR_RNDN); mpfr_div(r19366, r19365, r19352, MPFR_RNDN); mpfr_div(r19367, r19366, r19346, MPFR_RNDN); mpfr_add(r19368, r19363, r19367, MPFR_RNDN); mpfr_add(r19369, r19354, r19368, MPFR_RNDN); mpfr_div(r19370, r19344, r19341, MPFR_RNDN); mpfr_div(r19371, r19370, r19351, MPFR_RNDN); ; mpfr_mul(r19373, r19371, r19372, MPFR_RNDN); mpfr_add(r19374, r19344, r19373, MPFR_RNDN); ; mpfr_mul(r19376, r19375, r19344, MPFR_RNDN); mpfr_div(r19377, r19376, r19352, MPFR_RNDN); mpfr_add(r19378, r19374, r19377, MPFR_RNDN); ; mpfr_mul(r19380, r19344, r19379, MPFR_RNDN); mpfr_div(r19381, r19380, r19352, MPFR_RNDN); mpfr_div(r19382, r19381, r19346, MPFR_RNDN); mpfr_div(r19383, r19344, r19351, MPFR_RNDN); mpfr_mul(r19384, r19383, r19383, MPFR_RNDN); ; mpfr_pow(r19386, r19341, r19385, MPFR_RNDN); ; mpfr_div(r19388, r19386, r19387, MPFR_RNDN); mpfr_div(r19389, r19384, r19388, MPFR_RNDN); mpfr_add(r19390, r19382, r19389, MPFR_RNDN); mpfr_add(r19391, r19378, r19390, MPFR_RNDN); ; mpfr_div(r19393, r19392, r19346, MPFR_RNDN); ; mpfr_pow(r19395, r19341, r19394, MPFR_RNDN); ; mpfr_div(r19397, r19395, r19396, MPFR_RNDN); mpfr_div(r19398, r19384, r19397, MPFR_RNDN); mpfr_add(r19399, r19393, r19398, MPFR_RNDN); mpfr_add(r19400, r19391, r19399, MPFR_RNDN); mpfr_mul(r19401, r19400, r19341, MPFR_RNDN); mpfr_div(r19402, r19369, r19401, MPFR_RNDN); ; mpfr_set_si(r19404, mpfr_cmp(r19341, r19403) <= 0, MPFR_RNDN); mpfr_mul(r19405, r19351, r19351, MPFR_RNDN); mpfr_sqr(r19406, r19341, MPFR_RNDN); mpfr_mul(r19407, r19364, r19406, MPFR_RNDN); mpfr_add(r19408, r19349, r19407, MPFR_RNDN); mpfr_mul(r19409, r19405, r19408, MPFR_RNDN); mpfr_mul(r19410, r19341, r19355, MPFR_RNDN); mpfr_mul(r19411, r19410, r19351, MPFR_RNDN); mpfr_add(r19412, r19344, r19411, MPFR_RNDN); mpfr_mul(r19413, r19345, r19406, MPFR_RNDN); mpfr_add(r19414, r19412, r19413, MPFR_RNDN); mpfr_add(r19415, r19409, r19414, MPFR_RNDN); mpfr_sqr(r19416, r19406, MPFR_RNDN); mpfr_mul(r19417, r19405, r19416, MPFR_RNDN); mpfr_mul(r19418, r19417, r19359, MPFR_RNDN); mpfr_add(r19419, r19415, r19418, MPFR_RNDN); mpfr_mul(r19420, r19395, r19396, MPFR_RNDN); mpfr_mul(r19421, r19405, r19420, MPFR_RNDN); mpfr_add(r19422, r19344, r19421, MPFR_RNDN); mpfr_mul(r19423, r19379, r19341, MPFR_RNDN); mpfr_mul(r19424, r19341, r19423, MPFR_RNDN); mpfr_add(r19425, r19424, r19375, MPFR_RNDN); mpfr_mul(r19426, r19405, r19425, MPFR_RNDN); mpfr_add(r19427, r19422, r19426, MPFR_RNDN); mpfr_mul(r19428, r19386, r19387, MPFR_RNDN); mpfr_mul(r19429, r19428, r19405, MPFR_RNDN); mpfr_mul(r19430, r19372, r19416, MPFR_RNDN); mpfr_mul(r19431, r19341, r19392, MPFR_RNDN); mpfr_mul(r19432, r19431, r19341, MPFR_RNDN); mpfr_add(r19433, r19430, r19432, MPFR_RNDN); mpfr_add(r19434, r19429, r19433, MPFR_RNDN); mpfr_add(r19435, r19427, r19434, MPFR_RNDN); mpfr_div(r19436, r19341, r19435, MPFR_RNDN); mpfr_mul(r19437, r19419, r19436, MPFR_RNDN); if (mpfr_get_si(r19404, MPFR_RNDN)) { mpfr_set(r19438, r19437, MPFR_RNDN); } else { mpfr_set(r19438, r19402, MPFR_RNDN); }; if (mpfr_get_si(r19343, MPFR_RNDN)) { mpfr_set(r19439, r19402, MPFR_RNDN); } else { mpfr_set(r19439, r19438, MPFR_RNDN); }; return mpfr_get_d(r19439, MPFR_RNDN); }