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

char *name = "Data.Histogram.Bin.BinF:$cfromIndex from histogram-fill-0.8.4.1";

double f_if(float x, float y, float z) {
        float r56541 = x;
        float r56542 = 2.0;
        float r56543 = r56541 / r56542;
        float r56544 = y;
        float r56545 = r56544 * r56541;
        float r56546 = r56543 + r56545;
        float r56547 = z;
        float r56548 = r56546 + r56547;
        return r56548;
}

double f_id(double x, double y, double z) {
        double r56549 = x;
        double r56550 = 2.0;
        double r56551 = r56549 / r56550;
        double r56552 = y;
        double r56553 = r56552 * r56549;
        double r56554 = r56551 + r56553;
        double r56555 = z;
        double r56556 = r56554 + r56555;
        return r56556;
}


double f_of(float x, float y, float z) {
        float r56557 = x;
        float r56558 = 2.0;
        float r56559 = r56557 / r56558;
        float r56560 = y;
        float r56561 = r56560 * r56557;
        float r56562 = r56559 + r56561;
        float r56563 = z;
        float r56564 = r56562 + r56563;
        return r56564;
}

double f_od(double x, double y, double z) {
        double r56565 = x;
        double r56566 = 2.0;
        double r56567 = r56565 / r56566;
        double r56568 = y;
        double r56569 = r56568 * r56565;
        double r56570 = r56567 + r56569;
        double r56571 = z;
        double r56572 = r56570 + r56571;
        return r56572;
}

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 r56573, r56574, r56575, r56576, r56577, r56578, r56579, r56580;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56573);
        mpfr_init_set_str(r56574, "2.0", 10, MPFR_RNDN);
        mpfr_init(r56575);
        mpfr_init(r56576);
        mpfr_init(r56577);
        mpfr_init(r56578);
        mpfr_init(r56579);
        mpfr_init(r56580);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r56573, x, MPFR_RNDN);
        ;
        mpfr_div(r56575, r56573, r56574, MPFR_RNDN);
        mpfr_set_d(r56576, y, MPFR_RNDN);
        mpfr_mul(r56577, r56576, r56573, MPFR_RNDN);
        mpfr_add(r56578, r56575, r56577, MPFR_RNDN);
        mpfr_set_d(r56579, z, MPFR_RNDN);
        mpfr_add(r56580, r56578, r56579, MPFR_RNDN);
        return mpfr_get_d(r56580, MPFR_RNDN);
}

static mpfr_t r56581, r56582, r56583, r56584, r56585, r56586, r56587, r56588;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56581);
        mpfr_init_set_str(r56582, "2.0", 10, MPFR_RNDN);
        mpfr_init(r56583);
        mpfr_init(r56584);
        mpfr_init(r56585);
        mpfr_init(r56586);
        mpfr_init(r56587);
        mpfr_init(r56588);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r56581, x, MPFR_RNDN);
        ;
        mpfr_div(r56583, r56581, r56582, MPFR_RNDN);
        mpfr_set_d(r56584, y, MPFR_RNDN);
        mpfr_mul(r56585, r56584, r56581, MPFR_RNDN);
        mpfr_add(r56586, r56583, r56585, MPFR_RNDN);
        mpfr_set_d(r56587, z, MPFR_RNDN);
        mpfr_add(r56588, r56586, r56587, MPFR_RNDN);
        return mpfr_get_d(r56588, MPFR_RNDN);
}

static mpfr_t r56589, r56590, r56591, r56592, r56593, r56594, r56595, r56596;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56589);
        mpfr_init_set_str(r56590, "2.0", 10, MPFR_RNDN);
        mpfr_init(r56591);
        mpfr_init(r56592);
        mpfr_init(r56593);
        mpfr_init(r56594);
        mpfr_init(r56595);
        mpfr_init(r56596);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r56589, x, MPFR_RNDN);
        ;
        mpfr_div(r56591, r56589, r56590, MPFR_RNDN);
        mpfr_set_d(r56592, y, MPFR_RNDN);
        mpfr_mul(r56593, r56592, r56589, MPFR_RNDN);
        mpfr_add(r56594, r56591, r56593, MPFR_RNDN);
        mpfr_set_d(r56595, z, MPFR_RNDN);
        mpfr_add(r56596, r56594, r56595, MPFR_RNDN);
        return mpfr_get_d(r56596, MPFR_RNDN);
}

