#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 r50911 = a;
        float r50912 = b;
        float r50913 = c;
        float r50914 = r50912 + r50913;
        float r50915 = d;
        float r50916 = r50914 + r50915;
        float r50917 = r50911 * r50916;
        return r50917;
}

double f_id(double a, double b, double c, double d) {
        double r50918 = a;
        double r50919 = b;
        double r50920 = c;
        double r50921 = r50919 + r50920;
        double r50922 = d;
        double r50923 = r50921 + r50922;
        double r50924 = r50918 * r50923;
        return r50924;
}


double f_of(float a, float b, float c, float d) {
        float r50925 = a;
        float r50926 = b;
        float r50927 = c;
        float r50928 = r50926 + r50927;
        float r50929 = d;
        float r50930 = r50928 + r50929;
        float r50931 = r50925 * r50930;
        return r50931;
}

double f_od(double a, double b, double c, double d) {
        double r50932 = a;
        double r50933 = b;
        double r50934 = c;
        double r50935 = r50933 + r50934;
        double r50936 = d;
        double r50937 = r50935 + r50936;
        double r50938 = r50932 * r50937;
        return r50938;
}

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 r50939, r50940, r50941, r50942, r50943, r50944, r50945;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r50939);
        mpfr_init(r50940);
        mpfr_init(r50941);
        mpfr_init(r50942);
        mpfr_init(r50943);
        mpfr_init(r50944);
        mpfr_init(r50945);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r50939, a, MPFR_RNDN);
        mpfr_set_d(r50940, b, MPFR_RNDN);
        mpfr_set_d(r50941, c, MPFR_RNDN);
        mpfr_add(r50942, r50940, r50941, MPFR_RNDN);
        mpfr_set_d(r50943, d, MPFR_RNDN);
        mpfr_add(r50944, r50942, r50943, MPFR_RNDN);
        mpfr_mul(r50945, r50939, r50944, MPFR_RNDN);
        return mpfr_get_d(r50945, MPFR_RNDN);
}

static mpfr_t r50946, r50947, r50948, r50949, r50950, r50951, r50952;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50946);
        mpfr_init(r50947);
        mpfr_init(r50948);
        mpfr_init(r50949);
        mpfr_init(r50950);
        mpfr_init(r50951);
        mpfr_init(r50952);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r50946, a, MPFR_RNDN);
        mpfr_set_d(r50947, b, MPFR_RNDN);
        mpfr_set_d(r50948, c, MPFR_RNDN);
        mpfr_add(r50949, r50947, r50948, MPFR_RNDN);
        mpfr_set_d(r50950, d, MPFR_RNDN);
        mpfr_add(r50951, r50949, r50950, MPFR_RNDN);
        mpfr_mul(r50952, r50946, r50951, MPFR_RNDN);
        return mpfr_get_d(r50952, MPFR_RNDN);
}

static mpfr_t r50953, r50954, r50955, r50956, r50957, r50958, r50959;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50953);
        mpfr_init(r50954);
        mpfr_init(r50955);
        mpfr_init(r50956);
        mpfr_init(r50957);
        mpfr_init(r50958);
        mpfr_init(r50959);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r50953, a, MPFR_RNDN);
        mpfr_set_d(r50954, b, MPFR_RNDN);
        mpfr_set_d(r50955, c, MPFR_RNDN);
        mpfr_add(r50956, r50954, r50955, MPFR_RNDN);
        mpfr_set_d(r50957, d, MPFR_RNDN);
        mpfr_add(r50958, r50956, r50957, MPFR_RNDN);
        mpfr_mul(r50959, r50953, r50958, MPFR_RNDN);
        return mpfr_get_d(r50959, MPFR_RNDN);
}