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