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

char *name = "FastMath test3";

double f_if(float d1, float d2, float d3) {
        float r18742 = d1;
        float r18743 = 3.0f;
        float r18744 = r18742 * r18743;
        float r18745 = d2;
        float r18746 = r18742 * r18745;
        float r18747 = r18744 + r18746;
        float r18748 = d3;
        float r18749 = r18742 * r18748;
        float r18750 = r18747 + r18749;
        return r18750;
}

double f_id(double d1, double d2, double d3) {
        double r18751 = d1;
        double r18752 = 3.0;
        double r18753 = r18751 * r18752;
        double r18754 = d2;
        double r18755 = r18751 * r18754;
        double r18756 = r18753 + r18755;
        double r18757 = d3;
        double r18758 = r18751 * r18757;
        double r18759 = r18756 + r18758;
        return r18759;
}


double f_of(float d1, float d2, float d3) {
        float r18760 = d1;
        float r18761 = 3.0f;
        float r18762 = d2;
        float r18763 = r18761 + r18762;
        float r18764 = d3;
        float r18765 = r18763 + r18764;
        float r18766 = r18760 * r18765;
        return r18766;
}

double f_od(double d1, double d2, double d3) {
        double r18767 = d1;
        double r18768 = 3.0;
        double r18769 = d2;
        double r18770 = r18768 + r18769;
        double r18771 = d3;
        double r18772 = r18770 + r18771;
        double r18773 = r18767 * r18772;
        return r18773;
}

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 r18774, r18775, r18776, r18777, r18778, r18779, r18780, r18781, r18782;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18774);
        mpfr_init_set_str(r18775, "3", 10, MPFR_RNDN);
        mpfr_init(r18776);
        mpfr_init(r18777);
        mpfr_init(r18778);
        mpfr_init(r18779);
        mpfr_init(r18780);
        mpfr_init(r18781);
        mpfr_init(r18782);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r18774, d1, MPFR_RNDN);
        ;
        mpfr_mul(r18776, r18774, r18775, MPFR_RNDN);
        mpfr_set_d(r18777, d2, MPFR_RNDN);
        mpfr_mul(r18778, r18774, r18777, MPFR_RNDN);
        mpfr_add(r18779, r18776, r18778, MPFR_RNDN);
        mpfr_set_d(r18780, d3, MPFR_RNDN);
        mpfr_mul(r18781, r18774, r18780, MPFR_RNDN);
        mpfr_add(r18782, r18779, r18781, MPFR_RNDN);
        return mpfr_get_d(r18782, MPFR_RNDN);
}

static mpfr_t r18783, r18784, r18785, r18786, r18787, r18788, r18789;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18783);
        mpfr_init_set_str(r18784, "3", 10, MPFR_RNDN);
        mpfr_init(r18785);
        mpfr_init(r18786);
        mpfr_init(r18787);
        mpfr_init(r18788);
        mpfr_init(r18789);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r18783, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r18785, d2, MPFR_RNDN);
        mpfr_add(r18786, r18784, r18785, MPFR_RNDN);
        mpfr_set_d(r18787, d3, MPFR_RNDN);
        mpfr_add(r18788, r18786, r18787, MPFR_RNDN);
        mpfr_mul(r18789, r18783, r18788, MPFR_RNDN);
        return mpfr_get_d(r18789, MPFR_RNDN);
}

static mpfr_t r18790, r18791, r18792, r18793, r18794, r18795, r18796;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18790);
        mpfr_init_set_str(r18791, "3", 10, MPFR_RNDN);
        mpfr_init(r18792);
        mpfr_init(r18793);
        mpfr_init(r18794);
        mpfr_init(r18795);
        mpfr_init(r18796);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r18790, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r18792, d2, MPFR_RNDN);
        mpfr_add(r18793, r18791, r18792, MPFR_RNDN);
        mpfr_set_d(r18794, d3, MPFR_RNDN);
        mpfr_add(r18795, r18793, r18794, MPFR_RNDN);
        mpfr_mul(r18796, r18790, r18795, MPFR_RNDN);
        return mpfr_get_d(r18796, MPFR_RNDN);
}