#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 r17347 = x;
        float r17348 = 18.0f;
        float r17349 = r17347 * r17348;
        float r17350 = y;
        float r17351 = r17349 * r17350;
        float r17352 = z;
        float r17353 = r17351 * r17352;
        float r17354 = t;
        float r17355 = r17353 * r17354;
        float r17356 = a;
        float r17357 = 4.0f;
        float r17358 = r17356 * r17357;
        float r17359 = r17358 * r17354;
        float r17360 = r17355 - r17359;
        float r17361 = b;
        float r17362 = c;
        float r17363 = r17361 * r17362;
        float r17364 = r17360 + r17363;
        float r17365 = r17347 * r17357;
        float r17366 = i;
        float r17367 = r17365 * r17366;
        float r17368 = r17364 - r17367;
        float r17369 = j;
        float r17370 = 27.0f;
        float r17371 = r17369 * r17370;
        float r17372 = k;
        float r17373 = r17371 * r17372;
        float r17374 = r17368 - r17373;
        return r17374;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17375 = x;
        double r17376 = 18.0;
        double r17377 = r17375 * r17376;
        double r17378 = y;
        double r17379 = r17377 * r17378;
        double r17380 = z;
        double r17381 = r17379 * r17380;
        double r17382 = t;
        double r17383 = r17381 * r17382;
        double r17384 = a;
        double r17385 = 4.0;
        double r17386 = r17384 * r17385;
        double r17387 = r17386 * r17382;
        double r17388 = r17383 - r17387;
        double r17389 = b;
        double r17390 = c;
        double r17391 = r17389 * r17390;
        double r17392 = r17388 + r17391;
        double r17393 = r17375 * r17385;
        double r17394 = i;
        double r17395 = r17393 * r17394;
        double r17396 = r17392 - r17395;
        double r17397 = j;
        double r17398 = 27.0;
        double r17399 = r17397 * r17398;
        double r17400 = k;
        double r17401 = r17399 * r17400;
        double r17402 = r17396 - r17401;
        return r17402;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) {
        float r17403 = t;
        float r17404 = z;
        float r17405 = r17403 * r17404;
        float r17406 = 18.0f;
        float r17407 = y;
        float r17408 = x;
        float r17409 = r17407 * r17408;
        float r17410 = r17406 * r17409;
        float r17411 = b;
        float r17412 = c;
        float r17413 = r17411 * r17412;
        float r17414 = fma(r17405, r17410, r17413);
        float r17415 = 4.0f;
        float r17416 = i;
        float r17417 = a;
        float r17418 = r17403 * r17417;
        float r17419 = fma(r17416, r17408, r17418);
        float r17420 = 27.0f;
        float r17421 = j;
        float r17422 = k;
        float r17423 = r17421 * r17422;
        float r17424 = r17420 * r17423;
        float r17425 = fma(r17415, r17419, r17424);
        float r17426 = r17414 - r17425;
        return r17426;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r17427 = t;
        double r17428 = z;
        double r17429 = r17427 * r17428;
        double r17430 = 18.0;
        double r17431 = y;
        double r17432 = x;
        double r17433 = r17431 * r17432;
        double r17434 = r17430 * r17433;
        double r17435 = b;
        double r17436 = c;
        double r17437 = r17435 * r17436;
        double r17438 = fma(r17429, r17434, r17437);
        double r17439 = 4.0;
        double r17440 = i;
        double r17441 = a;
        double r17442 = r17427 * r17441;
        double r17443 = fma(r17440, r17432, r17442);
        double r17444 = 27.0;
        double r17445 = j;
        double r17446 = k;
        double r17447 = r17445 * r17446;
        double r17448 = r17444 * r17447;
        double r17449 = fma(r17439, r17443, r17448);
        double r17450 = r17438 - r17449;
        return r17450;
}

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 r17451, r17452, r17453, r17454, r17455, r17456, r17457, r17458, r17459, r17460, r17461, r17462, r17463, r17464, r17465, r17466, r17467, r17468, r17469, r17470, r17471, r17472, r17473, r17474, r17475, r17476, r17477, r17478;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17451);
        mpfr_init_set_str(r17452, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17453);
        mpfr_init(r17454);
        mpfr_init(r17455);
        mpfr_init(r17456);
        mpfr_init(r17457);
        mpfr_init(r17458);
        mpfr_init(r17459);
        mpfr_init(r17460);
        mpfr_init_set_str(r17461, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17462);
        mpfr_init(r17463);
        mpfr_init(r17464);
        mpfr_init(r17465);
        mpfr_init(r17466);
        mpfr_init(r17467);
        mpfr_init(r17468);
        mpfr_init(r17469);
        mpfr_init(r17470);
        mpfr_init(r17471);
        mpfr_init(r17472);
        mpfr_init(r17473);
        mpfr_init_set_str(r17474, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17475);
        mpfr_init(r17476);
        mpfr_init(r17477);
        mpfr_init(r17478);
}

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(r17451, x, MPFR_RNDN);
        ;
        mpfr_mul(r17453, r17451, r17452, MPFR_RNDN);
        mpfr_set_d(r17454, y, MPFR_RNDN);
        mpfr_mul(r17455, r17453, r17454, MPFR_RNDN);
        mpfr_set_d(r17456, z, MPFR_RNDN);
        mpfr_mul(r17457, r17455, r17456, MPFR_RNDN);
        mpfr_set_d(r17458, t, MPFR_RNDN);
        mpfr_mul(r17459, r17457, r17458, MPFR_RNDN);
        mpfr_set_d(r17460, a, MPFR_RNDN);
        ;
        mpfr_mul(r17462, r17460, r17461, MPFR_RNDN);
        mpfr_mul(r17463, r17462, r17458, MPFR_RNDN);
        mpfr_sub(r17464, r17459, r17463, MPFR_RNDN);
        mpfr_set_d(r17465, b, MPFR_RNDN);
        mpfr_set_d(r17466, c, MPFR_RNDN);
        mpfr_mul(r17467, r17465, r17466, MPFR_RNDN);
        mpfr_add(r17468, r17464, r17467, MPFR_RNDN);
        mpfr_mul(r17469, r17451, r17461, MPFR_RNDN);
        mpfr_set_d(r17470, i, MPFR_RNDN);
        mpfr_mul(r17471, r17469, r17470, MPFR_RNDN);
        mpfr_sub(r17472, r17468, r17471, MPFR_RNDN);
        mpfr_set_d(r17473, j, MPFR_RNDN);
        ;
        mpfr_mul(r17475, r17473, r17474, MPFR_RNDN);
        mpfr_set_d(r17476, k, MPFR_RNDN);
        mpfr_mul(r17477, r17475, r17476, MPFR_RNDN);
        mpfr_sub(r17478, r17472, r17477, MPFR_RNDN);
        return mpfr_get_d(r17478, MPFR_RNDN);
}

