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