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

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

double f_if(float x, float y) {
        float r56661 = x;
        float r56662 = y;
        float r56663 = r56661 * r56662;
        float r56664 = r56663 - r56661;
        return r56664;
}

double f_id(double x, double y) {
        double r56665 = x;
        double r56666 = y;
        double r56667 = r56665 * r56666;
        double r56668 = r56667 - r56665;
        return r56668;
}


double f_of(float x, float y) {
        float r56669 = x;
        float r56670 = y;
        float r56671 = r56669 * r56670;
        float r56672 = r56671 - r56669;
        return r56672;
}

double f_od(double x, double y) {
        double r56673 = x;
        double r56674 = y;
        double r56675 = r56673 * r56674;
        double r56676 = r56675 - r56673;
        return r56676;
}

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 r56677, r56678, r56679, r56680;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56677);
        mpfr_init(r56678);
        mpfr_init(r56679);
        mpfr_init(r56680);
}

double f_im(double x, double y) {
        mpfr_set_d(r56677, x, MPFR_RNDN);
        mpfr_set_d(r56678, y, MPFR_RNDN);
        mpfr_mul(r56679, r56677, r56678, MPFR_RNDN);
        mpfr_sub(r56680, r56679, r56677, MPFR_RNDN);
        return mpfr_get_d(r56680, MPFR_RNDN);
}

static mpfr_t r56681, r56682, r56683, r56684;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56681);
        mpfr_init(r56682);
        mpfr_init(r56683);
        mpfr_init(r56684);
}

double f_fm(double x, double y) {
        mpfr_set_d(r56681, x, MPFR_RNDN);
        mpfr_set_d(r56682, y, MPFR_RNDN);
        mpfr_mul(r56683, r56681, r56682, MPFR_RNDN);
        mpfr_sub(r56684, r56683, r56681, MPFR_RNDN);
        return mpfr_get_d(r56684, MPFR_RNDN);
}

static mpfr_t r56685, r56686, r56687, r56688;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56685);
        mpfr_init(r56686);
        mpfr_init(r56687);
        mpfr_init(r56688);
}

double f_dm(double x, double y) {
        mpfr_set_d(r56685, x, MPFR_RNDN);
        mpfr_set_d(r56686, y, MPFR_RNDN);
        mpfr_mul(r56687, r56685, r56686, MPFR_RNDN);
        mpfr_sub(r56688, r56687, r56685, MPFR_RNDN);
        return mpfr_get_d(r56688, MPFR_RNDN);
}