#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r24486 = 1; float r24487 = 0.1049934947; float r24488 = x; float r24489 = r24488 * r24488; float r24490 = r24487 * r24489; float r24491 = r24486 + r24490; float r24492 = 0.0424060604; float r24493 = r24489 * r24489; float r24494 = r24492 * r24493; float r24495 = r24491 + r24494; float r24496 = 0.0072644182; float r24497 = r24493 * r24489; float r24498 = r24496 * r24497; float r24499 = r24495 + r24498; float r24500 = 0.0005064034; float r24501 = r24497 * r24489; float r24502 = r24500 * r24501; float r24503 = r24499 + r24502; float r24504 = 0.0001789971; float r24505 = r24501 * r24489; float r24506 = r24504 * r24505; float r24507 = r24503 + r24506; float r24508 = 0.7715471019; float r24509 = r24508 * r24489; float r24510 = r24486 + r24509; float r24511 = 0.2909738639; float r24512 = r24511 * r24493; float r24513 = r24510 + r24512; float r24514 = 0.0694555761; float r24515 = r24514 * r24497; float r24516 = r24513 + r24515; float r24517 = 0.0140005442; float r24518 = r24517 * r24501; float r24519 = r24516 + r24518; float r24520 = 0.0008327945; float r24521 = r24520 * r24505; float r24522 = r24519 + r24521; float r24523 = 2; float r24524 = r24523 * r24504; float r24525 = r24505 * r24489; float r24526 = r24524 * r24525; float r24527 = r24522 + r24526; float r24528 = r24507 / r24527; float r24529 = r24528 * r24488; return r24529; } double f_id(double x) { double r24530 = 1; double r24531 = 0.1049934947; double r24532 = x; double r24533 = r24532 * r24532; double r24534 = r24531 * r24533; double r24535 = r24530 + r24534; double r24536 = 0.0424060604; double r24537 = r24533 * r24533; double r24538 = r24536 * r24537; double r24539 = r24535 + r24538; double r24540 = 0.0072644182; double r24541 = r24537 * r24533; double r24542 = r24540 * r24541; double r24543 = r24539 + r24542; double r24544 = 0.0005064034; double r24545 = r24541 * r24533; double r24546 = r24544 * r24545; double r24547 = r24543 + r24546; double r24548 = 0.0001789971; double r24549 = r24545 * r24533; double r24550 = r24548 * r24549; double r24551 = r24547 + r24550; double r24552 = 0.7715471019; double r24553 = r24552 * r24533; double r24554 = r24530 + r24553; double r24555 = 0.2909738639; double r24556 = r24555 * r24537; double r24557 = r24554 + r24556; double r24558 = 0.0694555761; double r24559 = r24558 * r24541; double r24560 = r24557 + r24559; double r24561 = 0.0140005442; double r24562 = r24561 * r24545; double r24563 = r24560 + r24562; double r24564 = 0.0008327945; double r24565 = r24564 * r24549; double r24566 = r24563 + r24565; double r24567 = 2; double r24568 = r24567 * r24548; double r24569 = r24549 * r24533; double r24570 = r24568 * r24569; double r24571 = r24566 + r24570; double r24572 = r24551 / r24571; double r24573 = r24572 * r24532; return r24573; } double f_of(float x) { float r24574 = 0.5; float r24575 = x; float r24576 = r24574 / r24575; float r24577 = 0.2514179000665375; float r24578 = 4; float r24579 = pow(r24575, r24578); float r24580 = r24577 / r24579; float r24581 = 0.15298196345929327; float r24582 = 6; float r24583 = pow(r24575, r24582); float r24584 = r24581 / r24583; float r24585 = r24580 + r24584; float r24586 = r24585 * r24575; float r24587 = r24576 + r24586; float r24588 = -5.4456471621350795e-08; bool r24589 = r24587 <= r24588; float r24590 = 1; float r24591 = 0.1049934947; float r24592 = r24575 * r24575; float r24593 = r24591 * r24592; float r24594 = r24590 + r24593; float r24595 = 0.0424060604; float r24596 = r24592 * r24592; float r24597 = r24595 * r24596; float r24598 = r24594 + r24597; float r24599 = 0.0072644182; float r24600 = r24596 * r24592; float r24601 = r24599 * r24600; float r24602 = r24598 + r24601; float r24603 = 0.0005064034; float r24604 = r24600 * r24592; float r24605 = r24603 * r24604; float r24606 = r24602 + r24605; float r24607 = 0.0001789971; float r24608 = r24604 * r24592; float r24609 = r24607 * r24608; float r24610 = r24606 + r24609; float r24611 = 0.7715471019; float r24612 = r24611 * r24592; float r24613 = r24590 + r24612; float r24614 = 0.2909738639; float r24615 = r24614 * r24596; float r24616 = r24613 + r24615; float r24617 = 0.0694555761; float r24618 = r24617 * r24600; float r24619 = r24616 + r24618; float r24620 = 0.0140005442; float r24621 = r24620 * r24604; float r24622 = r24619 + r24621; float r24623 = 0.0008327945; float r24624 = r24623 * r24608; float r24625 = r24622 + r24624; float r24626 = 2; float r24627 = r24626 * r24607; float r24628 = r24608 * r24592; float r24629 = r24627 * r24628; float r24630 = r24625 + r24629; float r24631 = r24610 / r24630; float r24632 = r24631 * r24575; float r24633 = 1.2852790717449512e-07; bool r24634 = r24587 <= r24633; float r24635 = 0.265709700396151; float r24636 = r24635 * r24579; float r24637 = r24590 + r24636; float r24638 = 0.6665536072; float r24639 = pow(r24575, r24626); float r24640 = r24638 * r24639; float r24641 = r24637 - r24640; float r24642 = r24641 * r24575; float r24643 = r24634 ? r24587 : r24642; float r24644 = r24589 ? r24632 : r24643; return r24644; } double f_od(double x) { double r24645 = 0.5; double r24646 = x; double r24647 = r24645 / r24646; double r24648 = 0.2514179000665375; double r24649 = 4; double r24650 = pow(r24646, r24649); double r24651 = r24648 / r24650; double r24652 = 0.15298196345929327; double r24653 = 6; double r24654 = pow(r24646, r24653); double r24655 = r24652 / r24654; double r24656 = r24651 + r24655; double r24657 = r24656 * r24646; double r24658 = r24647 + r24657; double r24659 = -5.4456471621350795e-08; bool r24660 = r24658 <= r24659; double r24661 = 1; double r24662 = 0.1049934947; double r24663 = r24646 * r24646; double r24664 = r24662 * r24663; double r24665 = r24661 + r24664; double r24666 = 0.0424060604; double r24667 = r24663 * r24663; double r24668 = r24666 * r24667; double r24669 = r24665 + r24668; double r24670 = 0.0072644182; double r24671 = r24667 * r24663; double r24672 = r24670 * r24671; double r24673 = r24669 + r24672; double r24674 = 0.0005064034; double r24675 = r24671 * r24663; double r24676 = r24674 * r24675; double r24677 = r24673 + r24676; double r24678 = 0.0001789971; double r24679 = r24675 * r24663; double r24680 = r24678 * r24679; double r24681 = r24677 + r24680; double r24682 = 0.7715471019; double r24683 = r24682 * r24663; double r24684 = r24661 + r24683; double r24685 = 0.2909738639; double r24686 = r24685 * r24667; double r24687 = r24684 + r24686; double r24688 = 0.0694555761; double r24689 = r24688 * r24671; double r24690 = r24687 + r24689; double r24691 = 0.0140005442; double r24692 = r24691 * r24675; double r24693 = r24690 + r24692; double r24694 = 0.0008327945; double r24695 = r24694 * r24679; double r24696 = r24693 + r24695; double r24697 = 2; double r24698 = r24697 * r24678; double r24699 = r24679 * r24663; double r24700 = r24698 * r24699; double r24701 = r24696 + r24700; double r24702 = r24681 / r24701; double r24703 = r24702 * r24646; double r24704 = 1.2852790717449512e-07; bool r24705 = r24658 <= r24704; double r24706 = 0.265709700396151; double r24707 = r24706 * r24650; double r24708 = r24661 + r24707; double r24709 = 0.6665536072; double r24710 = pow(r24646, r24697); double r24711 = r24709 * r24710; double r24712 = r24708 - r24711; double r24713 = r24712 * r24646; double r24714 = r24705 ? r24658 : r24713; double r24715 = r24660 ? r24703 : r24714; return r24715; } 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 r24716, r24717, r24718, r24719, r24720, r24721, r24722, r24723, r24724, r24725, r24726, r24727, r24728, r24729, r24730, r24731, r24732, r24733, r24734, r24735, r24736, r24737, r24738, r24739, r24740, r24741, r24742, r24743, r24744, r24745, r24746, r24747, r24748, r24749, r24750, r24751, r24752, r24753, r24754, r24755, r24756, r24757, r24758, r24759; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24716, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24717, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r24718); mpfr_init(r24719); mpfr_init(r24720); mpfr_init(r24721); mpfr_init_set_str(r24722, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r24723); mpfr_init(r24724); mpfr_init(r24725); mpfr_init_set_str(r24726, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r24727); mpfr_init(r24728); mpfr_init(r24729); mpfr_init_set_str(r24730, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r24731); mpfr_init(r24732); mpfr_init(r24733); mpfr_init_set_str(r24734, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r24735); mpfr_init(r24736); mpfr_init(r24737); mpfr_init_set_str(r24738, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r24739); mpfr_init(r24740); mpfr_init_set_str(r24741, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r24742); mpfr_init(r24743); mpfr_init_set_str(r24744, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r24745); mpfr_init(r24746); mpfr_init_set_str(r24747, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r24748); mpfr_init(r24749); mpfr_init_set_str(r24750, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r24751); mpfr_init(r24752); mpfr_init_set_str(r24753, "2", 10, MPFR_RNDN); mpfr_init(r24754); mpfr_init(r24755); mpfr_init(r24756); mpfr_init(r24757); mpfr_init(r24758); mpfr_init(r24759); } double f_im(double x) { ; ; mpfr_set_d(r24718, x, MPFR_RNDN); mpfr_mul(r24719, r24718, r24718, MPFR_RNDN); mpfr_mul(r24720, r24717, r24719, MPFR_RNDN); mpfr_add(r24721, r24716, r24720, MPFR_RNDN); ; mpfr_mul(r24723, r24719, r24719, MPFR_RNDN); mpfr_mul(r24724, r24722, r24723, MPFR_RNDN); mpfr_add(r24725, r24721, r24724, MPFR_RNDN); ; mpfr_mul(r24727, r24723, r24719, MPFR_RNDN); mpfr_mul(r24728, r24726, r24727, MPFR_RNDN); mpfr_add(r24729, r24725, r24728, MPFR_RNDN); ; mpfr_mul(r24731, r24727, r24719, MPFR_RNDN); mpfr_mul(r24732, r24730, r24731, MPFR_RNDN); mpfr_add(r24733, r24729, r24732, MPFR_RNDN); ; mpfr_mul(r24735, r24731, r24719, MPFR_RNDN); mpfr_mul(r24736, r24734, r24735, MPFR_RNDN); mpfr_add(r24737, r24733, r24736, MPFR_RNDN); ; mpfr_mul(r24739, r24738, r24719, MPFR_RNDN); mpfr_add(r24740, r24716, r24739, MPFR_RNDN); ; mpfr_mul(r24742, r24741, r24723, MPFR_RNDN); mpfr_add(r24743, r24740, r24742, MPFR_RNDN); ; mpfr_mul(r24745, r24744, r24727, MPFR_RNDN); mpfr_add(r24746, r24743, r24745, MPFR_RNDN); ; mpfr_mul(r24748, r24747, r24731, MPFR_RNDN); mpfr_add(r24749, r24746, r24748, MPFR_RNDN); ; mpfr_mul(r24751, r24750, r24735, MPFR_RNDN); mpfr_add(r24752, r24749, r24751, MPFR_RNDN); ; mpfr_mul(r24754, r24753, r24734, MPFR_RNDN); mpfr_mul(r24755, r24735, r24719, MPFR_RNDN); mpfr_mul(r24756, r24754, r24755, MPFR_RNDN); mpfr_add(r24757, r24752, r24756, MPFR_RNDN); mpfr_div(r24758, r24737, r24757, MPFR_RNDN); mpfr_mul(r24759, r24758, r24718, MPFR_RNDN); return mpfr_get_d(r24759, MPFR_RNDN); } static mpfr_t r24760, r24761, r24762, r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773, r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792, r24793, r24794, r24795, r24796, r24797, r24798, r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812, r24813, r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826, r24827, r24828, r24829, r24830; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24760, "0.5", 10, MPFR_RNDN); mpfr_init(r24761); mpfr_init(r24762); mpfr_init_set_str(r24763, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r24764, "4", 10, MPFR_RNDN); mpfr_init(r24765); mpfr_init(r24766); mpfr_init_set_str(r24767, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r24768, "6", 10, MPFR_RNDN); mpfr_init(r24769); mpfr_init(r24770); mpfr_init(r24771); mpfr_init(r24772); mpfr_init(r24773); mpfr_init_set_str(r24774, "-5.4456471621350795e-08", 10, MPFR_RNDN); mpfr_init(r24775); mpfr_init_set_str(r24776, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24777, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init_set_str(r24781, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r24782); mpfr_init(r24783); mpfr_init(r24784); mpfr_init_set_str(r24785, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init_set_str(r24789, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); mpfr_init_set_str(r24793, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r24794); mpfr_init(r24795); mpfr_init(r24796); mpfr_init_set_str(r24797, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r24798); mpfr_init(r24799); mpfr_init_set_str(r24800, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r24801); mpfr_init(r24802); mpfr_init_set_str(r24803, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r24804); mpfr_init(r24805); mpfr_init_set_str(r24806, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r24807); mpfr_init(r24808); mpfr_init_set_str(r24809, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r24810); mpfr_init(r24811); mpfr_init_set_str(r24812, "2", 10, MPFR_RNDN); mpfr_init(r24813); mpfr_init(r24814); mpfr_init(r24815); mpfr_init(r24816); mpfr_init(r24817); mpfr_init(r24818); mpfr_init_set_str(r24819, "1.2852790717449512e-07", 10, MPFR_RNDN); mpfr_init(r24820); mpfr_init_set_str(r24821, "0.265709700396151", 10, MPFR_RNDN); mpfr_init(r24822); mpfr_init(r24823); mpfr_init_set_str(r24824, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); mpfr_init(r24829); mpfr_init(r24830); } double f_fm(double x) { ; mpfr_set_d(r24761, x, MPFR_RNDN); mpfr_div(r24762, r24760, r24761, MPFR_RNDN); ; ; mpfr_pow(r24765, r24761, r24764, MPFR_RNDN); mpfr_div(r24766, r24763, r24765, MPFR_RNDN); ; ; mpfr_pow(r24769, r24761, r24768, MPFR_RNDN); mpfr_div(r24770, r24767, r24769, MPFR_RNDN); mpfr_add(r24771, r24766, r24770, MPFR_RNDN); mpfr_mul(r24772, r24771, r24761, MPFR_RNDN); mpfr_add(r24773, r24762, r24772, MPFR_RNDN); ; mpfr_set_si(r24775, mpfr_cmp(r24773, r24774) <= 0, MPFR_RNDN); ; ; mpfr_mul(r24778, r24761, r24761, MPFR_RNDN); mpfr_mul(r24779, r24777, r24778, MPFR_RNDN); mpfr_add(r24780, r24776, r24779, MPFR_RNDN); ; mpfr_mul(r24782, r24778, r24778, MPFR_RNDN); mpfr_mul(r24783, r24781, r24782, MPFR_RNDN); mpfr_add(r24784, r24780, r24783, MPFR_RNDN); ; mpfr_mul(r24786, r24782, r24778, MPFR_RNDN); mpfr_mul(r24787, r24785, r24786, MPFR_RNDN); mpfr_add(r24788, r24784, r24787, MPFR_RNDN); ; mpfr_mul(r24790, r24786, r24778, MPFR_RNDN); mpfr_mul(r24791, r24789, r24790, MPFR_RNDN); mpfr_add(r24792, r24788, r24791, MPFR_RNDN); ; mpfr_mul(r24794, r24790, r24778, MPFR_RNDN); mpfr_mul(r24795, r24793, r24794, MPFR_RNDN); mpfr_add(r24796, r24792, r24795, MPFR_RNDN); ; mpfr_mul(r24798, r24797, r24778, MPFR_RNDN); mpfr_add(r24799, r24776, r24798, MPFR_RNDN); ; mpfr_mul(r24801, r24800, r24782, MPFR_RNDN); mpfr_add(r24802, r24799, r24801, MPFR_RNDN); ; mpfr_mul(r24804, r24803, r24786, MPFR_RNDN); mpfr_add(r24805, r24802, r24804, MPFR_RNDN); ; mpfr_mul(r24807, r24806, r24790, MPFR_RNDN); mpfr_add(r24808, r24805, r24807, MPFR_RNDN); ; mpfr_mul(r24810, r24809, r24794, MPFR_RNDN); mpfr_add(r24811, r24808, r24810, MPFR_RNDN); ; mpfr_mul(r24813, r24812, r24793, MPFR_RNDN); mpfr_mul(r24814, r24794, r24778, MPFR_RNDN); mpfr_mul(r24815, r24813, r24814, MPFR_RNDN); mpfr_add(r24816, r24811, r24815, MPFR_RNDN); mpfr_div(r24817, r24796, r24816, MPFR_RNDN); mpfr_mul(r24818, r24817, r24761, MPFR_RNDN); ; mpfr_set_si(r24820, mpfr_cmp(r24773, r24819) <= 0, MPFR_RNDN); ; mpfr_mul(r24822, r24821, r24765, MPFR_RNDN); mpfr_add(r24823, r24776, r24822, MPFR_RNDN); ; mpfr_pow(r24825, r24761, r24812, MPFR_RNDN); mpfr_mul(r24826, r24824, r24825, MPFR_RNDN); mpfr_sub(r24827, r24823, r24826, MPFR_RNDN); mpfr_mul(r24828, r24827, r24761, MPFR_RNDN); if (mpfr_get_si(r24820, MPFR_RNDN)) { mpfr_set(r24829, r24773, MPFR_RNDN); } else { mpfr_set(r24829, r24828, MPFR_RNDN); }; if (mpfr_get_si(r24775, MPFR_RNDN)) { mpfr_set(r24830, r24818, MPFR_RNDN); } else { mpfr_set(r24830, r24829, MPFR_RNDN); }; return mpfr_get_d(r24830, MPFR_RNDN); } static mpfr_t r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840, r24841, r24842, r24843, r24844, r24845, r24846, r24847, r24848, r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868, r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888, r24889, r24890, r24891, r24892, r24893, r24894, r24895, r24896, r24897, r24898, r24899, r24900, r24901; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24831, "0.5", 10, MPFR_RNDN); mpfr_init(r24832); mpfr_init(r24833); mpfr_init_set_str(r24834, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r24835, "4", 10, MPFR_RNDN); mpfr_init(r24836); mpfr_init(r24837); mpfr_init_set_str(r24838, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r24839, "6", 10, MPFR_RNDN); mpfr_init(r24840); mpfr_init(r24841); mpfr_init(r24842); mpfr_init(r24843); mpfr_init(r24844); mpfr_init_set_str(r24845, "-5.4456471621350795e-08", 10, MPFR_RNDN); mpfr_init(r24846); mpfr_init_set_str(r24847, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24848, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r24849); mpfr_init(r24850); mpfr_init(r24851); mpfr_init_set_str(r24852, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r24853); mpfr_init(r24854); mpfr_init(r24855); mpfr_init_set_str(r24856, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r24857); mpfr_init(r24858); mpfr_init(r24859); mpfr_init_set_str(r24860, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r24861); mpfr_init(r24862); mpfr_init(r24863); mpfr_init_set_str(r24864, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r24865); mpfr_init(r24866); mpfr_init(r24867); mpfr_init_set_str(r24868, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r24869); mpfr_init(r24870); mpfr_init_set_str(r24871, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r24872); mpfr_init(r24873); mpfr_init_set_str(r24874, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r24875); mpfr_init(r24876); mpfr_init_set_str(r24877, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r24878); mpfr_init(r24879); mpfr_init_set_str(r24880, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r24881); mpfr_init(r24882); mpfr_init_set_str(r24883, "2", 10, MPFR_RNDN); mpfr_init(r24884); mpfr_init(r24885); mpfr_init(r24886); mpfr_init(r24887); mpfr_init(r24888); mpfr_init(r24889); mpfr_init_set_str(r24890, "1.2852790717449512e-07", 10, MPFR_RNDN); mpfr_init(r24891); mpfr_init_set_str(r24892, "0.265709700396151", 10, MPFR_RNDN); mpfr_init(r24893); mpfr_init(r24894); mpfr_init_set_str(r24895, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r24896); mpfr_init(r24897); mpfr_init(r24898); mpfr_init(r24899); mpfr_init(r24900); mpfr_init(r24901); } double f_dm(double x) { ; mpfr_set_d(r24832, x, MPFR_RNDN); mpfr_div(r24833, r24831, r24832, MPFR_RNDN); ; ; mpfr_pow(r24836, r24832, r24835, MPFR_RNDN); mpfr_div(r24837, r24834, r24836, MPFR_RNDN); ; ; mpfr_pow(r24840, r24832, r24839, MPFR_RNDN); mpfr_div(r24841, r24838, r24840, MPFR_RNDN); mpfr_add(r24842, r24837, r24841, MPFR_RNDN); mpfr_mul(r24843, r24842, r24832, MPFR_RNDN); mpfr_add(r24844, r24833, r24843, MPFR_RNDN); ; mpfr_set_si(r24846, mpfr_cmp(r24844, r24845) <= 0, MPFR_RNDN); ; ; mpfr_mul(r24849, r24832, r24832, MPFR_RNDN); mpfr_mul(r24850, r24848, r24849, MPFR_RNDN); mpfr_add(r24851, r24847, r24850, MPFR_RNDN); ; mpfr_mul(r24853, r24849, r24849, MPFR_RNDN); mpfr_mul(r24854, r24852, r24853, MPFR_RNDN); mpfr_add(r24855, r24851, r24854, MPFR_RNDN); ; mpfr_mul(r24857, r24853, r24849, MPFR_RNDN); mpfr_mul(r24858, r24856, r24857, MPFR_RNDN); mpfr_add(r24859, r24855, r24858, MPFR_RNDN); ; mpfr_mul(r24861, r24857, r24849, MPFR_RNDN); mpfr_mul(r24862, r24860, r24861, MPFR_RNDN); mpfr_add(r24863, r24859, r24862, MPFR_RNDN); ; mpfr_mul(r24865, r24861, r24849, MPFR_RNDN); mpfr_mul(r24866, r24864, r24865, MPFR_RNDN); mpfr_add(r24867, r24863, r24866, MPFR_RNDN); ; mpfr_mul(r24869, r24868, r24849, MPFR_RNDN); mpfr_add(r24870, r24847, r24869, MPFR_RNDN); ; mpfr_mul(r24872, r24871, r24853, MPFR_RNDN); mpfr_add(r24873, r24870, r24872, MPFR_RNDN); ; mpfr_mul(r24875, r24874, r24857, MPFR_RNDN); mpfr_add(r24876, r24873, r24875, MPFR_RNDN); ; mpfr_mul(r24878, r24877, r24861, MPFR_RNDN); mpfr_add(r24879, r24876, r24878, MPFR_RNDN); ; mpfr_mul(r24881, r24880, r24865, MPFR_RNDN); mpfr_add(r24882, r24879, r24881, MPFR_RNDN); ; mpfr_mul(r24884, r24883, r24864, MPFR_RNDN); mpfr_mul(r24885, r24865, r24849, MPFR_RNDN); mpfr_mul(r24886, r24884, r24885, MPFR_RNDN); mpfr_add(r24887, r24882, r24886, MPFR_RNDN); mpfr_div(r24888, r24867, r24887, MPFR_RNDN); mpfr_mul(r24889, r24888, r24832, MPFR_RNDN); ; mpfr_set_si(r24891, mpfr_cmp(r24844, r24890) <= 0, MPFR_RNDN); ; mpfr_mul(r24893, r24892, r24836, MPFR_RNDN); mpfr_add(r24894, r24847, r24893, MPFR_RNDN); ; mpfr_pow(r24896, r24832, r24883, MPFR_RNDN); mpfr_mul(r24897, r24895, r24896, MPFR_RNDN); mpfr_sub(r24898, r24894, r24897, MPFR_RNDN); mpfr_mul(r24899, r24898, r24832, MPFR_RNDN); if (mpfr_get_si(r24891, MPFR_RNDN)) { mpfr_set(r24900, r24844, MPFR_RNDN); } else { mpfr_set(r24900, r24899, MPFR_RNDN); }; if (mpfr_get_si(r24846, MPFR_RNDN)) { mpfr_set(r24901, r24889, MPFR_RNDN); } else { mpfr_set(r24901, r24900, MPFR_RNDN); }; return mpfr_get_d(r24901, MPFR_RNDN); }