#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19183 = 1.0f; float r19184 = 0.10499349236488342f; float r19185 = x; float r19186 = r19185 * r19185; float r19187 = r19184 * r19186; float r19188 = r19183 + r19187; float r19189 = 0.04240605980157852f; float r19190 = r19186 * r19186; float r19191 = r19189 * r19190; float r19192 = r19188 + r19191; float r19193 = 0.007264418061822653f; float r19194 = r19190 * r19186; float r19195 = r19193 * r19194; float r19196 = r19192 + r19195; float r19197 = 0.0005064033903181553f; float r19198 = r19194 * r19186; float r19199 = r19197 * r19198; float r19200 = r19196 + r19199; float r19201 = 0.00017899709928315133f; float r19202 = r19198 * r19186; float r19203 = r19201 * r19202; float r19204 = r19200 + r19203; float r19205 = 0.7715470790863037f; float r19206 = r19205 * r19186; float r19207 = r19183 + r19206; float r19208 = 0.29097387194633484f; float r19209 = r19208 * r19190; float r19210 = r19207 + r19209; float r19211 = 0.0694555789232254f; float r19212 = r19211 * r19194; float r19213 = r19210 + r19212; float r19214 = 0.01400054432451725f; float r19215 = r19214 * r19198; float r19216 = r19213 + r19215; float r19217 = 0.0008327945251949131f; float r19218 = r19217 * r19202; float r19219 = r19216 + r19218; float r19220 = 2.0f; float r19221 = r19220 * r19201; float r19222 = r19202 * r19186; float r19223 = r19221 * r19222; float r19224 = r19219 + r19223; float r19225 = r19204 / r19224; float r19226 = r19225 * r19185; return r19226; } double f_id(double x) { double r19227 = 1.0; double r19228 = 0.10499349236488342; double r19229 = x; double r19230 = r19229 * r19229; double r19231 = r19228 * r19230; double r19232 = r19227 + r19231; double r19233 = 0.04240605980157852; double r19234 = r19230 * r19230; double r19235 = r19233 * r19234; double r19236 = r19232 + r19235; double r19237 = 0.007264418061822653; double r19238 = r19234 * r19230; double r19239 = r19237 * r19238; double r19240 = r19236 + r19239; double r19241 = 0.0005064033903181553; double r19242 = r19238 * r19230; double r19243 = r19241 * r19242; double r19244 = r19240 + r19243; double r19245 = 0.00017899709928315133; double r19246 = r19242 * r19230; double r19247 = r19245 * r19246; double r19248 = r19244 + r19247; double r19249 = 0.7715470790863037; double r19250 = r19249 * r19230; double r19251 = r19227 + r19250; double r19252 = 0.29097387194633484; double r19253 = r19252 * r19234; double r19254 = r19251 + r19253; double r19255 = 0.0694555789232254; double r19256 = r19255 * r19238; double r19257 = r19254 + r19256; double r19258 = 0.01400054432451725; double r19259 = r19258 * r19242; double r19260 = r19257 + r19259; double r19261 = 0.0008327945251949131; double r19262 = r19261 * r19246; double r19263 = r19260 + r19262; double r19264 = 2.0; double r19265 = r19264 * r19245; double r19266 = r19246 * r19230; double r19267 = r19265 * r19266; double r19268 = r19263 + r19267; double r19269 = r19248 / r19268; double r19270 = r19269 * r19229; return r19270; } double f_of(float x) { float r19271 = x; float r19272 = 0.00017899709928315133f; float r19273 = r19272 * r19271; float r19274 = r19273 * r19271; float r19275 = r19271 * r19271; float r19276 = r19275 * r19275; float r19277 = r19276 * r19276; float r19278 = r19274 * r19277; float r19279 = 1.0f; float r19280 = 0.04240605980157852f; float r19281 = r19271 * r19280; float r19282 = r19271 * (r19271 * r19271); float r19283 = r19281 * r19282; float r19284 = 0.10499349236488342f; float r19285 = r19284 * r19271; float r19286 = r19285 * r19271; float r19287 = r19283 + r19286; float r19288 = r19279 + r19287; float r19289 = r19278 + r19288; float r19290 = 0.0005064033903181553f; float r19291 = r19275 * r19290; float r19292 = 0.007264418061822653f; float r19293 = r19291 + r19292; float r19294 = r19282 * r19282; float r19295 = r19293 * r19294; float r19296 = r19289 + r19295; float r19297 = 0.7715470790863037f; float r19298 = r19297 * r19275; float r19299 = 0.29097387194633484f; float r19300 = r19299 * r19271; float r19301 = r19300 * r19282; float r19302 = r19298 + r19301; float r19303 = r19279 + r19302; float r19304 = 2.0f; float r19305 = r19304 * r19272; float r19306 = r19305 * r19275; float r19307 = 0.0008327945251949131f; float r19308 = r19306 + r19307; float r19309 = r19276 * r19294; float r19310 = r19308 * r19309; float r19311 = r19303 + r19310; float r19312 = 0.0694555789232254f; float r19313 = 0.01400054432451725f; float r19314 = r19275 * r19313; float r19315 = r19312 + r19314; float r19316 = r19294 * r19315; float r19317 = r19311 + r19316; float r19318 = r19296 / r19317; float r19319 = r19271 * r19318; return r19319; } double f_od(double x) { double r19320 = x; double r19321 = 0.00017899709928315133; double r19322 = r19321 * r19320; double r19323 = r19322 * r19320; double r19324 = r19320 * r19320; double r19325 = r19324 * r19324; double r19326 = r19325 * r19325; double r19327 = r19323 * r19326; double r19328 = 1.0; double r19329 = 0.04240605980157852; double r19330 = r19320 * r19329; double r19331 = r19320 * (r19320 * r19320); double r19332 = r19330 * r19331; double r19333 = 0.10499349236488342; double r19334 = r19333 * r19320; double r19335 = r19334 * r19320; double r19336 = r19332 + r19335; double r19337 = r19328 + r19336; double r19338 = r19327 + r19337; double r19339 = 0.0005064033903181553; double r19340 = r19324 * r19339; double r19341 = 0.007264418061822653; double r19342 = r19340 + r19341; double r19343 = r19331 * r19331; double r19344 = r19342 * r19343; double r19345 = r19338 + r19344; double r19346 = 0.7715470790863037; double r19347 = r19346 * r19324; double r19348 = 0.29097387194633484; double r19349 = r19348 * r19320; double r19350 = r19349 * r19331; double r19351 = r19347 + r19350; double r19352 = r19328 + r19351; double r19353 = 2.0; double r19354 = r19353 * r19321; double r19355 = r19354 * r19324; double r19356 = 0.0008327945251949131; double r19357 = r19355 + r19356; double r19358 = r19325 * r19343; double r19359 = r19357 * r19358; double r19360 = r19352 + r19359; double r19361 = 0.0694555789232254; double r19362 = 0.01400054432451725; double r19363 = r19324 * r19362; double r19364 = r19361 + r19363; double r19365 = r19343 * r19364; double r19366 = r19360 + r19365; double r19367 = r19345 / r19366; double r19368 = r19320 * r19367; return r19368; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19369, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19370, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19371); mpfr_init(r19372); mpfr_init(r19373); mpfr_init(r19374); mpfr_init_set_str(r19375, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19376); mpfr_init(r19377); mpfr_init(r19378); mpfr_init_set_str(r19379, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init_set_str(r19383, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19384); mpfr_init(r19385); mpfr_init(r19386); mpfr_init_set_str(r19387, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init_set_str(r19391, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19392); mpfr_init(r19393); mpfr_init_set_str(r19394, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19395); mpfr_init(r19396); mpfr_init_set_str(r19397, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19398); mpfr_init(r19399); mpfr_init_set_str(r19400, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19401); mpfr_init(r19402); mpfr_init_set_str(r19403, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19404); mpfr_init(r19405); mpfr_init_set_str(r19406, "2", 10, MPFR_RNDN); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); } double f_im(double x) { ; ; mpfr_set_d(r19371, x, MPFR_RNDN); mpfr_mul(r19372, r19371, r19371, MPFR_RNDN); mpfr_mul(r19373, r19370, r19372, MPFR_RNDN); mpfr_add(r19374, r19369, r19373, MPFR_RNDN); ; mpfr_mul(r19376, r19372, r19372, MPFR_RNDN); mpfr_mul(r19377, r19375, r19376, MPFR_RNDN); mpfr_add(r19378, r19374, r19377, MPFR_RNDN); ; mpfr_mul(r19380, r19376, r19372, MPFR_RNDN); mpfr_mul(r19381, r19379, r19380, MPFR_RNDN); mpfr_add(r19382, r19378, r19381, MPFR_RNDN); ; mpfr_mul(r19384, r19380, r19372, MPFR_RNDN); mpfr_mul(r19385, r19383, r19384, MPFR_RNDN); mpfr_add(r19386, r19382, r19385, MPFR_RNDN); ; mpfr_mul(r19388, r19384, r19372, MPFR_RNDN); mpfr_mul(r19389, r19387, r19388, MPFR_RNDN); mpfr_add(r19390, r19386, r19389, MPFR_RNDN); ; mpfr_mul(r19392, r19391, r19372, MPFR_RNDN); mpfr_add(r19393, r19369, r19392, MPFR_RNDN); ; mpfr_mul(r19395, r19394, r19376, MPFR_RNDN); mpfr_add(r19396, r19393, r19395, MPFR_RNDN); ; mpfr_mul(r19398, r19397, r19380, MPFR_RNDN); mpfr_add(r19399, r19396, r19398, MPFR_RNDN); ; mpfr_mul(r19401, r19400, r19384, MPFR_RNDN); mpfr_add(r19402, r19399, r19401, MPFR_RNDN); ; mpfr_mul(r19404, r19403, r19388, MPFR_RNDN); mpfr_add(r19405, r19402, r19404, MPFR_RNDN); ; mpfr_mul(r19407, r19406, r19387, MPFR_RNDN); mpfr_mul(r19408, r19388, r19372, MPFR_RNDN); mpfr_mul(r19409, r19407, r19408, MPFR_RNDN); mpfr_add(r19410, r19405, r19409, MPFR_RNDN); mpfr_div(r19411, r19390, r19410, MPFR_RNDN); mpfr_mul(r19412, r19411, r19371, MPFR_RNDN); return mpfr_get_d(r19412, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19413); mpfr_init_set_str(r19414, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19415); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init_set_str(r19421, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19422, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19423); mpfr_init(r19424); mpfr_init(r19425); mpfr_init_set_str(r19426, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19427); mpfr_init(r19428); mpfr_init(r19429); mpfr_init(r19430); mpfr_init(r19431); mpfr_init_set_str(r19432, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19433); mpfr_init_set_str(r19434, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19435); mpfr_init(r19436); mpfr_init(r19437); mpfr_init(r19438); mpfr_init_set_str(r19439, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19440); mpfr_init_set_str(r19441, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19442); mpfr_init(r19443); mpfr_init(r19444); mpfr_init(r19445); mpfr_init_set_str(r19446, "2", 10, MPFR_RNDN); mpfr_init(r19447); mpfr_init(r19448); mpfr_init_set_str(r19449, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init_set_str(r19454, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19455, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19456); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init(r19460); mpfr_init(r19461); } double f_fm(double x) { mpfr_set_d(r19413, x, MPFR_RNDN); ; mpfr_mul(r19415, r19414, r19413, MPFR_RNDN); mpfr_mul(r19416, r19415, r19413, MPFR_RNDN); mpfr_mul(r19417, r19413, r19413, MPFR_RNDN); mpfr_mul(r19418, r19417, r19417, MPFR_RNDN); mpfr_mul(r19419, r19418, r19418, MPFR_RNDN); mpfr_mul(r19420, r19416, r19419, MPFR_RNDN); ; ; mpfr_mul(r19423, r19413, r19422, MPFR_RNDN); mpfr_mul(r19424, r19413, r19413, MPFR_RNDN); mpfr_mul(r19424, r19424, r19413, MPFR_RNDN); mpfr_mul(r19425, r19423, r19424, MPFR_RNDN); ; mpfr_mul(r19427, r19426, r19413, MPFR_RNDN); mpfr_mul(r19428, r19427, r19413, MPFR_RNDN); mpfr_add(r19429, r19425, r19428, MPFR_RNDN); mpfr_add(r19430, r19421, r19429, MPFR_RNDN); mpfr_add(r19431, r19420, r19430, MPFR_RNDN); ; mpfr_mul(r19433, r19417, r19432, MPFR_RNDN); ; mpfr_add(r19435, r19433, r19434, MPFR_RNDN); mpfr_mul(r19436, r19424, r19424, MPFR_RNDN); mpfr_mul(r19437, r19435, r19436, MPFR_RNDN); mpfr_add(r19438, r19431, r19437, MPFR_RNDN); ; mpfr_mul(r19440, r19439, r19417, MPFR_RNDN); ; mpfr_mul(r19442, r19441, r19413, MPFR_RNDN); mpfr_mul(r19443, r19442, r19424, MPFR_RNDN); mpfr_add(r19444, r19440, r19443, MPFR_RNDN); mpfr_add(r19445, r19421, r19444, MPFR_RNDN); ; mpfr_mul(r19447, r19446, r19414, MPFR_RNDN); mpfr_mul(r19448, r19447, r19417, MPFR_RNDN); ; mpfr_add(r19450, r19448, r19449, MPFR_RNDN); mpfr_mul(r19451, r19418, r19436, MPFR_RNDN); mpfr_mul(r19452, r19450, r19451, MPFR_RNDN); mpfr_add(r19453, r19445, r19452, MPFR_RNDN); ; ; mpfr_mul(r19456, r19417, r19455, MPFR_RNDN); mpfr_add(r19457, r19454, r19456, MPFR_RNDN); mpfr_mul(r19458, r19436, r19457, MPFR_RNDN); mpfr_add(r19459, r19453, r19458, MPFR_RNDN); mpfr_div(r19460, r19438, r19459, MPFR_RNDN); mpfr_mul(r19461, r19413, r19460, MPFR_RNDN); return mpfr_get_d(r19461, MPFR_RNDN); } static mpfr_t r19462, r19463, r19464, r19465, r19466, r19467, r19468, r19469, r19470, r19471, r19472, r19473, r19474, r19475, r19476, r19477, r19478, r19479, r19480, r19481, r19482, r19483, r19484, r19485, r19486, r19487, r19488, r19489, r19490, r19491, r19492, r19493, r19494, r19495, r19496, r19497, r19498, r19499, r19500, r19501, r19502, r19503, r19504, r19505, r19506, r19507, r19508, r19509, r19510; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19462); mpfr_init_set_str(r19463, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19464); mpfr_init(r19465); mpfr_init(r19466); mpfr_init(r19467); mpfr_init(r19468); mpfr_init(r19469); mpfr_init_set_str(r19470, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19471, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19472); mpfr_init(r19473); mpfr_init(r19474); mpfr_init_set_str(r19475, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19476); mpfr_init(r19477); mpfr_init(r19478); mpfr_init(r19479); mpfr_init(r19480); mpfr_init_set_str(r19481, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19482); mpfr_init_set_str(r19483, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19484); mpfr_init(r19485); mpfr_init(r19486); mpfr_init(r19487); mpfr_init_set_str(r19488, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19489); mpfr_init_set_str(r19490, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19491); mpfr_init(r19492); mpfr_init(r19493); mpfr_init(r19494); mpfr_init_set_str(r19495, "2", 10, MPFR_RNDN); mpfr_init(r19496); mpfr_init(r19497); mpfr_init_set_str(r19498, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19499); mpfr_init(r19500); mpfr_init(r19501); mpfr_init(r19502); mpfr_init_set_str(r19503, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19504, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19505); mpfr_init(r19506); mpfr_init(r19507); mpfr_init(r19508); mpfr_init(r19509); mpfr_init(r19510); } double f_dm(double x) { mpfr_set_d(r19462, x, MPFR_RNDN); ; mpfr_mul(r19464, r19463, r19462, MPFR_RNDN); mpfr_mul(r19465, r19464, r19462, MPFR_RNDN); mpfr_mul(r19466, r19462, r19462, MPFR_RNDN); mpfr_mul(r19467, r19466, r19466, MPFR_RNDN); mpfr_mul(r19468, r19467, r19467, MPFR_RNDN); mpfr_mul(r19469, r19465, r19468, MPFR_RNDN); ; ; mpfr_mul(r19472, r19462, r19471, MPFR_RNDN); mpfr_mul(r19473, r19462, r19462, MPFR_RNDN); mpfr_mul(r19473, r19473, r19462, MPFR_RNDN); mpfr_mul(r19474, r19472, r19473, MPFR_RNDN); ; mpfr_mul(r19476, r19475, r19462, MPFR_RNDN); mpfr_mul(r19477, r19476, r19462, MPFR_RNDN); mpfr_add(r19478, r19474, r19477, MPFR_RNDN); mpfr_add(r19479, r19470, r19478, MPFR_RNDN); mpfr_add(r19480, r19469, r19479, MPFR_RNDN); ; mpfr_mul(r19482, r19466, r19481, MPFR_RNDN); ; mpfr_add(r19484, r19482, r19483, MPFR_RNDN); mpfr_mul(r19485, r19473, r19473, MPFR_RNDN); mpfr_mul(r19486, r19484, r19485, MPFR_RNDN); mpfr_add(r19487, r19480, r19486, MPFR_RNDN); ; mpfr_mul(r19489, r19488, r19466, MPFR_RNDN); ; mpfr_mul(r19491, r19490, r19462, MPFR_RNDN); mpfr_mul(r19492, r19491, r19473, MPFR_RNDN); mpfr_add(r19493, r19489, r19492, MPFR_RNDN); mpfr_add(r19494, r19470, r19493, MPFR_RNDN); ; mpfr_mul(r19496, r19495, r19463, MPFR_RNDN); mpfr_mul(r19497, r19496, r19466, MPFR_RNDN); ; mpfr_add(r19499, r19497, r19498, MPFR_RNDN); mpfr_mul(r19500, r19467, r19485, MPFR_RNDN); mpfr_mul(r19501, r19499, r19500, MPFR_RNDN); mpfr_add(r19502, r19494, r19501, MPFR_RNDN); ; ; mpfr_mul(r19505, r19466, r19504, MPFR_RNDN); mpfr_add(r19506, r19503, r19505, MPFR_RNDN); mpfr_mul(r19507, r19485, r19506, MPFR_RNDN); mpfr_add(r19508, r19502, r19507, MPFR_RNDN); mpfr_div(r19509, r19487, r19508, MPFR_RNDN); mpfr_mul(r19510, r19462, r19509, MPFR_RNDN); return mpfr_get_d(r19510, MPFR_RNDN); }