#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 r51558 = a;
        float r51559 = b;
        float r51560 = c;
        float r51561 = r51559 + r51560;
        float r51562 = d;
        float r51563 = r51561 + r51562;
        float r51564 = r51558 * r51563;
        return r51564;
}

double f_id(double a, double b, double c, double d) {
        double r51565 = a;
        double r51566 = b;
        double r51567 = c;
        double r51568 = r51566 + r51567;
        double r51569 = d;
        double r51570 = r51568 + r51569;
        double r51571 = r51565 * r51570;
        return r51571;
}


double f_of(float a, float b, float c, float d) {
        float r51572 = a;
        float r51573 = b;
        float r51574 = c;
        float r51575 = r51573 + r51574;
        float r51576 = d;
        float r51577 = r51575 + r51576;
        float r51578 = r51572 * r51577;
        return r51578;
}

double f_od(double a, double b, double c, double d) {
        double r51579 = a;
        double r51580 = b;
        double r51581 = c;
        double r51582 = r51580 + r51581;
        double r51583 = d;
        double r51584 = r51582 + r51583;
        double r51585 = r51579 * r51584;
        return r51585;
}

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 r51586, r51587, r51588, r51589, r51590, r51591, r51592;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r51586);
        mpfr_init(r51587);
        mpfr_init(r51588);
        mpfr_init(r51589);
        mpfr_init(r51590);
        mpfr_init(r51591);
        mpfr_init(r51592);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r51586, a, MPFR_RNDN);
        mpfr_set_d(r51587, b, MPFR_RNDN);
        mpfr_set_d(r51588, c, MPFR_RNDN);
        mpfr_add(r51589, r51587, r51588, MPFR_RNDN);
        mpfr_set_d(r51590, d, MPFR_RNDN);
        mpfr_add(r51591, r51589, r51590, MPFR_RNDN);
        mpfr_mul(r51592, r51586, r51591, MPFR_RNDN);
        return mpfr_get_d(r51592, MPFR_RNDN);
}

static mpfr_t r51593, r51594, r51595, r51596, r51597, r51598, r51599;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51593);
        mpfr_init(r51594);
        mpfr_init(r51595);
        mpfr_init(r51596);
        mpfr_init(r51597);
        mpfr_init(r51598);
        mpfr_init(r51599);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r51593, a, MPFR_RNDN);
        mpfr_set_d(r51594, b, MPFR_RNDN);
        mpfr_set_d(r51595, c, MPFR_RNDN);
        mpfr_add(r51596, r51594, r51595, MPFR_RNDN);
        mpfr_set_d(r51597, d, MPFR_RNDN);
        mpfr_add(r51598, r51596, r51597, MPFR_RNDN);
        mpfr_mul(r51599, r51593, r51598, MPFR_RNDN);
        return mpfr_get_d(r51599, MPFR_RNDN);
}

static mpfr_t r51600, r51601, r51602, r51603, r51604, r51605, r51606;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51600);
        mpfr_init(r51601);
        mpfr_init(r51602);
        mpfr_init(r51603);
        mpfr_init(r51604);
        mpfr_init(r51605);
        mpfr_init(r51606);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r51600, a, MPFR_RNDN);
        mpfr_set_d(r51601, b, MPFR_RNDN);
        mpfr_set_d(r51602, c, MPFR_RNDN);
        mpfr_add(r51603, r51601, r51602, MPFR_RNDN);
        mpfr_set_d(r51604, d, MPFR_RNDN);
        mpfr_add(r51605, r51603, r51604, MPFR_RNDN);
        mpfr_mul(r51606, r51600, r51605, MPFR_RNDN);
        return mpfr_get_d(r51606, MPFR_RNDN);
}