#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 r17632 = x;
        float r17633 = 18.0f;
        float r17634 = r17632 * r17633;
        float r17635 = y;
        float r17636 = r17634 * r17635;
        float r17637 = z;
        float r17638 = r17636 * r17637;
        float r17639 = t;
        float r17640 = r17638 * r17639;
        float r17641 = a;
        float r17642 = 4.0f;
        float r17643 = r17641 * r17642;
        float r17644 = r17643 * r17639;
        float r17645 = r17640 - r17644;
        float r17646 = b;
        float r17647 = c;
        float r17648 = r17646 * r17647;
        float r17649 = r17645 + r17648;
        float r17650 = r17632 * r17642;
        float r17651 = i;
        float r17652 = r17650 * r17651;
        float r17653 = r17649 - r17652;
        float r17654 = j;
        float r17655 = 27.0f;
        float r17656 = r17654 * r17655;
        float r17657 = k;
        float r17658 = r17656 * r17657;
        float r17659 = r17653 - r17658;
        return r17659;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17660 = x;
        double r17661 = 18.0;
        double r17662 = r17660 * r17661;
        double r17663 = y;
        double r17664 = r17662 * r17663;
        double r17665 = z;
        double r17666 = r17664 * r17665;
        double r17667 = t;
        double r17668 = r17666 * r17667;
        double r17669 = a;
        double r17670 = 4.0;
        double r17671 = r17669 * r17670;
        double r17672 = r17671 * r17667;
        double r17673 = r17668 - r17672;
        double r17674 = b;
        double r17675 = c;
        double r17676 = r17674 * r17675;
        double r17677 = r17673 + r17676;
        double r17678 = r17660 * r17670;
        double r17679 = i;
        double r17680 = r17678 * r17679;
        double r17681 = r17677 - r17680;
        double r17682 = j;
        double r17683 = 27.0;
        double r17684 = r17682 * r17683;
        double r17685 = k;
        double r17686 = r17684 * r17685;
        double r17687 = r17681 - r17686;
        return r17687;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17688 = t;
        float r17689 = z;
        float r17690 = r17688 * r17689;
        float r17691 = 18.0f;
        float r17692 = x;
        float r17693 = r17691 * r17692;
        float r17694 = y;
        float r17695 = r17693 * r17694;
        float r17696 = c;
        float r17697 = b;
        float r17698 = r17696 * r17697;
        float r17699 = fma(r17690, r17695, r17698);
        float r17700 = 4.0f;
        float r17701 = i;
        float r17702 = a;
        float r17703 = r17688 * r17702;
        float r17704 = fma(r17701, r17692, r17703);
        float r17705 = j;
        float r17706 = 27.0f;
        float r17707 = r17705 * r17706;
        float r17708 = k;
        float r17709 = r17707 * r17708;
        float r17710 = fma(r17700, r17704, r17709);
        float r17711 = r17699 - r17710;
        return r17711;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17712 = t;
        double r17713 = z;
        double r17714 = r17712 * r17713;
        double r17715 = 18.0;
        double r17716 = x;
        double r17717 = r17715 * r17716;
        double r17718 = y;
        double r17719 = r17717 * r17718;
        double r17720 = c;
        double r17721 = b;
        double r17722 = r17720 * r17721;
        double r17723 = fma(r17714, r17719, r17722);
        double r17724 = 4.0;
        double r17725 = i;
        double r17726 = a;
        double r17727 = r17712 * r17726;
        double r17728 = fma(r17725, r17716, r17727);
        double r17729 = j;
        double r17730 = 27.0;
        double r17731 = r17729 * r17730;
        double r17732 = k;
        double r17733 = r17731 * r17732;
        double r17734 = fma(r17724, r17728, r17733);
        double r17735 = r17723 - r17734;
        return r17735;
}

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 r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743, r17744, r17745, r17746, r17747, r17748, r17749, r17750, r17751, r17752, r17753, r17754, r17755, r17756, r17757, r17758, r17759, r17760, r17761, r17762, r17763;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17736);
        mpfr_init_set_str(r17737, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17738);
        mpfr_init(r17739);
        mpfr_init(r17740);
        mpfr_init(r17741);
        mpfr_init(r17742);
        mpfr_init(r17743);
        mpfr_init(r17744);
        mpfr_init(r17745);
        mpfr_init_set_str(r17746, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17747);
        mpfr_init(r17748);
        mpfr_init(r17749);
        mpfr_init(r17750);
        mpfr_init(r17751);
        mpfr_init(r17752);
        mpfr_init(r17753);
        mpfr_init(r17754);
        mpfr_init(r17755);
        mpfr_init(r17756);
        mpfr_init(r17757);
        mpfr_init(r17758);
        mpfr_init_set_str(r17759, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17760);
        mpfr_init(r17761);
        mpfr_init(r17762);
        mpfr_init(r17763);
}

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(r17736, x, MPFR_RNDN);
        ;
        mpfr_mul(r17738, r17736, r17737, MPFR_RNDN);
        mpfr_set_d(r17739, y, MPFR_RNDN);
        mpfr_mul(r17740, r17738, r17739, MPFR_RNDN);
        mpfr_set_d(r17741, z, MPFR_RNDN);
        mpfr_mul(r17742, r17740, r17741, MPFR_RNDN);
        mpfr_set_d(r17743, t, MPFR_RNDN);
        mpfr_mul(r17744, r17742, r17743, MPFR_RNDN);
        mpfr_set_d(r17745, a, MPFR_RNDN);
        ;
        mpfr_mul(r17747, r17745, r17746, MPFR_RNDN);
        mpfr_mul(r17748, r17747, r17743, MPFR_RNDN);
        mpfr_sub(r17749, r17744, r17748, MPFR_RNDN);
        mpfr_set_d(r17750, b, MPFR_RNDN);
        mpfr_set_d(r17751, c, MPFR_RNDN);
        mpfr_mul(r17752, r17750, r17751, MPFR_RNDN);
        mpfr_add(r17753, r17749, r17752, MPFR_RNDN);
        mpfr_mul(r17754, r17736, r17746, MPFR_RNDN);
        mpfr_set_d(r17755, i, MPFR_RNDN);
        mpfr_mul(r17756, r17754, r17755, MPFR_RNDN);
        mpfr_sub(r17757, r17753, r17756, MPFR_RNDN);
        mpfr_set_d(r17758, j, MPFR_RNDN);
        ;
        mpfr_mul(r17760, r17758, r17759, MPFR_RNDN);
        mpfr_set_d(r17761, k, MPFR_RNDN);
        mpfr_mul(r17762, r17760, r17761, MPFR_RNDN);
        mpfr_sub(r17763, r17757, r17762, MPFR_RNDN);
        return mpfr_get_d(r17763, MPFR_RNDN);
}

