#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 r19535 = d1;
        float r19536 = d2;
        float r19537 = r19535 * r19536;
        float r19538 = d3;
        float r19539 = 5.0f;
        float r19540 = r19538 + r19539;
        float r19541 = r19540 * r19535;
        float r19542 = r19537 + r19541;
        float r19543 = 32.0f;
        float r19544 = r19535 * r19543;
        float r19545 = r19542 + r19544;
        return r19545;
}

double f_id(double d1, double d2, double d3) {
        double r19546 = d1;
        double r19547 = d2;
        double r19548 = r19546 * r19547;
        double r19549 = d3;
        double r19550 = 5.0;
        double r19551 = r19549 + r19550;
        double r19552 = r19551 * r19546;
        double r19553 = r19548 + r19552;
        double r19554 = 32.0;
        double r19555 = r19546 * r19554;
        double r19556 = r19553 + r19555;
        return r19556;
}


double f_of(float d1, float d2, float d3) {
        float r19557 = d1;
        float r19558 = d2;
        float r19559 = d3;
        float r19560 = r19558 + r19559;
        float r19561 = r19557 * r19560;
        float r19562 = 5.0f;
        float r19563 = 32.0f;
        float r19564 = r19562 + r19563;
        float r19565 = r19557 * r19564;
        float r19566 = r19561 + r19565;
        return r19566;
}

double f_od(double d1, double d2, double d3) {
        double r19567 = d1;
        double r19568 = d2;
        double r19569 = d3;
        double r19570 = r19568 + r19569;
        double r19571 = r19567 * r19570;
        double r19572 = 5.0;
        double r19573 = 32.0;
        double r19574 = r19572 + r19573;
        double r19575 = r19567 * r19574;
        double r19576 = r19571 + r19575;
        return r19576;
}

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 r19577, r19578, r19579, r19580, r19581, r19582, r19583, r19584, r19585, r19586, r19587;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r19577);
        mpfr_init(r19578);
        mpfr_init(r19579);
        mpfr_init(r19580);
        mpfr_init_set_str(r19581, "5", 10, MPFR_RNDN);
        mpfr_init(r19582);
        mpfr_init(r19583);
        mpfr_init(r19584);
        mpfr_init_set_str(r19585, "32", 10, MPFR_RNDN);
        mpfr_init(r19586);
        mpfr_init(r19587);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r19577, d1, MPFR_RNDN);
        mpfr_set_d(r19578, d2, MPFR_RNDN);
        mpfr_mul(r19579, r19577, r19578, MPFR_RNDN);
        mpfr_set_d(r19580, d3, MPFR_RNDN);
        ;
        mpfr_add(r19582, r19580, r19581, MPFR_RNDN);
        mpfr_mul(r19583, r19582, r19577, MPFR_RNDN);
        mpfr_add(r19584, r19579, r19583, MPFR_RNDN);
        ;
        mpfr_mul(r19586, r19577, r19585, MPFR_RNDN);
        mpfr_add(r19587, r19584, r19586, MPFR_RNDN);
        return mpfr_get_d(r19587, MPFR_RNDN);
}

static mpfr_t r19588, r19589, r19590, r19591, r19592, r19593, r19594, r19595, r19596, r19597;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19588);
        mpfr_init(r19589);
        mpfr_init(r19590);
        mpfr_init(r19591);
        mpfr_init(r19592);
        mpfr_init_set_str(r19593, "5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19594, "32", 10, MPFR_RNDN);
        mpfr_init(r19595);
        mpfr_init(r19596);
        mpfr_init(r19597);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r19588, d1, MPFR_RNDN);
        mpfr_set_d(r19589, d2, MPFR_RNDN);
        mpfr_set_d(r19590, d3, MPFR_RNDN);
        mpfr_add(r19591, r19589, r19590, MPFR_RNDN);
        mpfr_mul(r19592, r19588, r19591, MPFR_RNDN);
        ;
        ;
        mpfr_add(r19595, r19593, r19594, MPFR_RNDN);
        mpfr_mul(r19596, r19588, r19595, MPFR_RNDN);
        mpfr_add(r19597, r19592, r19596, MPFR_RNDN);
        return mpfr_get_d(r19597, MPFR_RNDN);
}

static mpfr_t r19598, r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19598);
        mpfr_init(r19599);
        mpfr_init(r19600);
        mpfr_init(r19601);
        mpfr_init(r19602);
        mpfr_init_set_str(r19603, "5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19604, "32", 10, MPFR_RNDN);
        mpfr_init(r19605);
        mpfr_init(r19606);
        mpfr_init(r19607);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r19598, d1, MPFR_RNDN);
        mpfr_set_d(r19599, d2, MPFR_RNDN);
        mpfr_set_d(r19600, d3, MPFR_RNDN);
        mpfr_add(r19601, r19599, r19600, MPFR_RNDN);
        mpfr_mul(r19602, r19598, r19601, MPFR_RNDN);
        ;
        ;
        mpfr_add(r19605, r19603, r19604, MPFR_RNDN);
        mpfr_mul(r19606, r19598, r19605, MPFR_RNDN);
        mpfr_add(r19607, r19602, r19606, MPFR_RNDN);
        return mpfr_get_d(r19607, MPFR_RNDN);
}