#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 r20788 = d1;
        float r20789 = d2;
        float r20790 = r20788 * r20789;
        float r20791 = d3;
        float r20792 = r20788 * r20791;
        float r20793 = r20790 + r20792;
        return r20793;
}

double f_id(double d1, double d2, double d3) {
        double r20794 = d1;
        double r20795 = d2;
        double r20796 = r20794 * r20795;
        double r20797 = d3;
        double r20798 = r20794 * r20797;
        double r20799 = r20796 + r20798;
        return r20799;
}


double f_of(float d1, float d2, float d3) {
        float r20800 = d1;
        float r20801 = d2;
        float r20802 = d3;
        float r20803 = r20801 + r20802;
        float r20804 = r20800 * r20803;
        return r20804;
}

double f_od(double d1, double d2, double d3) {
        double r20805 = d1;
        double r20806 = d2;
        double r20807 = d3;
        double r20808 = r20806 + r20807;
        double r20809 = r20805 * r20808;
        return r20809;
}

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 r20810, r20811, r20812, r20813, r20814, r20815;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20810);
        mpfr_init(r20811);
        mpfr_init(r20812);
        mpfr_init(r20813);
        mpfr_init(r20814);
        mpfr_init(r20815);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r20810, d1, MPFR_RNDN);
        mpfr_set_d(r20811, d2, MPFR_RNDN);
        mpfr_mul(r20812, r20810, r20811, MPFR_RNDN);
        mpfr_set_d(r20813, d3, MPFR_RNDN);
        mpfr_mul(r20814, r20810, r20813, MPFR_RNDN);
        mpfr_add(r20815, r20812, r20814, MPFR_RNDN);
        return mpfr_get_d(r20815, MPFR_RNDN);
}

static mpfr_t r20816, r20817, r20818, r20819, r20820;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20816);
        mpfr_init(r20817);
        mpfr_init(r20818);
        mpfr_init(r20819);
        mpfr_init(r20820);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r20816, d1, MPFR_RNDN);
        mpfr_set_d(r20817, d2, MPFR_RNDN);
        mpfr_set_d(r20818, d3, MPFR_RNDN);
        mpfr_add(r20819, r20817, r20818, MPFR_RNDN);
        mpfr_mul(r20820, r20816, r20819, MPFR_RNDN);
        return mpfr_get_d(r20820, MPFR_RNDN);
}

static mpfr_t r20821, r20822, r20823, r20824, r20825;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20821);
        mpfr_init(r20822);
        mpfr_init(r20823);
        mpfr_init(r20824);
        mpfr_init(r20825);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r20821, d1, MPFR_RNDN);
        mpfr_set_d(r20822, d2, MPFR_RNDN);
        mpfr_set_d(r20823, d3, MPFR_RNDN);
        mpfr_add(r20824, r20822, r20823, MPFR_RNDN);
        mpfr_mul(r20825, r20821, r20824, MPFR_RNDN);
        return mpfr_get_d(r20825, MPFR_RNDN);
}