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

char *name = "FastMath dist";

double f_if(float d1, float d2, float d3) {
        float r27372 = d1;
        float r27373 = d2;
        float r27374 = r27372 * r27373;
        float r27375 = d3;
        float r27376 = r27372 * r27375;
        float r27377 = r27374 + r27376;
        return r27377;
}

double f_id(double d1, double d2, double d3) {
        double r27378 = d1;
        double r27379 = d2;
        double r27380 = r27378 * r27379;
        double r27381 = d3;
        double r27382 = r27378 * r27381;
        double r27383 = r27380 + r27382;
        return r27383;
}


double f_of(float d1, float d2, float d3) {
        float r27384 = d1;
        float r27385 = d3;
        float r27386 = d2;
        float r27387 = r27385 + r27386;
        float r27388 = r27384 * r27387;
        return r27388;
}

double f_od(double d1, double d2, double d3) {
        double r27389 = d1;
        double r27390 = d3;
        double r27391 = d2;
        double r27392 = r27390 + r27391;
        double r27393 = r27389 * r27392;
        return r27393;
}

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 r27394, r27395, r27396, r27397, r27398, r27399;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27394);
        mpfr_init(r27395);
        mpfr_init(r27396);
        mpfr_init(r27397);
        mpfr_init(r27398);
        mpfr_init(r27399);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r27394, d1, MPFR_RNDN);
        mpfr_set_d(r27395, d2, MPFR_RNDN);
        mpfr_mul(r27396, r27394, r27395, MPFR_RNDN);
        mpfr_set_d(r27397, d3, MPFR_RNDN);
        mpfr_mul(r27398, r27394, r27397, MPFR_RNDN);
        mpfr_add(r27399, r27396, r27398, MPFR_RNDN);
        return mpfr_get_d(r27399, MPFR_RNDN);
}

static mpfr_t r27400, r27401, r27402, r27403, r27404;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27400);
        mpfr_init(r27401);
        mpfr_init(r27402);
        mpfr_init(r27403);
        mpfr_init(r27404);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r27400, d1, MPFR_RNDN);
        mpfr_set_d(r27401, d3, MPFR_RNDN);
        mpfr_set_d(r27402, d2, MPFR_RNDN);
        mpfr_add(r27403, r27401, r27402, MPFR_RNDN);
        mpfr_mul(r27404, r27400, r27403, MPFR_RNDN);
        return mpfr_get_d(r27404, MPFR_RNDN);
}

static mpfr_t r27405, r27406, r27407, r27408, r27409;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27405);
        mpfr_init(r27406);
        mpfr_init(r27407);
        mpfr_init(r27408);
        mpfr_init(r27409);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r27405, d1, MPFR_RNDN);
        mpfr_set_d(r27406, d3, MPFR_RNDN);
        mpfr_set_d(r27407, d2, MPFR_RNDN);
        mpfr_add(r27408, r27406, r27407, MPFR_RNDN);
        mpfr_mul(r27409, r27405, r27408, MPFR_RNDN);
        return mpfr_get_d(r27409, MPFR_RNDN);
}