#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

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 r17549 = x;
        float r17550 = 18.0f;
        float r17551 = r17549 * r17550;
        float r17552 = y;
        float r17553 = r17551 * r17552;
        float r17554 = z;
        float r17555 = r17553 * r17554;
        float r17556 = t;
        float r17557 = r17555 * r17556;
        float r17558 = a;
        float r17559 = 4.0f;
        float r17560 = r17558 * r17559;
        float r17561 = r17560 * r17556;
        float r17562 = r17557 - r17561;
        float r17563 = b;
        float r17564 = c;
        float r17565 = r17563 * r17564;
        float r17566 = r17562 + r17565;
        float r17567 = r17549 * r17559;
        float r17568 = i;
        float r17569 = r17567 * r17568;
        float r17570 = r17566 - r17569;
        float r17571 = j;
        float r17572 = 27.0f;
        float r17573 = r17571 * r17572;
        float r17574 = k;
        float r17575 = r17573 * r17574;
        float r17576 = r17570 - r17575;
        return r17576;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17577 = x;
        double r17578 = 18.0;
        double r17579 = r17577 * r17578;
        double r17580 = y;
        double r17581 = r17579 * r17580;
        double r17582 = z;
        double r17583 = r17581 * r17582;
        double r17584 = t;
        double r17585 = r17583 * r17584;
        double r17586 = a;
        double r17587 = 4.0;
        double r17588 = r17586 * r17587;
        double r17589 = r17588 * r17584;
        double r17590 = r17585 - r17589;
        double r17591 = b;
        double r17592 = c;
        double r17593 = r17591 * r17592;
        double r17594 = r17590 + r17593;
        double r17595 = r17577 * r17587;
        double r17596 = i;
        double r17597 = r17595 * r17596;
        double r17598 = r17594 - r17597;
        double r17599 = j;
        double r17600 = 27.0;
        double r17601 = r17599 * r17600;
        double r17602 = k;
        double r17603 = r17601 * r17602;
        double r17604 = r17598 - r17603;
        return r17604;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17605 = t;
        float r17606 = z;
        float r17607 = r17605 * r17606;
        float r17608 = 18.0f;
        float r17609 = x;
        float r17610 = r17608 * r17609;
        float r17611 = y;
        float r17612 = r17610 * r17611;
        float r17613 = r17607 * r17612;
        float r17614 = 4.0f;
        float r17615 = a;
        float r17616 = r17605 * r17615;
        float r17617 = i;
        float r17618 = r17617 * r17609;
        float r17619 = r17616 + r17618;
        float r17620 = r17614 * r17619;
        float r17621 = r17613 - r17620;
        float r17622 = c;
        float r17623 = b;
        float r17624 = r17622 * r17623;
        float r17625 = j;
        float r17626 = 27.0f;
        float r17627 = k;
        float r17628 = r17626 * r17627;
        float r17629 = r17625 * r17628;
        float r17630 = r17624 - r17629;
        float r17631 = r17621 + r17630;
        return r17631;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17632 = t;
        double r17633 = z;
        double r17634 = r17632 * r17633;
        double r17635 = 18.0;
        double r17636 = x;
        double r17637 = r17635 * r17636;
        double r17638 = y;
        double r17639 = r17637 * r17638;
        double r17640 = r17634 * r17639;
        double r17641 = 4.0;
        double r17642 = a;
        double r17643 = r17632 * r17642;
        double r17644 = i;
        double r17645 = r17644 * r17636;
        double r17646 = r17643 + r17645;
        double r17647 = r17641 * r17646;
        double r17648 = r17640 - r17647;
        double r17649 = c;
        double r17650 = b;
        double r17651 = r17649 * r17650;
        double r17652 = j;
        double r17653 = 27.0;
        double r17654 = k;
        double r17655 = r17653 * r17654;
        double r17656 = r17652 * r17655;
        double r17657 = r17651 - r17656;
        double r17658 = r17648 + r17657;
        return r17658;
}

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 r17659, r17660, r17661, r17662, r17663, r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17659);
        mpfr_init_set_str(r17660, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17661);
        mpfr_init(r17662);
        mpfr_init(r17663);
        mpfr_init(r17664);
        mpfr_init(r17665);
        mpfr_init(r17666);
        mpfr_init(r17667);
        mpfr_init(r17668);
        mpfr_init_set_str(r17669, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17670);
        mpfr_init(r17671);
        mpfr_init(r17672);
        mpfr_init(r17673);
        mpfr_init(r17674);
        mpfr_init(r17675);
        mpfr_init(r17676);
        mpfr_init(r17677);
        mpfr_init(r17678);
        mpfr_init(r17679);
        mpfr_init(r17680);
        mpfr_init(r17681);
        mpfr_init_set_str(r17682, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17683);
        mpfr_init(r17684);
        mpfr_init(r17685);
        mpfr_init(r17686);
}

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(r17659, x, MPFR_RNDN);
        ;
        mpfr_mul(r17661, r17659, r17660, MPFR_RNDN);
        mpfr_set_d(r17662, y, MPFR_RNDN);
        mpfr_mul(r17663, r17661, r17662, MPFR_RNDN);
        mpfr_set_d(r17664, z, MPFR_RNDN);
        mpfr_mul(r17665, r17663, r17664, MPFR_RNDN);
        mpfr_set_d(r17666, t, MPFR_RNDN);
        mpfr_mul(r17667, r17665, r17666, MPFR_RNDN);
        mpfr_set_d(r17668, a, MPFR_RNDN);
        ;
        mpfr_mul(r17670, r17668, r17669, MPFR_RNDN);
        mpfr_mul(r17671, r17670, r17666, MPFR_RNDN);
        mpfr_sub(r17672, r17667, r17671, MPFR_RNDN);
        mpfr_set_d(r17673, b, MPFR_RNDN);
        mpfr_set_d(r17674, c, MPFR_RNDN);
        mpfr_mul(r17675, r17673, r17674, MPFR_RNDN);
        mpfr_add(r17676, r17672, r17675, MPFR_RNDN);
        mpfr_mul(r17677, r17659, r17669, MPFR_RNDN);
        mpfr_set_d(r17678, i, MPFR_RNDN);
        mpfr_mul(r17679, r17677, r17678, MPFR_RNDN);
        mpfr_sub(r17680, r17676, r17679, MPFR_RNDN);
        mpfr_set_d(r17681, j, MPFR_RNDN);
        ;
        mpfr_mul(r17683, r17681, r17682, MPFR_RNDN);
        mpfr_set_d(r17684, k, MPFR_RNDN);
        mpfr_mul(r17685, r17683, r17684, MPFR_RNDN);
        mpfr_sub(r17686, r17680, r17685, MPFR_RNDN);
        return mpfr_get_d(r17686, MPFR_RNDN);
}

