#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 r17031 = x;
        float r17032 = 18.0f;
        float r17033 = r17031 * r17032;
        float r17034 = y;
        float r17035 = r17033 * r17034;
        float r17036 = z;
        float r17037 = r17035 * r17036;
        float r17038 = t;
        float r17039 = r17037 * r17038;
        float r17040 = a;
        float r17041 = 4.0f;
        float r17042 = r17040 * r17041;
        float r17043 = r17042 * r17038;
        float r17044 = r17039 - r17043;
        float r17045 = b;
        float r17046 = c;
        float r17047 = r17045 * r17046;
        float r17048 = r17044 + r17047;
        float r17049 = r17031 * r17041;
        float r17050 = i;
        float r17051 = r17049 * r17050;
        float r17052 = r17048 - r17051;
        float r17053 = j;
        float r17054 = 27.0f;
        float r17055 = r17053 * r17054;
        float r17056 = k;
        float r17057 = r17055 * r17056;
        float r17058 = r17052 - r17057;
        return r17058;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17059 = x;
        double r17060 = 18.0;
        double r17061 = r17059 * r17060;
        double r17062 = y;
        double r17063 = r17061 * r17062;
        double r17064 = z;
        double r17065 = r17063 * r17064;
        double r17066 = t;
        double r17067 = r17065 * r17066;
        double r17068 = a;
        double r17069 = 4.0;
        double r17070 = r17068 * r17069;
        double r17071 = r17070 * r17066;
        double r17072 = r17067 - r17071;
        double r17073 = b;
        double r17074 = c;
        double r17075 = r17073 * r17074;
        double r17076 = r17072 + r17075;
        double r17077 = r17059 * r17069;
        double r17078 = i;
        double r17079 = r17077 * r17078;
        double r17080 = r17076 - r17079;
        double r17081 = j;
        double r17082 = 27.0;
        double r17083 = r17081 * r17082;
        double r17084 = k;
        double r17085 = r17083 * r17084;
        double r17086 = r17080 - r17085;
        return r17086;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17087 = t;
        float r17088 = z;
        float r17089 = r17087 * r17088;
        float r17090 = 18.0f;
        float r17091 = x;
        float r17092 = y;
        float r17093 = r17091 * r17092;
        float r17094 = r17090 * r17093;
        float r17095 = c;
        float r17096 = b;
        float r17097 = r17095 * r17096;
        float r17098 = fma(r17089, r17094, r17097);
        float r17099 = 4.0f;
        float r17100 = i;
        float r17101 = a;
        float r17102 = r17087 * r17101;
        float r17103 = fma(r17100, r17091, r17102);
        float r17104 = j;
        float r17105 = 27.0f;
        float r17106 = r17104 * r17105;
        float r17107 = k;
        float r17108 = r17106 * r17107;
        float r17109 = fma(r17099, r17103, r17108);
        float r17110 = r17098 - r17109;
        return r17110;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17111 = t;
        double r17112 = z;
        double r17113 = r17111 * r17112;
        double r17114 = 18.0;
        double r17115 = x;
        double r17116 = y;
        double r17117 = r17115 * r17116;
        double r17118 = r17114 * r17117;
        double r17119 = c;
        double r17120 = b;
        double r17121 = r17119 * r17120;
        double r17122 = fma(r17113, r17118, r17121);
        double r17123 = 4.0;
        double r17124 = i;
        double r17125 = a;
        double r17126 = r17111 * r17125;
        double r17127 = fma(r17124, r17115, r17126);
        double r17128 = j;
        double r17129 = 27.0;
        double r17130 = r17128 * r17129;
        double r17131 = k;
        double r17132 = r17130 * r17131;
        double r17133 = fma(r17123, r17127, r17132);
        double r17134 = r17122 - r17133;
        return r17134;
}

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 r17135, r17136, r17137, r17138, r17139, r17140, r17141, r17142, r17143, r17144, r17145, r17146, r17147, r17148, r17149, r17150, r17151, r17152, r17153, r17154, r17155, r17156, r17157, r17158, r17159, r17160, r17161, r17162;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17135);
        mpfr_init_set_str(r17136, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17137);
        mpfr_init(r17138);
        mpfr_init(r17139);
        mpfr_init(r17140);
        mpfr_init(r17141);
        mpfr_init(r17142);
        mpfr_init(r17143);
        mpfr_init(r17144);
        mpfr_init_set_str(r17145, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17146);
        mpfr_init(r17147);
        mpfr_init(r17148);
        mpfr_init(r17149);
        mpfr_init(r17150);
        mpfr_init(r17151);
        mpfr_init(r17152);
        mpfr_init(r17153);
        mpfr_init(r17154);
        mpfr_init(r17155);
        mpfr_init(r17156);
        mpfr_init(r17157);
        mpfr_init_set_str(r17158, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17159);
        mpfr_init(r17160);
        mpfr_init(r17161);
        mpfr_init(r17162);
}

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(r17135, x, MPFR_RNDN);
        ;
        mpfr_mul(r17137, r17135, r17136, MPFR_RNDN);
        mpfr_set_d(r17138, y, MPFR_RNDN);
        mpfr_mul(r17139, r17137, r17138, MPFR_RNDN);
        mpfr_set_d(r17140, z, MPFR_RNDN);
        mpfr_mul(r17141, r17139, r17140, MPFR_RNDN);
        mpfr_set_d(r17142, t, MPFR_RNDN);
        mpfr_mul(r17143, r17141, r17142, MPFR_RNDN);
        mpfr_set_d(r17144, a, MPFR_RNDN);
        ;
        mpfr_mul(r17146, r17144, r17145, MPFR_RNDN);
        mpfr_mul(r17147, r17146, r17142, MPFR_RNDN);
        mpfr_sub(r17148, r17143, r17147, MPFR_RNDN);
        mpfr_set_d(r17149, b, MPFR_RNDN);
        mpfr_set_d(r17150, c, MPFR_RNDN);
        mpfr_mul(r17151, r17149, r17150, MPFR_RNDN);
        mpfr_add(r17152, r17148, r17151, MPFR_RNDN);
        mpfr_mul(r17153, r17135, r17145, MPFR_RNDN);
        mpfr_set_d(r17154, i, MPFR_RNDN);
        mpfr_mul(r17155, r17153, r17154, MPFR_RNDN);
        mpfr_sub(r17156, r17152, r17155, MPFR_RNDN);
        mpfr_set_d(r17157, j, MPFR_RNDN);
        ;
        mpfr_mul(r17159, r17157, r17158, MPFR_RNDN);
        mpfr_set_d(r17160, k, MPFR_RNDN);
        mpfr_mul(r17161, r17159, r17160, MPFR_RNDN);
        mpfr_sub(r17162, r17156, r17161, MPFR_RNDN);
        return mpfr_get_d(r17162, MPFR_RNDN);
}

