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

char *name = "Expression 4, p15";

double f_if(float a, float b) {
        float r10361 = a;
        float r10362 = b;
        float r10363 = r10361 + r10362;
        float r10364 = r10363 * r10363;
        return r10364;
}

double f_id(double a, double b) {
        double r10365 = a;
        double r10366 = b;
        double r10367 = r10365 + r10366;
        double r10368 = r10367 * r10367;
        return r10368;
}


double f_of(float a, float b) {
        float r10369 = a;
        float r10370 = b;
        float r10371 = r10369 + r10370;
        float r10372 = r10371 * r10371;
        return r10372;
}

double f_od(double a, double b) {
        double r10373 = a;
        double r10374 = b;
        double r10375 = r10373 + r10374;
        double r10376 = r10375 * r10375;
        return r10376;
}

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 r10377, r10378, r10379, r10380;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r10377);
        mpfr_init(r10378);
        mpfr_init(r10379);
        mpfr_init(r10380);
}

double f_im(double a, double b) {
        mpfr_set_d(r10377, a, MPFR_RNDN);
        mpfr_set_d(r10378, b, MPFR_RNDN);
        mpfr_add(r10379, r10377, r10378, MPFR_RNDN);
        mpfr_mul(r10380, r10379, r10379, MPFR_RNDN);
        return mpfr_get_d(r10380, MPFR_RNDN);
}

static mpfr_t r10381, r10382, r10383, r10384;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r10381);
        mpfr_init(r10382);
        mpfr_init(r10383);
        mpfr_init(r10384);
}

double f_fm(double a, double b) {
        mpfr_set_d(r10381, a, MPFR_RNDN);
        mpfr_set_d(r10382, b, MPFR_RNDN);
        mpfr_add(r10383, r10381, r10382, MPFR_RNDN);
        mpfr_mul(r10384, r10383, r10383, MPFR_RNDN);
        return mpfr_get_d(r10384, MPFR_RNDN);
}

static mpfr_t r10385, r10386, r10387, r10388;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r10385);
        mpfr_init(r10386);
        mpfr_init(r10387);
        mpfr_init(r10388);
}

double f_dm(double a, double b) {
        mpfr_set_d(r10385, a, MPFR_RNDN);
        mpfr_set_d(r10386, b, MPFR_RNDN);
        mpfr_add(r10387, r10385, r10386, MPFR_RNDN);
        mpfr_mul(r10388, r10387, r10387, MPFR_RNDN);
        return mpfr_get_d(r10388, MPFR_RNDN);
}