#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 r20216 = a;
        float r20217 = b;
        float r20218 = c;
        float r20219 = r20217 + r20218;
        float r20220 = d;
        float r20221 = r20219 + r20220;
        float r20222 = r20216 * r20221;
        return r20222;
}

double f_id(double a, double b, double c, double d) {
        double r20223 = a;
        double r20224 = b;
        double r20225 = c;
        double r20226 = r20224 + r20225;
        double r20227 = d;
        double r20228 = r20226 + r20227;
        double r20229 = r20223 * r20228;
        return r20229;
}


double f_of(float a, float b, float c, float d) {
        float r20230 = a;
        float r20231 = b;
        float r20232 = c;
        float r20233 = r20231 + r20232;
        float r20234 = d;
        float r20235 = r20233 + r20234;
        float r20236 = r20230 * r20235;
        return r20236;
}

double f_od(double a, double b, double c, double d) {
        double r20237 = a;
        double r20238 = b;
        double r20239 = c;
        double r20240 = r20238 + r20239;
        double r20241 = d;
        double r20242 = r20240 + r20241;
        double r20243 = r20237 * r20242;
        return r20243;
}

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 r20244, r20245, r20246, r20247, r20248, r20249, r20250;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20244);
        mpfr_init(r20245);
        mpfr_init(r20246);
        mpfr_init(r20247);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r20244, a, MPFR_RNDN);
        mpfr_set_d(r20245, b, MPFR_RNDN);
        mpfr_set_d(r20246, c, MPFR_RNDN);
        mpfr_add(r20247, r20245, r20246, MPFR_RNDN);
        mpfr_set_d(r20248, d, MPFR_RNDN);
        mpfr_add(r20249, r20247, r20248, MPFR_RNDN);
        mpfr_mul(r20250, r20244, r20249, MPFR_RNDN);
        return mpfr_get_d(r20250, MPFR_RNDN);
}

static mpfr_t r20251, r20252, r20253, r20254, r20255, r20256, r20257;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r20251, a, MPFR_RNDN);
        mpfr_set_d(r20252, b, MPFR_RNDN);
        mpfr_set_d(r20253, c, MPFR_RNDN);
        mpfr_add(r20254, r20252, r20253, MPFR_RNDN);
        mpfr_set_d(r20255, d, MPFR_RNDN);
        mpfr_add(r20256, r20254, r20255, MPFR_RNDN);
        mpfr_mul(r20257, r20251, r20256, MPFR_RNDN);
        return mpfr_get_d(r20257, MPFR_RNDN);
}

static mpfr_t r20258, r20259, r20260, r20261, r20262, r20263, r20264;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20258);
        mpfr_init(r20259);
        mpfr_init(r20260);
        mpfr_init(r20261);
        mpfr_init(r20262);
        mpfr_init(r20263);
        mpfr_init(r20264);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r20258, a, MPFR_RNDN);
        mpfr_set_d(r20259, b, MPFR_RNDN);
        mpfr_set_d(r20260, c, MPFR_RNDN);
        mpfr_add(r20261, r20259, r20260, MPFR_RNDN);
        mpfr_set_d(r20262, d, MPFR_RNDN);
        mpfr_add(r20263, r20261, r20262, MPFR_RNDN);
        mpfr_mul(r20264, r20258, r20263, MPFR_RNDN);
        return mpfr_get_d(r20264, MPFR_RNDN);
}