#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r11501 = 1; float r11502 = 0.1049934947; float r11503 = x; float r11504 = r11503 * r11503; float r11505 = r11502 * r11504; float r11506 = r11501 + r11505; float r11507 = 0.0424060604; float r11508 = r11504 * r11504; float r11509 = r11507 * r11508; float r11510 = r11506 + r11509; float r11511 = 0.0072644182; float r11512 = r11508 * r11504; float r11513 = r11511 * r11512; float r11514 = r11510 + r11513; float r11515 = 0.0005064034; float r11516 = r11512 * r11504; float r11517 = r11515 * r11516; float r11518 = r11514 + r11517; float r11519 = 0.0001789971; float r11520 = r11516 * r11504; float r11521 = r11519 * r11520; float r11522 = r11518 + r11521; float r11523 = 0.7715471019; float r11524 = r11523 * r11504; float r11525 = r11501 + r11524; float r11526 = 0.2909738639; float r11527 = r11526 * r11508; float r11528 = r11525 + r11527; float r11529 = 0.0694555761; float r11530 = r11529 * r11512; float r11531 = r11528 + r11530; float r11532 = 0.0140005442; float r11533 = r11532 * r11516; float r11534 = r11531 + r11533; float r11535 = 0.0008327945; float r11536 = r11535 * r11520; float r11537 = r11534 + r11536; float r11538 = 2; float r11539 = r11538 * r11519; float r11540 = r11520 * r11504; float r11541 = r11539 * r11540; float r11542 = r11537 + r11541; float r11543 = r11522 / r11542; float r11544 = r11543 * r11503; return r11544; } double f_id(double x) { double r11545 = 1; double r11546 = 0.1049934947; double r11547 = x; double r11548 = r11547 * r11547; double r11549 = r11546 * r11548; double r11550 = r11545 + r11549; double r11551 = 0.0424060604; double r11552 = r11548 * r11548; double r11553 = r11551 * r11552; double r11554 = r11550 + r11553; double r11555 = 0.0072644182; double r11556 = r11552 * r11548; double r11557 = r11555 * r11556; double r11558 = r11554 + r11557; double r11559 = 0.0005064034; double r11560 = r11556 * r11548; double r11561 = r11559 * r11560; double r11562 = r11558 + r11561; double r11563 = 0.0001789971; double r11564 = r11560 * r11548; double r11565 = r11563 * r11564; double r11566 = r11562 + r11565; double r11567 = 0.7715471019; double r11568 = r11567 * r11548; double r11569 = r11545 + r11568; double r11570 = 0.2909738639; double r11571 = r11570 * r11552; double r11572 = r11569 + r11571; double r11573 = 0.0694555761; double r11574 = r11573 * r11556; double r11575 = r11572 + r11574; double r11576 = 0.0140005442; double r11577 = r11576 * r11560; double r11578 = r11575 + r11577; double r11579 = 0.0008327945; double r11580 = r11579 * r11564; double r11581 = r11578 + r11580; double r11582 = 2; double r11583 = r11582 * r11563; double r11584 = r11564 * r11548; double r11585 = r11583 * r11584; double r11586 = r11581 + r11585; double r11587 = r11566 / r11586; double r11588 = r11587 * r11547; return r11588; } double f_of(float x) { float r11589 = x; float r11590 = -65822.2175525377; bool r11591 = r11589 <= r11590; float r11592 = 659.3418912140232; bool r11593 = r11589 <= r11592; float r11594 = !r11593; bool r11595 = r11591 || r11594; float r11596 = 0.15298196345929327; float r11597 = 6; float r11598 = pow(r11589, r11597); float r11599 = r11596 / r11598; float r11600 = 0.2514179000665375; float r11601 = 4; float r11602 = pow(r11589, r11601); float r11603 = r11600 / r11602; float r11604 = r11599 + r11603; float r11605 = 0.5; float r11606 = r11605 / r11589; float r11607 = fma(r11604, r11589, r11606); float r11608 = 0.0005064034; float r11609 = r11589 * r11589; float r11610 = r11609 * r11609; float r11611 = r11610 * r11610; float r11612 = 0.0424060604; float r11613 = 0.1049934947; float r11614 = 1; float r11615 = fma(r11613, r11609, r11614); float r11616 = fma(r11610, r11612, r11615); float r11617 = 0.0072644182; float r11618 = r11609 * r11617; float r11619 = r11610 * r11618; float r11620 = r11616 + r11619; float r11621 = fma(r11608, r11611, r11620); float r11622 = 0.0001789971; float r11623 = r11622 * r11609; float r11624 = r11611 * r11623; float r11625 = r11621 + r11624; float r11626 = 0.0008327945; float r11627 = r11609 * r11626; float r11628 = r11627 * r11611; float r11629 = 2; float r11630 = r11629 * r11622; float r11631 = r11610 * r11630; float r11632 = r11631 * r11611; float r11633 = r11628 + r11632; float r11634 = 0.0140005442; float r11635 = 0.2909738639; float r11636 = 0.7715471019; float r11637 = fma(r11609, r11636, r11614); float r11638 = fma(r11635, r11610, r11637); float r11639 = 0.0694555761; float r11640 = r11639 * r11609; float r11641 = r11640 * r11610; float r11642 = r11638 + r11641; float r11643 = fma(r11611, r11634, r11642); float r11644 = r11633 + r11643; float r11645 = r11625 / r11644; float r11646 = r11589 * r11645; float r11647 = r11595 ? r11607 : r11646; return r11647; } double f_od(double x) { double r11648 = x; double r11649 = -65822.2175525377; bool r11650 = r11648 <= r11649; double r11651 = 659.3418912140232; bool r11652 = r11648 <= r11651; double r11653 = !r11652; bool r11654 = r11650 || r11653; double r11655 = 0.15298196345929327; double r11656 = 6; double r11657 = pow(r11648, r11656); double r11658 = r11655 / r11657; double r11659 = 0.2514179000665375; double r11660 = 4; double r11661 = pow(r11648, r11660); double r11662 = r11659 / r11661; double r11663 = r11658 + r11662; double r11664 = 0.5; double r11665 = r11664 / r11648; double r11666 = fma(r11663, r11648, r11665); double r11667 = 0.0005064034; double r11668 = r11648 * r11648; double r11669 = r11668 * r11668; double r11670 = r11669 * r11669; double r11671 = 0.0424060604; double r11672 = 0.1049934947; double r11673 = 1; double r11674 = fma(r11672, r11668, r11673); double r11675 = fma(r11669, r11671, r11674); double r11676 = 0.0072644182; double r11677 = r11668 * r11676; double r11678 = r11669 * r11677; double r11679 = r11675 + r11678; double r11680 = fma(r11667, r11670, r11679); double r11681 = 0.0001789971; double r11682 = r11681 * r11668; double r11683 = r11670 * r11682; double r11684 = r11680 + r11683; double r11685 = 0.0008327945; double r11686 = r11668 * r11685; double r11687 = r11686 * r11670; double r11688 = 2; double r11689 = r11688 * r11681; double r11690 = r11669 * r11689; double r11691 = r11690 * r11670; double r11692 = r11687 + r11691; double r11693 = 0.0140005442; double r11694 = 0.2909738639; double r11695 = 0.7715471019; double r11696 = fma(r11668, r11695, r11673); double r11697 = fma(r11694, r11669, r11696); double r11698 = 0.0694555761; double r11699 = r11698 * r11668; double r11700 = r11699 * r11669; double r11701 = r11697 + r11700; double r11702 = fma(r11670, r11693, r11701); double r11703 = r11692 + r11702; double r11704 = r11684 / r11703; double r11705 = r11648 * r11704; double r11706 = r11654 ? r11666 : r11705; return r11706; } 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 r11707, r11708, r11709, r11710, r11711, r11712, r11713, r11714, r11715, r11716, r11717, r11718, r11719, r11720, r11721, r11722, r11723, r11724, r11725, r11726, r11727, r11728, r11729, r11730, r11731, r11732, r11733, r11734, r11735, r11736, r11737, r11738, r11739, r11740, r11741, r11742, r11743, r11744, r11745, r11746, r11747, r11748, r11749, r11750; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r11707, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11708, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r11709); mpfr_init(r11710); mpfr_init(r11711); mpfr_init(r11712); mpfr_init_set_str(r11713, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11714); mpfr_init(r11715); mpfr_init(r11716); mpfr_init_set_str(r11717, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11718); mpfr_init(r11719); mpfr_init(r11720); mpfr_init_set_str(r11721, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11722); mpfr_init(r11723); mpfr_init(r11724); mpfr_init_set_str(r11725, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11726); mpfr_init(r11727); mpfr_init(r11728); mpfr_init_set_str(r11729, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11730); mpfr_init(r11731); mpfr_init_set_str(r11732, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r11733); mpfr_init(r11734); mpfr_init_set_str(r11735, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11736); mpfr_init(r11737); mpfr_init_set_str(r11738, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r11739); mpfr_init(r11740); mpfr_init_set_str(r11741, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11742); mpfr_init(r11743); mpfr_init_set_str(r11744, "2", 10, MPFR_RNDN); mpfr_init(r11745); mpfr_init(r11746); mpfr_init(r11747); mpfr_init(r11748); mpfr_init(r11749); mpfr_init(r11750); } double f_im(double x) { ; ; mpfr_set_d(r11709, x, MPFR_RNDN); mpfr_mul(r11710, r11709, r11709, MPFR_RNDN); mpfr_mul(r11711, r11708, r11710, MPFR_RNDN); mpfr_add(r11712, r11707, r11711, MPFR_RNDN); ; mpfr_mul(r11714, r11710, r11710, MPFR_RNDN); mpfr_mul(r11715, r11713, r11714, MPFR_RNDN); mpfr_add(r11716, r11712, r11715, MPFR_RNDN); ; mpfr_mul(r11718, r11714, r11710, MPFR_RNDN); mpfr_mul(r11719, r11717, r11718, MPFR_RNDN); mpfr_add(r11720, r11716, r11719, MPFR_RNDN); ; mpfr_mul(r11722, r11718, r11710, MPFR_RNDN); mpfr_mul(r11723, r11721, r11722, MPFR_RNDN); mpfr_add(r11724, r11720, r11723, MPFR_RNDN); ; mpfr_mul(r11726, r11722, r11710, MPFR_RNDN); mpfr_mul(r11727, r11725, r11726, MPFR_RNDN); mpfr_add(r11728, r11724, r11727, MPFR_RNDN); ; mpfr_mul(r11730, r11729, r11710, MPFR_RNDN); mpfr_add(r11731, r11707, r11730, MPFR_RNDN); ; mpfr_mul(r11733, r11732, r11714, MPFR_RNDN); mpfr_add(r11734, r11731, r11733, MPFR_RNDN); ; mpfr_mul(r11736, r11735, r11718, MPFR_RNDN); mpfr_add(r11737, r11734, r11736, MPFR_RNDN); ; mpfr_mul(r11739, r11738, r11722, MPFR_RNDN); mpfr_add(r11740, r11737, r11739, MPFR_RNDN); ; mpfr_mul(r11742, r11741, r11726, MPFR_RNDN); mpfr_add(r11743, r11740, r11742, MPFR_RNDN); ; mpfr_mul(r11745, r11744, r11725, MPFR_RNDN); mpfr_mul(r11746, r11726, r11710, MPFR_RNDN); mpfr_mul(r11747, r11745, r11746, MPFR_RNDN); mpfr_add(r11748, r11743, r11747, MPFR_RNDN); mpfr_div(r11749, r11728, r11748, MPFR_RNDN); mpfr_mul(r11750, r11749, r11709, MPFR_RNDN); return mpfr_get_d(r11750, MPFR_RNDN); } static mpfr_t r11751, r11752, r11753, r11754, r11755, r11756, r11757, r11758, r11759, r11760, r11761, r11762, r11763, r11764, r11765, r11766, r11767, r11768, r11769, r11770, r11771, r11772, r11773, r11774, r11775, r11776, r11777, r11778, r11779, r11780, r11781, r11782, r11783, r11784, r11785, r11786, r11787, r11788, r11789, r11790, r11791, r11792, r11793, r11794, r11795, r11796, r11797, r11798, r11799, r11800, r11801, r11802, r11803, r11804, r11805, r11806, r11807, r11808, r11809; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r11751); mpfr_init_set_str(r11752, "-65822.2175525377", 10, MPFR_RNDN); mpfr_init(r11753); mpfr_init_set_str(r11754, "659.3418912140232", 10, MPFR_RNDN); mpfr_init(r11755); mpfr_init(r11756); mpfr_init(r11757); mpfr_init_set_str(r11758, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r11759, "6", 10, MPFR_RNDN); mpfr_init(r11760); mpfr_init(r11761); mpfr_init_set_str(r11762, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r11763, "4", 10, MPFR_RNDN); mpfr_init(r11764); mpfr_init(r11765); mpfr_init(r11766); mpfr_init_set_str(r11767, "0.5", 10, MPFR_RNDN); mpfr_init(r11768); mpfr_init(r11769); mpfr_init_set_str(r11770, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11771); mpfr_init(r11772); mpfr_init(r11773); mpfr_init_set_str(r11774, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r11775, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r11776, "1", 10, MPFR_RNDN); mpfr_init(r11777); mpfr_init(r11778); mpfr_init_set_str(r11779, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11780); mpfr_init(r11781); mpfr_init(r11782); mpfr_init(r11783); mpfr_init_set_str(r11784, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11785); mpfr_init(r11786); mpfr_init(r11787); mpfr_init_set_str(r11788, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11789); mpfr_init(r11790); mpfr_init_set_str(r11791, "2", 10, MPFR_RNDN); mpfr_init(r11792); mpfr_init(r11793); mpfr_init(r11794); mpfr_init(r11795); mpfr_init_set_str(r11796, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r11797, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r11798, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11799); mpfr_init(r11800); mpfr_init_set_str(r11801, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11802); mpfr_init(r11803); mpfr_init(r11804); mpfr_init(r11805); mpfr_init(r11806); mpfr_init(r11807); mpfr_init(r11808); mpfr_init(r11809); } double f_fm(double x) { mpfr_set_d(r11751, x, MPFR_RNDN); ; mpfr_set_si(r11753, mpfr_cmp(r11751, r11752) <= 0, MPFR_RNDN); ; mpfr_set_si(r11755, mpfr_cmp(r11751, r11754) <= 0, MPFR_RNDN); mpfr_set_si(r11756, !mpfr_get_si(r11755, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11757, mpfr_get_si(r11753, MPFR_RNDN) || mpfr_get_si(r11756, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r11760, r11751, r11759, MPFR_RNDN); mpfr_div(r11761, r11758, r11760, MPFR_RNDN); ; ; mpfr_pow(r11764, r11751, r11763, MPFR_RNDN); mpfr_div(r11765, r11762, r11764, MPFR_RNDN); mpfr_add(r11766, r11761, r11765, MPFR_RNDN); ; mpfr_div(r11768, r11767, r11751, MPFR_RNDN); mpfr_fma(r11769, r11766, r11751, r11768, MPFR_RNDN); ; mpfr_mul(r11771, r11751, r11751, MPFR_RNDN); mpfr_mul(r11772, r11771, r11771, MPFR_RNDN); mpfr_mul(r11773, r11772, r11772, MPFR_RNDN); ; ; ; mpfr_fma(r11777, r11775, r11771, r11776, MPFR_RNDN); mpfr_fma(r11778, r11772, r11774, r11777, MPFR_RNDN); ; mpfr_mul(r11780, r11771, r11779, MPFR_RNDN); mpfr_mul(r11781, r11772, r11780, MPFR_RNDN); mpfr_add(r11782, r11778, r11781, MPFR_RNDN); mpfr_fma(r11783, r11770, r11773, r11782, MPFR_RNDN); ; mpfr_mul(r11785, r11784, r11771, MPFR_RNDN); mpfr_mul(r11786, r11773, r11785, MPFR_RNDN); mpfr_add(r11787, r11783, r11786, MPFR_RNDN); ; mpfr_mul(r11789, r11771, r11788, MPFR_RNDN); mpfr_mul(r11790, r11789, r11773, MPFR_RNDN); ; mpfr_mul(r11792, r11791, r11784, MPFR_RNDN); mpfr_mul(r11793, r11772, r11792, MPFR_RNDN); mpfr_mul(r11794, r11793, r11773, MPFR_RNDN); mpfr_add(r11795, r11790, r11794, MPFR_RNDN); ; ; ; mpfr_fma(r11799, r11771, r11798, r11776, MPFR_RNDN); mpfr_fma(r11800, r11797, r11772, r11799, MPFR_RNDN); ; mpfr_mul(r11802, r11801, r11771, MPFR_RNDN); mpfr_mul(r11803, r11802, r11772, MPFR_RNDN); mpfr_add(r11804, r11800, r11803, MPFR_RNDN); mpfr_fma(r11805, r11773, r11796, r11804, MPFR_RNDN); mpfr_add(r11806, r11795, r11805, MPFR_RNDN); mpfr_div(r11807, r11787, r11806, MPFR_RNDN); mpfr_mul(r11808, r11751, r11807, MPFR_RNDN); if (mpfr_get_si(r11757, MPFR_RNDN)) { mpfr_set(r11809, r11769, MPFR_RNDN); } else { mpfr_set(r11809, r11808, MPFR_RNDN); }; return mpfr_get_d(r11809, MPFR_RNDN); } static mpfr_t r11810, r11811, r11812, r11813, r11814, r11815, r11816, r11817, r11818, r11819, r11820, r11821, r11822, r11823, r11824, r11825, r11826, r11827, r11828, r11829, r11830, r11831, r11832, r11833, r11834, r11835, r11836, r11837, r11838, r11839, r11840, r11841, r11842, r11843, r11844, r11845, r11846, r11847, r11848, r11849, r11850, r11851, r11852, r11853, r11854, r11855, r11856, r11857, r11858, r11859, r11860, r11861, r11862, r11863, r11864, r11865, r11866, r11867, r11868; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11810); mpfr_init_set_str(r11811, "-65822.2175525377", 10, MPFR_RNDN); mpfr_init(r11812); mpfr_init_set_str(r11813, "659.3418912140232", 10, MPFR_RNDN); mpfr_init(r11814); mpfr_init(r11815); mpfr_init(r11816); mpfr_init_set_str(r11817, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r11818, "6", 10, MPFR_RNDN); mpfr_init(r11819); mpfr_init(r11820); mpfr_init_set_str(r11821, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r11822, "4", 10, MPFR_RNDN); mpfr_init(r11823); mpfr_init(r11824); mpfr_init(r11825); mpfr_init_set_str(r11826, "0.5", 10, MPFR_RNDN); mpfr_init(r11827); mpfr_init(r11828); mpfr_init_set_str(r11829, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11830); mpfr_init(r11831); mpfr_init(r11832); mpfr_init_set_str(r11833, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r11834, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r11835, "1", 10, MPFR_RNDN); mpfr_init(r11836); mpfr_init(r11837); mpfr_init_set_str(r11838, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r11839); mpfr_init(r11840); mpfr_init(r11841); mpfr_init(r11842); mpfr_init_set_str(r11843, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r11844); mpfr_init(r11845); mpfr_init(r11846); mpfr_init_set_str(r11847, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r11848); mpfr_init(r11849); mpfr_init_set_str(r11850, "2", 10, MPFR_RNDN); mpfr_init(r11851); mpfr_init(r11852); mpfr_init(r11853); mpfr_init(r11854); mpfr_init_set_str(r11855, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r11856, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r11857, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r11858); mpfr_init(r11859); mpfr_init_set_str(r11860, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r11861); mpfr_init(r11862); mpfr_init(r11863); mpfr_init(r11864); mpfr_init(r11865); mpfr_init(r11866); mpfr_init(r11867); mpfr_init(r11868); } double f_dm(double x) { mpfr_set_d(r11810, x, MPFR_RNDN); ; mpfr_set_si(r11812, mpfr_cmp(r11810, r11811) <= 0, MPFR_RNDN); ; mpfr_set_si(r11814, mpfr_cmp(r11810, r11813) <= 0, MPFR_RNDN); mpfr_set_si(r11815, !mpfr_get_si(r11814, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11816, mpfr_get_si(r11812, MPFR_RNDN) || mpfr_get_si(r11815, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r11819, r11810, r11818, MPFR_RNDN); mpfr_div(r11820, r11817, r11819, MPFR_RNDN); ; ; mpfr_pow(r11823, r11810, r11822, MPFR_RNDN); mpfr_div(r11824, r11821, r11823, MPFR_RNDN); mpfr_add(r11825, r11820, r11824, MPFR_RNDN); ; mpfr_div(r11827, r11826, r11810, MPFR_RNDN); mpfr_fma(r11828, r11825, r11810, r11827, MPFR_RNDN); ; mpfr_mul(r11830, r11810, r11810, MPFR_RNDN); mpfr_mul(r11831, r11830, r11830, MPFR_RNDN); mpfr_mul(r11832, r11831, r11831, MPFR_RNDN); ; ; ; mpfr_fma(r11836, r11834, r11830, r11835, MPFR_RNDN); mpfr_fma(r11837, r11831, r11833, r11836, MPFR_RNDN); ; mpfr_mul(r11839, r11830, r11838, MPFR_RNDN); mpfr_mul(r11840, r11831, r11839, MPFR_RNDN); mpfr_add(r11841, r11837, r11840, MPFR_RNDN); mpfr_fma(r11842, r11829, r11832, r11841, MPFR_RNDN); ; mpfr_mul(r11844, r11843, r11830, MPFR_RNDN); mpfr_mul(r11845, r11832, r11844, MPFR_RNDN); mpfr_add(r11846, r11842, r11845, MPFR_RNDN); ; mpfr_mul(r11848, r11830, r11847, MPFR_RNDN); mpfr_mul(r11849, r11848, r11832, MPFR_RNDN); ; mpfr_mul(r11851, r11850, r11843, MPFR_RNDN); mpfr_mul(r11852, r11831, r11851, MPFR_RNDN); mpfr_mul(r11853, r11852, r11832, MPFR_RNDN); mpfr_add(r11854, r11849, r11853, MPFR_RNDN); ; ; ; mpfr_fma(r11858, r11830, r11857, r11835, MPFR_RNDN); mpfr_fma(r11859, r11856, r11831, r11858, MPFR_RNDN); ; mpfr_mul(r11861, r11860, r11830, MPFR_RNDN); mpfr_mul(r11862, r11861, r11831, MPFR_RNDN); mpfr_add(r11863, r11859, r11862, MPFR_RNDN); mpfr_fma(r11864, r11832, r11855, r11863, MPFR_RNDN); mpfr_add(r11865, r11854, r11864, MPFR_RNDN); mpfr_div(r11866, r11846, r11865, MPFR_RNDN); mpfr_mul(r11867, r11810, r11866, MPFR_RNDN); if (mpfr_get_si(r11816, MPFR_RNDN)) { mpfr_set(r11868, r11828, MPFR_RNDN); } else { mpfr_set(r11868, r11867, MPFR_RNDN); }; return mpfr_get_d(r11868, MPFR_RNDN); }