#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19127 = 1.0f; float r19128 = 0.10499349236488342f; float r19129 = x; float r19130 = r19129 * r19129; float r19131 = r19128 * r19130; float r19132 = r19127 + r19131; float r19133 = 0.04240605980157852f; float r19134 = r19130 * r19130; float r19135 = r19133 * r19134; float r19136 = r19132 + r19135; float r19137 = 0.007264418061822653f; float r19138 = r19134 * r19130; float r19139 = r19137 * r19138; float r19140 = r19136 + r19139; float r19141 = 0.0005064033903181553f; float r19142 = r19138 * r19130; float r19143 = r19141 * r19142; float r19144 = r19140 + r19143; float r19145 = 0.00017899709928315133f; float r19146 = r19142 * r19130; float r19147 = r19145 * r19146; float r19148 = r19144 + r19147; float r19149 = 0.7715470790863037f; float r19150 = r19149 * r19130; float r19151 = r19127 + r19150; float r19152 = 0.29097387194633484f; float r19153 = r19152 * r19134; float r19154 = r19151 + r19153; float r19155 = 0.0694555789232254f; float r19156 = r19155 * r19138; float r19157 = r19154 + r19156; float r19158 = 0.01400054432451725f; float r19159 = r19158 * r19142; float r19160 = r19157 + r19159; float r19161 = 0.0008327945251949131f; float r19162 = r19161 * r19146; float r19163 = r19160 + r19162; float r19164 = 2.0f; float r19165 = r19164 * r19145; float r19166 = r19146 * r19130; float r19167 = r19165 * r19166; float r19168 = r19163 + r19167; float r19169 = r19148 / r19168; float r19170 = r19169 * r19129; return r19170; } double f_id(double x) { double r19171 = 1.0; double r19172 = 0.10499349236488342; double r19173 = x; double r19174 = r19173 * r19173; double r19175 = r19172 * r19174; double r19176 = r19171 + r19175; double r19177 = 0.04240605980157852; double r19178 = r19174 * r19174; double r19179 = r19177 * r19178; double r19180 = r19176 + r19179; double r19181 = 0.007264418061822653; double r19182 = r19178 * r19174; double r19183 = r19181 * r19182; double r19184 = r19180 + r19183; double r19185 = 0.0005064033903181553; double r19186 = r19182 * r19174; double r19187 = r19185 * r19186; double r19188 = r19184 + r19187; double r19189 = 0.00017899709928315133; double r19190 = r19186 * r19174; double r19191 = r19189 * r19190; double r19192 = r19188 + r19191; double r19193 = 0.7715470790863037; double r19194 = r19193 * r19174; double r19195 = r19171 + r19194; double r19196 = 0.29097387194633484; double r19197 = r19196 * r19178; double r19198 = r19195 + r19197; double r19199 = 0.0694555789232254; double r19200 = r19199 * r19182; double r19201 = r19198 + r19200; double r19202 = 0.01400054432451725; double r19203 = r19202 * r19186; double r19204 = r19201 + r19203; double r19205 = 0.0008327945251949131; double r19206 = r19205 * r19190; double r19207 = r19204 + r19206; double r19208 = 2.0; double r19209 = r19208 * r19189; double r19210 = r19190 * r19174; double r19211 = r19209 * r19210; double r19212 = r19207 + r19211; double r19213 = r19192 / r19212; double r19214 = r19213 * r19173; return r19214; } double f_of(float x) { float r19215 = x; float r19216 = r19215 * (r19215 * r19215); float r19217 = r19216 * r19216; float r19218 = 0.0694555789232254f; float r19219 = 0.01400054432451725f; float r19220 = r19219 * r19215; float r19221 = r19215 * r19220; float r19222 = r19218 + r19221; float r19223 = r19217 * r19222; float r19224 = r19215 * r19215; float r19225 = r19216 * r19224; float r19226 = r19225 * r19225; float r19227 = 0.00017899709928315133f; float r19228 = r19227 * r19215; float r19229 = 2.0f; float r19230 = r19215 * r19229; float r19231 = r19228 * r19230; float r19232 = 0.0008327945251949131f; float r19233 = r19231 + r19232; float r19234 = r19226 * r19233; float r19235 = r19223 + r19234; float r19236 = 0.7715470790863037f; float r19237 = r19215 * r19236; float r19238 = r19237 * r19215; float r19239 = 0.29097387194633484f; float r19240 = r19215 * r19239; float r19241 = r19216 * r19240; float r19242 = r19238 + r19241; float r19243 = 1.0f; float r19244 = r19242 + r19243; float r19245 = r19235 + r19244; float r19246 = r19215 / r19245; float r19247 = 0.10499349236488342f; float r19248 = r19215 * r19247; float r19249 = r19215 * r19248; float r19250 = 0.04240605980157852f; float r19251 = r19250 * r19215; float r19252 = r19251 * r19216; float r19253 = r19252 + r19243; float r19254 = r19249 + r19253; float r19255 = r19224 * r19224; float r19256 = r19255 * r19255; float r19257 = r19224 * r19227; float r19258 = 0.0005064033903181553f; float r19259 = r19257 + r19258; float r19260 = r19256 * r19259; float r19261 = 0.007264418061822653f; float r19262 = r19261 * r19217; float r19263 = r19260 + r19262; float r19264 = r19254 + r19263; float r19265 = r19246 * r19264; return r19265; } double f_od(double x) { double r19266 = x; double r19267 = r19266 * (r19266 * r19266); double r19268 = r19267 * r19267; double r19269 = 0.0694555789232254; double r19270 = 0.01400054432451725; double r19271 = r19270 * r19266; double r19272 = r19266 * r19271; double r19273 = r19269 + r19272; double r19274 = r19268 * r19273; double r19275 = r19266 * r19266; double r19276 = r19267 * r19275; double r19277 = r19276 * r19276; double r19278 = 0.00017899709928315133; double r19279 = r19278 * r19266; double r19280 = 2.0; double r19281 = r19266 * r19280; double r19282 = r19279 * r19281; double r19283 = 0.0008327945251949131; double r19284 = r19282 + r19283; double r19285 = r19277 * r19284; double r19286 = r19274 + r19285; double r19287 = 0.7715470790863037; double r19288 = r19266 * r19287; double r19289 = r19288 * r19266; double r19290 = 0.29097387194633484; double r19291 = r19266 * r19290; double r19292 = r19267 * r19291; double r19293 = r19289 + r19292; double r19294 = 1.0; double r19295 = r19293 + r19294; double r19296 = r19286 + r19295; double r19297 = r19266 / r19296; double r19298 = 0.10499349236488342; double r19299 = r19266 * r19298; double r19300 = r19266 * r19299; double r19301 = 0.04240605980157852; double r19302 = r19301 * r19266; double r19303 = r19302 * r19267; double r19304 = r19303 + r19294; double r19305 = r19300 + r19304; double r19306 = r19275 * r19275; double r19307 = r19306 * r19306; double r19308 = r19275 * r19278; double r19309 = 0.0005064033903181553; double r19310 = r19308 + r19309; double r19311 = r19307 * r19310; double r19312 = 0.007264418061822653; double r19313 = r19312 * r19268; double r19314 = r19311 + r19313; double r19315 = r19305 + r19314; double r19316 = r19297 * r19315; return r19316; } 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 r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328, r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336, r19337, r19338, r19339, r19340, r19341, r19342, r19343, r19344, r19345, r19346, r19347, r19348, r19349, r19350, r19351, r19352, r19353, r19354, r19355, r19356, r19357, r19358, r19359, r19360; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19317, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19318, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19319); mpfr_init(r19320); mpfr_init(r19321); mpfr_init(r19322); mpfr_init_set_str(r19323, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init_set_str(r19327, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19328); mpfr_init(r19329); mpfr_init(r19330); mpfr_init_set_str(r19331, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19332); mpfr_init(r19333); mpfr_init(r19334); mpfr_init_set_str(r19335, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19336); mpfr_init(r19337); mpfr_init(r19338); mpfr_init_set_str(r19339, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19340); mpfr_init(r19341); mpfr_init_set_str(r19342, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19343); mpfr_init(r19344); mpfr_init_set_str(r19345, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19346); mpfr_init(r19347); mpfr_init_set_str(r19348, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19349); mpfr_init(r19350); mpfr_init_set_str(r19351, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19352); mpfr_init(r19353); mpfr_init_set_str(r19354, "2", 10, MPFR_RNDN); mpfr_init(r19355); mpfr_init(r19356); mpfr_init(r19357); mpfr_init(r19358); mpfr_init(r19359); mpfr_init(r19360); } double f_im(double x) { ; ; mpfr_set_d(r19319, x, MPFR_RNDN); mpfr_mul(r19320, r19319, r19319, MPFR_RNDN); mpfr_mul(r19321, r19318, r19320, MPFR_RNDN); mpfr_add(r19322, r19317, r19321, MPFR_RNDN); ; mpfr_mul(r19324, r19320, r19320, MPFR_RNDN); mpfr_mul(r19325, r19323, r19324, MPFR_RNDN); mpfr_add(r19326, r19322, r19325, MPFR_RNDN); ; mpfr_mul(r19328, r19324, r19320, MPFR_RNDN); mpfr_mul(r19329, r19327, r19328, MPFR_RNDN); mpfr_add(r19330, r19326, r19329, MPFR_RNDN); ; mpfr_mul(r19332, r19328, r19320, MPFR_RNDN); mpfr_mul(r19333, r19331, r19332, MPFR_RNDN); mpfr_add(r19334, r19330, r19333, MPFR_RNDN); ; mpfr_mul(r19336, r19332, r19320, MPFR_RNDN); mpfr_mul(r19337, r19335, r19336, MPFR_RNDN); mpfr_add(r19338, r19334, r19337, MPFR_RNDN); ; mpfr_mul(r19340, r19339, r19320, MPFR_RNDN); mpfr_add(r19341, r19317, r19340, MPFR_RNDN); ; mpfr_mul(r19343, r19342, r19324, MPFR_RNDN); mpfr_add(r19344, r19341, r19343, MPFR_RNDN); ; mpfr_mul(r19346, r19345, r19328, MPFR_RNDN); mpfr_add(r19347, r19344, r19346, MPFR_RNDN); ; mpfr_mul(r19349, r19348, r19332, MPFR_RNDN); mpfr_add(r19350, r19347, r19349, MPFR_RNDN); ; mpfr_mul(r19352, r19351, r19336, MPFR_RNDN); mpfr_add(r19353, r19350, r19352, MPFR_RNDN); ; mpfr_mul(r19355, r19354, r19335, MPFR_RNDN); mpfr_mul(r19356, r19336, r19320, MPFR_RNDN); mpfr_mul(r19357, r19355, r19356, MPFR_RNDN); mpfr_add(r19358, r19353, r19357, MPFR_RNDN); mpfr_div(r19359, r19338, r19358, MPFR_RNDN); mpfr_mul(r19360, r19359, r19319, MPFR_RNDN); return mpfr_get_d(r19360, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19361); mpfr_init(r19362); mpfr_init(r19363); mpfr_init_set_str(r19364, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19365, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19366); mpfr_init(r19367); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init(r19372); mpfr_init_set_str(r19373, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19374); mpfr_init_set_str(r19375, "2", 10, MPFR_RNDN); mpfr_init(r19376); mpfr_init(r19377); mpfr_init_set_str(r19378, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19379); mpfr_init(r19380); mpfr_init(r19381); mpfr_init_set_str(r19382, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19383); mpfr_init(r19384); mpfr_init_set_str(r19385, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init_set_str(r19389, "1", 10, MPFR_RNDN); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init_set_str(r19393, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19394); mpfr_init(r19395); mpfr_init_set_str(r19396, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init(r19403); mpfr_init_set_str(r19404, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19405); mpfr_init(r19406); mpfr_init_set_str(r19407, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); } double f_fm(double x) { mpfr_set_d(r19361, x, MPFR_RNDN); mpfr_mul(r19362, r19361, r19361, MPFR_RNDN); mpfr_mul(r19362, r19362, r19361, MPFR_RNDN); mpfr_mul(r19363, r19362, r19362, MPFR_RNDN); ; ; mpfr_mul(r19366, r19365, r19361, MPFR_RNDN); mpfr_mul(r19367, r19361, r19366, MPFR_RNDN); mpfr_add(r19368, r19364, r19367, MPFR_RNDN); mpfr_mul(r19369, r19363, r19368, MPFR_RNDN); mpfr_sqr(r19370, r19361, MPFR_RNDN); mpfr_mul(r19371, r19362, r19370, MPFR_RNDN); mpfr_mul(r19372, r19371, r19371, MPFR_RNDN); ; mpfr_mul(r19374, r19373, r19361, MPFR_RNDN); ; mpfr_mul(r19376, r19361, r19375, MPFR_RNDN); mpfr_mul(r19377, r19374, r19376, MPFR_RNDN); ; mpfr_add(r19379, r19377, r19378, MPFR_RNDN); mpfr_mul(r19380, r19372, r19379, MPFR_RNDN); mpfr_add(r19381, r19369, r19380, MPFR_RNDN); ; mpfr_mul(r19383, r19361, r19382, MPFR_RNDN); mpfr_mul(r19384, r19383, r19361, MPFR_RNDN); ; mpfr_mul(r19386, r19361, r19385, MPFR_RNDN); mpfr_mul(r19387, r19362, r19386, MPFR_RNDN); mpfr_add(r19388, r19384, r19387, MPFR_RNDN); ; mpfr_add(r19390, r19388, r19389, MPFR_RNDN); mpfr_add(r19391, r19381, r19390, MPFR_RNDN); mpfr_div(r19392, r19361, r19391, MPFR_RNDN); ; mpfr_mul(r19394, r19361, r19393, MPFR_RNDN); mpfr_mul(r19395, r19361, r19394, MPFR_RNDN); ; mpfr_mul(r19397, r19396, r19361, MPFR_RNDN); mpfr_mul(r19398, r19397, r19362, MPFR_RNDN); mpfr_add(r19399, r19398, r19389, MPFR_RNDN); mpfr_add(r19400, r19395, r19399, MPFR_RNDN); mpfr_mul(r19401, r19370, r19370, MPFR_RNDN); mpfr_sqr(r19402, r19401, MPFR_RNDN); mpfr_mul(r19403, r19370, r19373, MPFR_RNDN); ; mpfr_add(r19405, r19403, r19404, MPFR_RNDN); mpfr_mul(r19406, r19402, r19405, MPFR_RNDN); ; mpfr_mul(r19408, r19407, r19363, MPFR_RNDN); mpfr_add(r19409, r19406, r19408, MPFR_RNDN); mpfr_add(r19410, r19400, r19409, MPFR_RNDN); mpfr_mul(r19411, r19392, r19410, MPFR_RNDN); return mpfr_get_d(r19411, MPFR_RNDN); } static mpfr_t 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, r19440, r19441, r19442, r19443, r19444, r19445, r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459, r19460, r19461, r19462; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init_set_str(r19415, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19416, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); mpfr_init(r19422); mpfr_init(r19423); mpfr_init_set_str(r19424, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19425); mpfr_init_set_str(r19426, "2", 10, MPFR_RNDN); mpfr_init(r19427); mpfr_init(r19428); mpfr_init_set_str(r19429, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19430); mpfr_init(r19431); mpfr_init(r19432); mpfr_init_set_str(r19433, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19434); mpfr_init(r19435); mpfr_init_set_str(r19436, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19437); mpfr_init(r19438); mpfr_init(r19439); mpfr_init_set_str(r19440, "1", 10, MPFR_RNDN); mpfr_init(r19441); mpfr_init(r19442); mpfr_init(r19443); mpfr_init_set_str(r19444, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19445); mpfr_init(r19446); mpfr_init_set_str(r19447, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19448); mpfr_init(r19449); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init(r19454); mpfr_init_set_str(r19455, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19456); mpfr_init(r19457); mpfr_init_set_str(r19458, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19459); mpfr_init(r19460); mpfr_init(r19461); mpfr_init(r19462); } double f_dm(double x) { mpfr_set_d(r19412, x, MPFR_RNDN); mpfr_mul(r19413, r19412, r19412, MPFR_RNDN); mpfr_mul(r19413, r19413, r19412, MPFR_RNDN); mpfr_mul(r19414, r19413, r19413, MPFR_RNDN); ; ; mpfr_mul(r19417, r19416, r19412, MPFR_RNDN); mpfr_mul(r19418, r19412, r19417, MPFR_RNDN); mpfr_add(r19419, r19415, r19418, MPFR_RNDN); mpfr_mul(r19420, r19414, r19419, MPFR_RNDN); mpfr_sqr(r19421, r19412, MPFR_RNDN); mpfr_mul(r19422, r19413, r19421, MPFR_RNDN); mpfr_mul(r19423, r19422, r19422, MPFR_RNDN); ; mpfr_mul(r19425, r19424, r19412, MPFR_RNDN); ; mpfr_mul(r19427, r19412, r19426, MPFR_RNDN); mpfr_mul(r19428, r19425, r19427, MPFR_RNDN); ; mpfr_add(r19430, r19428, r19429, MPFR_RNDN); mpfr_mul(r19431, r19423, r19430, MPFR_RNDN); mpfr_add(r19432, r19420, r19431, MPFR_RNDN); ; mpfr_mul(r19434, r19412, r19433, MPFR_RNDN); mpfr_mul(r19435, r19434, r19412, MPFR_RNDN); ; mpfr_mul(r19437, r19412, r19436, MPFR_RNDN); mpfr_mul(r19438, r19413, r19437, MPFR_RNDN); mpfr_add(r19439, r19435, r19438, MPFR_RNDN); ; mpfr_add(r19441, r19439, r19440, MPFR_RNDN); mpfr_add(r19442, r19432, r19441, MPFR_RNDN); mpfr_div(r19443, r19412, r19442, MPFR_RNDN); ; mpfr_mul(r19445, r19412, r19444, MPFR_RNDN); mpfr_mul(r19446, r19412, r19445, MPFR_RNDN); ; mpfr_mul(r19448, r19447, r19412, MPFR_RNDN); mpfr_mul(r19449, r19448, r19413, MPFR_RNDN); mpfr_add(r19450, r19449, r19440, MPFR_RNDN); mpfr_add(r19451, r19446, r19450, MPFR_RNDN); mpfr_mul(r19452, r19421, r19421, MPFR_RNDN); mpfr_sqr(r19453, r19452, MPFR_RNDN); mpfr_mul(r19454, r19421, r19424, MPFR_RNDN); ; mpfr_add(r19456, r19454, r19455, MPFR_RNDN); mpfr_mul(r19457, r19453, r19456, MPFR_RNDN); ; mpfr_mul(r19459, r19458, r19414, MPFR_RNDN); mpfr_add(r19460, r19457, r19459, MPFR_RNDN); mpfr_add(r19461, r19451, r19460, MPFR_RNDN); mpfr_mul(r19462, r19443, r19461, MPFR_RNDN); return mpfr_get_d(r19462, MPFR_RNDN); }