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

char *name = "System.Random.MWC.Distributions:standard from mwc-random-0.13.3.2";

double f_if(float x, float y) {
        float r25357 = 0.5;
        float r25358 = x;
        float r25359 = r25358 * r25358;
        float r25360 = y;
        float r25361 = r25359 - r25360;
        float r25362 = r25357 * r25361;
        return r25362;
}

double f_id(double x, double y) {
        double r25363 = 0.5;
        double r25364 = x;
        double r25365 = r25364 * r25364;
        double r25366 = y;
        double r25367 = r25365 - r25366;
        double r25368 = r25363 * r25367;
        return r25368;
}


double f_of(float x, float y) {
        float r25369 = 0.5;
        float r25370 = x;
        float r25371 = r25370 * r25370;
        float r25372 = y;
        float r25373 = r25371 - r25372;
        float r25374 = r25369 * r25373;
        return r25374;
}

double f_od(double x, double y) {
        double r25375 = 0.5;
        double r25376 = x;
        double r25377 = r25376 * r25376;
        double r25378 = y;
        double r25379 = r25377 - r25378;
        double r25380 = r25375 * r25379;
        return r25380;
}

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 r25381, r25382, r25383, r25384, r25385, r25386;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25381, "0.5", 10, MPFR_RNDN);
        mpfr_init(r25382);
        mpfr_init(r25383);
        mpfr_init(r25384);
        mpfr_init(r25385);
        mpfr_init(r25386);
}

double f_im(double x, double y) {
        ;
        mpfr_set_d(r25382, x, MPFR_RNDN);
        mpfr_mul(r25383, r25382, r25382, MPFR_RNDN);
        mpfr_set_d(r25384, y, MPFR_RNDN);
        mpfr_sub(r25385, r25383, r25384, MPFR_RNDN);
        mpfr_mul(r25386, r25381, r25385, MPFR_RNDN);
        return mpfr_get_d(r25386, MPFR_RNDN);
}

static mpfr_t r25387, r25388, r25389, r25390, r25391, r25392;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25387, "0.5", 10, MPFR_RNDN);
        mpfr_init(r25388);
        mpfr_init(r25389);
        mpfr_init(r25390);
        mpfr_init(r25391);
        mpfr_init(r25392);
}

double f_fm(double x, double y) {
        ;
        mpfr_set_d(r25388, x, MPFR_RNDN);
        mpfr_mul(r25389, r25388, r25388, MPFR_RNDN);
        mpfr_set_d(r25390, y, MPFR_RNDN);
        mpfr_sub(r25391, r25389, r25390, MPFR_RNDN);
        mpfr_mul(r25392, r25387, r25391, MPFR_RNDN);
        return mpfr_get_d(r25392, MPFR_RNDN);
}

static mpfr_t r25393, r25394, r25395, r25396, r25397, r25398;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25393, "0.5", 10, MPFR_RNDN);
        mpfr_init(r25394);
        mpfr_init(r25395);
        mpfr_init(r25396);
        mpfr_init(r25397);
        mpfr_init(r25398);
}

double f_dm(double x, double y) {
        ;
        mpfr_set_d(r25394, x, MPFR_RNDN);
        mpfr_mul(r25395, r25394, r25394, MPFR_RNDN);
        mpfr_set_d(r25396, y, MPFR_RNDN);
        mpfr_sub(r25397, r25395, r25396, MPFR_RNDN);
        mpfr_mul(r25398, r25393, r25397, MPFR_RNDN);
        return mpfr_get_d(r25398, MPFR_RNDN);
}