#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19463 = 1.0f; float r19464 = 0.1049934947f; float r19465 = x; float r19466 = r19465 * r19465; float r19467 = r19464 * r19466; float r19468 = r19463 + r19467; float r19469 = 0.0424060604f; float r19470 = r19466 * r19466; float r19471 = r19469 * r19470; float r19472 = r19468 + r19471; float r19473 = 0.0072644182f; float r19474 = r19470 * r19466; float r19475 = r19473 * r19474; float r19476 = r19472 + r19475; float r19477 = 0.0005064034f; float r19478 = r19474 * r19466; float r19479 = r19477 * r19478; float r19480 = r19476 + r19479; float r19481 = 0.0001789971f; float r19482 = r19478 * r19466; float r19483 = r19481 * r19482; float r19484 = r19480 + r19483; float r19485 = 0.7715471019f; float r19486 = r19485 * r19466; float r19487 = r19463 + r19486; float r19488 = 0.2909738639f; float r19489 = r19488 * r19470; float r19490 = r19487 + r19489; float r19491 = 0.0694555761f; float r19492 = r19491 * r19474; float r19493 = r19490 + r19492; float r19494 = 0.0140005442f; float r19495 = r19494 * r19478; float r19496 = r19493 + r19495; float r19497 = 0.0008327945f; float r19498 = r19497 * r19482; float r19499 = r19496 + r19498; float r19500 = 2.0f; float r19501 = r19500 * r19481; float r19502 = r19482 * r19466; float r19503 = r19501 * r19502; float r19504 = r19499 + r19503; float r19505 = r19484 / r19504; float r19506 = r19505 * r19465; return r19506; } double f_id(double x) { double r19507 = 1.0; double r19508 = 0.1049934947; double r19509 = x; double r19510 = r19509 * r19509; double r19511 = r19508 * r19510; double r19512 = r19507 + r19511; double r19513 = 0.0424060604; double r19514 = r19510 * r19510; double r19515 = r19513 * r19514; double r19516 = r19512 + r19515; double r19517 = 0.0072644182; double r19518 = r19514 * r19510; double r19519 = r19517 * r19518; double r19520 = r19516 + r19519; double r19521 = 0.0005064034; double r19522 = r19518 * r19510; double r19523 = r19521 * r19522; double r19524 = r19520 + r19523; double r19525 = 0.0001789971; double r19526 = r19522 * r19510; double r19527 = r19525 * r19526; double r19528 = r19524 + r19527; double r19529 = 0.7715471019; double r19530 = r19529 * r19510; double r19531 = r19507 + r19530; double r19532 = 0.2909738639; double r19533 = r19532 * r19514; double r19534 = r19531 + r19533; double r19535 = 0.0694555761; double r19536 = r19535 * r19518; double r19537 = r19534 + r19536; double r19538 = 0.0140005442; double r19539 = r19538 * r19522; double r19540 = r19537 + r19539; double r19541 = 0.0008327945; double r19542 = r19541 * r19526; double r19543 = r19540 + r19542; double r19544 = 2.0; double r19545 = r19544 * r19525; double r19546 = r19526 * r19510; double r19547 = r19545 * r19546; double r19548 = r19543 + r19547; double r19549 = r19528 / r19548; double r19550 = r19549 * r19509; return r19550; } double f_of(float x) { float r19551 = x; float r19552 = r19551 * r19551; float r19553 = r19552 * r19552; float r19554 = r19551 * (r19551 * r19551); float r19555 = r19554 * r19554; float r19556 = r19553 * r19555; float r19557 = 2.0f; float r19558 = 0.0001789971f; float r19559 = r19557 * r19558; float r19560 = r19559 * r19552; float r19561 = 0.0008327945f; float r19562 = r19560 + r19561; float r19563 = r19556 * r19562; float r19564 = 1.0f; float r19565 = 0.2909738639f; float r19566 = r19551 * r19565; float r19567 = r19566 * r19554; float r19568 = r19564 + r19567; float r19569 = 0.7715471019f; float r19570 = r19569 * r19551; float r19571 = r19570 * r19551; float r19572 = r19568 + r19571; float r19573 = r19563 + r19572; float r19574 = 0.0694555761f; float r19575 = 0.0140005442f; float r19576 = r19575 * r19552; float r19577 = r19574 + r19576; float r19578 = r19553 * r19577; float r19579 = r19552 * r19578; float r19580 = r19573 + r19579; float r19581 = r19551 / r19580; float r19582 = 0.0424060604f; float r19583 = r19551 * r19582; float r19584 = r19583 * r19554; float r19585 = 0.1049934947f; float r19586 = r19552 * r19585; float r19587 = r19586 + r19564; float r19588 = r19584 + r19587; float r19589 = r19552 * r19558; float r19590 = r19554 * r19551; float r19591 = r19553 * r19590; float r19592 = r19589 * r19591; float r19593 = r19588 + r19592; float r19594 = 0.0072644182f; float r19595 = 0.0005064034f; float r19596 = r19595 * r19551; float r19597 = r19596 * r19551; float r19598 = r19594 + r19597; float r19599 = r19553 * r19598; float r19600 = r19552 * r19599; float r19601 = r19593 + r19600; float r19602 = r19581 * r19601; return r19602; } double f_od(double x) { double r19603 = x; double r19604 = r19603 * r19603; double r19605 = r19604 * r19604; double r19606 = r19603 * (r19603 * r19603); double r19607 = r19606 * r19606; double r19608 = r19605 * r19607; double r19609 = 2.0; double r19610 = 0.0001789971; double r19611 = r19609 * r19610; double r19612 = r19611 * r19604; double r19613 = 0.0008327945; double r19614 = r19612 + r19613; double r19615 = r19608 * r19614; double r19616 = 1.0; double r19617 = 0.2909738639; double r19618 = r19603 * r19617; double r19619 = r19618 * r19606; double r19620 = r19616 + r19619; double r19621 = 0.7715471019; double r19622 = r19621 * r19603; double r19623 = r19622 * r19603; double r19624 = r19620 + r19623; double r19625 = r19615 + r19624; double r19626 = 0.0694555761; double r19627 = 0.0140005442; double r19628 = r19627 * r19604; double r19629 = r19626 + r19628; double r19630 = r19605 * r19629; double r19631 = r19604 * r19630; double r19632 = r19625 + r19631; double r19633 = r19603 / r19632; double r19634 = 0.0424060604; double r19635 = r19603 * r19634; double r19636 = r19635 * r19606; double r19637 = 0.1049934947; double r19638 = r19604 * r19637; double r19639 = r19638 + r19616; double r19640 = r19636 + r19639; double r19641 = r19604 * r19610; double r19642 = r19606 * r19603; double r19643 = r19605 * r19642; double r19644 = r19641 * r19643; double r19645 = r19640 + r19644; double r19646 = 0.0072644182; double r19647 = 0.0005064034; double r19648 = r19647 * r19603; double r19649 = r19648 * r19603; double r19650 = r19646 + r19649; double r19651 = r19605 * r19650; double r19652 = r19604 * r19651; double r19653 = r19645 + r19652; double r19654 = r19633 * r19653; return r19654; } 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 r19655, r19656, r19657, r19658, r19659, r19660, r19661, r19662, r19663, r19664, r19665, r19666, r19667, r19668, r19669, r19670, r19671, r19672, r19673, r19674, r19675, r19676, r19677, r19678, r19679, r19680, r19681, r19682, r19683, r19684, r19685, r19686, r19687, r19688, r19689, r19690, r19691, r19692, r19693, r19694, r19695, r19696, r19697, r19698; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19655, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19656, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19657); mpfr_init(r19658); mpfr_init(r19659); mpfr_init(r19660); mpfr_init_set_str(r19661, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19662); mpfr_init(r19663); mpfr_init(r19664); mpfr_init_set_str(r19665, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19666); mpfr_init(r19667); mpfr_init(r19668); mpfr_init_set_str(r19669, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19670); mpfr_init(r19671); mpfr_init(r19672); mpfr_init_set_str(r19673, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19674); mpfr_init(r19675); mpfr_init(r19676); mpfr_init_set_str(r19677, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19678); mpfr_init(r19679); mpfr_init_set_str(r19680, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19681); mpfr_init(r19682); mpfr_init_set_str(r19683, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19684); mpfr_init(r19685); mpfr_init_set_str(r19686, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19687); mpfr_init(r19688); mpfr_init_set_str(r19689, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19690); mpfr_init(r19691); mpfr_init_set_str(r19692, "2", 10, MPFR_RNDN); mpfr_init(r19693); mpfr_init(r19694); mpfr_init(r19695); mpfr_init(r19696); mpfr_init(r19697); mpfr_init(r19698); } double f_im(double x) { ; ; mpfr_set_d(r19657, x, MPFR_RNDN); mpfr_mul(r19658, r19657, r19657, MPFR_RNDN); mpfr_mul(r19659, r19656, r19658, MPFR_RNDN); mpfr_add(r19660, r19655, r19659, MPFR_RNDN); ; mpfr_mul(r19662, r19658, r19658, MPFR_RNDN); mpfr_mul(r19663, r19661, r19662, MPFR_RNDN); mpfr_add(r19664, r19660, r19663, MPFR_RNDN); ; mpfr_mul(r19666, r19662, r19658, MPFR_RNDN); mpfr_mul(r19667, r19665, r19666, MPFR_RNDN); mpfr_add(r19668, r19664, r19667, MPFR_RNDN); ; mpfr_mul(r19670, r19666, r19658, MPFR_RNDN); mpfr_mul(r19671, r19669, r19670, MPFR_RNDN); mpfr_add(r19672, r19668, r19671, MPFR_RNDN); ; mpfr_mul(r19674, r19670, r19658, MPFR_RNDN); mpfr_mul(r19675, r19673, r19674, MPFR_RNDN); mpfr_add(r19676, r19672, r19675, MPFR_RNDN); ; mpfr_mul(r19678, r19677, r19658, MPFR_RNDN); mpfr_add(r19679, r19655, r19678, MPFR_RNDN); ; mpfr_mul(r19681, r19680, r19662, MPFR_RNDN); mpfr_add(r19682, r19679, r19681, MPFR_RNDN); ; mpfr_mul(r19684, r19683, r19666, MPFR_RNDN); mpfr_add(r19685, r19682, r19684, MPFR_RNDN); ; mpfr_mul(r19687, r19686, r19670, MPFR_RNDN); mpfr_add(r19688, r19685, r19687, MPFR_RNDN); ; mpfr_mul(r19690, r19689, r19674, MPFR_RNDN); mpfr_add(r19691, r19688, r19690, MPFR_RNDN); ; mpfr_mul(r19693, r19692, r19673, MPFR_RNDN); mpfr_mul(r19694, r19674, r19658, MPFR_RNDN); mpfr_mul(r19695, r19693, r19694, MPFR_RNDN); mpfr_add(r19696, r19691, r19695, MPFR_RNDN); mpfr_div(r19697, r19676, r19696, MPFR_RNDN); mpfr_mul(r19698, r19697, r19657, MPFR_RNDN); return mpfr_get_d(r19698, MPFR_RNDN); } static mpfr_t r19699, r19700, r19701, r19702, r19703, r19704, r19705, r19706, r19707, r19708, r19709, r19710, r19711, r19712, r19713, r19714, r19715, r19716, r19717, r19718, r19719, r19720, r19721, r19722, r19723, r19724, r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19699); mpfr_init(r19700); mpfr_init(r19701); mpfr_init(r19702); mpfr_init(r19703); mpfr_init(r19704); mpfr_init_set_str(r19705, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19706, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19707); mpfr_init(r19708); mpfr_init_set_str(r19709, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19710); mpfr_init(r19711); mpfr_init_set_str(r19712, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19713, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19714); mpfr_init(r19715); mpfr_init(r19716); mpfr_init_set_str(r19717, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19718); mpfr_init(r19719); mpfr_init(r19720); mpfr_init(r19721); mpfr_init_set_str(r19722, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19723, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19724); mpfr_init(r19725); mpfr_init(r19726); mpfr_init(r19727); mpfr_init(r19728); mpfr_init(r19729); mpfr_init_set_str(r19730, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19731); mpfr_init(r19732); mpfr_init_set_str(r19733, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19734); mpfr_init(r19735); mpfr_init(r19736); mpfr_init(r19737); mpfr_init(r19738); mpfr_init(r19739); mpfr_init(r19740); mpfr_init(r19741); mpfr_init_set_str(r19742, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19743, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19744); mpfr_init(r19745); mpfr_init(r19746); mpfr_init(r19747); mpfr_init(r19748); mpfr_init(r19749); mpfr_init(r19750); } double f_fm(double x) { mpfr_set_d(r19699, x, MPFR_RNDN); mpfr_mul(r19700, r19699, r19699, MPFR_RNDN); mpfr_mul(r19701, r19700, r19700, MPFR_RNDN); mpfr_mul(r19702, r19699, r19699, MPFR_RNDN); mpfr_mul(r19702, r19702, r19699, MPFR_RNDN); mpfr_mul(r19703, r19702, r19702, MPFR_RNDN); mpfr_mul(r19704, r19701, r19703, MPFR_RNDN); ; ; mpfr_mul(r19707, r19705, r19706, MPFR_RNDN); mpfr_mul(r19708, r19707, r19700, MPFR_RNDN); ; mpfr_add(r19710, r19708, r19709, MPFR_RNDN); mpfr_mul(r19711, r19704, r19710, MPFR_RNDN); ; ; mpfr_mul(r19714, r19699, r19713, MPFR_RNDN); mpfr_mul(r19715, r19714, r19702, MPFR_RNDN); mpfr_add(r19716, r19712, r19715, MPFR_RNDN); ; mpfr_mul(r19718, r19717, r19699, MPFR_RNDN); mpfr_mul(r19719, r19718, r19699, MPFR_RNDN); mpfr_add(r19720, r19716, r19719, MPFR_RNDN); mpfr_add(r19721, r19711, r19720, MPFR_RNDN); ; ; mpfr_mul(r19724, r19723, r19700, MPFR_RNDN); mpfr_add(r19725, r19722, r19724, MPFR_RNDN); mpfr_mul(r19726, r19701, r19725, MPFR_RNDN); mpfr_mul(r19727, r19700, r19726, MPFR_RNDN); mpfr_add(r19728, r19721, r19727, MPFR_RNDN); mpfr_div(r19729, r19699, r19728, MPFR_RNDN); ; mpfr_mul(r19731, r19699, r19730, MPFR_RNDN); mpfr_mul(r19732, r19731, r19702, MPFR_RNDN); ; mpfr_mul(r19734, r19700, r19733, MPFR_RNDN); mpfr_add(r19735, r19734, r19712, MPFR_RNDN); mpfr_add(r19736, r19732, r19735, MPFR_RNDN); mpfr_mul(r19737, r19700, r19706, MPFR_RNDN); mpfr_mul(r19738, r19702, r19699, MPFR_RNDN); mpfr_mul(r19739, r19701, r19738, MPFR_RNDN); mpfr_mul(r19740, r19737, r19739, MPFR_RNDN); mpfr_add(r19741, r19736, r19740, MPFR_RNDN); ; ; mpfr_mul(r19744, r19743, r19699, MPFR_RNDN); mpfr_mul(r19745, r19744, r19699, MPFR_RNDN); mpfr_add(r19746, r19742, r19745, MPFR_RNDN); mpfr_mul(r19747, r19701, r19746, MPFR_RNDN); mpfr_mul(r19748, r19700, r19747, MPFR_RNDN); mpfr_add(r19749, r19741, r19748, MPFR_RNDN); mpfr_mul(r19750, r19729, r19749, MPFR_RNDN); return mpfr_get_d(r19750, MPFR_RNDN); } static mpfr_t r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19751); mpfr_init(r19752); mpfr_init(r19753); mpfr_init(r19754); mpfr_init(r19755); mpfr_init(r19756); mpfr_init_set_str(r19757, "2", 10, MPFR_RNDN); mpfr_init_set_str(r19758, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19759); mpfr_init(r19760); mpfr_init_set_str(r19761, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19762); mpfr_init(r19763); mpfr_init_set_str(r19764, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19765, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19766); mpfr_init(r19767); mpfr_init(r19768); mpfr_init_set_str(r19769, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19770); mpfr_init(r19771); mpfr_init(r19772); mpfr_init(r19773); mpfr_init_set_str(r19774, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19775, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19776); mpfr_init(r19777); mpfr_init(r19778); mpfr_init(r19779); mpfr_init(r19780); mpfr_init(r19781); mpfr_init_set_str(r19782, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19783); mpfr_init(r19784); mpfr_init_set_str(r19785, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19786); mpfr_init(r19787); mpfr_init(r19788); mpfr_init(r19789); mpfr_init(r19790); mpfr_init(r19791); mpfr_init(r19792); mpfr_init(r19793); mpfr_init_set_str(r19794, "0.0072644182", 10, MPFR_RNDN); mpfr_init_set_str(r19795, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19796); mpfr_init(r19797); mpfr_init(r19798); mpfr_init(r19799); mpfr_init(r19800); mpfr_init(r19801); mpfr_init(r19802); } double f_dm(double x) { mpfr_set_d(r19751, x, MPFR_RNDN); mpfr_mul(r19752, r19751, r19751, MPFR_RNDN); mpfr_mul(r19753, r19752, r19752, MPFR_RNDN); mpfr_mul(r19754, r19751, r19751, MPFR_RNDN); mpfr_mul(r19754, r19754, r19751, MPFR_RNDN); mpfr_mul(r19755, r19754, r19754, MPFR_RNDN); mpfr_mul(r19756, r19753, r19755, MPFR_RNDN); ; ; mpfr_mul(r19759, r19757, r19758, MPFR_RNDN); mpfr_mul(r19760, r19759, r19752, MPFR_RNDN); ; mpfr_add(r19762, r19760, r19761, MPFR_RNDN); mpfr_mul(r19763, r19756, r19762, MPFR_RNDN); ; ; mpfr_mul(r19766, r19751, r19765, MPFR_RNDN); mpfr_mul(r19767, r19766, r19754, MPFR_RNDN); mpfr_add(r19768, r19764, r19767, MPFR_RNDN); ; mpfr_mul(r19770, r19769, r19751, MPFR_RNDN); mpfr_mul(r19771, r19770, r19751, MPFR_RNDN); mpfr_add(r19772, r19768, r19771, MPFR_RNDN); mpfr_add(r19773, r19763, r19772, MPFR_RNDN); ; ; mpfr_mul(r19776, r19775, r19752, MPFR_RNDN); mpfr_add(r19777, r19774, r19776, MPFR_RNDN); mpfr_mul(r19778, r19753, r19777, MPFR_RNDN); mpfr_mul(r19779, r19752, r19778, MPFR_RNDN); mpfr_add(r19780, r19773, r19779, MPFR_RNDN); mpfr_div(r19781, r19751, r19780, MPFR_RNDN); ; mpfr_mul(r19783, r19751, r19782, MPFR_RNDN); mpfr_mul(r19784, r19783, r19754, MPFR_RNDN); ; mpfr_mul(r19786, r19752, r19785, MPFR_RNDN); mpfr_add(r19787, r19786, r19764, MPFR_RNDN); mpfr_add(r19788, r19784, r19787, MPFR_RNDN); mpfr_mul(r19789, r19752, r19758, MPFR_RNDN); mpfr_mul(r19790, r19754, r19751, MPFR_RNDN); mpfr_mul(r19791, r19753, r19790, MPFR_RNDN); mpfr_mul(r19792, r19789, r19791, MPFR_RNDN); mpfr_add(r19793, r19788, r19792, MPFR_RNDN); ; ; mpfr_mul(r19796, r19795, r19751, MPFR_RNDN); mpfr_mul(r19797, r19796, r19751, MPFR_RNDN); mpfr_add(r19798, r19794, r19797, MPFR_RNDN); mpfr_mul(r19799, r19753, r19798, MPFR_RNDN); mpfr_mul(r19800, r19752, r19799, MPFR_RNDN); mpfr_add(r19801, r19793, r19800, MPFR_RNDN); mpfr_mul(r19802, r19781, r19801, MPFR_RNDN); return mpfr_get_d(r19802, MPFR_RNDN); }