#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 r26093 = d1;
        float r26094 = d2;
        float r26095 = r26093 * r26094;
        float r26096 = d3;
        float r26097 = 5;
        float r26098 = r26096 + r26097;
        float r26099 = r26098 * r26093;
        float r26100 = r26095 + r26099;
        float r26101 = 32;
        float r26102 = r26093 * r26101;
        float r26103 = r26100 + r26102;
        return r26103;
}

double f_id(double d1, double d2, double d3) {
        double r26104 = d1;
        double r26105 = d2;
        double r26106 = r26104 * r26105;
        double r26107 = d3;
        double r26108 = 5;
        double r26109 = r26107 + r26108;
        double r26110 = r26109 * r26104;
        double r26111 = r26106 + r26110;
        double r26112 = 32;
        double r26113 = r26104 * r26112;
        double r26114 = r26111 + r26113;
        return r26114;
}


double f_of(float d1, float d2, float d3) {
        float r26115 = d1;
        float r26116 = d3;
        float r26117 = 5;
        float r26118 = r26116 + r26117;
        float r26119 = 32;
        float r26120 = d2;
        float r26121 = r26119 + r26120;
        float r26122 = r26118 + r26121;
        float r26123 = r26115 * r26122;
        return r26123;
}

double f_od(double d1, double d2, double d3) {
        double r26124 = d1;
        double r26125 = d3;
        double r26126 = 5;
        double r26127 = r26125 + r26126;
        double r26128 = 32;
        double r26129 = d2;
        double r26130 = r26128 + r26129;
        double r26131 = r26127 + r26130;
        double r26132 = r26124 * r26131;
        return r26132;
}

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 r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140, r26141, r26142, r26143;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init(r26133);
        mpfr_init(r26134);
        mpfr_init(r26135);
        mpfr_init(r26136);
        mpfr_init_set_str(r26137, "5", 10, MPFR_RNDN);
        mpfr_init(r26138);
        mpfr_init(r26139);
        mpfr_init(r26140);
        mpfr_init_set_str(r26141, "32", 10, MPFR_RNDN);
        mpfr_init(r26142);
        mpfr_init(r26143);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r26133, d1, MPFR_RNDN);
        mpfr_set_d(r26134, d2, MPFR_RNDN);
        mpfr_mul(r26135, r26133, r26134, MPFR_RNDN);
        mpfr_set_d(r26136, d3, MPFR_RNDN);
        ;
        mpfr_add(r26138, r26136, r26137, MPFR_RNDN);
        mpfr_mul(r26139, r26138, r26133, MPFR_RNDN);
        mpfr_add(r26140, r26135, r26139, MPFR_RNDN);
        ;
        mpfr_mul(r26142, r26133, r26141, MPFR_RNDN);
        mpfr_add(r26143, r26140, r26142, MPFR_RNDN);
        return mpfr_get_d(r26143, MPFR_RNDN);
}

static mpfr_t r26144, r26145, r26146, r26147, r26148, r26149, r26150, r26151, r26152;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init(r26144);
        mpfr_init(r26145);
        mpfr_init_set_str(r26146, "5", 10, MPFR_RNDN);
        mpfr_init(r26147);
        mpfr_init_set_str(r26148, "32", 10, MPFR_RNDN);
        mpfr_init(r26149);
        mpfr_init(r26150);
        mpfr_init(r26151);
        mpfr_init(r26152);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r26144, d1, MPFR_RNDN);
        mpfr_set_d(r26145, d3, MPFR_RNDN);
        ;
        mpfr_add(r26147, r26145, r26146, MPFR_RNDN);
        ;
        mpfr_set_d(r26149, d2, MPFR_RNDN);
        mpfr_add(r26150, r26148, r26149, MPFR_RNDN);
        mpfr_add(r26151, r26147, r26150, MPFR_RNDN);
        mpfr_mul(r26152, r26144, r26151, MPFR_RNDN);
        return mpfr_get_d(r26152, MPFR_RNDN);
}

static mpfr_t r26153, r26154, r26155, r26156, r26157, r26158, r26159, r26160, r26161;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init(r26153);
        mpfr_init(r26154);
        mpfr_init_set_str(r26155, "5", 10, MPFR_RNDN);
        mpfr_init(r26156);
        mpfr_init_set_str(r26157, "32", 10, MPFR_RNDN);
        mpfr_init(r26158);
        mpfr_init(r26159);
        mpfr_init(r26160);
        mpfr_init(r26161);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r26153, d1, MPFR_RNDN);
        mpfr_set_d(r26154, d3, MPFR_RNDN);
        ;
        mpfr_add(r26156, r26154, r26155, MPFR_RNDN);
        ;
        mpfr_set_d(r26158, d2, MPFR_RNDN);
        mpfr_add(r26159, r26157, r26158, MPFR_RNDN);
        mpfr_add(r26160, r26156, r26159, MPFR_RNDN);
        mpfr_mul(r26161, r26153, r26160, MPFR_RNDN);
        return mpfr_get_d(r26161, MPFR_RNDN);
}