#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 r10560 = x; float r10561 = 18.0; float r10562 = r10560 * r10561; float r10563 = y; float r10564 = r10562 * r10563; float r10565 = z; float r10566 = r10564 * r10565; float r10567 = t; float r10568 = r10566 * r10567; float r10569 = a; float r10570 = 4.0; float r10571 = r10569 * r10570; float r10572 = r10571 * r10567; float r10573 = r10568 - r10572; float r10574 = b; float r10575 = c; float r10576 = r10574 * r10575; float r10577 = r10573 + r10576; float r10578 = r10560 * r10570; float r10579 = i; float r10580 = r10578 * r10579; float r10581 = r10577 - r10580; float r10582 = j; float r10583 = 27.0; float r10584 = r10582 * r10583; float r10585 = k; float r10586 = r10584 * r10585; float r10587 = r10581 - r10586; return r10587; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r10588 = x; double r10589 = 18.0; double r10590 = r10588 * r10589; double r10591 = y; double r10592 = r10590 * r10591; double r10593 = z; double r10594 = r10592 * r10593; double r10595 = t; double r10596 = r10594 * r10595; double r10597 = a; double r10598 = 4.0; double r10599 = r10597 * r10598; double r10600 = r10599 * r10595; double r10601 = r10596 - r10600; double r10602 = b; double r10603 = c; double r10604 = r10602 * r10603; double r10605 = r10601 + r10604; double r10606 = r10588 * r10598; double r10607 = i; double r10608 = r10606 * r10607; double r10609 = r10605 - r10608; double r10610 = j; double r10611 = 27.0; double r10612 = r10610 * r10611; double r10613 = k; double r10614 = r10612 * r10613; double r10615 = r10609 - r10614; return r10615; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r10616 = c; float r10617 = b; float r10618 = r10616 * r10617; float r10619 = 27.0; float r10620 = j; float r10621 = k; float r10622 = r10620 * r10621; float r10623 = r10619 * r10622; float r10624 = r10618 - r10623; float r10625 = x; float r10626 = i; float r10627 = 4.0; float r10628 = r10626 * r10627; float r10629 = z; float r10630 = t; float r10631 = r10629 * r10630; float r10632 = y; float r10633 = 18.0; float r10634 = r10632 * r10633; float r10635 = r10631 * r10634; float r10636 = r10628 - r10635; float r10637 = r10625 * r10636; float r10638 = cbrt(r10637); float r10639 = r10638 * r10638; float r10640 = r10639 * r10638; float r10641 = r10624 - r10640; float r10642 = a; float r10643 = r10642 * r10627; float r10644 = -r10643; float r10645 = r10630 * r10644; float r10646 = r10641 + r10645; float r10647 = -1.103406987904582e+307; bool r10648 = r10646 <= r10647; float r10649 = r10625 * r10633; float r10650 = r10649 * r10632; float r10651 = r10650 * r10629; float r10652 = r10651 * r10630; float r10653 = r10643 * r10630; float r10654 = r10652 - r10653; float r10655 = r10617 * r10616; float r10656 = r10654 + r10655; float r10657 = r10625 * r10627; float r10658 = r10657 * r10626; float r10659 = r10656 - r10658; float r10660 = r10620 * r10619; float r10661 = r10660 * r10621; float r10662 = r10659 - r10661; float r10663 = 1.3127925574492975e+300; bool r10664 = r10646 <= r10663; float r10665 = r10631 * r10632; float r10666 = r10665 * r10633; float r10667 = r10628 - r10666; float r10668 = r10625 * r10667; float r10669 = r10624 - r10668; float r10670 = r10669 + r10645; float r10671 = r10664 ? r10670 : r10662; float r10672 = r10648 ? r10662 : r10671; return r10672; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r10673 = c; double r10674 = b; double r10675 = r10673 * r10674; double r10676 = 27.0; double r10677 = j; double r10678 = k; double r10679 = r10677 * r10678; double r10680 = r10676 * r10679; double r10681 = r10675 - r10680; double r10682 = x; double r10683 = i; double r10684 = 4.0; double r10685 = r10683 * r10684; double r10686 = z; double r10687 = t; double r10688 = r10686 * r10687; double r10689 = y; double r10690 = 18.0; double r10691 = r10689 * r10690; double r10692 = r10688 * r10691; double r10693 = r10685 - r10692; double r10694 = r10682 * r10693; double r10695 = cbrt(r10694); double r10696 = r10695 * r10695; double r10697 = r10696 * r10695; double r10698 = r10681 - r10697; double r10699 = a; double r10700 = r10699 * r10684; double r10701 = -r10700; double r10702 = r10687 * r10701; double r10703 = r10698 + r10702; double r10704 = -1.103406987904582e+307; bool r10705 = r10703 <= r10704; double r10706 = r10682 * r10690; double r10707 = r10706 * r10689; double r10708 = r10707 * r10686; double r10709 = r10708 * r10687; double r10710 = r10700 * r10687; double r10711 = r10709 - r10710; double r10712 = r10674 * r10673; double r10713 = r10711 + r10712; double r10714 = r10682 * r10684; double r10715 = r10714 * r10683; double r10716 = r10713 - r10715; double r10717 = r10677 * r10676; double r10718 = r10717 * r10678; double r10719 = r10716 - r10718; double r10720 = 1.3127925574492975e+300; bool r10721 = r10703 <= r10720; double r10722 = r10688 * r10689; double r10723 = r10722 * r10690; double r10724 = r10685 - r10723; double r10725 = r10682 * r10724; double r10726 = r10681 - r10725; double r10727 = r10726 + r10702; double r10728 = r10721 ? r10727 : r10719; double r10729 = r10705 ? r10719 : r10728; return r10729; } 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 r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739, r10740, r10741, r10742, r10743, r10744, r10745, r10746, r10747, r10748, r10749, r10750, r10751, r10752, r10753, r10754, r10755, r10756, r10757; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10730); mpfr_init_set_str(r10731, "18.0", 10, MPFR_RNDN); mpfr_init(r10732); mpfr_init(r10733); mpfr_init(r10734); mpfr_init(r10735); mpfr_init(r10736); mpfr_init(r10737); mpfr_init(r10738); mpfr_init(r10739); mpfr_init_set_str(r10740, "4.0", 10, MPFR_RNDN); mpfr_init(r10741); mpfr_init(r10742); mpfr_init(r10743); mpfr_init(r10744); mpfr_init(r10745); mpfr_init(r10746); mpfr_init(r10747); mpfr_init(r10748); mpfr_init(r10749); mpfr_init(r10750); mpfr_init(r10751); mpfr_init(r10752); mpfr_init_set_str(r10753, "27.0", 10, MPFR_RNDN); mpfr_init(r10754); mpfr_init(r10755); mpfr_init(r10756); mpfr_init(r10757); } 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(r10730, x, MPFR_RNDN); ; mpfr_mul(r10732, r10730, r10731, MPFR_RNDN); mpfr_set_d(r10733, y, MPFR_RNDN); mpfr_mul(r10734, r10732, r10733, MPFR_RNDN); mpfr_set_d(r10735, z, MPFR_RNDN); mpfr_mul(r10736, r10734, r10735, MPFR_RNDN); mpfr_set_d(r10737, t, MPFR_RNDN); mpfr_mul(r10738, r10736, r10737, MPFR_RNDN); mpfr_set_d(r10739, a, MPFR_RNDN); ; mpfr_mul(r10741, r10739, r10740, MPFR_RNDN); mpfr_mul(r10742, r10741, r10737, MPFR_RNDN); mpfr_sub(r10743, r10738, r10742, MPFR_RNDN); mpfr_set_d(r10744, b, MPFR_RNDN); mpfr_set_d(r10745, c, MPFR_RNDN); mpfr_mul(r10746, r10744, r10745, MPFR_RNDN); mpfr_add(r10747, r10743, r10746, MPFR_RNDN); mpfr_mul(r10748, r10730, r10740, MPFR_RNDN); mpfr_set_d(r10749, i, MPFR_RNDN); mpfr_mul(r10750, r10748, r10749, MPFR_RNDN); mpfr_sub(r10751, r10747, r10750, MPFR_RNDN); mpfr_set_d(r10752, j, MPFR_RNDN); ; mpfr_mul(r10754, r10752, r10753, MPFR_RNDN); mpfr_set_d(r10755, k, MPFR_RNDN); mpfr_mul(r10756, r10754, r10755, MPFR_RNDN); mpfr_sub(r10757, r10751, r10756, MPFR_RNDN); return mpfr_get_d(r10757, MPFR_RNDN); } static mpfr_t r10758, r10759, r10760, r10761, r10762, r10763, r10764, r10765, r10766, r10767, r10768, r10769, r10770, r10771, r10772, r10773, r10774, r10775, r10776, r10777, r10778, r10779, r10780, r10781, r10782, r10783, r10784, r10785, r10786, r10787, r10788, r10789, r10790, r10791, r10792, r10793, r10794, r10795, r10796, r10797, r10798, r10799, r10800, r10801, r10802, r10803, r10804, r10805, r10806, r10807, r10808, r10809, r10810, r10811, r10812, r10813, r10814; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10758); mpfr_init(r10759); mpfr_init(r10760); mpfr_init_set_str(r10761, "27.0", 10, MPFR_RNDN); mpfr_init(r10762); mpfr_init(r10763); mpfr_init(r10764); mpfr_init(r10765); mpfr_init(r10766); mpfr_init(r10767); mpfr_init(r10768); mpfr_init_set_str(r10769, "4.0", 10, MPFR_RNDN); mpfr_init(r10770); mpfr_init(r10771); mpfr_init(r10772); mpfr_init(r10773); mpfr_init(r10774); mpfr_init_set_str(r10775, "18.0", 10, MPFR_RNDN); mpfr_init(r10776); mpfr_init(r10777); mpfr_init(r10778); mpfr_init(r10779); mpfr_init(r10780); mpfr_init(r10781); mpfr_init(r10782); mpfr_init(r10783); mpfr_init(r10784); mpfr_init(r10785); mpfr_init(r10786); mpfr_init(r10787); mpfr_init(r10788); mpfr_init_set_str(r10789, "-1.103406987904582e+307", 10, MPFR_RNDN); mpfr_init(r10790); mpfr_init(r10791); mpfr_init(r10792); mpfr_init(r10793); mpfr_init(r10794); mpfr_init(r10795); mpfr_init(r10796); mpfr_init(r10797); mpfr_init(r10798); mpfr_init(r10799); mpfr_init(r10800); mpfr_init(r10801); mpfr_init(r10802); mpfr_init(r10803); mpfr_init(r10804); mpfr_init_set_str(r10805, "1.3127925574492975e+300", 10, MPFR_RNDN); mpfr_init(r10806); mpfr_init(r10807); mpfr_init(r10808); mpfr_init(r10809); mpfr_init(r10810); mpfr_init(r10811); mpfr_init(r10812); mpfr_init(r10813); mpfr_init(r10814); } 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(r10758, c, MPFR_RNDN); mpfr_set_d(r10759, b, MPFR_RNDN); mpfr_mul(r10760, r10758, r10759, MPFR_RNDN); ; mpfr_set_d(r10762, j, MPFR_RNDN); mpfr_set_d(r10763, k, MPFR_RNDN); mpfr_mul(r10764, r10762, r10763, MPFR_RNDN); mpfr_mul(r10765, r10761, r10764, MPFR_RNDN); mpfr_sub(r10766, r10760, r10765, MPFR_RNDN); mpfr_set_d(r10767, x, MPFR_RNDN); mpfr_set_d(r10768, i, MPFR_RNDN); ; mpfr_mul(r10770, r10768, r10769, MPFR_RNDN); mpfr_set_d(r10771, z, MPFR_RNDN); mpfr_set_d(r10772, t, MPFR_RNDN); mpfr_mul(r10773, r10771, r10772, MPFR_RNDN); mpfr_set_d(r10774, y, MPFR_RNDN); ; mpfr_mul(r10776, r10774, r10775, MPFR_RNDN); mpfr_mul(r10777, r10773, r10776, MPFR_RNDN); mpfr_sub(r10778, r10770, r10777, MPFR_RNDN); mpfr_mul(r10779, r10767, r10778, MPFR_RNDN); mpfr_cbrt(r10780, r10779, MPFR_RNDN); mpfr_mul(r10781, r10780, r10780, MPFR_RNDN); mpfr_mul(r10782, r10781, r10780, MPFR_RNDN); mpfr_sub(r10783, r10766, r10782, MPFR_RNDN); mpfr_set_d(r10784, a, MPFR_RNDN); mpfr_mul(r10785, r10784, r10769, MPFR_RNDN); mpfr_neg(r10786, r10785, MPFR_RNDN); mpfr_mul(r10787, r10772, r10786, MPFR_RNDN); mpfr_add(r10788, r10783, r10787, MPFR_RNDN); ; mpfr_set_si(r10790, mpfr_cmp(r10788, r10789) <= 0, MPFR_RNDN); mpfr_mul(r10791, r10767, r10775, MPFR_RNDN); mpfr_mul(r10792, r10791, r10774, MPFR_RNDN); mpfr_mul(r10793, r10792, r10771, MPFR_RNDN); mpfr_mul(r10794, r10793, r10772, MPFR_RNDN); mpfr_mul(r10795, r10785, r10772, MPFR_RNDN); mpfr_sub(r10796, r10794, r10795, MPFR_RNDN); mpfr_mul(r10797, r10759, r10758, MPFR_RNDN); mpfr_add(r10798, r10796, r10797, MPFR_RNDN); mpfr_mul(r10799, r10767, r10769, MPFR_RNDN); mpfr_mul(r10800, r10799, r10768, MPFR_RNDN); mpfr_sub(r10801, r10798, r10800, MPFR_RNDN); mpfr_mul(r10802, r10762, r10761, MPFR_RNDN); mpfr_mul(r10803, r10802, r10763, MPFR_RNDN); mpfr_sub(r10804, r10801, r10803, MPFR_RNDN); ; mpfr_set_si(r10806, mpfr_cmp(r10788, r10805) <= 0, MPFR_RNDN); mpfr_mul(r10807, r10773, r10774, MPFR_RNDN); mpfr_mul(r10808, r10807, r10775, MPFR_RNDN); mpfr_sub(r10809, r10770, r10808, MPFR_RNDN); mpfr_mul(r10810, r10767, r10809, MPFR_RNDN); mpfr_sub(r10811, r10766, r10810, MPFR_RNDN); mpfr_add(r10812, r10811, r10787, MPFR_RNDN); if (mpfr_get_si(r10806, MPFR_RNDN)) { mpfr_set(r10813, r10812, MPFR_RNDN); } else { mpfr_set(r10813, r10804, MPFR_RNDN); }; if (mpfr_get_si(r10790, MPFR_RNDN)) { mpfr_set(r10814, r10804, MPFR_RNDN); } else { mpfr_set(r10814, r10813, MPFR_RNDN); }; return mpfr_get_d(r10814, MPFR_RNDN); } static mpfr_t r10815, r10816, r10817, r10818, r10819, r10820, r10821, r10822, r10823, r10824, r10825, r10826, r10827, r10828, r10829, r10830, r10831, r10832, r10833, r10834, r10835, r10836, r10837, r10838, r10839, r10840, r10841, r10842, r10843, r10844, r10845, r10846, r10847, r10848, r10849, r10850, r10851, r10852, r10853, r10854, r10855, r10856, r10857, r10858, r10859, r10860, r10861, r10862, r10863, r10864, r10865, r10866, r10867, r10868, r10869, r10870, r10871; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10815); mpfr_init(r10816); mpfr_init(r10817); mpfr_init_set_str(r10818, "27.0", 10, MPFR_RNDN); mpfr_init(r10819); mpfr_init(r10820); mpfr_init(r10821); mpfr_init(r10822); mpfr_init(r10823); mpfr_init(r10824); mpfr_init(r10825); mpfr_init_set_str(r10826, "4.0", 10, MPFR_RNDN); mpfr_init(r10827); mpfr_init(r10828); mpfr_init(r10829); mpfr_init(r10830); mpfr_init(r10831); mpfr_init_set_str(r10832, "18.0", 10, MPFR_RNDN); mpfr_init(r10833); mpfr_init(r10834); mpfr_init(r10835); mpfr_init(r10836); mpfr_init(r10837); mpfr_init(r10838); mpfr_init(r10839); mpfr_init(r10840); mpfr_init(r10841); mpfr_init(r10842); mpfr_init(r10843); mpfr_init(r10844); mpfr_init(r10845); mpfr_init_set_str(r10846, "-1.103406987904582e+307", 10, MPFR_RNDN); mpfr_init(r10847); mpfr_init(r10848); mpfr_init(r10849); mpfr_init(r10850); mpfr_init(r10851); mpfr_init(r10852); mpfr_init(r10853); mpfr_init(r10854); mpfr_init(r10855); mpfr_init(r10856); mpfr_init(r10857); mpfr_init(r10858); mpfr_init(r10859); mpfr_init(r10860); mpfr_init(r10861); mpfr_init_set_str(r10862, "1.3127925574492975e+300", 10, MPFR_RNDN); mpfr_init(r10863); mpfr_init(r10864); mpfr_init(r10865); mpfr_init(r10866); mpfr_init(r10867); mpfr_init(r10868); mpfr_init(r10869); mpfr_init(r10870); mpfr_init(r10871); } 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(r10815, c, MPFR_RNDN); mpfr_set_d(r10816, b, MPFR_RNDN); mpfr_mul(r10817, r10815, r10816, MPFR_RNDN); ; mpfr_set_d(r10819, j, MPFR_RNDN); mpfr_set_d(r10820, k, MPFR_RNDN); mpfr_mul(r10821, r10819, r10820, MPFR_RNDN); mpfr_mul(r10822, r10818, r10821, MPFR_RNDN); mpfr_sub(r10823, r10817, r10822, MPFR_RNDN); mpfr_set_d(r10824, x, MPFR_RNDN); mpfr_set_d(r10825, i, MPFR_RNDN); ; mpfr_mul(r10827, r10825, r10826, MPFR_RNDN); mpfr_set_d(r10828, z, MPFR_RNDN); mpfr_set_d(r10829, t, MPFR_RNDN); mpfr_mul(r10830, r10828, r10829, MPFR_RNDN); mpfr_set_d(r10831, y, MPFR_RNDN); ; mpfr_mul(r10833, r10831, r10832, MPFR_RNDN); mpfr_mul(r10834, r10830, r10833, MPFR_RNDN); mpfr_sub(r10835, r10827, r10834, MPFR_RNDN); mpfr_mul(r10836, r10824, r10835, MPFR_RNDN); mpfr_cbrt(r10837, r10836, MPFR_RNDN); mpfr_mul(r10838, r10837, r10837, MPFR_RNDN); mpfr_mul(r10839, r10838, r10837, MPFR_RNDN); mpfr_sub(r10840, r10823, r10839, MPFR_RNDN); mpfr_set_d(r10841, a, MPFR_RNDN); mpfr_mul(r10842, r10841, r10826, MPFR_RNDN); mpfr_neg(r10843, r10842, MPFR_RNDN); mpfr_mul(r10844, r10829, r10843, MPFR_RNDN); mpfr_add(r10845, r10840, r10844, MPFR_RNDN); ; mpfr_set_si(r10847, mpfr_cmp(r10845, r10846) <= 0, MPFR_RNDN); mpfr_mul(r10848, r10824, r10832, MPFR_RNDN); mpfr_mul(r10849, r10848, r10831, MPFR_RNDN); mpfr_mul(r10850, r10849, r10828, MPFR_RNDN); mpfr_mul(r10851, r10850, r10829, MPFR_RNDN); mpfr_mul(r10852, r10842, r10829, MPFR_RNDN); mpfr_sub(r10853, r10851, r10852, MPFR_RNDN); mpfr_mul(r10854, r10816, r10815, MPFR_RNDN); mpfr_add(r10855, r10853, r10854, MPFR_RNDN); mpfr_mul(r10856, r10824, r10826, MPFR_RNDN); mpfr_mul(r10857, r10856, r10825, MPFR_RNDN); mpfr_sub(r10858, r10855, r10857, MPFR_RNDN); mpfr_mul(r10859, r10819, r10818, MPFR_RNDN); mpfr_mul(r10860, r10859, r10820, MPFR_RNDN); mpfr_sub(r10861, r10858, r10860, MPFR_RNDN); ; mpfr_set_si(r10863, mpfr_cmp(r10845, r10862) <= 0, MPFR_RNDN); mpfr_mul(r10864, r10830, r10831, MPFR_RNDN); mpfr_mul(r10865, r10864, r10832, MPFR_RNDN); mpfr_sub(r10866, r10827, r10865, MPFR_RNDN); mpfr_mul(r10867, r10824, r10866, MPFR_RNDN); mpfr_sub(r10868, r10823, r10867, MPFR_RNDN); mpfr_add(r10869, r10868, r10844, MPFR_RNDN); if (mpfr_get_si(r10863, MPFR_RNDN)) { mpfr_set(r10870, r10869, MPFR_RNDN); } else { mpfr_set(r10870, r10861, MPFR_RNDN); }; if (mpfr_get_si(r10847, MPFR_RNDN)) { mpfr_set(r10871, r10861, MPFR_RNDN); } else { mpfr_set(r10871, r10870, MPFR_RNDN); }; return mpfr_get_d(r10871, MPFR_RNDN); }