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

char *name = "Data.Metrics.Snapshot:quantile from metrics-0.3.0.2";

double f_if(float x, float y, float z, float t) {
        float r55029 = x;
        float r55030 = y;
        float r55031 = z;
        float r55032 = r55030 - r55031;
        float r55033 = t;
        float r55034 = r55033 - r55029;
        float r55035 = r55032 * r55034;
        float r55036 = r55029 + r55035;
        return r55036;
}

double f_id(double x, double y, double z, double t) {
        double r55037 = x;
        double r55038 = y;
        double r55039 = z;
        double r55040 = r55038 - r55039;
        double r55041 = t;
        double r55042 = r55041 - r55037;
        double r55043 = r55040 * r55042;
        double r55044 = r55037 + r55043;
        return r55044;
}


double f_of(float x, float y, float z, float t) {
        float r55045 = x;
        float r55046 = y;
        float r55047 = z;
        float r55048 = r55046 - r55047;
        float r55049 = t;
        float r55050 = r55049 - r55045;
        float r55051 = r55048 * r55050;
        float r55052 = r55045 + r55051;
        return r55052;
}

double f_od(double x, double y, double z, double t) {
        double r55053 = x;
        double r55054 = y;
        double r55055 = z;
        double r55056 = r55054 - r55055;
        double r55057 = t;
        double r55058 = r55057 - r55053;
        double r55059 = r55056 * r55058;
        double r55060 = r55053 + r55059;
        return r55060;
}

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 r55061, r55062, r55063, r55064, r55065, r55066, r55067, r55068;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r55061);
        mpfr_init(r55062);
        mpfr_init(r55063);
        mpfr_init(r55064);
        mpfr_init(r55065);
        mpfr_init(r55066);
        mpfr_init(r55067);
        mpfr_init(r55068);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r55061, x, MPFR_RNDN);
        mpfr_set_d(r55062, y, MPFR_RNDN);
        mpfr_set_d(r55063, z, MPFR_RNDN);
        mpfr_sub(r55064, r55062, r55063, MPFR_RNDN);
        mpfr_set_d(r55065, t, MPFR_RNDN);
        mpfr_sub(r55066, r55065, r55061, MPFR_RNDN);
        mpfr_mul(r55067, r55064, r55066, MPFR_RNDN);
        mpfr_add(r55068, r55061, r55067, MPFR_RNDN);
        return mpfr_get_d(r55068, MPFR_RNDN);
}

static mpfr_t r55069, r55070, r55071, r55072, r55073, r55074, r55075, r55076;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55069);
        mpfr_init(r55070);
        mpfr_init(r55071);
        mpfr_init(r55072);
        mpfr_init(r55073);
        mpfr_init(r55074);
        mpfr_init(r55075);
        mpfr_init(r55076);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r55069, x, MPFR_RNDN);
        mpfr_set_d(r55070, y, MPFR_RNDN);
        mpfr_set_d(r55071, z, MPFR_RNDN);
        mpfr_sub(r55072, r55070, r55071, MPFR_RNDN);
        mpfr_set_d(r55073, t, MPFR_RNDN);
        mpfr_sub(r55074, r55073, r55069, MPFR_RNDN);
        mpfr_mul(r55075, r55072, r55074, MPFR_RNDN);
        mpfr_add(r55076, r55069, r55075, MPFR_RNDN);
        return mpfr_get_d(r55076, MPFR_RNDN);
}

static mpfr_t r55077, r55078, r55079, r55080, r55081, r55082, r55083, r55084;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55077);
        mpfr_init(r55078);
        mpfr_init(r55079);
        mpfr_init(r55080);
        mpfr_init(r55081);
        mpfr_init(r55082);
        mpfr_init(r55083);
        mpfr_init(r55084);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r55077, x, MPFR_RNDN);
        mpfr_set_d(r55078, y, MPFR_RNDN);
        mpfr_set_d(r55079, z, MPFR_RNDN);
        mpfr_sub(r55080, r55078, r55079, MPFR_RNDN);
        mpfr_set_d(r55081, t, MPFR_RNDN);
        mpfr_sub(r55082, r55081, r55077, MPFR_RNDN);
        mpfr_mul(r55083, r55080, r55082, MPFR_RNDN);
        mpfr_add(r55084, r55077, r55083, MPFR_RNDN);
        return mpfr_get_d(r55084, MPFR_RNDN);
}