static mpfr_t r17479, r17480, r17481, r17482, r17483, r17484, r17485, r17486, r17487, r17488, r17489, r17490, r17491, r17492, r17493, r17494, r17495, r17496, r17497, r17498, r17499, r17500, r17501, r17502;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17479);
        mpfr_init(r17480);
        mpfr_init(r17481);
        mpfr_init_set_str(r17482, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17483);
        mpfr_init(r17484);
        mpfr_init(r17485);
        mpfr_init(r17486);
        mpfr_init(r17487);
        mpfr_init(r17488);
        mpfr_init(r17489);
        mpfr_init(r17490);
        mpfr_init_set_str(r17491, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17492);
        mpfr_init(r17493);
        mpfr_init(r17494);
        mpfr_init(r17495);
        mpfr_init_set_str(r17496, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17497);
        mpfr_init(r17498);
        mpfr_init(r17499);
        mpfr_init(r17500);
        mpfr_init(r17501);
        mpfr_init(r17502);
}

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(r17479, t, MPFR_RNDN);
        mpfr_set_d(r17480, z, MPFR_RNDN);
        mpfr_mul(r17481, r17479, r17480, MPFR_RNDN);
        ;
        mpfr_set_d(r17483, y, MPFR_RNDN);
        mpfr_set_d(r17484, x, MPFR_RNDN);
        mpfr_mul(r17485, r17483, r17484, MPFR_RNDN);
        mpfr_mul(r17486, r17482, r17485, MPFR_RNDN);
        mpfr_set_d(r17487, b, MPFR_RNDN);
        mpfr_set_d(r17488, c, MPFR_RNDN);
        mpfr_mul(r17489, r17487, r17488, MPFR_RNDN);
        mpfr_fma(r17490, r17481, r17486, r17489, MPFR_RNDN);
        ;
        mpfr_set_d(r17492, i, MPFR_RNDN);
        mpfr_set_d(r17493, a, MPFR_RNDN);
        mpfr_mul(r17494, r17479, r17493, MPFR_RNDN);
        mpfr_fma(r17495, r17492, r17484, r17494, MPFR_RNDN);
        ;
        mpfr_set_d(r17497, j, MPFR_RNDN);
        mpfr_set_d(r17498, k, MPFR_RNDN);
        mpfr_mul(r17499, r17497, r17498, MPFR_RNDN);
        mpfr_mul(r17500, r17496, r17499, MPFR_RNDN);
        mpfr_fma(r17501, r17491, r17495, r17500, MPFR_RNDN);
        mpfr_sub(r17502, r17490, r17501, MPFR_RNDN);
        return mpfr_get_d(r17502, MPFR_RNDN);
}

