#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 r26009 = a;
        float r26010 = b;
        float r26011 = c;
        float r26012 = r26010 + r26011;
        float r26013 = d;
        float r26014 = r26012 + r26013;
        float r26015 = r26009 * r26014;
        return r26015;
}

double f_id(double a, double b, double c, double d) {
        double r26016 = a;
        double r26017 = b;
        double r26018 = c;
        double r26019 = r26017 + r26018;
        double r26020 = d;
        double r26021 = r26019 + r26020;
        double r26022 = r26016 * r26021;
        return r26022;
}


double f_of(float a, float b, float c, float d) {
        float r26023 = a;
        float r26024 = b;
        float r26025 = c;
        float r26026 = r26024 + r26025;
        float r26027 = d;
        float r26028 = r26026 + r26027;
        float r26029 = r26023 * r26028;
        return r26029;
}

double f_od(double a, double b, double c, double d) {
        double r26030 = a;
        double r26031 = b;
        double r26032 = c;
        double r26033 = r26031 + r26032;
        double r26034 = d;
        double r26035 = r26033 + r26034;
        double r26036 = r26030 * r26035;
        return r26036;
}

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 r26037, r26038, r26039, r26040, r26041, r26042, r26043;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26037);
        mpfr_init(r26038);
        mpfr_init(r26039);
        mpfr_init(r26040);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26037, a, MPFR_RNDN);
        mpfr_set_d(r26038, b, MPFR_RNDN);
        mpfr_set_d(r26039, c, MPFR_RNDN);
        mpfr_add(r26040, r26038, r26039, MPFR_RNDN);
        mpfr_set_d(r26041, d, MPFR_RNDN);
        mpfr_add(r26042, r26040, r26041, MPFR_RNDN);
        mpfr_mul(r26043, r26037, r26042, MPFR_RNDN);
        return mpfr_get_d(r26043, MPFR_RNDN);
}

static mpfr_t r26044, r26045, r26046, r26047, r26048, r26049, r26050;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26044);
        mpfr_init(r26045);
        mpfr_init(r26046);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26044, a, MPFR_RNDN);
        mpfr_set_d(r26045, b, MPFR_RNDN);
        mpfr_set_d(r26046, c, MPFR_RNDN);
        mpfr_add(r26047, r26045, r26046, MPFR_RNDN);
        mpfr_set_d(r26048, d, MPFR_RNDN);
        mpfr_add(r26049, r26047, r26048, MPFR_RNDN);
        mpfr_mul(r26050, r26044, r26049, MPFR_RNDN);
        return mpfr_get_d(r26050, MPFR_RNDN);
}

static mpfr_t r26051, r26052, r26053, r26054, r26055, r26056, r26057;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
        mpfr_init(r26056);
        mpfr_init(r26057);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26051, a, MPFR_RNDN);
        mpfr_set_d(r26052, b, MPFR_RNDN);
        mpfr_set_d(r26053, c, MPFR_RNDN);
        mpfr_add(r26054, r26052, r26053, MPFR_RNDN);
        mpfr_set_d(r26055, d, MPFR_RNDN);
        mpfr_add(r26056, r26054, r26055, MPFR_RNDN);
        mpfr_mul(r26057, r26051, r26056, MPFR_RNDN);
        return mpfr_get_d(r26057, MPFR_RNDN);
}