#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 r17511 = x;
        float r17512 = 18.0f;
        float r17513 = r17511 * r17512;
        float r17514 = y;
        float r17515 = r17513 * r17514;
        float r17516 = z;
        float r17517 = r17515 * r17516;
        float r17518 = t;
        float r17519 = r17517 * r17518;
        float r17520 = a;
        float r17521 = 4.0f;
        float r17522 = r17520 * r17521;
        float r17523 = r17522 * r17518;
        float r17524 = r17519 - r17523;
        float r17525 = b;
        float r17526 = c;
        float r17527 = r17525 * r17526;
        float r17528 = r17524 + r17527;
        float r17529 = r17511 * r17521;
        float r17530 = i;
        float r17531 = r17529 * r17530;
        float r17532 = r17528 - r17531;
        float r17533 = j;
        float r17534 = 27.0f;
        float r17535 = r17533 * r17534;
        float r17536 = k;
        float r17537 = r17535 * r17536;
        float r17538 = r17532 - r17537;
        return r17538;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17539 = x;
        double r17540 = 18.0;
        double r17541 = r17539 * r17540;
        double r17542 = y;
        double r17543 = r17541 * r17542;
        double r17544 = z;
        double r17545 = r17543 * r17544;
        double r17546 = t;
        double r17547 = r17545 * r17546;
        double r17548 = a;
        double r17549 = 4.0;
        double r17550 = r17548 * r17549;
        double r17551 = r17550 * r17546;
        double r17552 = r17547 - r17551;
        double r17553 = b;
        double r17554 = c;
        double r17555 = r17553 * r17554;
        double r17556 = r17552 + r17555;
        double r17557 = r17539 * r17549;
        double r17558 = i;
        double r17559 = r17557 * r17558;
        double r17560 = r17556 - r17559;
        double r17561 = j;
        double r17562 = 27.0;
        double r17563 = r17561 * r17562;
        double r17564 = k;
        double r17565 = r17563 * r17564;
        double r17566 = r17560 - r17565;
        return r17566;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17567 = t;
        float r17568 = z;
        float r17569 = r17567 * r17568;
        float r17570 = 18.0f;
        float r17571 = x;
        float r17572 = y;
        float r17573 = r17571 * r17572;
        float r17574 = r17570 * r17573;
        float r17575 = c;
        float r17576 = b;
        float r17577 = r17575 * r17576;
        float r17578 = fma(r17569, r17574, r17577);
        float r17579 = 4.0f;
        float r17580 = i;
        float r17581 = a;
        float r17582 = r17567 * r17581;
        float r17583 = fma(r17580, r17571, r17582);
        float r17584 = j;
        float r17585 = 27.0f;
        float r17586 = r17584 * r17585;
        float r17587 = k;
        float r17588 = r17586 * r17587;
        float r17589 = fma(r17579, r17583, r17588);
        float r17590 = r17578 - r17589;
        return r17590;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17591 = t;
        double r17592 = z;
        double r17593 = r17591 * r17592;
        double r17594 = 18.0;
        double r17595 = x;
        double r17596 = y;
        double r17597 = r17595 * r17596;
        double r17598 = r17594 * r17597;
        double r17599 = c;
        double r17600 = b;
        double r17601 = r17599 * r17600;
        double r17602 = fma(r17593, r17598, r17601);
        double r17603 = 4.0;
        double r17604 = i;
        double r17605 = a;
        double r17606 = r17591 * r17605;
        double r17607 = fma(r17604, r17595, r17606);
        double r17608 = j;
        double r17609 = 27.0;
        double r17610 = r17608 * r17609;
        double r17611 = k;
        double r17612 = r17610 * r17611;
        double r17613 = fma(r17603, r17607, r17612);
        double r17614 = r17602 - r17613;
        return r17614;
}

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 r17615, r17616, r17617, r17618, r17619, r17620, r17621, r17622, r17623, r17624, r17625, r17626, r17627, r17628, r17629, r17630, r17631, r17632, r17633, r17634, r17635, r17636, r17637, r17638, r17639, r17640, r17641, r17642;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17615);
        mpfr_init_set_str(r17616, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17617);
        mpfr_init(r17618);
        mpfr_init(r17619);
        mpfr_init(r17620);
        mpfr_init(r17621);
        mpfr_init(r17622);
        mpfr_init(r17623);
        mpfr_init(r17624);
        mpfr_init_set_str(r17625, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17626);
        mpfr_init(r17627);
        mpfr_init(r17628);
        mpfr_init(r17629);
        mpfr_init(r17630);
        mpfr_init(r17631);
        mpfr_init(r17632);
        mpfr_init(r17633);
        mpfr_init(r17634);
        mpfr_init(r17635);
        mpfr_init(r17636);
        mpfr_init(r17637);
        mpfr_init_set_str(r17638, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17639);
        mpfr_init(r17640);
        mpfr_init(r17641);
        mpfr_init(r17642);
}

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(r17615, x, MPFR_RNDN);
        ;
        mpfr_mul(r17617, r17615, r17616, MPFR_RNDN);
        mpfr_set_d(r17618, y, MPFR_RNDN);
        mpfr_mul(r17619, r17617, r17618, MPFR_RNDN);
        mpfr_set_d(r17620, z, MPFR_RNDN);
        mpfr_mul(r17621, r17619, r17620, MPFR_RNDN);
        mpfr_set_d(r17622, t, MPFR_RNDN);
        mpfr_mul(r17623, r17621, r17622, MPFR_RNDN);
        mpfr_set_d(r17624, a, MPFR_RNDN);
        ;
        mpfr_mul(r17626, r17624, r17625, MPFR_RNDN);
        mpfr_mul(r17627, r17626, r17622, MPFR_RNDN);
        mpfr_sub(r17628, r17623, r17627, MPFR_RNDN);
        mpfr_set_d(r17629, b, MPFR_RNDN);
        mpfr_set_d(r17630, c, MPFR_RNDN);
        mpfr_mul(r17631, r17629, r17630, MPFR_RNDN);
        mpfr_add(r17632, r17628, r17631, MPFR_RNDN);
        mpfr_mul(r17633, r17615, r17625, MPFR_RNDN);
        mpfr_set_d(r17634, i, MPFR_RNDN);
        mpfr_mul(r17635, r17633, r17634, MPFR_RNDN);
        mpfr_sub(r17636, r17632, r17635, MPFR_RNDN);
        mpfr_set_d(r17637, j, MPFR_RNDN);
        ;
        mpfr_mul(r17639, r17637, r17638, MPFR_RNDN);
        mpfr_set_d(r17640, k, MPFR_RNDN);
        mpfr_mul(r17641, r17639, r17640, MPFR_RNDN);
        mpfr_sub(r17642, r17636, r17641, MPFR_RNDN);
        return mpfr_get_d(r17642, MPFR_RNDN);
}

