#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 r20041 = d1;
        float r20042 = d2;
        float r20043 = r20041 * r20042;
        float r20044 = d3;
        float r20045 = r20041 * r20044;
        float r20046 = r20043 + r20045;
        return r20046;
}

double f_id(double d1, double d2, double d3) {
        double r20047 = d1;
        double r20048 = d2;
        double r20049 = r20047 * r20048;
        double r20050 = d3;
        double r20051 = r20047 * r20050;
        double r20052 = r20049 + r20051;
        return r20052;
}


double f_of(float d1, float d2, float d3) {
        float r20053 = d1;
        float r20054 = d2;
        float r20055 = r20053 * r20054;
        float r20056 = d3;
        float r20057 = r20053 * r20056;
        float r20058 = r20055 + r20057;
        return r20058;
}

double f_od(double d1, double d2, double d3) {
        double r20059 = d1;
        double r20060 = d2;
        double r20061 = r20059 * r20060;
        double r20062 = d3;
        double r20063 = r20059 * r20062;
        double r20064 = r20061 + r20063;
        return r20064;
}

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 r20065, r20066, r20067, r20068, r20069, r20070;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20065);
        mpfr_init(r20066);
        mpfr_init(r20067);
        mpfr_init(r20068);
        mpfr_init(r20069);
        mpfr_init(r20070);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r20065, d1, MPFR_RNDN);
        mpfr_set_d(r20066, d2, MPFR_RNDN);
        mpfr_mul(r20067, r20065, r20066, MPFR_RNDN);
        mpfr_set_d(r20068, d3, MPFR_RNDN);
        mpfr_mul(r20069, r20065, r20068, MPFR_RNDN);
        mpfr_add(r20070, r20067, r20069, MPFR_RNDN);
        return mpfr_get_d(r20070, MPFR_RNDN);
}

static mpfr_t r20071, r20072, r20073, r20074, r20075, r20076;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20071);
        mpfr_init(r20072);
        mpfr_init(r20073);
        mpfr_init(r20074);
        mpfr_init(r20075);
        mpfr_init(r20076);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r20071, d1, MPFR_RNDN);
        mpfr_set_d(r20072, d2, MPFR_RNDN);
        mpfr_mul(r20073, r20071, r20072, MPFR_RNDN);
        mpfr_set_d(r20074, d3, MPFR_RNDN);
        mpfr_mul(r20075, r20071, r20074, MPFR_RNDN);
        mpfr_add(r20076, r20073, r20075, MPFR_RNDN);
        return mpfr_get_d(r20076, MPFR_RNDN);
}

static mpfr_t r20077, r20078, r20079, r20080, r20081, r20082;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20077);
        mpfr_init(r20078);
        mpfr_init(r20079);
        mpfr_init(r20080);
        mpfr_init(r20081);
        mpfr_init(r20082);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r20077, d1, MPFR_RNDN);
        mpfr_set_d(r20078, d2, MPFR_RNDN);
        mpfr_mul(r20079, r20077, r20078, MPFR_RNDN);
        mpfr_set_d(r20080, d3, MPFR_RNDN);
        mpfr_mul(r20081, r20077, r20080, MPFR_RNDN);
        mpfr_add(r20082, r20079, r20081, MPFR_RNDN);
        return mpfr_get_d(r20082, MPFR_RNDN);
}