#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 r54985 = x;
        float r54986 = y;
        float r54987 = z;
        float r54988 = r54986 - r54987;
        float r54989 = t;
        float r54990 = r54989 - r54985;
        float r54991 = r54988 * r54990;
        float r54992 = r54985 + r54991;
        return r54992;
}

double f_id(double x, double y, double z, double t) {
        double r54993 = x;
        double r54994 = y;
        double r54995 = z;
        double r54996 = r54994 - r54995;
        double r54997 = t;
        double r54998 = r54997 - r54993;
        double r54999 = r54996 * r54998;
        double r55000 = r54993 + r54999;
        return r55000;
}


double f_of(float x, float y, float z, float t) {
        float r55001 = x;
        float r55002 = y;
        float r55003 = z;
        float r55004 = r55002 - r55003;
        float r55005 = t;
        float r55006 = r55004 * r55005;
        float r55007 = -r55001;
        float r55008 = r55004 * r55007;
        float r55009 = r55006 + r55008;
        float r55010 = r55001 + r55009;
        return r55010;
}

double f_od(double x, double y, double z, double t) {
        double r55011 = x;
        double r55012 = y;
        double r55013 = z;
        double r55014 = r55012 - r55013;
        double r55015 = t;
        double r55016 = r55014 * r55015;
        double r55017 = -r55011;
        double r55018 = r55014 * r55017;
        double r55019 = r55016 + r55018;
        double r55020 = r55011 + r55019;
        return r55020;
}

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 r55021, r55022, r55023, r55024, r55025, r55026, r55027, r55028;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r55021);
        mpfr_init(r55022);
        mpfr_init(r55023);
        mpfr_init(r55024);
        mpfr_init(r55025);
        mpfr_init(r55026);
        mpfr_init(r55027);
        mpfr_init(r55028);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r55021, x, MPFR_RNDN);
        mpfr_set_d(r55022, y, MPFR_RNDN);
        mpfr_set_d(r55023, z, MPFR_RNDN);
        mpfr_sub(r55024, r55022, r55023, MPFR_RNDN);
        mpfr_set_d(r55025, t, MPFR_RNDN);
        mpfr_sub(r55026, r55025, r55021, MPFR_RNDN);
        mpfr_mul(r55027, r55024, r55026, MPFR_RNDN);
        mpfr_add(r55028, r55021, r55027, MPFR_RNDN);
        return mpfr_get_d(r55028, MPFR_RNDN);
}

static mpfr_t r55029, r55030, r55031, r55032, r55033, r55034, r55035, r55036, r55037, r55038;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r55029);
        mpfr_init(r55030);
        mpfr_init(r55031);
        mpfr_init(r55032);
        mpfr_init(r55033);
        mpfr_init(r55034);
        mpfr_init(r55035);
        mpfr_init(r55036);
        mpfr_init(r55037);
        mpfr_init(r55038);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r55029, x, MPFR_RNDN);
        mpfr_set_d(r55030, y, MPFR_RNDN);
        mpfr_set_d(r55031, z, MPFR_RNDN);
        mpfr_sub(r55032, r55030, r55031, MPFR_RNDN);
        mpfr_set_d(r55033, t, MPFR_RNDN);
        mpfr_mul(r55034, r55032, r55033, MPFR_RNDN);
        mpfr_neg(r55035, r55029, MPFR_RNDN);
        mpfr_mul(r55036, r55032, r55035, MPFR_RNDN);
        mpfr_add(r55037, r55034, r55036, MPFR_RNDN);
        mpfr_add(r55038, r55029, r55037, MPFR_RNDN);
        return mpfr_get_d(r55038, MPFR_RNDN);
}

static mpfr_t r55039, r55040, r55041, r55042, r55043, r55044, r55045, r55046, r55047, r55048;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r55039);
        mpfr_init(r55040);
        mpfr_init(r55041);
        mpfr_init(r55042);
        mpfr_init(r55043);
        mpfr_init(r55044);
        mpfr_init(r55045);
        mpfr_init(r55046);
        mpfr_init(r55047);
        mpfr_init(r55048);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r55039, x, MPFR_RNDN);
        mpfr_set_d(r55040, y, MPFR_RNDN);
        mpfr_set_d(r55041, z, MPFR_RNDN);
        mpfr_sub(r55042, r55040, r55041, MPFR_RNDN);
        mpfr_set_d(r55043, t, MPFR_RNDN);
        mpfr_mul(r55044, r55042, r55043, MPFR_RNDN);
        mpfr_neg(r55045, r55039, MPFR_RNDN);
        mpfr_mul(r55046, r55042, r55045, MPFR_RNDN);
        mpfr_add(r55047, r55044, r55046, MPFR_RNDN);
        mpfr_add(r55048, r55039, r55047, MPFR_RNDN);
        return mpfr_get_d(r55048, MPFR_RNDN);
}