static mpfr_t r17643, r17644, r17645, r17646, r17647, r17648, r17649, r17650, r17651, r17652, r17653, r17654, r17655, r17656, r17657, r17658, r17659, r17660, r17661, r17662, r17663, r17664, r17665, r17666;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17643);
        mpfr_init(r17644);
        mpfr_init(r17645);
        mpfr_init_set_str(r17646, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17647);
        mpfr_init(r17648);
        mpfr_init(r17649);
        mpfr_init(r17650);
        mpfr_init(r17651);
        mpfr_init(r17652);
        mpfr_init(r17653);
        mpfr_init(r17654);
        mpfr_init_set_str(r17655, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17656);
        mpfr_init(r17657);
        mpfr_init(r17658);
        mpfr_init(r17659);
        mpfr_init(r17660);
        mpfr_init_set_str(r17661, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17662);
        mpfr_init(r17663);
        mpfr_init(r17664);
        mpfr_init(r17665);
        mpfr_init(r17666);
}

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(r17643, t, MPFR_RNDN);
        mpfr_set_d(r17644, z, MPFR_RNDN);
        mpfr_mul(r17645, r17643, r17644, MPFR_RNDN);
        ;
        mpfr_set_d(r17647, x, MPFR_RNDN);
        mpfr_set_d(r17648, y, MPFR_RNDN);
        mpfr_mul(r17649, r17647, r17648, MPFR_RNDN);
        mpfr_mul(r17650, r17646, r17649, MPFR_RNDN);
        mpfr_set_d(r17651, c, MPFR_RNDN);
        mpfr_set_d(r17652, b, MPFR_RNDN);
        mpfr_mul(r17653, r17651, r17652, MPFR_RNDN);
        mpfr_fma(r17654, r17645, r17650, r17653, MPFR_RNDN);
        ;
        mpfr_set_d(r17656, i, MPFR_RNDN);
        mpfr_set_d(r17657, a, MPFR_RNDN);
        mpfr_mul(r17658, r17643, r17657, MPFR_RNDN);
        mpfr_fma(r17659, r17656, r17647, r17658, MPFR_RNDN);
        mpfr_set_d(r17660, j, MPFR_RNDN);
        ;
        mpfr_mul(r17662, r17660, r17661, MPFR_RNDN);
        mpfr_set_d(r17663, k, MPFR_RNDN);
        mpfr_mul(r17664, r17662, r17663, MPFR_RNDN);
        mpfr_fma(r17665, r17655, r17659, r17664, MPFR_RNDN);
        mpfr_sub(r17666, r17654, r17665, MPFR_RNDN);
        return mpfr_get_d(r17666, MPFR_RNDN);
}