static mpfr_t r17687, r17688, r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703, r17704, r17705, r17706, r17707, r17708, r17709, r17710, r17711, r17712, r17713;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17687);
        mpfr_init(r17688);
        mpfr_init(r17689);
        mpfr_init_set_str(r17690, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17691);
        mpfr_init(r17692);
        mpfr_init(r17693);
        mpfr_init(r17694);
        mpfr_init(r17695);
        mpfr_init_set_str(r17696, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17697);
        mpfr_init(r17698);
        mpfr_init(r17699);
        mpfr_init(r17700);
        mpfr_init(r17701);
        mpfr_init(r17702);
        mpfr_init(r17703);
        mpfr_init(r17704);
        mpfr_init(r17705);
        mpfr_init(r17706);
        mpfr_init(r17707);
        mpfr_init_set_str(r17708, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17709);
        mpfr_init(r17710);
        mpfr_init(r17711);
        mpfr_init(r17712);
        mpfr_init(r17713);
}

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(r17687, t, MPFR_RNDN);
        mpfr_set_d(r17688, z, MPFR_RNDN);
        mpfr_mul(r17689, r17687, r17688, MPFR_RNDN);
        ;
        mpfr_set_d(r17691, x, MPFR_RNDN);
        mpfr_mul(r17692, r17690, r17691, MPFR_RNDN);
        mpfr_set_d(r17693, y, MPFR_RNDN);
        mpfr_mul(r17694, r17692, r17693, MPFR_RNDN);
        mpfr_mul(r17695, r17689, r17694, MPFR_RNDN);
        ;
        mpfr_set_d(r17697, a, MPFR_RNDN);
        mpfr_mul(r17698, r17687, r17697, MPFR_RNDN);
        mpfr_set_d(r17699, i, MPFR_RNDN);
        mpfr_mul(r17700, r17699, r17691, MPFR_RNDN);
        mpfr_add(r17701, r17698, r17700, MPFR_RNDN);
        mpfr_mul(r17702, r17696, r17701, MPFR_RNDN);
        mpfr_sub(r17703, r17695, r17702, MPFR_RNDN);
        mpfr_set_d(r17704, c, MPFR_RNDN);
        mpfr_set_d(r17705, b, MPFR_RNDN);
        mpfr_mul(r17706, r17704, r17705, MPFR_RNDN);
        mpfr_set_d(r17707, j, MPFR_RNDN);
        ;
        mpfr_set_d(r17709, k, MPFR_RNDN);
        mpfr_mul(r17710, r17708, r17709, MPFR_RNDN);
        mpfr_mul(r17711, r17707, r17710, MPFR_RNDN);
        mpfr_sub(r17712, r17706, r17711, MPFR_RNDN);
        mpfr_add(r17713, r17703, r17712, MPFR_RNDN);
        return mpfr_get_d(r17713, MPFR_RNDN);
}