static mpfr_t r17163, r17164, r17165, r17166, r17167, r17168, r17169, r17170, r17171, r17172, r17173, r17174, r17175, r17176, r17177, r17178, r17179, r17180, r17181, r17182, r17183, r17184, r17185, r17186;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17163);
        mpfr_init(r17164);
        mpfr_init(r17165);
        mpfr_init_set_str(r17166, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17167);
        mpfr_init(r17168);
        mpfr_init(r17169);
        mpfr_init(r17170);
        mpfr_init(r17171);
        mpfr_init(r17172);
        mpfr_init(r17173);
        mpfr_init(r17174);
        mpfr_init_set_str(r17175, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17176);
        mpfr_init(r17177);
        mpfr_init(r17178);
        mpfr_init(r17179);
        mpfr_init(r17180);
        mpfr_init_set_str(r17181, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17182);
        mpfr_init(r17183);
        mpfr_init(r17184);
        mpfr_init(r17185);
        mpfr_init(r17186);
}

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(r17163, t, MPFR_RNDN);
        mpfr_set_d(r17164, z, MPFR_RNDN);
        mpfr_mul(r17165, r17163, r17164, MPFR_RNDN);
        ;
        mpfr_set_d(r17167, x, MPFR_RNDN);
        mpfr_set_d(r17168, y, MPFR_RNDN);
        mpfr_mul(r17169, r17167, r17168, MPFR_RNDN);
        mpfr_mul(r17170, r17166, r17169, MPFR_RNDN);
        mpfr_set_d(r17171, c, MPFR_RNDN);
        mpfr_set_d(r17172, b, MPFR_RNDN);
        mpfr_mul(r17173, r17171, r17172, MPFR_RNDN);
        mpfr_fma(r17174, r17165, r17170, r17173, MPFR_RNDN);
        ;
        mpfr_set_d(r17176, i, MPFR_RNDN);
        mpfr_set_d(r17177, a, MPFR_RNDN);
        mpfr_mul(r17178, r17163, r17177, MPFR_RNDN);
        mpfr_fma(r17179, r17176, r17167, r17178, MPFR_RNDN);
        mpfr_set_d(r17180, j, MPFR_RNDN);
        ;
        mpfr_mul(r17182, r17180, r17181, MPFR_RNDN);
        mpfr_set_d(r17183, k, MPFR_RNDN);
        mpfr_mul(r17184, r17182, r17183, MPFR_RNDN);
        mpfr_fma(r17185, r17175, r17179, r17184, MPFR_RNDN);
        mpfr_sub(r17186, r17174, r17185, MPFR_RNDN);
        return mpfr_get_d(r17186, MPFR_RNDN);
}

