#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 r17050 = x;
        float r17051 = 18.0f;
        float r17052 = r17050 * r17051;
        float r17053 = y;
        float r17054 = r17052 * r17053;
        float r17055 = z;
        float r17056 = r17054 * r17055;
        float r17057 = t;
        float r17058 = r17056 * r17057;
        float r17059 = a;
        float r17060 = 4.0f;
        float r17061 = r17059 * r17060;
        float r17062 = r17061 * r17057;
        float r17063 = r17058 - r17062;
        float r17064 = b;
        float r17065 = c;
        float r17066 = r17064 * r17065;
        float r17067 = r17063 + r17066;
        float r17068 = r17050 * r17060;
        float r17069 = i;
        float r17070 = r17068 * r17069;
        float r17071 = r17067 - r17070;
        float r17072 = j;
        float r17073 = 27.0f;
        float r17074 = r17072 * r17073;
        float r17075 = k;
        float r17076 = r17074 * r17075;
        float r17077 = r17071 - r17076;
        return r17077;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17078 = x;
        double r17079 = 18.0;
        double r17080 = r17078 * r17079;
        double r17081 = y;
        double r17082 = r17080 * r17081;
        double r17083 = z;
        double r17084 = r17082 * r17083;
        double r17085 = t;
        double r17086 = r17084 * r17085;
        double r17087 = a;
        double r17088 = 4.0;
        double r17089 = r17087 * r17088;
        double r17090 = r17089 * r17085;
        double r17091 = r17086 - r17090;
        double r17092 = b;
        double r17093 = c;
        double r17094 = r17092 * r17093;
        double r17095 = r17091 + r17094;
        double r17096 = r17078 * r17088;
        double r17097 = i;
        double r17098 = r17096 * r17097;
        double r17099 = r17095 - r17098;
        double r17100 = j;
        double r17101 = 27.0;
        double r17102 = r17100 * r17101;
        double r17103 = k;
        double r17104 = r17102 * r17103;
        double r17105 = r17099 - r17104;
        return r17105;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17106 = t;
        float r17107 = z;
        float r17108 = r17106 * r17107;
        float r17109 = 18.0f;
        float r17110 = x;
        float r17111 = y;
        float r17112 = r17110 * r17111;
        float r17113 = r17109 * r17112;
        float r17114 = c;
        float r17115 = b;
        float r17116 = r17114 * r17115;
        float r17117 = fma(r17108, r17113, r17116);
        float r17118 = 4.0f;
        float r17119 = i;
        float r17120 = a;
        float r17121 = r17106 * r17120;
        float r17122 = fma(r17119, r17110, r17121);
        float r17123 = j;
        float r17124 = 27.0f;
        float r17125 = r17123 * r17124;
        float r17126 = k;
        float r17127 = r17125 * r17126;
        float r17128 = fma(r17118, r17122, r17127);
        float r17129 = r17117 - r17128;
        return r17129;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17130 = t;
        double r17131 = z;
        double r17132 = r17130 * r17131;
        double r17133 = 18.0;
        double r17134 = x;
        double r17135 = y;
        double r17136 = r17134 * r17135;
        double r17137 = r17133 * r17136;
        double r17138 = c;
        double r17139 = b;
        double r17140 = r17138 * r17139;
        double r17141 = fma(r17132, r17137, r17140);
        double r17142 = 4.0;
        double r17143 = i;
        double r17144 = a;
        double r17145 = r17130 * r17144;
        double r17146 = fma(r17143, r17134, r17145);
        double r17147 = j;
        double r17148 = 27.0;
        double r17149 = r17147 * r17148;
        double r17150 = k;
        double r17151 = r17149 * r17150;
        double r17152 = fma(r17142, r17146, r17151);
        double r17153 = r17141 - r17152;
        return r17153;
}

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 r17154, r17155, r17156, r17157, r17158, r17159, r17160, r17161, r17162, r17163, r17164, r17165, r17166, r17167, r17168, r17169, r17170, r17171, r17172, r17173, r17174, r17175, r17176, r17177, r17178, r17179, r17180, r17181;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17154);
        mpfr_init_set_str(r17155, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17156);
        mpfr_init(r17157);
        mpfr_init(r17158);
        mpfr_init(r17159);
        mpfr_init(r17160);
        mpfr_init(r17161);
        mpfr_init(r17162);
        mpfr_init(r17163);
        mpfr_init_set_str(r17164, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17165);
        mpfr_init(r17166);
        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(r17175);
        mpfr_init(r17176);
        mpfr_init_set_str(r17177, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17178);
        mpfr_init(r17179);
        mpfr_init(r17180);
        mpfr_init(r17181);
}

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(r17154, x, MPFR_RNDN);
        ;
        mpfr_mul(r17156, r17154, r17155, MPFR_RNDN);
        mpfr_set_d(r17157, y, MPFR_RNDN);
        mpfr_mul(r17158, r17156, r17157, MPFR_RNDN);
        mpfr_set_d(r17159, z, MPFR_RNDN);
        mpfr_mul(r17160, r17158, r17159, MPFR_RNDN);
        mpfr_set_d(r17161, t, MPFR_RNDN);
        mpfr_mul(r17162, r17160, r17161, MPFR_RNDN);
        mpfr_set_d(r17163, a, MPFR_RNDN);
        ;
        mpfr_mul(r17165, r17163, r17164, MPFR_RNDN);
        mpfr_mul(r17166, r17165, r17161, MPFR_RNDN);
        mpfr_sub(r17167, r17162, r17166, MPFR_RNDN);
        mpfr_set_d(r17168, b, MPFR_RNDN);
        mpfr_set_d(r17169, c, MPFR_RNDN);
        mpfr_mul(r17170, r17168, r17169, MPFR_RNDN);
        mpfr_add(r17171, r17167, r17170, MPFR_RNDN);
        mpfr_mul(r17172, r17154, r17164, MPFR_RNDN);
        mpfr_set_d(r17173, i, MPFR_RNDN);
        mpfr_mul(r17174, r17172, r17173, MPFR_RNDN);
        mpfr_sub(r17175, r17171, r17174, MPFR_RNDN);
        mpfr_set_d(r17176, j, MPFR_RNDN);
        ;
        mpfr_mul(r17178, r17176, r17177, MPFR_RNDN);
        mpfr_set_d(r17179, k, MPFR_RNDN);
        mpfr_mul(r17180, r17178, r17179, MPFR_RNDN);
        mpfr_sub(r17181, r17175, r17180, MPFR_RNDN);
        return mpfr_get_d(r17181, MPFR_RNDN);
}

