#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 r17954 = x;
        float r17955 = 18.0f;
        float r17956 = r17954 * r17955;
        float r17957 = y;
        float r17958 = r17956 * r17957;
        float r17959 = z;
        float r17960 = r17958 * r17959;
        float r17961 = t;
        float r17962 = r17960 * r17961;
        float r17963 = a;
        float r17964 = 4.0f;
        float r17965 = r17963 * r17964;
        float r17966 = r17965 * r17961;
        float r17967 = r17962 - r17966;
        float r17968 = b;
        float r17969 = c;
        float r17970 = r17968 * r17969;
        float r17971 = r17967 + r17970;
        float r17972 = r17954 * r17964;
        float r17973 = i;
        float r17974 = r17972 * r17973;
        float r17975 = r17971 - r17974;
        float r17976 = j;
        float r17977 = 27.0f;
        float r17978 = r17976 * r17977;
        float r17979 = k;
        float r17980 = r17978 * r17979;
        float r17981 = r17975 - r17980;
        return r17981;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17982 = x;
        double r17983 = 18.0;
        double r17984 = r17982 * r17983;
        double r17985 = y;
        double r17986 = r17984 * r17985;
        double r17987 = z;
        double r17988 = r17986 * r17987;
        double r17989 = t;
        double r17990 = r17988 * r17989;
        double r17991 = a;
        double r17992 = 4.0;
        double r17993 = r17991 * r17992;
        double r17994 = r17993 * r17989;
        double r17995 = r17990 - r17994;
        double r17996 = b;
        double r17997 = c;
        double r17998 = r17996 * r17997;
        double r17999 = r17995 + r17998;
        double r18000 = r17982 * r17992;
        double r18001 = i;
        double r18002 = r18000 * r18001;
        double r18003 = r17999 - r18002;
        double r18004 = j;
        double r18005 = 27.0;
        double r18006 = r18004 * r18005;
        double r18007 = k;
        double r18008 = r18006 * r18007;
        double r18009 = r18003 - r18008;
        return r18009;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r18010 = z;
        float r18011 = t;
        float r18012 = r18010 * r18011;
        float r18013 = y;
        float r18014 = x;
        float r18015 = r18013 * r18014;
        float r18016 = 18.0f;
        float r18017 = r18015 * r18016;
        float r18018 = b;
        float r18019 = c;
        float r18020 = r18018 * r18019;
        float r18021 = fma(r18012, r18017, r18020);
        float r18022 = 4.0f;
        float r18023 = i;
        float r18024 = a;
        float r18025 = r18011 * r18024;
        float r18026 = fma(r18023, r18014, r18025);
        float r18027 = 27.0f;
        float r18028 = j;
        float r18029 = r18027 * r18028;
        float r18030 = k;
        float r18031 = r18029 * r18030;
        float r18032 = fma(r18022, r18026, r18031);
        float r18033 = r18021 - r18032;
        return r18033;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r18034 = z;
        double r18035 = t;
        double r18036 = r18034 * r18035;
        double r18037 = y;
        double r18038 = x;
        double r18039 = r18037 * r18038;
        double r18040 = 18.0;
        double r18041 = r18039 * r18040;
        double r18042 = b;
        double r18043 = c;
        double r18044 = r18042 * r18043;
        double r18045 = fma(r18036, r18041, r18044);
        double r18046 = 4.0;
        double r18047 = i;
        double r18048 = a;
        double r18049 = r18035 * r18048;
        double r18050 = fma(r18047, r18038, r18049);
        double r18051 = 27.0;
        double r18052 = j;
        double r18053 = r18051 * r18052;
        double r18054 = k;
        double r18055 = r18053 * r18054;
        double r18056 = fma(r18046, r18050, r18055);
        double r18057 = r18045 - r18056;
        return r18057;
}

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 r18058, r18059, r18060, r18061, r18062, r18063, r18064, r18065, r18066, r18067, r18068, r18069, r18070, r18071, r18072, r18073, r18074, r18075, r18076, r18077, r18078, r18079, r18080, r18081, r18082, r18083, r18084, r18085;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18058);
        mpfr_init_set_str(r18059, "18.0", 10, MPFR_RNDN);
        mpfr_init(r18060);
        mpfr_init(r18061);
        mpfr_init(r18062);
        mpfr_init(r18063);
        mpfr_init(r18064);
        mpfr_init(r18065);
        mpfr_init(r18066);
        mpfr_init(r18067);
        mpfr_init_set_str(r18068, "4.0", 10, MPFR_RNDN);
        mpfr_init(r18069);
        mpfr_init(r18070);
        mpfr_init(r18071);
        mpfr_init(r18072);
        mpfr_init(r18073);
        mpfr_init(r18074);
        mpfr_init(r18075);
        mpfr_init(r18076);
        mpfr_init(r18077);
        mpfr_init(r18078);
        mpfr_init(r18079);
        mpfr_init(r18080);
        mpfr_init_set_str(r18081, "27.0", 10, MPFR_RNDN);
        mpfr_init(r18082);
        mpfr_init(r18083);
        mpfr_init(r18084);
        mpfr_init(r18085);
}

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(r18058, x, MPFR_RNDN);
        ;
        mpfr_mul(r18060, r18058, r18059, MPFR_RNDN);
        mpfr_set_d(r18061, y, MPFR_RNDN);
        mpfr_mul(r18062, r18060, r18061, MPFR_RNDN);
        mpfr_set_d(r18063, z, MPFR_RNDN);
        mpfr_mul(r18064, r18062, r18063, MPFR_RNDN);
        mpfr_set_d(r18065, t, MPFR_RNDN);
        mpfr_mul(r18066, r18064, r18065, MPFR_RNDN);
        mpfr_set_d(r18067, a, MPFR_RNDN);
        ;
        mpfr_mul(r18069, r18067, r18068, MPFR_RNDN);
        mpfr_mul(r18070, r18069, r18065, MPFR_RNDN);
        mpfr_sub(r18071, r18066, r18070, MPFR_RNDN);
        mpfr_set_d(r18072, b, MPFR_RNDN);
        mpfr_set_d(r18073, c, MPFR_RNDN);
        mpfr_mul(r18074, r18072, r18073, MPFR_RNDN);
        mpfr_add(r18075, r18071, r18074, MPFR_RNDN);
        mpfr_mul(r18076, r18058, r18068, MPFR_RNDN);
        mpfr_set_d(r18077, i, MPFR_RNDN);
        mpfr_mul(r18078, r18076, r18077, MPFR_RNDN);
        mpfr_sub(r18079, r18075, r18078, MPFR_RNDN);
        mpfr_set_d(r18080, j, MPFR_RNDN);
        ;
        mpfr_mul(r18082, r18080, r18081, MPFR_RNDN);
        mpfr_set_d(r18083, k, MPFR_RNDN);
        mpfr_mul(r18084, r18082, r18083, MPFR_RNDN);
        mpfr_sub(r18085, r18079, r18084, MPFR_RNDN);
        return mpfr_get_d(r18085, MPFR_RNDN);
}

