#include #include #include #include #include char *name = "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r6491 = x; float r6492 = 18.0; float r6493 = r6491 * r6492; float r6494 = y; float r6495 = r6493 * r6494; float r6496 = z; float r6497 = r6495 * r6496; float r6498 = t; float r6499 = r6497 * r6498; float r6500 = a; float r6501 = 4.0; float r6502 = r6500 * r6501; float r6503 = r6502 * r6498; float r6504 = r6499 - r6503; float r6505 = b; float r6506 = c; float r6507 = r6505 * r6506; float r6508 = r6504 + r6507; float r6509 = r6491 * r6501; float r6510 = i; float r6511 = r6509 * r6510; float r6512 = r6508 - r6511; float r6513 = j; float r6514 = 27.0; float r6515 = r6513 * r6514; float r6516 = k; float r6517 = r6515 * r6516; float r6518 = r6512 - r6517; return r6518; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r6519 = x; double r6520 = 18.0; double r6521 = r6519 * r6520; double r6522 = y; double r6523 = r6521 * r6522; double r6524 = z; double r6525 = r6523 * r6524; double r6526 = t; double r6527 = r6525 * r6526; double r6528 = a; double r6529 = 4.0; double r6530 = r6528 * r6529; double r6531 = r6530 * r6526; double r6532 = r6527 - r6531; double r6533 = b; double r6534 = c; double r6535 = r6533 * r6534; double r6536 = r6532 + r6535; double r6537 = r6519 * r6529; double r6538 = i; double r6539 = r6537 * r6538; double r6540 = r6536 - r6539; double r6541 = j; double r6542 = 27.0; double r6543 = r6541 * r6542; double r6544 = k; double r6545 = r6543 * r6544; double r6546 = r6540 - r6545; return r6546; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r6547 = x; float r6548 = 18.0; float r6549 = r6547 * r6548; float r6550 = y; float r6551 = r6549 * r6550; float r6552 = z; float r6553 = r6551 * r6552; float r6554 = t; float r6555 = r6553 * r6554; float r6556 = a; float r6557 = 4.0; float r6558 = r6556 * r6557; float r6559 = r6558 * r6554; float r6560 = r6555 - r6559; float r6561 = b; float r6562 = c; float r6563 = r6561 * r6562; float r6564 = r6560 + r6563; float r6565 = r6547 * r6557; float r6566 = i; float r6567 = r6565 * r6566; float r6568 = r6564 - r6567; float r6569 = -inf.0; bool r6570 = r6568 <= r6569; float r6571 = r6554 * r6547; float r6572 = r6550 * r6571; float r6573 = r6552 * r6572; float r6574 = r6548 * r6573; float r6575 = 27.0; float r6576 = j; float r6577 = k; float r6578 = r6576 * r6577; float r6579 = r6575 * r6578; float r6580 = r6556 * r6554; float r6581 = r6557 * r6580; float r6582 = r6579 + r6581; float r6583 = r6574 - r6582; float r6584 = 3.7026513727065744e+300; bool r6585 = r6568 <= r6584; float r6586 = r6566 * r6565; float r6587 = r6563 - r6586; float r6588 = r6575 * r6577; float r6589 = r6576 * r6588; float r6590 = r6587 - r6589; float r6591 = r6550 * r6547; float r6592 = r6552 * r6591; float r6593 = r6548 * r6592; float r6594 = r6557 * r6556; float r6595 = r6593 - r6594; float r6596 = r6554 * r6595; float r6597 = r6590 + r6596; float r6598 = r6575 * r6576; float r6599 = r6577 * r6598; float r6600 = r6563 - r6599; float r6601 = r6566 * r6557; float r6602 = r6552 * r6554; float r6603 = r6548 * r6550; float r6604 = r6602 * r6603; float r6605 = r6601 - r6604; float r6606 = r6547 * r6605; float r6607 = r6600 - r6606; float r6608 = -r6594; float r6609 = r6608 * r6554; float r6610 = r6607 + r6609; float r6611 = r6585 ? r6597 : r6610; float r6612 = r6570 ? r6583 : r6611; return r6612; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r6613 = x; double r6614 = 18.0; double r6615 = r6613 * r6614; double r6616 = y; double r6617 = r6615 * r6616; double r6618 = z; double r6619 = r6617 * r6618; double r6620 = t; double r6621 = r6619 * r6620; double r6622 = a; double r6623 = 4.0; double r6624 = r6622 * r6623; double r6625 = r6624 * r6620; double r6626 = r6621 - r6625; double r6627 = b; double r6628 = c; double r6629 = r6627 * r6628; double r6630 = r6626 + r6629; double r6631 = r6613 * r6623; double r6632 = i; double r6633 = r6631 * r6632; double r6634 = r6630 - r6633; double r6635 = -inf.0; bool r6636 = r6634 <= r6635; double r6637 = r6620 * r6613; double r6638 = r6616 * r6637; double r6639 = r6618 * r6638; double r6640 = r6614 * r6639; double r6641 = 27.0; double r6642 = j; double r6643 = k; double r6644 = r6642 * r6643; double r6645 = r6641 * r6644; double r6646 = r6622 * r6620; double r6647 = r6623 * r6646; double r6648 = r6645 + r6647; double r6649 = r6640 - r6648; double r6650 = 3.7026513727065744e+300; bool r6651 = r6634 <= r6650; double r6652 = r6632 * r6631; double r6653 = r6629 - r6652; double r6654 = r6641 * r6643; double r6655 = r6642 * r6654; double r6656 = r6653 - r6655; double r6657 = r6616 * r6613; double r6658 = r6618 * r6657; double r6659 = r6614 * r6658; double r6660 = r6623 * r6622; double r6661 = r6659 - r6660; double r6662 = r6620 * r6661; double r6663 = r6656 + r6662; double r6664 = r6641 * r6642; double r6665 = r6643 * r6664; double r6666 = r6629 - r6665; double r6667 = r6632 * r6623; double r6668 = r6618 * r6620; double r6669 = r6614 * r6616; double r6670 = r6668 * r6669; double r6671 = r6667 - r6670; double r6672 = r6613 * r6671; double r6673 = r6666 - r6672; double r6674 = -r6660; double r6675 = r6674 * r6620; double r6676 = r6673 + r6675; double r6677 = r6651 ? r6663 : r6676; double r6678 = r6636 ? r6649 : r6677; return r6678; } 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 r6679, r6680, r6681, r6682, r6683, r6684, r6685, r6686, r6687, r6688, r6689, r6690, r6691, r6692, r6693, r6694, r6695, r6696, r6697, r6698, r6699, r6700, r6701, r6702, r6703, r6704, r6705, r6706; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r6679); mpfr_init_set_str(r6680, "18.0", 10, MPFR_RNDN); mpfr_init(r6681); mpfr_init(r6682); mpfr_init(r6683); mpfr_init(r6684); mpfr_init(r6685); mpfr_init(r6686); mpfr_init(r6687); mpfr_init(r6688); mpfr_init_set_str(r6689, "4.0", 10, MPFR_RNDN); mpfr_init(r6690); mpfr_init(r6691); mpfr_init(r6692); mpfr_init(r6693); mpfr_init(r6694); mpfr_init(r6695); mpfr_init(r6696); mpfr_init(r6697); mpfr_init(r6698); mpfr_init(r6699); mpfr_init(r6700); mpfr_init(r6701); mpfr_init_set_str(r6702, "27.0", 10, MPFR_RNDN); mpfr_init(r6703); mpfr_init(r6704); mpfr_init(r6705); mpfr_init(r6706); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r6679, x, MPFR_RNDN); ; mpfr_mul(r6681, r6679, r6680, MPFR_RNDN); mpfr_set_d(r6682, y, MPFR_RNDN); mpfr_mul(r6683, r6681, r6682, MPFR_RNDN); mpfr_set_d(r6684, z, MPFR_RNDN); mpfr_mul(r6685, r6683, r6684, MPFR_RNDN); mpfr_set_d(r6686, t, MPFR_RNDN); mpfr_mul(r6687, r6685, r6686, MPFR_RNDN); mpfr_set_d(r6688, a, MPFR_RNDN); ; mpfr_mul(r6690, r6688, r6689, MPFR_RNDN); mpfr_mul(r6691, r6690, r6686, MPFR_RNDN); mpfr_sub(r6692, r6687, r6691, MPFR_RNDN); mpfr_set_d(r6693, b, MPFR_RNDN); mpfr_set_d(r6694, c, MPFR_RNDN); mpfr_mul(r6695, r6693, r6694, MPFR_RNDN); mpfr_add(r6696, r6692, r6695, MPFR_RNDN); mpfr_mul(r6697, r6679, r6689, MPFR_RNDN); mpfr_set_d(r6698, i, MPFR_RNDN); mpfr_mul(r6699, r6697, r6698, MPFR_RNDN); mpfr_sub(r6700, r6696, r6699, MPFR_RNDN); mpfr_set_d(r6701, j, MPFR_RNDN); ; mpfr_mul(r6703, r6701, r6702, MPFR_RNDN); mpfr_set_d(r6704, k, MPFR_RNDN); mpfr_mul(r6705, r6703, r6704, MPFR_RNDN); mpfr_sub(r6706, r6700, r6705, MPFR_RNDN); return mpfr_get_d(r6706, MPFR_RNDN); } static mpfr_t r6707, r6708, r6709, r6710, r6711, r6712, r6713, r6714, r6715, r6716, r6717, r6718, r6719, r6720, r6721, r6722, r6723, r6724, r6725, r6726, r6727, r6728, r6729, r6730, r6731, r6732, r6733, r6734, r6735, r6736, r6737, r6738, r6739, r6740, r6741, r6742, r6743, r6744, r6745, r6746, r6747, r6748, r6749, r6750, r6751, r6752, r6753, r6754, r6755, r6756, r6757, r6758, r6759, r6760, r6761, r6762, r6763, r6764, r6765, r6766, r6767, r6768, r6769, r6770, r6771, r6772; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r6707); mpfr_init_set_str(r6708, "18.0", 10, MPFR_RNDN); mpfr_init(r6709); mpfr_init(r6710); mpfr_init(r6711); mpfr_init(r6712); mpfr_init(r6713); mpfr_init(r6714); mpfr_init(r6715); mpfr_init(r6716); mpfr_init_set_str(r6717, "4.0", 10, MPFR_RNDN); mpfr_init(r6718); mpfr_init(r6719); mpfr_init(r6720); mpfr_init(r6721); mpfr_init(r6722); mpfr_init(r6723); mpfr_init(r6724); mpfr_init(r6725); mpfr_init(r6726); mpfr_init(r6727); mpfr_init(r6728); mpfr_init_set_str(r6729, "-inf.0", 10, MPFR_RNDN); mpfr_init(r6730); mpfr_init(r6731); mpfr_init(r6732); mpfr_init(r6733); mpfr_init(r6734); mpfr_init_set_str(r6735, "27.0", 10, MPFR_RNDN); mpfr_init(r6736); mpfr_init(r6737); mpfr_init(r6738); mpfr_init(r6739); mpfr_init(r6740); mpfr_init(r6741); mpfr_init(r6742); mpfr_init(r6743); mpfr_init_set_str(r6744, "3.7026513727065744e+300", 10, MPFR_RNDN); mpfr_init(r6745); mpfr_init(r6746); mpfr_init(r6747); mpfr_init(r6748); mpfr_init(r6749); mpfr_init(r6750); mpfr_init(r6751); mpfr_init(r6752); mpfr_init(r6753); mpfr_init(r6754); mpfr_init(r6755); mpfr_init(r6756); mpfr_init(r6757); mpfr_init(r6758); mpfr_init(r6759); mpfr_init(r6760); mpfr_init(r6761); mpfr_init(r6762); mpfr_init(r6763); mpfr_init(r6764); mpfr_init(r6765); mpfr_init(r6766); mpfr_init(r6767); mpfr_init(r6768); mpfr_init(r6769); mpfr_init(r6770); mpfr_init(r6771); mpfr_init(r6772); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r6707, x, MPFR_RNDN); ; mpfr_mul(r6709, r6707, r6708, MPFR_RNDN); mpfr_set_d(r6710, y, MPFR_RNDN); mpfr_mul(r6711, r6709, r6710, MPFR_RNDN); mpfr_set_d(r6712, z, MPFR_RNDN); mpfr_mul(r6713, r6711, r6712, MPFR_RNDN); mpfr_set_d(r6714, t, MPFR_RNDN); mpfr_mul(r6715, r6713, r6714, MPFR_RNDN); mpfr_set_d(r6716, a, MPFR_RNDN); ; mpfr_mul(r6718, r6716, r6717, MPFR_RNDN); mpfr_mul(r6719, r6718, r6714, MPFR_RNDN); mpfr_sub(r6720, r6715, r6719, MPFR_RNDN); mpfr_set_d(r6721, b, MPFR_RNDN); mpfr_set_d(r6722, c, MPFR_RNDN); mpfr_mul(r6723, r6721, r6722, MPFR_RNDN); mpfr_add(r6724, r6720, r6723, MPFR_RNDN); mpfr_mul(r6725, r6707, r6717, MPFR_RNDN); mpfr_set_d(r6726, i, MPFR_RNDN); mpfr_mul(r6727, r6725, r6726, MPFR_RNDN); mpfr_sub(r6728, r6724, r6727, MPFR_RNDN); ; mpfr_set_si(r6730, mpfr_cmp(r6728, r6729) <= 0, MPFR_RNDN); mpfr_mul(r6731, r6714, r6707, MPFR_RNDN); mpfr_mul(r6732, r6710, r6731, MPFR_RNDN); mpfr_mul(r6733, r6712, r6732, MPFR_RNDN); mpfr_mul(r6734, r6708, r6733, MPFR_RNDN); ; mpfr_set_d(r6736, j, MPFR_RNDN); mpfr_set_d(r6737, k, MPFR_RNDN); mpfr_mul(r6738, r6736, r6737, MPFR_RNDN); mpfr_mul(r6739, r6735, r6738, MPFR_RNDN); mpfr_mul(r6740, r6716, r6714, MPFR_RNDN); mpfr_mul(r6741, r6717, r6740, MPFR_RNDN); mpfr_add(r6742, r6739, r6741, MPFR_RNDN); mpfr_sub(r6743, r6734, r6742, MPFR_RNDN); ; mpfr_set_si(r6745, mpfr_cmp(r6728, r6744) <= 0, MPFR_RNDN); mpfr_mul(r6746, r6726, r6725, MPFR_RNDN); mpfr_sub(r6747, r6723, r6746, MPFR_RNDN); mpfr_mul(r6748, r6735, r6737, MPFR_RNDN); mpfr_mul(r6749, r6736, r6748, MPFR_RNDN); mpfr_sub(r6750, r6747, r6749, MPFR_RNDN); mpfr_mul(r6751, r6710, r6707, MPFR_RNDN); mpfr_mul(r6752, r6712, r6751, MPFR_RNDN); mpfr_mul(r6753, r6708, r6752, MPFR_RNDN); mpfr_mul(r6754, r6717, r6716, MPFR_RNDN); mpfr_sub(r6755, r6753, r6754, MPFR_RNDN); mpfr_mul(r6756, r6714, r6755, MPFR_RNDN); mpfr_add(r6757, r6750, r6756, MPFR_RNDN); mpfr_mul(r6758, r6735, r6736, MPFR_RNDN); mpfr_mul(r6759, r6737, r6758, MPFR_RNDN); mpfr_sub(r6760, r6723, r6759, MPFR_RNDN); mpfr_mul(r6761, r6726, r6717, MPFR_RNDN); mpfr_mul(r6762, r6712, r6714, MPFR_RNDN); mpfr_mul(r6763, r6708, r6710, MPFR_RNDN); mpfr_mul(r6764, r6762, r6763, MPFR_RNDN); mpfr_sub(r6765, r6761, r6764, MPFR_RNDN); mpfr_mul(r6766, r6707, r6765, MPFR_RNDN); mpfr_sub(r6767, r6760, r6766, MPFR_RNDN); mpfr_neg(r6768, r6754, MPFR_RNDN); mpfr_mul(r6769, r6768, r6714, MPFR_RNDN); mpfr_add(r6770, r6767, r6769, MPFR_RNDN); if (mpfr_get_si(r6745, MPFR_RNDN)) { mpfr_set(r6771, r6757, MPFR_RNDN); } else { mpfr_set(r6771, r6770, MPFR_RNDN); }; if (mpfr_get_si(r6730, MPFR_RNDN)) { mpfr_set(r6772, r6743, MPFR_RNDN); } else { mpfr_set(r6772, r6771, MPFR_RNDN); }; return mpfr_get_d(r6772, MPFR_RNDN); } static mpfr_t r6773, r6774, r6775, r6776, r6777, r6778, r6779, r6780, r6781, r6782, r6783, r6784, r6785, r6786, r6787, r6788, r6789, r6790, r6791, r6792, r6793, r6794, r6795, r6796, r6797, r6798, r6799, r6800, r6801, r6802, r6803, r6804, r6805, r6806, r6807, r6808, r6809, r6810, r6811, r6812, r6813, r6814, r6815, r6816, r6817, r6818, r6819, r6820, r6821, r6822, r6823, r6824, r6825, r6826, r6827, r6828, r6829, r6830, r6831, r6832, r6833, r6834, r6835, r6836, r6837, r6838; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r6773); mpfr_init_set_str(r6774, "18.0", 10, MPFR_RNDN); mpfr_init(r6775); mpfr_init(r6776); mpfr_init(r6777); mpfr_init(r6778); mpfr_init(r6779); mpfr_init(r6780); mpfr_init(r6781); mpfr_init(r6782); mpfr_init_set_str(r6783, "4.0", 10, MPFR_RNDN); mpfr_init(r6784); mpfr_init(r6785); mpfr_init(r6786); mpfr_init(r6787); mpfr_init(r6788); mpfr_init(r6789); mpfr_init(r6790); mpfr_init(r6791); mpfr_init(r6792); mpfr_init(r6793); mpfr_init(r6794); mpfr_init_set_str(r6795, "-inf.0", 10, MPFR_RNDN); mpfr_init(r6796); mpfr_init(r6797); mpfr_init(r6798); mpfr_init(r6799); mpfr_init(r6800); mpfr_init_set_str(r6801, "27.0", 10, MPFR_RNDN); mpfr_init(r6802); mpfr_init(r6803); mpfr_init(r6804); mpfr_init(r6805); mpfr_init(r6806); mpfr_init(r6807); mpfr_init(r6808); mpfr_init(r6809); mpfr_init_set_str(r6810, "3.7026513727065744e+300", 10, MPFR_RNDN); mpfr_init(r6811); mpfr_init(r6812); mpfr_init(r6813); mpfr_init(r6814); mpfr_init(r6815); mpfr_init(r6816); mpfr_init(r6817); mpfr_init(r6818); mpfr_init(r6819); mpfr_init(r6820); mpfr_init(r6821); mpfr_init(r6822); mpfr_init(r6823); mpfr_init(r6824); mpfr_init(r6825); mpfr_init(r6826); mpfr_init(r6827); mpfr_init(r6828); mpfr_init(r6829); mpfr_init(r6830); mpfr_init(r6831); mpfr_init(r6832); mpfr_init(r6833); mpfr_init(r6834); mpfr_init(r6835); mpfr_init(r6836); mpfr_init(r6837); mpfr_init(r6838); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r6773, x, MPFR_RNDN); ; mpfr_mul(r6775, r6773, r6774, MPFR_RNDN); mpfr_set_d(r6776, y, MPFR_RNDN); mpfr_mul(r6777, r6775, r6776, MPFR_RNDN); mpfr_set_d(r6778, z, MPFR_RNDN); mpfr_mul(r6779, r6777, r6778, MPFR_RNDN); mpfr_set_d(r6780, t, MPFR_RNDN); mpfr_mul(r6781, r6779, r6780, MPFR_RNDN); mpfr_set_d(r6782, a, MPFR_RNDN); ; mpfr_mul(r6784, r6782, r6783, MPFR_RNDN); mpfr_mul(r6785, r6784, r6780, MPFR_RNDN); mpfr_sub(r6786, r6781, r6785, MPFR_RNDN); mpfr_set_d(r6787, b, MPFR_RNDN); mpfr_set_d(r6788, c, MPFR_RNDN); mpfr_mul(r6789, r6787, r6788, MPFR_RNDN); mpfr_add(r6790, r6786, r6789, MPFR_RNDN); mpfr_mul(r6791, r6773, r6783, MPFR_RNDN); mpfr_set_d(r6792, i, MPFR_RNDN); mpfr_mul(r6793, r6791, r6792, MPFR_RNDN); mpfr_sub(r6794, r6790, r6793, MPFR_RNDN); ; mpfr_set_si(r6796, mpfr_cmp(r6794, r6795) <= 0, MPFR_RNDN); mpfr_mul(r6797, r6780, r6773, MPFR_RNDN); mpfr_mul(r6798, r6776, r6797, MPFR_RNDN); mpfr_mul(r6799, r6778, r6798, MPFR_RNDN); mpfr_mul(r6800, r6774, r6799, MPFR_RNDN); ; mpfr_set_d(r6802, j, MPFR_RNDN); mpfr_set_d(r6803, k, MPFR_RNDN); mpfr_mul(r6804, r6802, r6803, MPFR_RNDN); mpfr_mul(r6805, r6801, r6804, MPFR_RNDN); mpfr_mul(r6806, r6782, r6780, MPFR_RNDN); mpfr_mul(r6807, r6783, r6806, MPFR_RNDN); mpfr_add(r6808, r6805, r6807, MPFR_RNDN); mpfr_sub(r6809, r6800, r6808, MPFR_RNDN); ; mpfr_set_si(r6811, mpfr_cmp(r6794, r6810) <= 0, MPFR_RNDN); mpfr_mul(r6812, r6792, r6791, MPFR_RNDN); mpfr_sub(r6813, r6789, r6812, MPFR_RNDN); mpfr_mul(r6814, r6801, r6803, MPFR_RNDN); mpfr_mul(r6815, r6802, r6814, MPFR_RNDN); mpfr_sub(r6816, r6813, r6815, MPFR_RNDN); mpfr_mul(r6817, r6776, r6773, MPFR_RNDN); mpfr_mul(r6818, r6778, r6817, MPFR_RNDN); mpfr_mul(r6819, r6774, r6818, MPFR_RNDN); mpfr_mul(r6820, r6783, r6782, MPFR_RNDN); mpfr_sub(r6821, r6819, r6820, MPFR_RNDN); mpfr_mul(r6822, r6780, r6821, MPFR_RNDN); mpfr_add(r6823, r6816, r6822, MPFR_RNDN); mpfr_mul(r6824, r6801, r6802, MPFR_RNDN); mpfr_mul(r6825, r6803, r6824, MPFR_RNDN); mpfr_sub(r6826, r6789, r6825, MPFR_RNDN); mpfr_mul(r6827, r6792, r6783, MPFR_RNDN); mpfr_mul(r6828, r6778, r6780, MPFR_RNDN); mpfr_mul(r6829, r6774, r6776, MPFR_RNDN); mpfr_mul(r6830, r6828, r6829, MPFR_RNDN); mpfr_sub(r6831, r6827, r6830, MPFR_RNDN); mpfr_mul(r6832, r6773, r6831, MPFR_RNDN); mpfr_sub(r6833, r6826, r6832, MPFR_RNDN); mpfr_neg(r6834, r6820, MPFR_RNDN); mpfr_mul(r6835, r6834, r6780, MPFR_RNDN); mpfr_add(r6836, r6833, r6835, MPFR_RNDN); if (mpfr_get_si(r6811, MPFR_RNDN)) { mpfr_set(r6837, r6823, MPFR_RNDN); } else { mpfr_set(r6837, r6836, MPFR_RNDN); }; if (mpfr_get_si(r6796, MPFR_RNDN)) { mpfr_set(r6838, r6809, MPFR_RNDN); } else { mpfr_set(r6838, r6837, MPFR_RNDN); }; return mpfr_get_d(r6838, MPFR_RNDN); }