static mpfr_t r17182, r17183, r17184, r17185, r17186, r17187, r17188, r17189, r17190, r17191, r17192, r17193, r17194, r17195, r17196, r17197, r17198, r17199, r17200, r17201, r17202, r17203, r17204, r17205;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17182);
        mpfr_init(r17183);
        mpfr_init(r17184);
        mpfr_init_set_str(r17185, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17186);
        mpfr_init(r17187);
        mpfr_init(r17188);
        mpfr_init(r17189);
        mpfr_init(r17190);
        mpfr_init(r17191);
        mpfr_init(r17192);
        mpfr_init(r17193);
        mpfr_init_set_str(r17194, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17195);
        mpfr_init(r17196);
        mpfr_init(r17197);
        mpfr_init(r17198);
        mpfr_init(r17199);
        mpfr_init_set_str(r17200, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17201);
        mpfr_init(r17202);
        mpfr_init(r17203);
        mpfr_init(r17204);
        mpfr_init(r17205);
}

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

static mpfr_t r17206, r17207, r17208, r17209, r17210, r17211, r17212, r17213, r17214, r17215, r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17206);
        mpfr_init(r17207);
        mpfr_init(r17208);
        mpfr_init_set_str(r17209, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17210);
        mpfr_init(r17211);
        mpfr_init(r17212);
        mpfr_init(r17213);
        mpfr_init(r17214);
        mpfr_init(r17215);
        mpfr_init(r17216);
        mpfr_init(r17217);
        mpfr_init_set_str(r17218, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17219);
        mpfr_init(r17220);
        mpfr_init(r17221);
        mpfr_init(r17222);
        mpfr_init(r17223);
        mpfr_init_set_str(r17224, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17225);
        mpfr_init(r17226);
        mpfr_init(r17227);
        mpfr_init(r17228);
        mpfr_init(r17229);
}

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(r17206, t, MPFR_RNDN);
        mpfr_set_d(r17207, z, MPFR_RNDN);
        mpfr_mul(r17208, r17206, r17207, MPFR_RNDN);
        ;
        mpfr_set_d(r17210, x, MPFR_RNDN);
        mpfr_set_d(r17211, y, MPFR_RNDN);
        mpfr_mul(r17212, r17210, r17211, MPFR_RNDN);
        mpfr_mul(r17213, r17209, r17212, MPFR_RNDN);
        mpfr_set_d(r17214, c, MPFR_RNDN);
        mpfr_set_d(r17215, b, MPFR_RNDN);
        mpfr_mul(r17216, r17214, r17215, MPFR_RNDN);
        mpfr_fma(r17217, r17208, r17213, r17216, MPFR_RNDN);
        ;
        mpfr_set_d(r17219, i, MPFR_RNDN);
        mpfr_set_d(r17220, a, MPFR_RNDN);
        mpfr_mul(r17221, r17206, r17220, MPFR_RNDN);
        mpfr_fma(r17222, r17219, r17210, r17221, MPFR_RNDN);
        mpfr_set_d(r17223, j, MPFR_RNDN);
        ;
        mpfr_mul(r17225, r17223, r17224, MPFR_RNDN);
        mpfr_set_d(r17226, k, MPFR_RNDN);
        mpfr_mul(r17227, r17225, r17226, MPFR_RNDN);
        mpfr_fma(r17228, r17218, r17222, r17227, MPFR_RNDN);
        mpfr_sub(r17229, r17217, r17228, MPFR_RNDN);
        return mpfr_get_d(r17229, MPFR_RNDN);
}