static mpfr_t r17187, r17188, r17189, r17190, r17191, r17192, r17193, r17194, r17195, r17196, r17197, r17198, r17199, r17200, r17201, r17202, r17203, r17204, r17205, r17206, r17207, r17208, r17209, r17210;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17187);
        mpfr_init(r17188);
        mpfr_init(r17189);
        mpfr_init_set_str(r17190, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17191);
        mpfr_init(r17192);
        mpfr_init(r17193);
        mpfr_init(r17194);
        mpfr_init(r17195);
        mpfr_init(r17196);
        mpfr_init(r17197);
        mpfr_init(r17198);
        mpfr_init_set_str(r17199, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17200);
        mpfr_init(r17201);
        mpfr_init(r17202);
        mpfr_init(r17203);
        mpfr_init(r17204);
        mpfr_init_set_str(r17205, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17206);
        mpfr_init(r17207);
        mpfr_init(r17208);
        mpfr_init(r17209);
        mpfr_init(r17210);
}

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(r17187, t, MPFR_RNDN);
        mpfr_set_d(r17188, z, MPFR_RNDN);
        mpfr_mul(r17189, r17187, r17188, MPFR_RNDN);
        ;
        mpfr_set_d(r17191, x, MPFR_RNDN);
        mpfr_set_d(r17192, y, MPFR_RNDN);
        mpfr_mul(r17193, r17191, r17192, MPFR_RNDN);
        mpfr_mul(r17194, r17190, r17193, MPFR_RNDN);
        mpfr_set_d(r17195, c, MPFR_RNDN);
        mpfr_set_d(r17196, b, MPFR_RNDN);
        mpfr_mul(r17197, r17195, r17196, MPFR_RNDN);
        mpfr_fma(r17198, r17189, r17194, r17197, MPFR_RNDN);
        ;
        mpfr_set_d(r17200, i, MPFR_RNDN);
        mpfr_set_d(r17201, a, MPFR_RNDN);
        mpfr_mul(r17202, r17187, r17201, MPFR_RNDN);
        mpfr_fma(r17203, r17200, r17191, r17202, MPFR_RNDN);
        mpfr_set_d(r17204, j, MPFR_RNDN);
        ;
        mpfr_mul(r17206, r17204, r17205, MPFR_RNDN);
        mpfr_set_d(r17207, k, MPFR_RNDN);
        mpfr_mul(r17208, r17206, r17207, MPFR_RNDN);
        mpfr_fma(r17209, r17199, r17203, r17208, MPFR_RNDN);
        mpfr_sub(r17210, r17198, r17209, MPFR_RNDN);
        return mpfr_get_d(r17210, MPFR_RNDN);
}