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

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

double f_if(float x, float y) {
        float r34626 = x;
        float r34627 = 3.0;
        float r34628 = 8.0;
        float r34629 = r34627 / r34628;
        float r34630 = y;
        float r34631 = r34629 * r34630;
        float r34632 = r34626 - r34631;
        return r34632;
}

double f_id(double x, double y) {
        double r34633 = x;
        double r34634 = 3.0;
        double r34635 = 8.0;
        double r34636 = r34634 / r34635;
        double r34637 = y;
        double r34638 = r34636 * r34637;
        double r34639 = r34633 - r34638;
        return r34639;
}


double f_of(float x, float y) {
        float r34640 = x;
        float r34641 = 3.0;
        float r34642 = 8.0;
        float r34643 = r34641 / r34642;
        float r34644 = y;
        float r34645 = r34643 * r34644;
        float r34646 = r34640 - r34645;
        return r34646;
}

double f_od(double x, double y) {
        double r34647 = x;
        double r34648 = 3.0;
        double r34649 = 8.0;
        double r34650 = r34648 / r34649;
        double r34651 = y;
        double r34652 = r34650 * r34651;
        double r34653 = r34647 - r34652;
        return r34653;
}

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 r34654, r34655, r34656, r34657, r34658, r34659, r34660;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r34654);
        mpfr_init_set_str(r34655, "3.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34656, "8.0", 10, MPFR_RNDN);
        mpfr_init(r34657);
        mpfr_init(r34658);
        mpfr_init(r34659);
        mpfr_init(r34660);
}

double f_im(double x, double y) {
        mpfr_set_d(r34654, x, MPFR_RNDN);
        ;
        ;
        mpfr_div(r34657, r34655, r34656, MPFR_RNDN);
        mpfr_set_d(r34658, y, MPFR_RNDN);
        mpfr_mul(r34659, r34657, r34658, MPFR_RNDN);
        mpfr_sub(r34660, r34654, r34659, MPFR_RNDN);
        return mpfr_get_d(r34660, MPFR_RNDN);
}

static mpfr_t r34661, r34662, r34663, r34664, r34665, r34666, r34667;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r34661);
        mpfr_init_set_str(r34662, "3.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34663, "8.0", 10, MPFR_RNDN);
        mpfr_init(r34664);
        mpfr_init(r34665);
        mpfr_init(r34666);
        mpfr_init(r34667);
}

double f_fm(double x, double y) {
        mpfr_set_d(r34661, x, MPFR_RNDN);
        ;
        ;
        mpfr_div(r34664, r34662, r34663, MPFR_RNDN);
        mpfr_set_d(r34665, y, MPFR_RNDN);
        mpfr_mul(r34666, r34664, r34665, MPFR_RNDN);
        mpfr_sub(r34667, r34661, r34666, MPFR_RNDN);
        return mpfr_get_d(r34667, MPFR_RNDN);
}

static mpfr_t r34668, r34669, r34670, r34671, r34672, r34673, r34674;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r34668);
        mpfr_init_set_str(r34669, "3.0", 10, MPFR_RNDN);
        mpfr_init_set_str(r34670, "8.0", 10, MPFR_RNDN);
        mpfr_init(r34671);
        mpfr_init(r34672);
        mpfr_init(r34673);
        mpfr_init(r34674);
}

double f_dm(double x, double y) {
        mpfr_set_d(r34668, x, MPFR_RNDN);
        ;
        ;
        mpfr_div(r34671, r34669, r34670, MPFR_RNDN);
        mpfr_set_d(r34672, y, MPFR_RNDN);
        mpfr_mul(r34673, r34671, r34672, MPFR_RNDN);
        mpfr_sub(r34674, r34668, r34673, MPFR_RNDN);
        return mpfr_get_d(r34674, MPFR_RNDN);
}