#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 r11499 = x; float r11500 = 18.0; float r11501 = r11499 * r11500; float r11502 = y; float r11503 = r11501 * r11502; float r11504 = z; float r11505 = r11503 * r11504; float r11506 = t; float r11507 = r11505 * r11506; float r11508 = a; float r11509 = 4.0; float r11510 = r11508 * r11509; float r11511 = r11510 * r11506; float r11512 = r11507 - r11511; float r11513 = b; float r11514 = c; float r11515 = r11513 * r11514; float r11516 = r11512 + r11515; float r11517 = r11499 * r11509; float r11518 = i; float r11519 = r11517 * r11518; float r11520 = r11516 - r11519; float r11521 = j; float r11522 = 27.0; float r11523 = r11521 * r11522; float r11524 = k; float r11525 = r11523 * r11524; float r11526 = r11520 - r11525; return r11526; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r11527 = x; double r11528 = 18.0; double r11529 = r11527 * r11528; double r11530 = y; double r11531 = r11529 * r11530; double r11532 = z; double r11533 = r11531 * r11532; double r11534 = t; double r11535 = r11533 * r11534; double r11536 = a; double r11537 = 4.0; double r11538 = r11536 * r11537; double r11539 = r11538 * r11534; double r11540 = r11535 - r11539; double r11541 = b; double r11542 = c; double r11543 = r11541 * r11542; double r11544 = r11540 + r11543; double r11545 = r11527 * r11537; double r11546 = i; double r11547 = r11545 * r11546; double r11548 = r11544 - r11547; double r11549 = j; double r11550 = 27.0; double r11551 = r11549 * r11550; double r11552 = k; double r11553 = r11551 * r11552; double r11554 = r11548 - r11553; return r11554; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r11555 = z; float r11556 = y; float r11557 = r11555 * r11556; float r11558 = 18.0; float r11559 = r11557 * r11558; float r11560 = -3.239769753391218e+287; bool r11561 = r11559 <= r11560; float r11562 = t; float r11563 = r11555 * r11562; float r11564 = x; float r11565 = r11558 * r11564; float r11566 = r11565 * r11556; float r11567 = r11563 * r11566; float r11568 = 4.0; float r11569 = a; float r11570 = r11568 * r11569; float r11571 = r11562 * r11570; float r11572 = r11567 - r11571; float r11573 = b; float r11574 = c; float r11575 = r11573 * r11574; float r11576 = r11572 + r11575; float r11577 = i; float r11578 = r11568 * r11564; float r11579 = r11577 * r11578; float r11580 = r11576 - r11579; float r11581 = k; float r11582 = 27.0; float r11583 = j; float r11584 = r11582 * r11583; float r11585 = r11581 * r11584; float r11586 = r11580 - r11585; float r11587 = -1.324853790717036e-196; bool r11588 = r11559 <= r11587; float r11589 = 1.2766226871705525e-285; bool r11590 = r11559 <= r11589; float r11591 = 6.301398134605551e+265; bool r11592 = r11559 <= r11591; float r11593 = !r11592; bool r11594 = r11590 || r11593; float r11595 = !r11594; bool r11596 = r11588 || r11595; float r11597 = r11564 * r11562; float r11598 = r11559 * r11597; float r11599 = r11575 + r11598; float r11600 = r11581 * r11582; float r11601 = r11583 * r11600; float r11602 = r11577 * r11564; float r11603 = r11562 * r11569; float r11604 = r11602 + r11603; float r11605 = r11604 * r11568; float r11606 = r11601 + r11605; float r11607 = r11599 - r11606; float r11608 = r11596 ? r11607 : r11586; float r11609 = r11561 ? r11586 : r11608; return r11609; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r11610 = z; double r11611 = y; double r11612 = r11610 * r11611; double r11613 = 18.0; double r11614 = r11612 * r11613; double r11615 = -3.239769753391218e+287; bool r11616 = r11614 <= r11615; double r11617 = t; double r11618 = r11610 * r11617; double r11619 = x; double r11620 = r11613 * r11619; double r11621 = r11620 * r11611; double r11622 = r11618 * r11621; double r11623 = 4.0; double r11624 = a; double r11625 = r11623 * r11624; double r11626 = r11617 * r11625; double r11627 = r11622 - r11626; double r11628 = b; double r11629 = c; double r11630 = r11628 * r11629; double r11631 = r11627 + r11630; double r11632 = i; double r11633 = r11623 * r11619; double r11634 = r11632 * r11633; double r11635 = r11631 - r11634; double r11636 = k; double r11637 = 27.0; double r11638 = j; double r11639 = r11637 * r11638; double r11640 = r11636 * r11639; double r11641 = r11635 - r11640; double r11642 = -1.324853790717036e-196; bool r11643 = r11614 <= r11642; double r11644 = 1.2766226871705525e-285; bool r11645 = r11614 <= r11644; double r11646 = 6.301398134605551e+265; bool r11647 = r11614 <= r11646; double r11648 = !r11647; bool r11649 = r11645 || r11648; double r11650 = !r11649; bool r11651 = r11643 || r11650; double r11652 = r11619 * r11617; double r11653 = r11614 * r11652; double r11654 = r11630 + r11653; double r11655 = r11636 * r11637; double r11656 = r11638 * r11655; double r11657 = r11632 * r11619; double r11658 = r11617 * r11624; double r11659 = r11657 + r11658; double r11660 = r11659 * r11623; double r11661 = r11656 + r11660; double r11662 = r11654 - r11661; double r11663 = r11651 ? r11662 : r11641; double r11664 = r11616 ? r11641 : r11663; return r11664; } 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 r11665, r11666, r11667, r11668, r11669, r11670, r11671, r11672, r11673, r11674, r11675, r11676, r11677, r11678, r11679, r11680, r11681, r11682, r11683, r11684, r11685, r11686, r11687, r11688, r11689, r11690, r11691, r11692; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r11665); mpfr_init_set_str(r11666, "18.0", 10, MPFR_RNDN); mpfr_init(r11667); mpfr_init(r11668); mpfr_init(r11669); mpfr_init(r11670); mpfr_init(r11671); mpfr_init(r11672); mpfr_init(r11673); mpfr_init(r11674); mpfr_init_set_str(r11675, "4.0", 10, MPFR_RNDN); mpfr_init(r11676); mpfr_init(r11677); mpfr_init(r11678); mpfr_init(r11679); mpfr_init(r11680); mpfr_init(r11681); mpfr_init(r11682); mpfr_init(r11683); mpfr_init(r11684); mpfr_init(r11685); mpfr_init(r11686); mpfr_init(r11687); mpfr_init_set_str(r11688, "27.0", 10, MPFR_RNDN); mpfr_init(r11689); mpfr_init(r11690); mpfr_init(r11691); mpfr_init(r11692); } 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(r11665, x, MPFR_RNDN); ; mpfr_mul(r11667, r11665, r11666, MPFR_RNDN); mpfr_set_d(r11668, y, MPFR_RNDN); mpfr_mul(r11669, r11667, r11668, MPFR_RNDN); mpfr_set_d(r11670, z, MPFR_RNDN); mpfr_mul(r11671, r11669, r11670, MPFR_RNDN); mpfr_set_d(r11672, t, MPFR_RNDN); mpfr_mul(r11673, r11671, r11672, MPFR_RNDN); mpfr_set_d(r11674, a, MPFR_RNDN); ; mpfr_mul(r11676, r11674, r11675, MPFR_RNDN); mpfr_mul(r11677, r11676, r11672, MPFR_RNDN); mpfr_sub(r11678, r11673, r11677, MPFR_RNDN); mpfr_set_d(r11679, b, MPFR_RNDN); mpfr_set_d(r11680, c, MPFR_RNDN); mpfr_mul(r11681, r11679, r11680, MPFR_RNDN); mpfr_add(r11682, r11678, r11681, MPFR_RNDN); mpfr_mul(r11683, r11665, r11675, MPFR_RNDN); mpfr_set_d(r11684, i, MPFR_RNDN); mpfr_mul(r11685, r11683, r11684, MPFR_RNDN); mpfr_sub(r11686, r11682, r11685, MPFR_RNDN); mpfr_set_d(r11687, j, MPFR_RNDN); ; mpfr_mul(r11689, r11687, r11688, MPFR_RNDN); mpfr_set_d(r11690, k, MPFR_RNDN); mpfr_mul(r11691, r11689, r11690, MPFR_RNDN); mpfr_sub(r11692, r11686, r11691, MPFR_RNDN); return mpfr_get_d(r11692, MPFR_RNDN); } static mpfr_t r11693, r11694, r11695, r11696, r11697, r11698, r11699, r11700, r11701, r11702, r11703, r11704, r11705, r11706, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r11693); mpfr_init(r11694); mpfr_init(r11695); mpfr_init_set_str(r11696, "18.0", 10, MPFR_RNDN); mpfr_init(r11697); mpfr_init_set_str(r11698, "-3.239769753391218e+287", 10, MPFR_RNDN); mpfr_init(r11699); mpfr_init(r11700); mpfr_init(r11701); mpfr_init(r11702); mpfr_init(r11703); mpfr_init(r11704); mpfr_init(r11705); mpfr_init_set_str(r11706, "4.0", 10, MPFR_RNDN); mpfr_init(r11707); mpfr_init(r11708); mpfr_init(r11709); mpfr_init(r11710); mpfr_init(r11711); mpfr_init(r11712); mpfr_init(r11713); mpfr_init(r11714); mpfr_init(r11715); mpfr_init(r11716); mpfr_init(r11717); mpfr_init(r11718); mpfr_init(r11719); mpfr_init_set_str(r11720, "27.0", 10, MPFR_RNDN); mpfr_init(r11721); mpfr_init(r11722); mpfr_init(r11723); mpfr_init(r11724); mpfr_init_set_str(r11725, "-1.324853790717036e-196", 10, MPFR_RNDN); mpfr_init(r11726); mpfr_init_set_str(r11727, "1.2766226871705525e-285", 10, MPFR_RNDN); mpfr_init(r11728); mpfr_init_set_str(r11729, "6.301398134605551e+265", 10, MPFR_RNDN); mpfr_init(r11730); mpfr_init(r11731); mpfr_init(r11732); mpfr_init(r11733); mpfr_init(r11734); mpfr_init(r11735); mpfr_init(r11736); mpfr_init(r11737); mpfr_init(r11738); mpfr_init(r11739); mpfr_init(r11740); mpfr_init(r11741); mpfr_init(r11742); mpfr_init(r11743); mpfr_init(r11744); mpfr_init(r11745); mpfr_init(r11746); mpfr_init(r11747); } 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(r11693, z, MPFR_RNDN); mpfr_set_d(r11694, y, MPFR_RNDN); mpfr_mul(r11695, r11693, r11694, MPFR_RNDN); ; mpfr_mul(r11697, r11695, r11696, MPFR_RNDN); ; mpfr_set_si(r11699, mpfr_cmp(r11697, r11698) <= 0, MPFR_RNDN); mpfr_set_d(r11700, t, MPFR_RNDN); mpfr_mul(r11701, r11693, r11700, MPFR_RNDN); mpfr_set_d(r11702, x, MPFR_RNDN); mpfr_mul(r11703, r11696, r11702, MPFR_RNDN); mpfr_mul(r11704, r11703, r11694, MPFR_RNDN); mpfr_mul(r11705, r11701, r11704, MPFR_RNDN); ; mpfr_set_d(r11707, a, MPFR_RNDN); mpfr_mul(r11708, r11706, r11707, MPFR_RNDN); mpfr_mul(r11709, r11700, r11708, MPFR_RNDN); mpfr_sub(r11710, r11705, r11709, MPFR_RNDN); mpfr_set_d(r11711, b, MPFR_RNDN); mpfr_set_d(r11712, c, MPFR_RNDN); mpfr_mul(r11713, r11711, r11712, MPFR_RNDN); mpfr_add(r11714, r11710, r11713, MPFR_RNDN); mpfr_set_d(r11715, i, MPFR_RNDN); mpfr_mul(r11716, r11706, r11702, MPFR_RNDN); mpfr_mul(r11717, r11715, r11716, MPFR_RNDN); mpfr_sub(r11718, r11714, r11717, MPFR_RNDN); mpfr_set_d(r11719, k, MPFR_RNDN); ; mpfr_set_d(r11721, j, MPFR_RNDN); mpfr_mul(r11722, r11720, r11721, MPFR_RNDN); mpfr_mul(r11723, r11719, r11722, MPFR_RNDN); mpfr_sub(r11724, r11718, r11723, MPFR_RNDN); ; mpfr_set_si(r11726, mpfr_cmp(r11697, r11725) <= 0, MPFR_RNDN); ; mpfr_set_si(r11728, mpfr_cmp(r11697, r11727) <= 0, MPFR_RNDN); ; mpfr_set_si(r11730, mpfr_cmp(r11697, r11729) <= 0, MPFR_RNDN); mpfr_set_si(r11731, !mpfr_get_si(r11730, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11732, mpfr_get_si(r11728, MPFR_RNDN) || mpfr_get_si(r11731, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11733, !mpfr_get_si(r11732, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11734, mpfr_get_si(r11726, MPFR_RNDN) || mpfr_get_si(r11733, MPFR_RNDN), MPFR_RNDN); mpfr_mul(r11735, r11702, r11700, MPFR_RNDN); mpfr_mul(r11736, r11697, r11735, MPFR_RNDN); mpfr_add(r11737, r11713, r11736, MPFR_RNDN); mpfr_mul(r11738, r11719, r11720, MPFR_RNDN); mpfr_mul(r11739, r11721, r11738, MPFR_RNDN); mpfr_mul(r11740, r11715, r11702, MPFR_RNDN); mpfr_mul(r11741, r11700, r11707, MPFR_RNDN); mpfr_add(r11742, r11740, r11741, MPFR_RNDN); mpfr_mul(r11743, r11742, r11706, MPFR_RNDN); mpfr_add(r11744, r11739, r11743, MPFR_RNDN); mpfr_sub(r11745, r11737, r11744, MPFR_RNDN); if (mpfr_get_si(r11734, MPFR_RNDN)) { mpfr_set(r11746, r11745, MPFR_RNDN); } else { mpfr_set(r11746, r11724, MPFR_RNDN); }; if (mpfr_get_si(r11699, MPFR_RNDN)) { mpfr_set(r11747, r11724, MPFR_RNDN); } else { mpfr_set(r11747, r11746, MPFR_RNDN); }; return mpfr_get_d(r11747, MPFR_RNDN); } static mpfr_t r11748, r11749, r11750, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11748); mpfr_init(r11749); mpfr_init(r11750); mpfr_init_set_str(r11751, "18.0", 10, MPFR_RNDN); mpfr_init(r11752); mpfr_init_set_str(r11753, "-3.239769753391218e+287", 10, MPFR_RNDN); mpfr_init(r11754); mpfr_init(r11755); mpfr_init(r11756); mpfr_init(r11757); mpfr_init(r11758); mpfr_init(r11759); mpfr_init(r11760); mpfr_init_set_str(r11761, "4.0", 10, MPFR_RNDN); mpfr_init(r11762); mpfr_init(r11763); mpfr_init(r11764); mpfr_init(r11765); mpfr_init(r11766); mpfr_init(r11767); mpfr_init(r11768); mpfr_init(r11769); mpfr_init(r11770); mpfr_init(r11771); mpfr_init(r11772); mpfr_init(r11773); mpfr_init(r11774); mpfr_init_set_str(r11775, "27.0", 10, MPFR_RNDN); mpfr_init(r11776); mpfr_init(r11777); mpfr_init(r11778); mpfr_init(r11779); mpfr_init_set_str(r11780, "-1.324853790717036e-196", 10, MPFR_RNDN); mpfr_init(r11781); mpfr_init_set_str(r11782, "1.2766226871705525e-285", 10, MPFR_RNDN); mpfr_init(r11783); mpfr_init_set_str(r11784, "6.301398134605551e+265", 10, MPFR_RNDN); mpfr_init(r11785); mpfr_init(r11786); mpfr_init(r11787); mpfr_init(r11788); mpfr_init(r11789); mpfr_init(r11790); mpfr_init(r11791); mpfr_init(r11792); mpfr_init(r11793); mpfr_init(r11794); mpfr_init(r11795); mpfr_init(r11796); mpfr_init(r11797); mpfr_init(r11798); mpfr_init(r11799); mpfr_init(r11800); mpfr_init(r11801); mpfr_init(r11802); } 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(r11748, z, MPFR_RNDN); mpfr_set_d(r11749, y, MPFR_RNDN); mpfr_mul(r11750, r11748, r11749, MPFR_RNDN); ; mpfr_mul(r11752, r11750, r11751, MPFR_RNDN); ; mpfr_set_si(r11754, mpfr_cmp(r11752, r11753) <= 0, MPFR_RNDN); mpfr_set_d(r11755, t, MPFR_RNDN); mpfr_mul(r11756, r11748, r11755, MPFR_RNDN); mpfr_set_d(r11757, x, MPFR_RNDN); mpfr_mul(r11758, r11751, r11757, MPFR_RNDN); mpfr_mul(r11759, r11758, r11749, MPFR_RNDN); mpfr_mul(r11760, r11756, r11759, MPFR_RNDN); ; mpfr_set_d(r11762, a, MPFR_RNDN); mpfr_mul(r11763, r11761, r11762, MPFR_RNDN); mpfr_mul(r11764, r11755, r11763, MPFR_RNDN); mpfr_sub(r11765, r11760, r11764, MPFR_RNDN); mpfr_set_d(r11766, b, MPFR_RNDN); mpfr_set_d(r11767, c, MPFR_RNDN); mpfr_mul(r11768, r11766, r11767, MPFR_RNDN); mpfr_add(r11769, r11765, r11768, MPFR_RNDN); mpfr_set_d(r11770, i, MPFR_RNDN); mpfr_mul(r11771, r11761, r11757, MPFR_RNDN); mpfr_mul(r11772, r11770, r11771, MPFR_RNDN); mpfr_sub(r11773, r11769, r11772, MPFR_RNDN); mpfr_set_d(r11774, k, MPFR_RNDN); ; mpfr_set_d(r11776, j, MPFR_RNDN); mpfr_mul(r11777, r11775, r11776, MPFR_RNDN); mpfr_mul(r11778, r11774, r11777, MPFR_RNDN); mpfr_sub(r11779, r11773, r11778, MPFR_RNDN); ; mpfr_set_si(r11781, mpfr_cmp(r11752, r11780) <= 0, MPFR_RNDN); ; mpfr_set_si(r11783, mpfr_cmp(r11752, r11782) <= 0, MPFR_RNDN); ; mpfr_set_si(r11785, mpfr_cmp(r11752, r11784) <= 0, MPFR_RNDN); mpfr_set_si(r11786, !mpfr_get_si(r11785, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11787, mpfr_get_si(r11783, MPFR_RNDN) || mpfr_get_si(r11786, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11788, !mpfr_get_si(r11787, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r11789, mpfr_get_si(r11781, MPFR_RNDN) || mpfr_get_si(r11788, MPFR_RNDN), MPFR_RNDN); mpfr_mul(r11790, r11757, r11755, MPFR_RNDN); mpfr_mul(r11791, r11752, r11790, MPFR_RNDN); mpfr_add(r11792, r11768, r11791, MPFR_RNDN); mpfr_mul(r11793, r11774, r11775, MPFR_RNDN); mpfr_mul(r11794, r11776, r11793, MPFR_RNDN); mpfr_mul(r11795, r11770, r11757, MPFR_RNDN); mpfr_mul(r11796, r11755, r11762, MPFR_RNDN); mpfr_add(r11797, r11795, r11796, MPFR_RNDN); mpfr_mul(r11798, r11797, r11761, MPFR_RNDN); mpfr_add(r11799, r11794, r11798, MPFR_RNDN); mpfr_sub(r11800, r11792, r11799, MPFR_RNDN); if (mpfr_get_si(r11789, MPFR_RNDN)) { mpfr_set(r11801, r11800, MPFR_RNDN); } else { mpfr_set(r11801, r11779, MPFR_RNDN); }; if (mpfr_get_si(r11754, MPFR_RNDN)) { mpfr_set(r11802, r11779, MPFR_RNDN); } else { mpfr_set(r11802, r11801, MPFR_RNDN); }; return mpfr_get_d(r11802, MPFR_RNDN); }