#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19311 = 1.0f; float r19312 = 0.1049934947f; float r19313 = x; float r19314 = r19313 * r19313; float r19315 = r19312 * r19314; float r19316 = r19311 + r19315; float r19317 = 0.0424060604f; float r19318 = r19314 * r19314; float r19319 = r19317 * r19318; float r19320 = r19316 + r19319; float r19321 = 0.0072644182f; float r19322 = r19318 * r19314; float r19323 = r19321 * r19322; float r19324 = r19320 + r19323; float r19325 = 0.0005064034f; float r19326 = r19322 * r19314; float r19327 = r19325 * r19326; float r19328 = r19324 + r19327; float r19329 = 0.0001789971f; float r19330 = r19326 * r19314; float r19331 = r19329 * r19330; float r19332 = r19328 + r19331; float r19333 = 0.7715471019f; float r19334 = r19333 * r19314; float r19335 = r19311 + r19334; float r19336 = 0.2909738639f; float r19337 = r19336 * r19318; float r19338 = r19335 + r19337; float r19339 = 0.0694555761f; float r19340 = r19339 * r19322; float r19341 = r19338 + r19340; float r19342 = 0.0140005442f; float r19343 = r19342 * r19326; float r19344 = r19341 + r19343; float r19345 = 0.0008327945f; float r19346 = r19345 * r19330; float r19347 = r19344 + r19346; float r19348 = 2.0f; float r19349 = r19348 * r19329; float r19350 = r19330 * r19314; float r19351 = r19349 * r19350; float r19352 = r19347 + r19351; float r19353 = r19332 / r19352; float r19354 = r19353 * r19313; return r19354; } double f_id(double x) { double r19355 = 1.0; double r19356 = 0.1049934947; double r19357 = x; double r19358 = r19357 * r19357; double r19359 = r19356 * r19358; double r19360 = r19355 + r19359; double r19361 = 0.0424060604; double r19362 = r19358 * r19358; double r19363 = r19361 * r19362; double r19364 = r19360 + r19363; double r19365 = 0.0072644182; double r19366 = r19362 * r19358; double r19367 = r19365 * r19366; double r19368 = r19364 + r19367; double r19369 = 0.0005064034; double r19370 = r19366 * r19358; double r19371 = r19369 * r19370; double r19372 = r19368 + r19371; double r19373 = 0.0001789971; double r19374 = r19370 * r19358; double r19375 = r19373 * r19374; double r19376 = r19372 + r19375; double r19377 = 0.7715471019; double r19378 = r19377 * r19358; double r19379 = r19355 + r19378; double r19380 = 0.2909738639; double r19381 = r19380 * r19362; double r19382 = r19379 + r19381; double r19383 = 0.0694555761; double r19384 = r19383 * r19366; double r19385 = r19382 + r19384; double r19386 = 0.0140005442; double r19387 = r19386 * r19370; double r19388 = r19385 + r19387; double r19389 = 0.0008327945; double r19390 = r19389 * r19374; double r19391 = r19388 + r19390; double r19392 = 2.0; double r19393 = r19392 * r19373; double r19394 = r19374 * r19358; double r19395 = r19393 * r19394; double r19396 = r19391 + r19395; double r19397 = r19376 / r19396; double r19398 = r19397 * r19357; return r19398; } double f_of(float x) { float r19399 = x; float r19400 = r19399 * r19399; float r19401 = r19400 * r19400; float r19402 = r19399 * (r19399 * r19399); float r19403 = r19402 * r19402; float r19404 = r19401 * r19403; float r19405 = 2.0f; float r19406 = 0.0001789971f; float r19407 = r19405 * r19406; float r19408 = r19407 * r19400; float r19409 = 0.0008327945f; float r19410 = r19408 + r19409; float r19411 = r19404 * r19410; float r19412 = 1.0f; float r19413 = 0.2909738639f; float r19414 = r19399 * r19413; float r19415 = r19414 * r19402; float r19416 = r19412 + r19415; float r19417 = 0.7715471019f; float r19418 = r19417 * r19399; float r19419 = r19418 * r19399; float r19420 = r19416 + r19419; float r19421 = r19411 + r19420; float r19422 = 0.0694555761f; float r19423 = 0.0140005442f; float r19424 = r19423 * r19400; float r19425 = r19422 + r19424; float r19426 = r19401 * r19425; float r19427 = r19400 * r19426; float r19428 = r19421 + r19427; float r19429 = r19399 / r19428; float r19430 = 0.0424060604f; float r19431 = r19399 * r19430; float r19432 = r19431 * r19402; float r19433 = 0.1049934947f; float r19434 = r19400 * r19433; float r19435 = r19434 + r19412; float r19436 = r19432 + r19435; float r19437 = r19400 * r19406; float r19438 = r19402 * r19399; float r19439 = r19401 * r19438; float r19440 = r19437 * r19439; float r19441 = r19436 + r19440; float r19442 = 0.0072644182f; float r19443 = 0.0005064034f; float r19444 = r19443 * r19399; float r19445 = r19444 * r19399; float r19446 = r19442 + r19445; float r19447 = r19401 * r19446; float r19448 = r19400 * r19447; float r19449 = r19441 + r19448; float r19450 = r19429 * r19449; return r19450; } double f_od(double x) { double r19451 = x; double r19452 = r19451 * r19451; double r19453 = r19452 * r19452; double r19454 = r19451 * (r19451 * r19451); double r19455 = r19454 * r19454; double r19456 = r19453 * r19455; double r19457 = 2.0; double r19458 = 0.0001789971; double r19459 = r19457 * r19458; double r19460 = r19459 * r19452; double r19461 = 0.0008327945; double r19462 = r19460 + r19461; double r19463 = r19456 * r19462; double r19464 = 1.0; double r19465 = 0.2909738639; double r19466 = r19451 * r19465; double r19467 = r19466 * r19454; double r19468 = r19464 + r19467; double r19469 = 0.7715471019; double r19470 = r19469 * r19451; double r19471 = r19470 * r19451; double r19472 = r19468 + r19471; double r19473 = r19463 + r19472; double r19474 = 0.0694555761; double r19475 = 0.0140005442; double r19476 = r19475 * r19452; double r19477 = r19474 + r19476; double r19478 = r19453 * r19477; double r19479 = r19452 * r19478; double r19480 = r19473 + r19479; double r19481 = r19451 / r19480; double r19482 = 0.0424060604; double r19483 = r19451 * r19482; double r19484 = r19483 * r19454; double r19485 = 0.1049934947; double r19486 = r19452 * r19485; double r19487 = r19486 + r19464; double r19488 = r19484 + r19487; double r19489 = r19452 * r19458; double r19490 = r19454 * r19451; double r19491 = r19453 * r19490; double r19492 = r19489 * r19491; double r19493 = r19488 + r19492; double r19494 = 0.0072644182; double r19495 = 0.0005064034; double r19496 = r19495 * r19451; double r19497 = r19496 * r19451; double r19498 = r19494 + r19497; double r19499 = r19453 * r19498; double r19500 = r19452 * r19499; double r19501 = r19493 + r19500; double r19502 = r19481 * r19501; return r19502; } 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 r19503, r19504, r19505, r19506, r19507, r19508, r19509, r19510, r19511, r19512, r19513, r19514, r19515, r19516, r19517, r19518, r19519, r19520, r19521, r19522, r19523, r19524, r19525, r19526, r19527, r19528, r19529, r19530, r19531, r19532, r19533, r19534, r19535, r19536, r19537, r19538, r19539, r19540, r19541, r19542, r19543, r19544, r19545, r19546; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19503, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19504, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19505); mpfr_init(r19506); mpfr_init(r19507); mpfr_init(r19508); mpfr_init_set_str(r19509, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19510); mpfr_init(r19511); mpfr_init(r19512); mpfr_init_set_str(r19513, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19514); mpfr_init(r19515); mpfr_init(r19516); mpfr_init_set_str(r19517, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19518); mpfr_init(r19519); mpfr_init(r19520); mpfr_init_set_str(r19521, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19522); mpfr_init(r19523); mpfr_init(r19524); mpfr_init_set_str(r19525, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19526); mpfr_init(r19527); mpfr_init_set_str(r19528, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19529); mpfr_init(r19530); mpfr_init_set_str(r19531, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19532); mpfr_init(r19533); mpfr_init_set_str(r19534, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19535); mpfr_init(r19536); mpfr_init_set_str(r19537, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19538); mpfr_init(r19539); mpfr_init_set_str(r19540, "2", 10, MPFR_RNDN); mpfr_init(r19541); mpfr_init(r19542); mpfr_init(r19543); mpfr_init(r19544); mpfr_init(r19545); mpfr_init(r19546); } double f_im(double x) { ; ; mpfr_set_d(r19505, x, MPFR_RNDN); mpfr_mul(r19506, r19505, r19505, MPFR_RNDN); mpfr_mul(r19507, r19504, r19506, MPFR_RNDN); mpfr_add(r19508, r19503, r19507, MPFR_RNDN); ; mpfr_mul(r19510, r19506, r19506, MPFR_RNDN); mpfr_mul(r19511, r19509, r19510, MPFR_RNDN); mpfr_add(r19512, r19508, r19511, MPFR_RNDN); ; mpfr_mul(r19514, r19510, r19506, MPFR_RNDN); mpfr_mul(r19515, r19513, r19514, MPFR_RNDN); mpfr_add(r19516, r19512, r19515, MPFR_RNDN); ; mpfr_mul(r19518, r19514, r19506, MPFR_RNDN); mpfr_mul(r19519, r19517, r19518, MPFR_RNDN); mpfr_add(r19520, r19516, r19519, MPFR_RNDN); ; mpfr_mul(r19522, r19518, r19506, MPFR_RNDN); mpfr_mul(r19523, r19521, r19522, MPFR_RNDN); mpfr_add(r19524, r19520, r19523, MPFR_RNDN); ; mpfr_mul(r19526, r19525, r19506, MPFR_RNDN); mpfr_add(r19527, r19503, r19526, MPFR_RNDN); ; mpfr_mul(r19529, r19528, r19510, MPFR_RNDN); mpfr_add(r19530, r19527, r19529, MPFR_RNDN); ; mpfr_mul(r19532, r19531, r19514, MPFR_RNDN); mpfr_add(r19533, r19530, r19532, MPFR_RNDN); ; mpfr_mul(r19535, r19534, r19518, MPFR_RNDN); mpfr_add(r19536, r19533, r19535, MPFR_RNDN); ; mpfr_mul(r19538, r19537, r19522, MPFR_RNDN); mpfr_add(r19539, r19536, r19538, MPFR_RNDN); ; mpfr_mul(r19541, r19540, r19521, MPFR_RNDN); mpfr_mul(r19542, r19522, r19506, MPFR_RNDN); mpfr_mul(r19543, r19541, r19542, MPFR_RNDN); mpfr_add(r19544, r19539, r19543, MPFR_RNDN); mpfr_div(r19545, r19524, r19544, MPFR_RNDN); mpfr_mul(r19546, r19545, r19505, MPFR_RNDN); return mpfr_get_d(r19546, MPFR_RNDN); } static mpfr_t r19547, r19548, r19549, r19550, r19551, r19552, r19553, r19554, r19555, r19556, r19557, r19558, r19559, r19560, r19561, r19562, r19563, r19564, r19565, r19566, r19567, r19568, r19569, r19570, r19571, r19572, r19573, r19574, r19575, r19576, r19577, r19578, r19579, r19580, r19581, r19582, r19583, r19584, r19585, r19586, r19587, r19588, r19589, r19590, r19591, r19592, r19593, r19594, r19595, r19596, r19597, r19598; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19547); mpfr_init(r19548); mpfr_init(r19549); mpfr_init(r19550); mpfr_init(r19551); mpfr_init(r19552); mpfr_init_set_str(r19553, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19554, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19555); mpfr_init(r19556); mpfr_init_set_str(r19557, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19558); mpfr_init(r19559); mpfr_init_set_str(r19560, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19561, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19562); mpfr_init(r19563); mpfr_init(r19564); mpfr_init_set_str(r19565, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19566); mpfr_init(r19567); mpfr_init(r19568); mpfr_init(r19569); mpfr_init_set_str(r19570, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19571, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19572); mpfr_init(r19573); mpfr_init(r19574); mpfr_init(r19575); mpfr_init(r19576); mpfr_init(r19577); mpfr_init_set_str(r19578, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19579); mpfr_init(r19580); mpfr_init_set_str(r19581, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19582); mpfr_init(r19583); mpfr_init(r19584); mpfr_init(r19585); mpfr_init(r19586); mpfr_init(r19587); mpfr_init(r19588); mpfr_init(r19589); mpfr_init_set_str(r19590, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19591, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19592); mpfr_init(r19593); mpfr_init(r19594); mpfr_init(r19595); mpfr_init(r19596); mpfr_init(r19597); mpfr_init(r19598); } double f_fm(double x) { mpfr_set_d(r19547, x, MPFR_RNDN); mpfr_mul(r19548, r19547, r19547, MPFR_RNDN); mpfr_mul(r19549, r19548, r19548, MPFR_RNDN); mpfr_mul(r19550, r19547, r19547, MPFR_RNDN); mpfr_mul(r19550, r19550, r19547, MPFR_RNDN); mpfr_mul(r19551, r19550, r19550, MPFR_RNDN); mpfr_mul(r19552, r19549, r19551, MPFR_RNDN); ; ; mpfr_mul(r19555, r19553, r19554, MPFR_RNDN); mpfr_mul(r19556, r19555, r19548, MPFR_RNDN); ; mpfr_add(r19558, r19556, r19557, MPFR_RNDN); mpfr_mul(r19559, r19552, r19558, MPFR_RNDN); ; ; mpfr_mul(r19562, r19547, r19561, MPFR_RNDN); mpfr_mul(r19563, r19562, r19550, MPFR_RNDN); mpfr_add(r19564, r19560, r19563, MPFR_RNDN); ; mpfr_mul(r19566, r19565, r19547, MPFR_RNDN); mpfr_mul(r19567, r19566, r19547, MPFR_RNDN); mpfr_add(r19568, r19564, r19567, MPFR_RNDN); mpfr_add(r19569, r19559, r19568, MPFR_RNDN); ; ; mpfr_mul(r19572, r19571, r19548, MPFR_RNDN); mpfr_add(r19573, r19570, r19572, MPFR_RNDN); mpfr_mul(r19574, r19549, r19573, MPFR_RNDN); mpfr_mul(r19575, r19548, r19574, MPFR_RNDN); mpfr_add(r19576, r19569, r19575, MPFR_RNDN); mpfr_div(r19577, r19547, r19576, MPFR_RNDN); ; mpfr_mul(r19579, r19547, r19578, MPFR_RNDN); mpfr_mul(r19580, r19579, r19550, MPFR_RNDN); ; mpfr_mul(r19582, r19548, r19581, MPFR_RNDN); mpfr_add(r19583, r19582, r19560, MPFR_RNDN); mpfr_add(r19584, r19580, r19583, MPFR_RNDN); mpfr_mul(r19585, r19548, r19554, MPFR_RNDN); mpfr_mul(r19586, r19550, r19547, MPFR_RNDN); mpfr_mul(r19587, r19549, r19586, MPFR_RNDN); mpfr_mul(r19588, r19585, r19587, MPFR_RNDN); mpfr_add(r19589, r19584, r19588, MPFR_RNDN); ; ; mpfr_mul(r19592, r19591, r19547, MPFR_RNDN); mpfr_mul(r19593, r19592, r19547, MPFR_RNDN); mpfr_add(r19594, r19590, r19593, MPFR_RNDN); mpfr_mul(r19595, r19549, r19594, MPFR_RNDN); mpfr_mul(r19596, r19548, r19595, MPFR_RNDN); mpfr_add(r19597, r19589, r19596, MPFR_RNDN); mpfr_mul(r19598, r19577, r19597, MPFR_RNDN); return mpfr_get_d(r19598, MPFR_RNDN); } static mpfr_t r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607, r19608, r19609, r19610, r19611, r19612, r19613, r19614, r19615, r19616, r19617, r19618, r19619, r19620, r19621, r19622, r19623, r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19599); mpfr_init(r19600); mpfr_init(r19601); mpfr_init(r19602); mpfr_init(r19603); mpfr_init(r19604); mpfr_init_set_str(r19605, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19606, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19607); mpfr_init(r19608); mpfr_init_set_str(r19609, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19610); mpfr_init(r19611); mpfr_init_set_str(r19612, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19613, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19614); mpfr_init(r19615); mpfr_init(r19616); mpfr_init_set_str(r19617, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19618); mpfr_init(r19619); mpfr_init(r19620); mpfr_init(r19621); mpfr_init_set_str(r19622, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19623, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19624); mpfr_init(r19625); mpfr_init(r19626); mpfr_init(r19627); mpfr_init(r19628); mpfr_init(r19629); mpfr_init_set_str(r19630, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19631); mpfr_init(r19632); mpfr_init_set_str(r19633, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); mpfr_init(r19637); mpfr_init(r19638); mpfr_init(r19639); mpfr_init(r19640); mpfr_init(r19641); mpfr_init_set_str(r19642, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19643, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19644); mpfr_init(r19645); mpfr_init(r19646); mpfr_init(r19647); mpfr_init(r19648); mpfr_init(r19649); mpfr_init(r19650); } double f_dm(double x) { mpfr_set_d(r19599, x, MPFR_RNDN); mpfr_mul(r19600, r19599, r19599, MPFR_RNDN); mpfr_mul(r19601, r19600, r19600, MPFR_RNDN); mpfr_mul(r19602, r19599, r19599, MPFR_RNDN); mpfr_mul(r19602, r19602, r19599, MPFR_RNDN); mpfr_mul(r19603, r19602, r19602, MPFR_RNDN); mpfr_mul(r19604, r19601, r19603, MPFR_RNDN); ; ; mpfr_mul(r19607, r19605, r19606, MPFR_RNDN); mpfr_mul(r19608, r19607, r19600, MPFR_RNDN); ; mpfr_add(r19610, r19608, r19609, MPFR_RNDN); mpfr_mul(r19611, r19604, r19610, MPFR_RNDN); ; ; mpfr_mul(r19614, r19599, r19613, MPFR_RNDN); mpfr_mul(r19615, r19614, r19602, MPFR_RNDN); mpfr_add(r19616, r19612, r19615, MPFR_RNDN); ; mpfr_mul(r19618, r19617, r19599, MPFR_RNDN); mpfr_mul(r19619, r19618, r19599, MPFR_RNDN); mpfr_add(r19620, r19616, r19619, MPFR_RNDN); mpfr_add(r19621, r19611, r19620, MPFR_RNDN); ; ; mpfr_mul(r19624, r19623, r19600, MPFR_RNDN); mpfr_add(r19625, r19622, r19624, MPFR_RNDN); mpfr_mul(r19626, r19601, r19625, MPFR_RNDN); mpfr_mul(r19627, r19600, r19626, MPFR_RNDN); mpfr_add(r19628, r19621, r19627, MPFR_RNDN); mpfr_div(r19629, r19599, r19628, MPFR_RNDN); ; mpfr_mul(r19631, r19599, r19630, MPFR_RNDN); mpfr_mul(r19632, r19631, r19602, MPFR_RNDN); ; mpfr_mul(r19634, r19600, r19633, MPFR_RNDN); mpfr_add(r19635, r19634, r19612, MPFR_RNDN); mpfr_add(r19636, r19632, r19635, MPFR_RNDN); mpfr_mul(r19637, r19600, r19606, MPFR_RNDN); mpfr_mul(r19638, r19602, r19599, MPFR_RNDN); mpfr_mul(r19639, r19601, r19638, MPFR_RNDN); mpfr_mul(r19640, r19637, r19639, MPFR_RNDN); mpfr_add(r19641, r19636, r19640, MPFR_RNDN); ; ; mpfr_mul(r19644, r19643, r19599, MPFR_RNDN); mpfr_mul(r19645, r19644, r19599, MPFR_RNDN); mpfr_add(r19646, r19642, r19645, MPFR_RNDN); mpfr_mul(r19647, r19601, r19646, MPFR_RNDN); mpfr_mul(r19648, r19600, r19647, MPFR_RNDN); mpfr_add(r19649, r19641, r19648, MPFR_RNDN); mpfr_mul(r19650, r19629, r19649, MPFR_RNDN); return mpfr_get_d(r19650, MPFR_RNDN); }