#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25510 = 1; float r25511 = 0.1049934947; float r25512 = x; float r25513 = r25512 * r25512; float r25514 = r25511 * r25513; float r25515 = r25510 + r25514; float r25516 = 0.0424060604; float r25517 = r25513 * r25513; float r25518 = r25516 * r25517; float r25519 = r25515 + r25518; float r25520 = 0.0072644182; float r25521 = r25517 * r25513; float r25522 = r25520 * r25521; float r25523 = r25519 + r25522; float r25524 = 0.0005064034; float r25525 = r25521 * r25513; float r25526 = r25524 * r25525; float r25527 = r25523 + r25526; float r25528 = 0.0001789971; float r25529 = r25525 * r25513; float r25530 = r25528 * r25529; float r25531 = r25527 + r25530; float r25532 = 0.7715471019; float r25533 = r25532 * r25513; float r25534 = r25510 + r25533; float r25535 = 0.2909738639; float r25536 = r25535 * r25517; float r25537 = r25534 + r25536; float r25538 = 0.0694555761; float r25539 = r25538 * r25521; float r25540 = r25537 + r25539; float r25541 = 0.0140005442; float r25542 = r25541 * r25525; float r25543 = r25540 + r25542; float r25544 = 0.0008327945; float r25545 = r25544 * r25529; float r25546 = r25543 + r25545; float r25547 = 2; float r25548 = r25547 * r25528; float r25549 = r25529 * r25513; float r25550 = r25548 * r25549; float r25551 = r25546 + r25550; float r25552 = r25531 / r25551; float r25553 = r25552 * r25512; return r25553; } double f_id(double x) { double r25554 = 1; double r25555 = 0.1049934947; double r25556 = x; double r25557 = r25556 * r25556; double r25558 = r25555 * r25557; double r25559 = r25554 + r25558; double r25560 = 0.0424060604; double r25561 = r25557 * r25557; double r25562 = r25560 * r25561; double r25563 = r25559 + r25562; double r25564 = 0.0072644182; double r25565 = r25561 * r25557; double r25566 = r25564 * r25565; double r25567 = r25563 + r25566; double r25568 = 0.0005064034; double r25569 = r25565 * r25557; double r25570 = r25568 * r25569; double r25571 = r25567 + r25570; double r25572 = 0.0001789971; double r25573 = r25569 * r25557; double r25574 = r25572 * r25573; double r25575 = r25571 + r25574; double r25576 = 0.7715471019; double r25577 = r25576 * r25557; double r25578 = r25554 + r25577; double r25579 = 0.2909738639; double r25580 = r25579 * r25561; double r25581 = r25578 + r25580; double r25582 = 0.0694555761; double r25583 = r25582 * r25565; double r25584 = r25581 + r25583; double r25585 = 0.0140005442; double r25586 = r25585 * r25569; double r25587 = r25584 + r25586; double r25588 = 0.0008327945; double r25589 = r25588 * r25573; double r25590 = r25587 + r25589; double r25591 = 2; double r25592 = r25591 * r25572; double r25593 = r25573 * r25557; double r25594 = r25592 * r25593; double r25595 = r25590 + r25594; double r25596 = r25575 / r25595; double r25597 = r25596 * r25556; return r25597; } double f_of(float x) { float r25598 = 0.5; float r25599 = x; float r25600 = r25598 / r25599; float r25601 = 0.2514179000665375; float r25602 = r25601 * r25599; float r25603 = 4; float r25604 = pow(r25599, r25603); float r25605 = r25602 / r25604; float r25606 = 0.15298196345929327; float r25607 = r25606 * r25599; float r25608 = 6; float r25609 = pow(r25599, r25608); float r25610 = r25607 / r25609; float r25611 = r25605 + r25610; float r25612 = r25600 + r25611; float r25613 = -1.429618249372218e-10; bool r25614 = r25612 <= r25613; float r25615 = 1; float r25616 = 0.1049934947; float r25617 = r25599 * r25599; float r25618 = r25616 * r25617; float r25619 = r25615 + r25618; float r25620 = 0.0424060604; float r25621 = r25617 * r25617; float r25622 = r25620 * r25621; float r25623 = r25619 + r25622; float r25624 = 0.0072644182; float r25625 = r25621 * r25617; float r25626 = r25624 * r25625; float r25627 = r25623 + r25626; float r25628 = 0.0005064034; float r25629 = r25625 * r25617; float r25630 = r25628 * r25629; float r25631 = r25627 + r25630; float r25632 = 0.0001789971; float r25633 = r25629 * r25617; float r25634 = r25632 * r25633; float r25635 = r25631 + r25634; float r25636 = 0.7715471019; float r25637 = r25636 * r25617; float r25638 = r25615 + r25637; float r25639 = 0.2909738639; float r25640 = r25639 * r25621; float r25641 = r25638 + r25640; float r25642 = 0.0694555761; float r25643 = r25642 * r25625; float r25644 = r25641 + r25643; float r25645 = 0.0140005442; float r25646 = r25645 * r25629; float r25647 = r25644 + r25646; float r25648 = 0.0008327945; float r25649 = r25648 * r25633; float r25650 = r25647 + r25649; float r25651 = 2; float r25652 = r25651 * r25632; float r25653 = r25633 * r25617; float r25654 = r25652 * r25653; float r25655 = r25650 + r25654; float r25656 = r25635 / r25655; float r25657 = r25656 * r25599; float r25658 = 5.4302095356491195e-11; bool r25659 = r25612 <= r25658; float r25660 = 0.265709700396151; float r25661 = r25660 * r25604; float r25662 = r25615 + r25661; float r25663 = 0.6665536072; float r25664 = pow(r25599, r25651); float r25665 = r25663 * r25664; float r25666 = r25662 - r25665; float r25667 = r25666 * r25599; float r25668 = r25659 ? r25612 : r25667; float r25669 = r25614 ? r25657 : r25668; return r25669; } double f_od(double x) { double r25670 = 0.5; double r25671 = x; double r25672 = r25670 / r25671; double r25673 = 0.2514179000665375; double r25674 = r25673 * r25671; double r25675 = 4; double r25676 = pow(r25671, r25675); double r25677 = r25674 / r25676; double r25678 = 0.15298196345929327; double r25679 = r25678 * r25671; double r25680 = 6; double r25681 = pow(r25671, r25680); double r25682 = r25679 / r25681; double r25683 = r25677 + r25682; double r25684 = r25672 + r25683; double r25685 = -1.429618249372218e-10; bool r25686 = r25684 <= r25685; double r25687 = 1; double r25688 = 0.1049934947; double r25689 = r25671 * r25671; double r25690 = r25688 * r25689; double r25691 = r25687 + r25690; double r25692 = 0.0424060604; double r25693 = r25689 * r25689; double r25694 = r25692 * r25693; double r25695 = r25691 + r25694; double r25696 = 0.0072644182; double r25697 = r25693 * r25689; double r25698 = r25696 * r25697; double r25699 = r25695 + r25698; double r25700 = 0.0005064034; double r25701 = r25697 * r25689; double r25702 = r25700 * r25701; double r25703 = r25699 + r25702; double r25704 = 0.0001789971; double r25705 = r25701 * r25689; double r25706 = r25704 * r25705; double r25707 = r25703 + r25706; double r25708 = 0.7715471019; double r25709 = r25708 * r25689; double r25710 = r25687 + r25709; double r25711 = 0.2909738639; double r25712 = r25711 * r25693; double r25713 = r25710 + r25712; double r25714 = 0.0694555761; double r25715 = r25714 * r25697; double r25716 = r25713 + r25715; double r25717 = 0.0140005442; double r25718 = r25717 * r25701; double r25719 = r25716 + r25718; double r25720 = 0.0008327945; double r25721 = r25720 * r25705; double r25722 = r25719 + r25721; double r25723 = 2; double r25724 = r25723 * r25704; double r25725 = r25705 * r25689; double r25726 = r25724 * r25725; double r25727 = r25722 + r25726; double r25728 = r25707 / r25727; double r25729 = r25728 * r25671; double r25730 = 5.4302095356491195e-11; bool r25731 = r25684 <= r25730; double r25732 = 0.265709700396151; double r25733 = r25732 * r25676; double r25734 = r25687 + r25733; double r25735 = 0.6665536072; double r25736 = pow(r25671, r25723); double r25737 = r25735 * r25736; double r25738 = r25734 - r25737; double r25739 = r25738 * r25671; double r25740 = r25731 ? r25684 : r25739; double r25741 = r25686 ? r25729 : r25740; return r25741; } 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 r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755, r25756, r25757, r25758, r25759, r25760, r25761, r25762, r25763, r25764, r25765, r25766, r25767, r25768, r25769, r25770, r25771, r25772, r25773, r25774, r25775, r25776, r25777, r25778, r25779, r25780, r25781, r25782, r25783, r25784, r25785; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r25742, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25743, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init(r25747); mpfr_init_set_str(r25748, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25749); mpfr_init(r25750); mpfr_init(r25751); mpfr_init_set_str(r25752, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init_set_str(r25756, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25757); mpfr_init(r25758); mpfr_init(r25759); mpfr_init_set_str(r25760, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25761); mpfr_init(r25762); mpfr_init(r25763); mpfr_init_set_str(r25764, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25765); mpfr_init(r25766); mpfr_init_set_str(r25767, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25768); mpfr_init(r25769); mpfr_init_set_str(r25770, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25771); mpfr_init(r25772); mpfr_init_set_str(r25773, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25774); mpfr_init(r25775); mpfr_init_set_str(r25776, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25777); mpfr_init(r25778); mpfr_init_set_str(r25779, "2", 10, MPFR_RNDN); mpfr_init(r25780); mpfr_init(r25781); mpfr_init(r25782); mpfr_init(r25783); mpfr_init(r25784); mpfr_init(r25785); } double f_im(double x) { ; ; mpfr_set_d(r25744, x, MPFR_RNDN); mpfr_mul(r25745, r25744, r25744, MPFR_RNDN); mpfr_mul(r25746, r25743, r25745, MPFR_RNDN); mpfr_add(r25747, r25742, r25746, MPFR_RNDN); ; mpfr_mul(r25749, r25745, r25745, MPFR_RNDN); mpfr_mul(r25750, r25748, r25749, MPFR_RNDN); mpfr_add(r25751, r25747, r25750, MPFR_RNDN); ; mpfr_mul(r25753, r25749, r25745, MPFR_RNDN); mpfr_mul(r25754, r25752, r25753, MPFR_RNDN); mpfr_add(r25755, r25751, r25754, MPFR_RNDN); ; mpfr_mul(r25757, r25753, r25745, MPFR_RNDN); mpfr_mul(r25758, r25756, r25757, MPFR_RNDN); mpfr_add(r25759, r25755, r25758, MPFR_RNDN); ; mpfr_mul(r25761, r25757, r25745, MPFR_RNDN); mpfr_mul(r25762, r25760, r25761, MPFR_RNDN); mpfr_add(r25763, r25759, r25762, MPFR_RNDN); ; mpfr_mul(r25765, r25764, r25745, MPFR_RNDN); mpfr_add(r25766, r25742, r25765, MPFR_RNDN); ; mpfr_mul(r25768, r25767, r25749, MPFR_RNDN); mpfr_add(r25769, r25766, r25768, MPFR_RNDN); ; mpfr_mul(r25771, r25770, r25753, MPFR_RNDN); mpfr_add(r25772, r25769, r25771, MPFR_RNDN); ; mpfr_mul(r25774, r25773, r25757, MPFR_RNDN); mpfr_add(r25775, r25772, r25774, MPFR_RNDN); ; mpfr_mul(r25777, r25776, r25761, MPFR_RNDN); mpfr_add(r25778, r25775, r25777, MPFR_RNDN); ; mpfr_mul(r25780, r25779, r25760, MPFR_RNDN); mpfr_mul(r25781, r25761, r25745, MPFR_RNDN); mpfr_mul(r25782, r25780, r25781, MPFR_RNDN); mpfr_add(r25783, r25778, r25782, MPFR_RNDN); mpfr_div(r25784, r25763, r25783, MPFR_RNDN); mpfr_mul(r25785, r25784, r25744, MPFR_RNDN); return mpfr_get_d(r25785, MPFR_RNDN); } static mpfr_t r25786, r25787, r25788, r25789, r25790, r25791, r25792, r25793, r25794, r25795, r25796, r25797, r25798, r25799, r25800, r25801, r25802, r25803, r25804, r25805, r25806, r25807, r25808, r25809, r25810, r25811, r25812, r25813, r25814, r25815, r25816, r25817, r25818, r25819, r25820, r25821, r25822, r25823, r25824, r25825, r25826, r25827, r25828, r25829, r25830, r25831, r25832, r25833, r25834, r25835, r25836, r25837, r25838, r25839, r25840, r25841, r25842, r25843, r25844, r25845, r25846, r25847, r25848, r25849, r25850, r25851, r25852, r25853, r25854, r25855, r25856, r25857; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25786, "0.5", 10, MPFR_RNDN); mpfr_init(r25787); mpfr_init(r25788); mpfr_init_set_str(r25789, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r25790); mpfr_init_set_str(r25791, "4", 10, MPFR_RNDN); mpfr_init(r25792); mpfr_init(r25793); mpfr_init_set_str(r25794, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r25795); mpfr_init_set_str(r25796, "6", 10, MPFR_RNDN); mpfr_init(r25797); mpfr_init(r25798); mpfr_init(r25799); mpfr_init(r25800); mpfr_init_set_str(r25801, "-1.429618249372218e-10", 10, MPFR_RNDN); mpfr_init(r25802); mpfr_init_set_str(r25803, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25804, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25805); mpfr_init(r25806); mpfr_init(r25807); mpfr_init_set_str(r25808, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25809); mpfr_init(r25810); mpfr_init(r25811); mpfr_init_set_str(r25812, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25813); mpfr_init(r25814); mpfr_init(r25815); mpfr_init_set_str(r25816, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25817); mpfr_init(r25818); mpfr_init(r25819); mpfr_init_set_str(r25820, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25821); mpfr_init(r25822); mpfr_init(r25823); mpfr_init_set_str(r25824, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25825); mpfr_init(r25826); mpfr_init_set_str(r25827, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25828); mpfr_init(r25829); mpfr_init_set_str(r25830, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25831); mpfr_init(r25832); mpfr_init_set_str(r25833, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25834); mpfr_init(r25835); mpfr_init_set_str(r25836, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25837); mpfr_init(r25838); mpfr_init_set_str(r25839, "2", 10, MPFR_RNDN); mpfr_init(r25840); mpfr_init(r25841); mpfr_init(r25842); mpfr_init(r25843); mpfr_init(r25844); mpfr_init(r25845); mpfr_init_set_str(r25846, "5.4302095356491195e-11", 10, MPFR_RNDN); mpfr_init(r25847); mpfr_init_set_str(r25848, "0.265709700396151", 10, MPFR_RNDN); mpfr_init(r25849); mpfr_init(r25850); mpfr_init_set_str(r25851, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r25852); mpfr_init(r25853); mpfr_init(r25854); mpfr_init(r25855); mpfr_init(r25856); mpfr_init(r25857); } double f_fm(double x) { ; mpfr_set_d(r25787, x, MPFR_RNDN); mpfr_div(r25788, r25786, r25787, MPFR_RNDN); ; mpfr_mul(r25790, r25789, r25787, MPFR_RNDN); ; mpfr_pow(r25792, r25787, r25791, MPFR_RNDN); mpfr_div(r25793, r25790, r25792, MPFR_RNDN); ; mpfr_mul(r25795, r25794, r25787, MPFR_RNDN); ; mpfr_pow(r25797, r25787, r25796, MPFR_RNDN); mpfr_div(r25798, r25795, r25797, MPFR_RNDN); mpfr_add(r25799, r25793, r25798, MPFR_RNDN); mpfr_add(r25800, r25788, r25799, MPFR_RNDN); ; mpfr_set_si(r25802, mpfr_cmp(r25800, r25801) <= 0, MPFR_RNDN); ; ; mpfr_mul(r25805, r25787, r25787, MPFR_RNDN); mpfr_mul(r25806, r25804, r25805, MPFR_RNDN); mpfr_add(r25807, r25803, r25806, MPFR_RNDN); ; mpfr_mul(r25809, r25805, r25805, MPFR_RNDN); mpfr_mul(r25810, r25808, r25809, MPFR_RNDN); mpfr_add(r25811, r25807, r25810, MPFR_RNDN); ; mpfr_mul(r25813, r25809, r25805, MPFR_RNDN); mpfr_mul(r25814, r25812, r25813, MPFR_RNDN); mpfr_add(r25815, r25811, r25814, MPFR_RNDN); ; mpfr_mul(r25817, r25813, r25805, MPFR_RNDN); mpfr_mul(r25818, r25816, r25817, MPFR_RNDN); mpfr_add(r25819, r25815, r25818, MPFR_RNDN); ; mpfr_mul(r25821, r25817, r25805, MPFR_RNDN); mpfr_mul(r25822, r25820, r25821, MPFR_RNDN); mpfr_add(r25823, r25819, r25822, MPFR_RNDN); ; mpfr_mul(r25825, r25824, r25805, MPFR_RNDN); mpfr_add(r25826, r25803, r25825, MPFR_RNDN); ; mpfr_mul(r25828, r25827, r25809, MPFR_RNDN); mpfr_add(r25829, r25826, r25828, MPFR_RNDN); ; mpfr_mul(r25831, r25830, r25813, MPFR_RNDN); mpfr_add(r25832, r25829, r25831, MPFR_RNDN); ; mpfr_mul(r25834, r25833, r25817, MPFR_RNDN); mpfr_add(r25835, r25832, r25834, MPFR_RNDN); ; mpfr_mul(r25837, r25836, r25821, MPFR_RNDN); mpfr_add(r25838, r25835, r25837, MPFR_RNDN); ; mpfr_mul(r25840, r25839, r25820, MPFR_RNDN); mpfr_mul(r25841, r25821, r25805, MPFR_RNDN); mpfr_mul(r25842, r25840, r25841, MPFR_RNDN); mpfr_add(r25843, r25838, r25842, MPFR_RNDN); mpfr_div(r25844, r25823, r25843, MPFR_RNDN); mpfr_mul(r25845, r25844, r25787, MPFR_RNDN); ; mpfr_set_si(r25847, mpfr_cmp(r25800, r25846) <= 0, MPFR_RNDN); ; mpfr_mul(r25849, r25848, r25792, MPFR_RNDN); mpfr_add(r25850, r25803, r25849, MPFR_RNDN); ; mpfr_pow(r25852, r25787, r25839, MPFR_RNDN); mpfr_mul(r25853, r25851, r25852, MPFR_RNDN); mpfr_sub(r25854, r25850, r25853, MPFR_RNDN); mpfr_mul(r25855, r25854, r25787, MPFR_RNDN); if (mpfr_get_si(r25847, MPFR_RNDN)) { mpfr_set(r25856, r25800, MPFR_RNDN); } else { mpfr_set(r25856, r25855, MPFR_RNDN); }; if (mpfr_get_si(r25802, MPFR_RNDN)) { mpfr_set(r25857, r25845, MPFR_RNDN); } else { mpfr_set(r25857, r25856, MPFR_RNDN); }; return mpfr_get_d(r25857, MPFR_RNDN); } static mpfr_t r25858, r25859, r25860, r25861, r25862, r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871, r25872, r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889, r25890, r25891, r25892, r25893, r25894, r25895, r25896, r25897, r25898, r25899, r25900, r25901, r25902, r25903, r25904, r25905, r25906, r25907, r25908, r25909, r25910, r25911, r25912, r25913, r25914, r25915, r25916, r25917, r25918, r25919, r25920, r25921, r25922, r25923, r25924, r25925, r25926, r25927, r25928, r25929; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25858, "0.5", 10, MPFR_RNDN); mpfr_init(r25859); mpfr_init(r25860); mpfr_init_set_str(r25861, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r25862); mpfr_init_set_str(r25863, "4", 10, MPFR_RNDN); mpfr_init(r25864); mpfr_init(r25865); mpfr_init_set_str(r25866, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r25867); mpfr_init_set_str(r25868, "6", 10, MPFR_RNDN); mpfr_init(r25869); mpfr_init(r25870); mpfr_init(r25871); mpfr_init(r25872); mpfr_init_set_str(r25873, "-1.429618249372218e-10", 10, MPFR_RNDN); mpfr_init(r25874); mpfr_init_set_str(r25875, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25876, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25877); mpfr_init(r25878); mpfr_init(r25879); mpfr_init_set_str(r25880, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25881); mpfr_init(r25882); mpfr_init(r25883); mpfr_init_set_str(r25884, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25885); mpfr_init(r25886); mpfr_init(r25887); mpfr_init_set_str(r25888, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25889); mpfr_init(r25890); mpfr_init(r25891); mpfr_init_set_str(r25892, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25893); mpfr_init(r25894); mpfr_init(r25895); mpfr_init_set_str(r25896, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25897); mpfr_init(r25898); mpfr_init_set_str(r25899, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25900); mpfr_init(r25901); mpfr_init_set_str(r25902, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25903); mpfr_init(r25904); mpfr_init_set_str(r25905, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25906); mpfr_init(r25907); mpfr_init_set_str(r25908, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25909); mpfr_init(r25910); mpfr_init_set_str(r25911, "2", 10, MPFR_RNDN); mpfr_init(r25912); mpfr_init(r25913); mpfr_init(r25914); mpfr_init(r25915); mpfr_init(r25916); mpfr_init(r25917); mpfr_init_set_str(r25918, "5.4302095356491195e-11", 10, MPFR_RNDN); mpfr_init(r25919); mpfr_init_set_str(r25920, "0.265709700396151", 10, MPFR_RNDN); mpfr_init(r25921); mpfr_init(r25922); mpfr_init_set_str(r25923, "0.6665536072", 10, MPFR_RNDN); mpfr_init(r25924); mpfr_init(r25925); mpfr_init(r25926); mpfr_init(r25927); mpfr_init(r25928); mpfr_init(r25929); } double f_dm(double x) { ; mpfr_set_d(r25859, x, MPFR_RNDN); mpfr_div(r25860, r25858, r25859, MPFR_RNDN); ; mpfr_mul(r25862, r25861, r25859, MPFR_RNDN); ; mpfr_pow(r25864, r25859, r25863, MPFR_RNDN); mpfr_div(r25865, r25862, r25864, MPFR_RNDN); ; mpfr_mul(r25867, r25866, r25859, MPFR_RNDN); ; mpfr_pow(r25869, r25859, r25868, MPFR_RNDN); mpfr_div(r25870, r25867, r25869, MPFR_RNDN); mpfr_add(r25871, r25865, r25870, MPFR_RNDN); mpfr_add(r25872, r25860, r25871, MPFR_RNDN); ; mpfr_set_si(r25874, mpfr_cmp(r25872, r25873) <= 0, MPFR_RNDN); ; ; mpfr_mul(r25877, r25859, r25859, MPFR_RNDN); mpfr_mul(r25878, r25876, r25877, MPFR_RNDN); mpfr_add(r25879, r25875, r25878, MPFR_RNDN); ; mpfr_mul(r25881, r25877, r25877, MPFR_RNDN); mpfr_mul(r25882, r25880, r25881, MPFR_RNDN); mpfr_add(r25883, r25879, r25882, MPFR_RNDN); ; mpfr_mul(r25885, r25881, r25877, MPFR_RNDN); mpfr_mul(r25886, r25884, r25885, MPFR_RNDN); mpfr_add(r25887, r25883, r25886, MPFR_RNDN); ; mpfr_mul(r25889, r25885, r25877, MPFR_RNDN); mpfr_mul(r25890, r25888, r25889, MPFR_RNDN); mpfr_add(r25891, r25887, r25890, MPFR_RNDN); ; mpfr_mul(r25893, r25889, r25877, MPFR_RNDN); mpfr_mul(r25894, r25892, r25893, MPFR_RNDN); mpfr_add(r25895, r25891, r25894, MPFR_RNDN); ; mpfr_mul(r25897, r25896, r25877, MPFR_RNDN); mpfr_add(r25898, r25875, r25897, MPFR_RNDN); ; mpfr_mul(r25900, r25899, r25881, MPFR_RNDN); mpfr_add(r25901, r25898, r25900, MPFR_RNDN); ; mpfr_mul(r25903, r25902, r25885, MPFR_RNDN); mpfr_add(r25904, r25901, r25903, MPFR_RNDN); ; mpfr_mul(r25906, r25905, r25889, MPFR_RNDN); mpfr_add(r25907, r25904, r25906, MPFR_RNDN); ; mpfr_mul(r25909, r25908, r25893, MPFR_RNDN); mpfr_add(r25910, r25907, r25909, MPFR_RNDN); ; mpfr_mul(r25912, r25911, r25892, MPFR_RNDN); mpfr_mul(r25913, r25893, r25877, MPFR_RNDN); mpfr_mul(r25914, r25912, r25913, MPFR_RNDN); mpfr_add(r25915, r25910, r25914, MPFR_RNDN); mpfr_div(r25916, r25895, r25915, MPFR_RNDN); mpfr_mul(r25917, r25916, r25859, MPFR_RNDN); ; mpfr_set_si(r25919, mpfr_cmp(r25872, r25918) <= 0, MPFR_RNDN); ; mpfr_mul(r25921, r25920, r25864, MPFR_RNDN); mpfr_add(r25922, r25875, r25921, MPFR_RNDN); ; mpfr_pow(r25924, r25859, r25911, MPFR_RNDN); mpfr_mul(r25925, r25923, r25924, MPFR_RNDN); mpfr_sub(r25926, r25922, r25925, MPFR_RNDN); mpfr_mul(r25927, r25926, r25859, MPFR_RNDN); if (mpfr_get_si(r25919, MPFR_RNDN)) { mpfr_set(r25928, r25872, MPFR_RNDN); } else { mpfr_set(r25928, r25927, MPFR_RNDN); }; if (mpfr_get_si(r25874, MPFR_RNDN)) { mpfr_set(r25929, r25917, MPFR_RNDN); } else { mpfr_set(r25929, r25928, MPFR_RNDN); }; return mpfr_get_d(r25929, MPFR_RNDN); }