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

char *name = "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B";

double f_if(float x, float y, float z, float t) {
        float r55690 = x;
        float r55691 = y;
        float r55692 = z;
        float r55693 = r55691 - r55692;
        float r55694 = t;
        float r55695 = r55694 - r55692;
        float r55696 = r55693 * r55695;
        float r55697 = r55690 / r55696;
        return r55697;
}

double f_id(double x, double y, double z, double t) {
        double r55698 = x;
        double r55699 = y;
        double r55700 = z;
        double r55701 = r55699 - r55700;
        double r55702 = t;
        double r55703 = r55702 - r55700;
        double r55704 = r55701 * r55703;
        double r55705 = r55698 / r55704;
        return r55705;
}


double f_of(float x, float y, float z, float t) {
        float r55706 = x;
        float r55707 = y;
        float r55708 = z;
        float r55709 = r55707 - r55708;
        float r55710 = r55706 / r55709;
        float r55711 = t;
        float r55712 = r55711 - r55708;
        float r55713 = r55710 / r55712;
        return r55713;
}

double f_od(double x, double y, double z, double t) {
        double r55714 = x;
        double r55715 = y;
        double r55716 = z;
        double r55717 = r55715 - r55716;
        double r55718 = r55714 / r55717;
        double r55719 = t;
        double r55720 = r55719 - r55716;
        double r55721 = r55718 / r55720;
        return r55721;
}

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 r55722, r55723, r55724, r55725, r55726, r55727, r55728, r55729;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r55722);
        mpfr_init(r55723);
        mpfr_init(r55724);
        mpfr_init(r55725);
        mpfr_init(r55726);
        mpfr_init(r55727);
        mpfr_init(r55728);
        mpfr_init(r55729);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r55722, x, MPFR_RNDN);
        mpfr_set_d(r55723, y, MPFR_RNDN);
        mpfr_set_d(r55724, z, MPFR_RNDN);
        mpfr_sub(r55725, r55723, r55724, MPFR_RNDN);
        mpfr_set_d(r55726, t, MPFR_RNDN);
        mpfr_sub(r55727, r55726, r55724, MPFR_RNDN);
        mpfr_mul(r55728, r55725, r55727, MPFR_RNDN);
        mpfr_div(r55729, r55722, r55728, MPFR_RNDN);
        return mpfr_get_d(r55729, MPFR_RNDN);
}

static mpfr_t r55730, r55731, r55732, r55733, r55734, r55735, r55736, r55737;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55730);
        mpfr_init(r55731);
        mpfr_init(r55732);
        mpfr_init(r55733);
        mpfr_init(r55734);
        mpfr_init(r55735);
        mpfr_init(r55736);
        mpfr_init(r55737);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r55730, x, MPFR_RNDN);
        mpfr_set_d(r55731, y, MPFR_RNDN);
        mpfr_set_d(r55732, z, MPFR_RNDN);
        mpfr_sub(r55733, r55731, r55732, MPFR_RNDN);
        mpfr_div(r55734, r55730, r55733, MPFR_RNDN);
        mpfr_set_d(r55735, t, MPFR_RNDN);
        mpfr_sub(r55736, r55735, r55732, MPFR_RNDN);
        mpfr_div(r55737, r55734, r55736, MPFR_RNDN);
        return mpfr_get_d(r55737, MPFR_RNDN);
}

static mpfr_t r55738, r55739, r55740, r55741, r55742, r55743, r55744, r55745;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55738);
        mpfr_init(r55739);
        mpfr_init(r55740);
        mpfr_init(r55741);
        mpfr_init(r55742);
        mpfr_init(r55743);
        mpfr_init(r55744);
        mpfr_init(r55745);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r55738, x, MPFR_RNDN);
        mpfr_set_d(r55739, y, MPFR_RNDN);
        mpfr_set_d(r55740, z, MPFR_RNDN);
        mpfr_sub(r55741, r55739, r55740, MPFR_RNDN);
        mpfr_div(r55742, r55738, r55741, MPFR_RNDN);
        mpfr_set_d(r55743, t, MPFR_RNDN);
        mpfr_sub(r55744, r55743, r55740, MPFR_RNDN);
        mpfr_div(r55745, r55742, r55744, MPFR_RNDN);
        return mpfr_get_d(r55745, MPFR_RNDN);
}