static mpfr_t r18086, r18087, r18088, r18089, r18090, r18091, r18092, r18093, r18094, r18095, r18096, r18097, r18098, r18099, r18100, r18101, r18102, r18103, r18104, r18105, r18106, r18107, r18108, r18109;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18086);
        mpfr_init(r18087);
        mpfr_init(r18088);
        mpfr_init(r18089);
        mpfr_init(r18090);
        mpfr_init(r18091);
        mpfr_init_set_str(r18092, "18.0", 10, MPFR_RNDN);
        mpfr_init(r18093);
        mpfr_init(r18094);
        mpfr_init(r18095);
        mpfr_init(r18096);
        mpfr_init(r18097);
        mpfr_init_set_str(r18098, "4.0", 10, MPFR_RNDN);
        mpfr_init(r18099);
        mpfr_init(r18100);
        mpfr_init(r18101);
        mpfr_init(r18102);
        mpfr_init_set_str(r18103, "27.0", 10, MPFR_RNDN);
        mpfr_init(r18104);
        mpfr_init(r18105);
        mpfr_init(r18106);
        mpfr_init(r18107);
        mpfr_init(r18108);
        mpfr_init(r18109);
}

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(r18086, z, MPFR_RNDN);
        mpfr_set_d(r18087, t, MPFR_RNDN);
        mpfr_mul(r18088, r18086, r18087, MPFR_RNDN);
        mpfr_set_d(r18089, y, MPFR_RNDN);
        mpfr_set_d(r18090, x, MPFR_RNDN);
        mpfr_mul(r18091, r18089, r18090, MPFR_RNDN);
        ;
        mpfr_mul(r18093, r18091, r18092, MPFR_RNDN);
        mpfr_set_d(r18094, b, MPFR_RNDN);
        mpfr_set_d(r18095, c, MPFR_RNDN);
        mpfr_mul(r18096, r18094, r18095, MPFR_RNDN);
        mpfr_fma(r18097, r18088, r18093, r18096, MPFR_RNDN);
        ;
        mpfr_set_d(r18099, i, MPFR_RNDN);
        mpfr_set_d(r18100, a, MPFR_RNDN);
        mpfr_mul(r18101, r18087, r18100, MPFR_RNDN);
        mpfr_fma(r18102, r18099, r18090, r18101, MPFR_RNDN);
        ;
        mpfr_set_d(r18104, j, MPFR_RNDN);
        mpfr_mul(r18105, r18103, r18104, MPFR_RNDN);
        mpfr_set_d(r18106, k, MPFR_RNDN);
        mpfr_mul(r18107, r18105, r18106, MPFR_RNDN);
        mpfr_fma(r18108, r18098, r18102, r18107, MPFR_RNDN);
        mpfr_sub(r18109, r18097, r18108, MPFR_RNDN);
        return mpfr_get_d(r18109, MPFR_RNDN);
}