static mpfr_t r17764, r17765, r17766, r17767, r17768, r17769, r17770, r17771, r17772, r17773, r17774, r17775, r17776, r17777, r17778, r17779, r17780, r17781, r17782, r17783, r17784, r17785, r17786, r17787;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17764);
        mpfr_init(r17765);
        mpfr_init(r17766);
        mpfr_init_set_str(r17767, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17768);
        mpfr_init(r17769);
        mpfr_init(r17770);
        mpfr_init(r17771);
        mpfr_init(r17772);
        mpfr_init(r17773);
        mpfr_init(r17774);
        mpfr_init(r17775);
        mpfr_init_set_str(r17776, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17777);
        mpfr_init(r17778);
        mpfr_init(r17779);
        mpfr_init(r17780);
        mpfr_init(r17781);
        mpfr_init_set_str(r17782, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17783);
        mpfr_init(r17784);
        mpfr_init(r17785);
        mpfr_init(r17786);
        mpfr_init(r17787);
}

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(r17764, t, MPFR_RNDN);
        mpfr_set_d(r17765, z, MPFR_RNDN);
        mpfr_mul(r17766, r17764, r17765, MPFR_RNDN);
        ;
        mpfr_set_d(r17768, x, MPFR_RNDN);
        mpfr_mul(r17769, r17767, r17768, MPFR_RNDN);
        mpfr_set_d(r17770, y, MPFR_RNDN);
        mpfr_mul(r17771, r17769, r17770, MPFR_RNDN);
        mpfr_set_d(r17772, c, MPFR_RNDN);
        mpfr_set_d(r17773, b, MPFR_RNDN);
        mpfr_mul(r17774, r17772, r17773, MPFR_RNDN);
        mpfr_fma(r17775, r17766, r17771, r17774, MPFR_RNDN);
        ;
        mpfr_set_d(r17777, i, MPFR_RNDN);
        mpfr_set_d(r17778, a, MPFR_RNDN);
        mpfr_mul(r17779, r17764, r17778, MPFR_RNDN);
        mpfr_fma(r17780, r17777, r17768, r17779, MPFR_RNDN);
        mpfr_set_d(r17781, j, MPFR_RNDN);
        ;
        mpfr_mul(r17783, r17781, r17782, MPFR_RNDN);
        mpfr_set_d(r17784, k, MPFR_RNDN);
        mpfr_mul(r17785, r17783, r17784, MPFR_RNDN);
        mpfr_fma(r17786, r17776, r17780, r17785, MPFR_RNDN);
        mpfr_sub(r17787, r17775, r17786, MPFR_RNDN);
        return mpfr_get_d(r17787, MPFR_RNDN);
}

