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

char *name = "Expression, p6";

double f_if(float a, float b, float c, float d) {
        float r51551 = a;
        float r51552 = b;
        float r51553 = c;
        float r51554 = d;
        float r51555 = r51553 + r51554;
        float r51556 = r51552 + r51555;
        float r51557 = r51551 + r51556;
        float r51558 = 2;
        float r51559 = r51557 * r51558;
        return r51559;
}

double f_id(double a, double b, double c, double d) {
        double r51560 = a;
        double r51561 = b;
        double r51562 = c;
        double r51563 = d;
        double r51564 = r51562 + r51563;
        double r51565 = r51561 + r51564;
        double r51566 = r51560 + r51565;
        double r51567 = 2;
        double r51568 = r51566 * r51567;
        return r51568;
}


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

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

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

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r51587);
        mpfr_init(r51588);
        mpfr_init(r51589);
        mpfr_init(r51590);
        mpfr_init(r51591);
        mpfr_init(r51592);
        mpfr_init(r51593);
        mpfr_init_set_str(r51594, "2", 10, MPFR_RNDN);
        mpfr_init(r51595);
}

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

static mpfr_t r51596, r51597, r51598, r51599, r51600, r51601, r51602, r51603, r51604;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51596);
        mpfr_init(r51597);
        mpfr_init(r51598);
        mpfr_init(r51599);
        mpfr_init(r51600);
        mpfr_init(r51601);
        mpfr_init(r51602);
        mpfr_init_set_str(r51603, "2", 10, MPFR_RNDN);
        mpfr_init(r51604);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r51596, a, MPFR_RNDN);
        mpfr_set_d(r51597, b, MPFR_RNDN);
        mpfr_set_d(r51598, c, MPFR_RNDN);
        mpfr_set_d(r51599, d, MPFR_RNDN);
        mpfr_add(r51600, r51598, r51599, MPFR_RNDN);
        mpfr_add(r51601, r51597, r51600, MPFR_RNDN);
        mpfr_add(r51602, r51596, r51601, MPFR_RNDN);
        ;
        mpfr_mul(r51604, r51602, r51603, MPFR_RNDN);
        return mpfr_get_d(r51604, MPFR_RNDN);
}

static mpfr_t r51605, r51606, r51607, r51608, r51609, r51610, r51611, r51612, r51613;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51605);
        mpfr_init(r51606);
        mpfr_init(r51607);
        mpfr_init(r51608);
        mpfr_init(r51609);
        mpfr_init(r51610);
        mpfr_init(r51611);
        mpfr_init_set_str(r51612, "2", 10, MPFR_RNDN);
        mpfr_init(r51613);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r51605, a, MPFR_RNDN);
        mpfr_set_d(r51606, b, MPFR_RNDN);
        mpfr_set_d(r51607, c, MPFR_RNDN);
        mpfr_set_d(r51608, d, MPFR_RNDN);
        mpfr_add(r51609, r51607, r51608, MPFR_RNDN);
        mpfr_add(r51610, r51606, r51609, MPFR_RNDN);
        mpfr_add(r51611, r51605, r51610, MPFR_RNDN);
        ;
        mpfr_mul(r51613, r51611, r51612, MPFR_RNDN);
        return mpfr_get_d(r51613, MPFR_RNDN);
}