static mpfr_t r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17667);
        mpfr_init(r17668);
        mpfr_init(r17669);
        mpfr_init_set_str(r17670, "18.0", 10, MPFR_RNDN);
        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_set_str(r17679, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17680);
        mpfr_init(r17681);
        mpfr_init(r17682);
        mpfr_init(r17683);
        mpfr_init(r17684);
        mpfr_init_set_str(r17685, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17686);
        mpfr_init(r17687);
        mpfr_init(r17688);
        mpfr_init(r17689);
        mpfr_init(r17690);
}

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(r17667, t, MPFR_RNDN);
        mpfr_set_d(r17668, z, MPFR_RNDN);
        mpfr_mul(r17669, r17667, r17668, MPFR_RNDN);
        ;
        mpfr_set_d(r17671, x, MPFR_RNDN);
        mpfr_set_d(r17672, y, MPFR_RNDN);
        mpfr_mul(r17673, r17671, r17672, MPFR_RNDN);
        mpfr_mul(r17674, r17670, r17673, MPFR_RNDN);
        mpfr_set_d(r17675, c, MPFR_RNDN);
        mpfr_set_d(r17676, b, MPFR_RNDN);
        mpfr_mul(r17677, r17675, r17676, MPFR_RNDN);
        mpfr_fma(r17678, r17669, r17674, r17677, MPFR_RNDN);
        ;
        mpfr_set_d(r17680, i, MPFR_RNDN);
        mpfr_set_d(r17681, a, MPFR_RNDN);
        mpfr_mul(r17682, r17667, r17681, MPFR_RNDN);
        mpfr_fma(r17683, r17680, r17671, r17682, MPFR_RNDN);
        mpfr_set_d(r17684, j, MPFR_RNDN);
        ;
        mpfr_mul(r17686, r17684, r17685, MPFR_RNDN);
        mpfr_set_d(r17687, k, MPFR_RNDN);
        mpfr_mul(r17688, r17686, r17687, MPFR_RNDN);
        mpfr_fma(r17689, r17679, r17683, r17688, MPFR_RNDN);
        mpfr_sub(r17690, r17678, r17689, MPFR_RNDN);
        return mpfr_get_d(r17690, MPFR_RNDN);
}