#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 r18594 = d1;
        float r18595 = 3.0f;
        float r18596 = r18594 * r18595;
        float r18597 = d2;
        float r18598 = r18594 * r18597;
        float r18599 = r18596 + r18598;
        float r18600 = d3;
        float r18601 = r18594 * r18600;
        float r18602 = r18599 + r18601;
        return r18602;
}

double f_id(double d1, double d2, double d3) {
        double r18603 = d1;
        double r18604 = 3.0;
        double r18605 = r18603 * r18604;
        double r18606 = d2;
        double r18607 = r18603 * r18606;
        double r18608 = r18605 + r18607;
        double r18609 = d3;
        double r18610 = r18603 * r18609;
        double r18611 = r18608 + r18610;
        return r18611;
}


double f_of(float d1, float d2, float d3) {
        float r18612 = d1;
        float r18613 = 3.0f;
        float r18614 = d2;
        float r18615 = r18613 + r18614;
        float r18616 = d3;
        float r18617 = r18615 + r18616;
        float r18618 = r18612 * r18617;
        return r18618;
}

double f_od(double d1, double d2, double d3) {
        double r18619 = d1;
        double r18620 = 3.0;
        double r18621 = d2;
        double r18622 = r18620 + r18621;
        double r18623 = d3;
        double r18624 = r18622 + r18623;
        double r18625 = r18619 * r18624;
        return r18625;
}

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 r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18626);
        mpfr_init_set_str(r18627, "3", 10, MPFR_RNDN);
        mpfr_init(r18628);
        mpfr_init(r18629);
        mpfr_init(r18630);
        mpfr_init(r18631);
        mpfr_init(r18632);
        mpfr_init(r18633);
        mpfr_init(r18634);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r18626, d1, MPFR_RNDN);
        ;
        mpfr_mul(r18628, r18626, r18627, MPFR_RNDN);
        mpfr_set_d(r18629, d2, MPFR_RNDN);
        mpfr_mul(r18630, r18626, r18629, MPFR_RNDN);
        mpfr_add(r18631, r18628, r18630, MPFR_RNDN);
        mpfr_set_d(r18632, d3, MPFR_RNDN);
        mpfr_mul(r18633, r18626, r18632, MPFR_RNDN);
        mpfr_add(r18634, r18631, r18633, MPFR_RNDN);
        return mpfr_get_d(r18634, MPFR_RNDN);
}

static mpfr_t r18635, r18636, r18637, r18638, r18639, r18640, r18641;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18635);
        mpfr_init_set_str(r18636, "3", 10, MPFR_RNDN);
        mpfr_init(r18637);
        mpfr_init(r18638);
        mpfr_init(r18639);
        mpfr_init(r18640);
        mpfr_init(r18641);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r18635, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r18637, d2, MPFR_RNDN);
        mpfr_add(r18638, r18636, r18637, MPFR_RNDN);
        mpfr_set_d(r18639, d3, MPFR_RNDN);
        mpfr_add(r18640, r18638, r18639, MPFR_RNDN);
        mpfr_mul(r18641, r18635, r18640, MPFR_RNDN);
        return mpfr_get_d(r18641, MPFR_RNDN);
}

static mpfr_t r18642, r18643, r18644, r18645, r18646, r18647, r18648;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18642);
        mpfr_init_set_str(r18643, "3", 10, MPFR_RNDN);
        mpfr_init(r18644);
        mpfr_init(r18645);
        mpfr_init(r18646);
        mpfr_init(r18647);
        mpfr_init(r18648);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r18642, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r18644, d2, MPFR_RNDN);
        mpfr_add(r18645, r18643, r18644, MPFR_RNDN);
        mpfr_set_d(r18646, d3, MPFR_RNDN);
        mpfr_add(r18647, r18645, r18646, MPFR_RNDN);
        mpfr_mul(r18648, r18642, r18647, MPFR_RNDN);
        return mpfr_get_d(r18648, MPFR_RNDN);
}