static mpfr_t r17788, r17789, r17790, r17791, r17792, r17793, r17794, r17795, r17796, r17797, r17798, r17799, r17800, r17801, r17802, r17803, r17804, r17805, r17806, r17807, r17808, r17809, r17810, r17811;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17788);
        mpfr_init(r17789);
        mpfr_init(r17790);
        mpfr_init_set_str(r17791, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17792);
        mpfr_init(r17793);
        mpfr_init(r17794);
        mpfr_init(r17795);
        mpfr_init(r17796);
        mpfr_init(r17797);
        mpfr_init(r17798);
        mpfr_init(r17799);
        mpfr_init_set_str(r17800, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17801);
        mpfr_init(r17802);
        mpfr_init(r17803);
        mpfr_init(r17804);
        mpfr_init(r17805);
        mpfr_init_set_str(r17806, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17807);
        mpfr_init(r17808);
        mpfr_init(r17809);
        mpfr_init(r17810);
        mpfr_init(r17811);
}

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(r17788, t, MPFR_RNDN);
        mpfr_set_d(r17789, z, MPFR_RNDN);
        mpfr_mul(r17790, r17788, r17789, MPFR_RNDN);
        ;
        mpfr_set_d(r17792, x, MPFR_RNDN);
        mpfr_mul(r17793, r17791, r17792, MPFR_RNDN);
        mpfr_set_d(r17794, y, MPFR_RNDN);
        mpfr_mul(r17795, r17793, r17794, MPFR_RNDN);
        mpfr_set_d(r17796, c, MPFR_RNDN);
        mpfr_set_d(r17797, b, MPFR_RNDN);
        mpfr_mul(r17798, r17796, r17797, MPFR_RNDN);
        mpfr_fma(r17799, r17790, r17795, r17798, MPFR_RNDN);
        ;
        mpfr_set_d(r17801, i, MPFR_RNDN);
        mpfr_set_d(r17802, a, MPFR_RNDN);
        mpfr_mul(r17803, r17788, r17802, MPFR_RNDN);
        mpfr_fma(r17804, r17801, r17792, r17803, MPFR_RNDN);
        mpfr_set_d(r17805, j, MPFR_RNDN);
        ;
        mpfr_mul(r17807, r17805, r17806, MPFR_RNDN);
        mpfr_set_d(r17808, k, MPFR_RNDN);
        mpfr_mul(r17809, r17807, r17808, MPFR_RNDN);
        mpfr_fma(r17810, r17800, r17804, r17809, MPFR_RNDN);
        mpfr_sub(r17811, r17799, r17810, MPFR_RNDN);
        return mpfr_get_d(r17811, MPFR_RNDN);
}