static mpfr_t r18110, r18111, r18112, r18113, r18114, r18115, r18116, r18117, r18118, r18119, r18120, r18121, r18122, r18123, r18124, r18125, r18126, r18127, r18128, r18129, r18130, r18131, r18132, r18133;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18110);
        mpfr_init(r18111);
        mpfr_init(r18112);
        mpfr_init(r18113);
        mpfr_init(r18114);
        mpfr_init(r18115);
        mpfr_init_set_str(r18116, "18.0", 10, MPFR_RNDN);
        mpfr_init(r18117);
        mpfr_init(r18118);
        mpfr_init(r18119);
        mpfr_init(r18120);
        mpfr_init(r18121);
        mpfr_init_set_str(r18122, "4.0", 10, MPFR_RNDN);
        mpfr_init(r18123);
        mpfr_init(r18124);
        mpfr_init(r18125);
        mpfr_init(r18126);
        mpfr_init_set_str(r18127, "27.0", 10, MPFR_RNDN);
        mpfr_init(r18128);
        mpfr_init(r18129);
        mpfr_init(r18130);
        mpfr_init(r18131);
        mpfr_init(r18132);
        mpfr_init(r18133);
}

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(r18110, z, MPFR_RNDN);
        mpfr_set_d(r18111, t, MPFR_RNDN);
        mpfr_mul(r18112, r18110, r18111, MPFR_RNDN);
        mpfr_set_d(r18113, y, MPFR_RNDN);
        mpfr_set_d(r18114, x, MPFR_RNDN);
        mpfr_mul(r18115, r18113, r18114, MPFR_RNDN);
        ;
        mpfr_mul(r18117, r18115, r18116, MPFR_RNDN);
        mpfr_set_d(r18118, b, MPFR_RNDN);
        mpfr_set_d(r18119, c, MPFR_RNDN);
        mpfr_mul(r18120, r18118, r18119, MPFR_RNDN);
        mpfr_fma(r18121, r18112, r18117, r18120, MPFR_RNDN);
        ;
        mpfr_set_d(r18123, i, MPFR_RNDN);
        mpfr_set_d(r18124, a, MPFR_RNDN);
        mpfr_mul(r18125, r18111, r18124, MPFR_RNDN);
        mpfr_fma(r18126, r18123, r18114, r18125, MPFR_RNDN);
        ;
        mpfr_set_d(r18128, j, MPFR_RNDN);
        mpfr_mul(r18129, r18127, r18128, MPFR_RNDN);
        mpfr_set_d(r18130, k, MPFR_RNDN);
        mpfr_mul(r18131, r18129, r18130, MPFR_RNDN);
        mpfr_fma(r18132, r18122, r18126, r18131, MPFR_RNDN);
        mpfr_sub(r18133, r18121, r18132, MPFR_RNDN);
        return mpfr_get_d(r18133, MPFR_RNDN);
}