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

char *name = "FastMath dist3";

double f_if(float d1, float d2, float d3) {
        float r19735 = d1;
        float r19736 = d2;
        float r19737 = r19735 * r19736;
        float r19738 = d3;
        float r19739 = 5.0f;
        float r19740 = r19738 + r19739;
        float r19741 = r19740 * r19735;
        float r19742 = r19737 + r19741;
        float r19743 = 32.0f;
        float r19744 = r19735 * r19743;
        float r19745 = r19742 + r19744;
        return r19745;
}

double f_id(double d1, double d2, double d3) {
        double r19746 = d1;
        double r19747 = d2;
        double r19748 = r19746 * r19747;
        double r19749 = d3;
        double r19750 = 5.0;
        double r19751 = r19749 + r19750;
        double r19752 = r19751 * r19746;
        double r19753 = r19748 + r19752;
        double r19754 = 32.0;
        double r19755 = r19746 * r19754;
        double r19756 = r19753 + r19755;
        return r19756;
}


double f_of(float d1, float d2, float d3) {
        float r19757 = d3;
        float r19758 = 37.0f;
        float r19759 = r19757 + r19758;
        float r19760 = d1;
        float r19761 = d2;
        float r19762 = r19760 * r19761;
        float r19763 = fma(r19759, r19760, r19762);
        return r19763;
}

double f_od(double d1, double d2, double d3) {
        double r19764 = d3;
        double r19765 = 37.0;
        double r19766 = r19764 + r19765;
        double r19767 = d1;
        double r19768 = d2;
        double r19769 = r19767 * r19768;
        double r19770 = fma(r19766, r19767, r19769);
        return r19770;
}

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 r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r19771);
        mpfr_init(r19772);
        mpfr_init(r19773);
        mpfr_init(r19774);
        mpfr_init_set_str(r19775, "5", 10, MPFR_RNDN);
        mpfr_init(r19776);
        mpfr_init(r19777);
        mpfr_init(r19778);
        mpfr_init_set_str(r19779, "32", 10, MPFR_RNDN);
        mpfr_init(r19780);
        mpfr_init(r19781);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r19771, d1, MPFR_RNDN);
        mpfr_set_d(r19772, d2, MPFR_RNDN);
        mpfr_mul(r19773, r19771, r19772, MPFR_RNDN);
        mpfr_set_d(r19774, d3, MPFR_RNDN);
        ;
        mpfr_add(r19776, r19774, r19775, MPFR_RNDN);
        mpfr_mul(r19777, r19776, r19771, MPFR_RNDN);
        mpfr_add(r19778, r19773, r19777, MPFR_RNDN);
        ;
        mpfr_mul(r19780, r19771, r19779, MPFR_RNDN);
        mpfr_add(r19781, r19778, r19780, MPFR_RNDN);
        return mpfr_get_d(r19781, MPFR_RNDN);
}

static mpfr_t r19782, r19783, r19784, r19785, r19786, r19787, r19788;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19782);
        mpfr_init_set_str(r19783, "37", 10, MPFR_RNDN);
        mpfr_init(r19784);
        mpfr_init(r19785);
        mpfr_init(r19786);
        mpfr_init(r19787);
        mpfr_init(r19788);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r19782, d3, MPFR_RNDN);
        ;
        mpfr_add(r19784, r19782, r19783, MPFR_RNDN);
        mpfr_set_d(r19785, d1, MPFR_RNDN);
        mpfr_set_d(r19786, d2, MPFR_RNDN);
        mpfr_mul(r19787, r19785, r19786, MPFR_RNDN);
        mpfr_fma(r19788, r19784, r19785, r19787, MPFR_RNDN);
        return mpfr_get_d(r19788, MPFR_RNDN);
}

static mpfr_t r19789, r19790, r19791, r19792, r19793, r19794, r19795;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19789);
        mpfr_init_set_str(r19790, "37", 10, MPFR_RNDN);
        mpfr_init(r19791);
        mpfr_init(r19792);
        mpfr_init(r19793);
        mpfr_init(r19794);
        mpfr_init(r19795);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r19789, d3, MPFR_RNDN);
        ;
        mpfr_add(r19791, r19789, r19790, MPFR_RNDN);
        mpfr_set_d(r19792, d1, MPFR_RNDN);
        mpfr_set_d(r19793, d2, MPFR_RNDN);
        mpfr_mul(r19794, r19792, r19793, MPFR_RNDN);
        mpfr_fma(r19795, r19791, r19792, r19794, MPFR_RNDN);
        return mpfr_get_d(r19795, MPFR_RNDN);
}