static mpfr_t r17503, r17504, r17505, r17506, r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515, r17516, r17517, r17518, r17519, r17520, r17521, r17522, r17523, r17524, r17525, r17526;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17503);
        mpfr_init(r17504);
        mpfr_init(r17505);
        mpfr_init_set_str(r17506, "18.0", 10, MPFR_RNDN);
        mpfr_init(r17507);
        mpfr_init(r17508);
        mpfr_init(r17509);
        mpfr_init(r17510);
        mpfr_init(r17511);
        mpfr_init(r17512);
        mpfr_init(r17513);
        mpfr_init(r17514);
        mpfr_init_set_str(r17515, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17516);
        mpfr_init(r17517);
        mpfr_init(r17518);
        mpfr_init(r17519);
        mpfr_init_set_str(r17520, "27.0", 10, MPFR_RNDN);
        mpfr_init(r17521);
        mpfr_init(r17522);
        mpfr_init(r17523);
        mpfr_init(r17524);
        mpfr_init(r17525);
        mpfr_init(r17526);
}

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(r17503, t, MPFR_RNDN);
        mpfr_set_d(r17504, z, MPFR_RNDN);
        mpfr_mul(r17505, r17503, r17504, MPFR_RNDN);
        ;
        mpfr_set_d(r17507, y, MPFR_RNDN);
        mpfr_set_d(r17508, x, MPFR_RNDN);
        mpfr_mul(r17509, r17507, r17508, MPFR_RNDN);
        mpfr_mul(r17510, r17506, r17509, MPFR_RNDN);
        mpfr_set_d(r17511, b, MPFR_RNDN);
        mpfr_set_d(r17512, c, MPFR_RNDN);
        mpfr_mul(r17513, r17511, r17512, MPFR_RNDN);
        mpfr_fma(r17514, r17505, r17510, r17513, MPFR_RNDN);
        ;
        mpfr_set_d(r17516, i, MPFR_RNDN);
        mpfr_set_d(r17517, a, MPFR_RNDN);
        mpfr_mul(r17518, r17503, r17517, MPFR_RNDN);
        mpfr_fma(r17519, r17516, r17508, r17518, MPFR_RNDN);
        ;
        mpfr_set_d(r17521, j, MPFR_RNDN);
        mpfr_set_d(r17522, k, MPFR_RNDN);
        mpfr_mul(r17523, r17521, r17522, MPFR_RNDN);
        mpfr_mul(r17524, r17520, r17523, MPFR_RNDN);
        mpfr_fma(r17525, r17515, r17519, r17524, MPFR_RNDN);
        mpfr_sub(r17526, r17514, r17525, MPFR_RNDN);
        return mpfr_get_d(r17526, MPFR_RNDN);
}