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

char *name = "Expression 1, p15";

double f_if(float a, float b, float c, float d, float e) {
        float r22294 = e;
        float r22295 = d;
        float r22296 = r22294 + r22295;
        float r22297 = c;
        float r22298 = r22296 + r22297;
        float r22299 = b;
        float r22300 = r22298 + r22299;
        float r22301 = a;
        float r22302 = r22300 + r22301;
        return r22302;
}

double f_id(double a, double b, double c, double d, double e) {
        double r22303 = e;
        double r22304 = d;
        double r22305 = r22303 + r22304;
        double r22306 = c;
        double r22307 = r22305 + r22306;
        double r22308 = b;
        double r22309 = r22307 + r22308;
        double r22310 = a;
        double r22311 = r22309 + r22310;
        return r22311;
}


double f_of(float a, float b, float c, float d, float e) {
        float r22312 = e;
        float r22313 = d;
        float r22314 = r22312 + r22313;
        float r22315 = c;
        float r22316 = r22314 + r22315;
        float r22317 = b;
        float r22318 = r22316 + r22317;
        float r22319 = a;
        float r22320 = r22318 + r22319;
        return r22320;
}

double f_od(double a, double b, double c, double d, double e) {
        double r22321 = e;
        double r22322 = d;
        double r22323 = r22321 + r22322;
        double r22324 = c;
        double r22325 = r22323 + r22324;
        double r22326 = b;
        double r22327 = r22325 + r22326;
        double r22328 = a;
        double r22329 = r22327 + r22328;
        return r22329;
}

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 r22330, r22331, r22332, r22333, r22334, r22335, r22336, r22337, r22338;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r22330);
        mpfr_init(r22331);
        mpfr_init(r22332);
        mpfr_init(r22333);
        mpfr_init(r22334);
        mpfr_init(r22335);
        mpfr_init(r22336);
        mpfr_init(r22337);
        mpfr_init(r22338);
}

double f_im(double a, double b, double c, double d, double e) {
        mpfr_set_d(r22330, e, MPFR_RNDN);
        mpfr_set_d(r22331, d, MPFR_RNDN);
        mpfr_add(r22332, r22330, r22331, MPFR_RNDN);
        mpfr_set_d(r22333, c, MPFR_RNDN);
        mpfr_add(r22334, r22332, r22333, MPFR_RNDN);
        mpfr_set_d(r22335, b, MPFR_RNDN);
        mpfr_add(r22336, r22334, r22335, MPFR_RNDN);
        mpfr_set_d(r22337, a, MPFR_RNDN);
        mpfr_add(r22338, r22336, r22337, MPFR_RNDN);
        return mpfr_get_d(r22338, MPFR_RNDN);
}

static mpfr_t r22339, r22340, r22341, r22342, r22343, r22344, r22345, r22346, r22347;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22339);
        mpfr_init(r22340);
        mpfr_init(r22341);
        mpfr_init(r22342);
        mpfr_init(r22343);
        mpfr_init(r22344);
        mpfr_init(r22345);
        mpfr_init(r22346);
        mpfr_init(r22347);
}

double f_fm(double a, double b, double c, double d, double e) {
        mpfr_set_d(r22339, e, MPFR_RNDN);
        mpfr_set_d(r22340, d, MPFR_RNDN);
        mpfr_add(r22341, r22339, r22340, MPFR_RNDN);
        mpfr_set_d(r22342, c, MPFR_RNDN);
        mpfr_add(r22343, r22341, r22342, MPFR_RNDN);
        mpfr_set_d(r22344, b, MPFR_RNDN);
        mpfr_add(r22345, r22343, r22344, MPFR_RNDN);
        mpfr_set_d(r22346, a, MPFR_RNDN);
        mpfr_add(r22347, r22345, r22346, MPFR_RNDN);
        return mpfr_get_d(r22347, MPFR_RNDN);
}

static mpfr_t r22348, r22349, r22350, r22351, r22352, r22353, r22354, r22355, r22356;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22348);
        mpfr_init(r22349);
        mpfr_init(r22350);
        mpfr_init(r22351);
        mpfr_init(r22352);
        mpfr_init(r22353);
        mpfr_init(r22354);
        mpfr_init(r22355);
        mpfr_init(r22356);
}

double f_dm(double a, double b, double c, double d, double e) {
        mpfr_set_d(r22348, e, MPFR_RNDN);
        mpfr_set_d(r22349, d, MPFR_RNDN);
        mpfr_add(r22350, r22348, r22349, MPFR_RNDN);
        mpfr_set_d(r22351, c, MPFR_RNDN);
        mpfr_add(r22352, r22350, r22351, MPFR_RNDN);
        mpfr_set_d(r22353, b, MPFR_RNDN);
        mpfr_add(r22354, r22352, r22353, MPFR_RNDN);
        mpfr_set_d(r22355, a, MPFR_RNDN);
        mpfr_add(r22356, r22354, r22355, MPFR_RNDN);
        return mpfr_get_d(r22356, MPFR_RNDN);
}

