#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Diagrams.Solve.Polynomial:quadForm from diagrams-solve-0.1, B";

double f_if(float x, float y, float z) {
        float r34346 = 1.0;
        float r34347 = 2.0;
        float r34348 = r34346 / r34347;
        float r34349 = x;
        float r34350 = y;
        float r34351 = z;
        float r34352 = sqrt(r34351);
        float r34353 = r34350 * r34352;
        float r34354 = r34349 + r34353;
        float r34355 = r34348 * r34354;
        return r34355;
}

double f_id(double x, double y, double z) {
        double r34356 = 1.0;
        double r34357 = 2.0;
        double r34358 = r34356 / r34357;
        double r34359 = x;
        double r34360 = y;
        double r34361 = z;
        double r34362 = sqrt(r34361);
        double r34363 = r34360 * r34362;
        double r34364 = r34359 + r34363;
        double r34365 = r34358 * r34364;
        return r34365;
}


double f_of(float x, float y, float z) {
        float r34366 = 1.0;
        float r34367 = 2.0;
        float r34368 = r34366 / r34367;
        float r34369 = x;
        float r34370 = y;
        float r34371 = z;
        float r34372 = sqrt(r34371);
        float r34373 = r34370 * r34372;
        float r34374 = r34369 + r34373;
        float r34375 = r34368 * r34374;
        return r34375;
}

double f_od(double x, double y, double z) {
        double r34376 = 1.0;
        double r34377 = 2.0;
        double r34378 = r34376 / r34377;
        double r34379 = x;
        double r34380 = y;
        double r34381 = z;
        double r34382 = sqrt(r34381);
        double r34383 = r34380 * r34382;
        double r34384 = r34379 + r34383;
        double r34385 = r34378 * r34384;
        return r34385;
}

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 r34386, r34387, r34388, r34389, r34390, r34391, r34392, r34393, r34394, r34395;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r34386, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34387, "2.0", 10, MPFR_RNDN);
        mpfr_init(r34388);
        mpfr_init(r34389);
        mpfr_init(r34390);
        mpfr_init(r34391);
        mpfr_init(r34392);
        mpfr_init(r34393);
        mpfr_init(r34394);
        mpfr_init(r34395);
}

double f_im(double x, double y, double z) {
        ;
        ;
        mpfr_div(r34388, r34386, r34387, MPFR_RNDN);
        mpfr_set_d(r34389, x, MPFR_RNDN);
        mpfr_set_d(r34390, y, MPFR_RNDN);
        mpfr_set_d(r34391, z, MPFR_RNDN);
        mpfr_sqrt(r34392, r34391, MPFR_RNDN);
        mpfr_mul(r34393, r34390, r34392, MPFR_RNDN);
        mpfr_add(r34394, r34389, r34393, MPFR_RNDN);
        mpfr_mul(r34395, r34388, r34394, MPFR_RNDN);
        return mpfr_get_d(r34395, MPFR_RNDN);
}

static mpfr_t r34396, r34397, r34398, r34399, r34400, r34401, r34402, r34403, r34404, r34405;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r34396, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34397, "2.0", 10, MPFR_RNDN);
        mpfr_init(r34398);
        mpfr_init(r34399);
        mpfr_init(r34400);
        mpfr_init(r34401);
        mpfr_init(r34402);
        mpfr_init(r34403);
        mpfr_init(r34404);
        mpfr_init(r34405);
}

double f_fm(double x, double y, double z) {
        ;
        ;
        mpfr_div(r34398, r34396, r34397, MPFR_RNDN);
        mpfr_set_d(r34399, x, MPFR_RNDN);
        mpfr_set_d(r34400, y, MPFR_RNDN);
        mpfr_set_d(r34401, z, MPFR_RNDN);
        mpfr_sqrt(r34402, r34401, MPFR_RNDN);
        mpfr_mul(r34403, r34400, r34402, MPFR_RNDN);
        mpfr_add(r34404, r34399, r34403, MPFR_RNDN);
        mpfr_mul(r34405, r34398, r34404, MPFR_RNDN);
        return mpfr_get_d(r34405, MPFR_RNDN);
}

static mpfr_t r34406, r34407, r34408, r34409, r34410, r34411, r34412, r34413, r34414, r34415;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r34406, "1.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34407, "2.0", 10, MPFR_RNDN);
        mpfr_init(r34408);
        mpfr_init(r34409);
        mpfr_init(r34410);
        mpfr_init(r34411);
        mpfr_init(r34412);
        mpfr_init(r34413);
        mpfr_init(r34414);
        mpfr_init(r34415);
}

double f_dm(double x, double y, double z) {
        ;
        ;
        mpfr_div(r34408, r34406, r34407, MPFR_RNDN);
        mpfr_set_d(r34409, x, MPFR_RNDN);
        mpfr_set_d(r34410, y, MPFR_RNDN);
        mpfr_set_d(r34411, z, MPFR_RNDN);
        mpfr_sqrt(r34412, r34411, MPFR_RNDN);
        mpfr_mul(r34413, r34410, r34412, MPFR_RNDN);
        mpfr_add(r34414, r34409, r34413, MPFR_RNDN);
        mpfr_mul(r34415, r34408, r34414, MPFR_RNDN);
        return mpfr_get_d(r34415, MPFR_RNDN);
}