#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 r56385 = x;
        float r56386 = y;
        float r56387 = r56385 * r56386;
        float r56388 = r56387 - r56385;
        return r56388;
}

double f_id(double x, double y) {
        double r56389 = x;
        double r56390 = y;
        double r56391 = r56389 * r56390;
        double r56392 = r56391 - r56389;
        return r56392;
}


double f_of(float x, float y) {
        float r56393 = x;
        float r56394 = y;
        float r56395 = r56393 * r56394;
        float r56396 = r56395 - r56393;
        return r56396;
}

double f_od(double x, double y) {
        double r56397 = x;
        double r56398 = y;
        double r56399 = r56397 * r56398;
        double r56400 = r56399 - r56397;
        return r56400;
}

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 r56401, r56402, r56403, r56404;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56401);
        mpfr_init(r56402);
        mpfr_init(r56403);
        mpfr_init(r56404);
}

double f_im(double x, double y) {
        mpfr_set_d(r56401, x, MPFR_RNDN);
        mpfr_set_d(r56402, y, MPFR_RNDN);
        mpfr_mul(r56403, r56401, r56402, MPFR_RNDN);
        mpfr_sub(r56404, r56403, r56401, MPFR_RNDN);
        return mpfr_get_d(r56404, MPFR_RNDN);
}

static mpfr_t r56405, r56406, r56407, r56408;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56405);
        mpfr_init(r56406);
        mpfr_init(r56407);
        mpfr_init(r56408);
}

double f_fm(double x, double y) {
        mpfr_set_d(r56405, x, MPFR_RNDN);
        mpfr_set_d(r56406, y, MPFR_RNDN);
        mpfr_mul(r56407, r56405, r56406, MPFR_RNDN);
        mpfr_sub(r56408, r56407, r56405, MPFR_RNDN);
        return mpfr_get_d(r56408, MPFR_RNDN);
}

static mpfr_t r56409, r56410, r56411, r56412;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56409);
        mpfr_init(r56410);
        mpfr_init(r56411);
        mpfr_init(r56412);
}

double f_dm(double x, double y) {
        mpfr_set_d(r56409, x, MPFR_RNDN);
        mpfr_set_d(r56410, y, MPFR_RNDN);
        mpfr_mul(r56411, r56409, r56410, MPFR_RNDN);
        mpfr_sub(r56412, r56411, r56409, MPFR_RNDN);
        return mpfr_get_d(r56412, MPFR_RNDN);
}