static mpfr_t r17714, r17715, r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732, r17733, r17734, r17735, r17736, r17737, r17738, r17739, r17740;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17714);
        mpfr_init(r17715);
        mpfr_init(r17716);
        mpfr_init_set_str(r17717, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17718);
        mpfr_init(r17719);
        mpfr_init(r17720);
        mpfr_init(r17721);
        mpfr_init(r17722);
        mpfr_init_set_str(r17723, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17724);
        mpfr_init(r17725);
        mpfr_init(r17726);
        mpfr_init(r17727);
        mpfr_init(r17728);
        mpfr_init(r17729);
        mpfr_init(r17730);
        mpfr_init(r17731);
        mpfr_init(r17732);
        mpfr_init(r17733);
        mpfr_init(r17734);
        mpfr_init_set_str(r17735, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17736);
        mpfr_init(r17737);
        mpfr_init(r17738);
        mpfr_init(r17739);
        mpfr_init(r17740);
}

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(r17714, t, MPFR_RNDN);
        mpfr_set_d(r17715, z, MPFR_RNDN);
        mpfr_mul(r17716, r17714, r17715, MPFR_RNDN);
        ;
        mpfr_set_d(r17718, x, MPFR_RNDN);
        mpfr_mul(r17719, r17717, r17718, MPFR_RNDN);
        mpfr_set_d(r17720, y, MPFR_RNDN);
        mpfr_mul(r17721, r17719, r17720, MPFR_RNDN);
        mpfr_mul(r17722, r17716, r17721, MPFR_RNDN);
        ;
        mpfr_set_d(r17724, a, MPFR_RNDN);
        mpfr_mul(r17725, r17714, r17724, MPFR_RNDN);
        mpfr_set_d(r17726, i, MPFR_RNDN);
        mpfr_mul(r17727, r17726, r17718, MPFR_RNDN);
        mpfr_add(r17728, r17725, r17727, MPFR_RNDN);
        mpfr_mul(r17729, r17723, r17728, MPFR_RNDN);
        mpfr_sub(r17730, r17722, r17729, MPFR_RNDN);
        mpfr_set_d(r17731, c, MPFR_RNDN);
        mpfr_set_d(r17732, b, MPFR_RNDN);
        mpfr_mul(r17733, r17731, r17732, MPFR_RNDN);
        mpfr_set_d(r17734, j, MPFR_RNDN);
        ;
        mpfr_set_d(r17736, k, MPFR_RNDN);
        mpfr_mul(r17737, r17735, r17736, MPFR_RNDN);
        mpfr_mul(r17738, r17734, r17737, MPFR_RNDN);
        mpfr_sub(r17739, r17733, r17738, MPFR_RNDN);
        mpfr_add(r17740, r17730, r17739, MPFR_RNDN);
        return mpfr_get_d(r17740, MPFR_RNDN);
}