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

char *name = "Expression 2, p15";

double f_if(float x) {
        float r51719 = x;
        float r51720 = r51719 * r51719;
        float r51721 = r51719 + r51720;
        return r51721;
}

double f_id(double x) {
        double r51722 = x;
        double r51723 = r51722 * r51722;
        double r51724 = r51722 + r51723;
        return r51724;
}


double f_of(float x) {
        float r51725 = x;
        float r51726 = r51725 * r51725;
        float r51727 = r51725 + r51726;
        return r51727;
}

double f_od(double x) {
        double r51728 = x;
        double r51729 = r51728 * r51728;
        double r51730 = r51728 + r51729;
        return r51730;
}

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 r51731, r51732, r51733;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r51731);
        mpfr_init(r51732);
        mpfr_init(r51733);
}

double f_im(double x) {
        mpfr_set_d(r51731, x, MPFR_RNDN);
        mpfr_mul(r51732, r51731, r51731, MPFR_RNDN);
        mpfr_add(r51733, r51731, r51732, MPFR_RNDN);
        return mpfr_get_d(r51733, MPFR_RNDN);
}

static mpfr_t r51734, r51735, r51736;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51734);
        mpfr_init(r51735);
        mpfr_init(r51736);
}

double f_fm(double x) {
        mpfr_set_d(r51734, x, MPFR_RNDN);
        mpfr_mul(r51735, r51734, r51734, MPFR_RNDN);
        mpfr_add(r51736, r51734, r51735, MPFR_RNDN);
        return mpfr_get_d(r51736, MPFR_RNDN);
}

static mpfr_t r51737, r51738, r51739;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51737);
        mpfr_init(r51738);
        mpfr_init(r51739);
}

double f_dm(double x) {
        mpfr_set_d(r51737, x, MPFR_RNDN);
        mpfr_mul(r51738, r51737, r51737, MPFR_RNDN);
        mpfr_add(r51739, r51737, r51738, MPFR_RNDN);
        return mpfr_get_d(r51739, MPFR_RNDN);
}