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

char *name = "Expression, p14";

double f_if(float a, float b, float c, float d) {
        float r23254 = a;
        float r23255 = b;
        float r23256 = c;
        float r23257 = r23255 + r23256;
        float r23258 = d;
        float r23259 = r23257 + r23258;
        float r23260 = r23254 * r23259;
        return r23260;
}

double f_id(double a, double b, double c, double d) {
        double r23261 = a;
        double r23262 = b;
        double r23263 = c;
        double r23264 = r23262 + r23263;
        double r23265 = d;
        double r23266 = r23264 + r23265;
        double r23267 = r23261 * r23266;
        return r23267;
}


double f_of(float a, float b, float c, float d) {
        float r23268 = a;
        float r23269 = b;
        float r23270 = c;
        float r23271 = r23269 + r23270;
        float r23272 = d;
        float r23273 = r23271 + r23272;
        float r23274 = r23268 * r23273;
        return r23274;
}

double f_od(double a, double b, double c, double d) {
        double r23275 = a;
        double r23276 = b;
        double r23277 = c;
        double r23278 = r23276 + r23277;
        double r23279 = d;
        double r23280 = r23278 + r23279;
        double r23281 = r23275 * r23280;
        return r23281;
}

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 r23282, r23283, r23284, r23285, r23286, r23287, r23288;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r23282);
        mpfr_init(r23283);
        mpfr_init(r23284);
        mpfr_init(r23285);
        mpfr_init(r23286);
        mpfr_init(r23287);
        mpfr_init(r23288);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r23282, a, MPFR_RNDN);
        mpfr_set_d(r23283, b, MPFR_RNDN);
        mpfr_set_d(r23284, c, MPFR_RNDN);
        mpfr_add(r23285, r23283, r23284, MPFR_RNDN);
        mpfr_set_d(r23286, d, MPFR_RNDN);
        mpfr_add(r23287, r23285, r23286, MPFR_RNDN);
        mpfr_mul(r23288, r23282, r23287, MPFR_RNDN);
        return mpfr_get_d(r23288, MPFR_RNDN);
}

static mpfr_t r23289, r23290, r23291, r23292, r23293, r23294, r23295;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23289);
        mpfr_init(r23290);
        mpfr_init(r23291);
        mpfr_init(r23292);
        mpfr_init(r23293);
        mpfr_init(r23294);
        mpfr_init(r23295);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r23289, a, MPFR_RNDN);
        mpfr_set_d(r23290, b, MPFR_RNDN);
        mpfr_set_d(r23291, c, MPFR_RNDN);
        mpfr_add(r23292, r23290, r23291, MPFR_RNDN);
        mpfr_set_d(r23293, d, MPFR_RNDN);
        mpfr_add(r23294, r23292, r23293, MPFR_RNDN);
        mpfr_mul(r23295, r23289, r23294, MPFR_RNDN);
        return mpfr_get_d(r23295, MPFR_RNDN);
}

static mpfr_t r23296, r23297, r23298, r23299, r23300, r23301, r23302;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r23296);
        mpfr_init(r23297);
        mpfr_init(r23298);
        mpfr_init(r23299);
        mpfr_init(r23300);
        mpfr_init(r23301);
        mpfr_init(r23302);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r23296, a, MPFR_RNDN);
        mpfr_set_d(r23297, b, MPFR_RNDN);
        mpfr_set_d(r23298, c, MPFR_RNDN);
        mpfr_add(r23299, r23297, r23298, MPFR_RNDN);
        mpfr_set_d(r23300, d, MPFR_RNDN);
        mpfr_add(r23301, r23299, r23300, MPFR_RNDN);
        mpfr_mul(r23302, r23296, r23301, MPFR_RNDN);
        return mpfr_get_d(r23302, MPFR_RNDN);
}