#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 r51016 = x;
        float r51017 = r51016 * r51016;
        float r51018 = r51016 + r51017;
        return r51018;
}

double f_id(double x) {
        double r51019 = x;
        double r51020 = r51019 * r51019;
        double r51021 = r51019 + r51020;
        return r51021;
}


double f_of(float x) {
        float r51022 = x;
        float r51023 = r51022 * r51022;
        float r51024 = r51022 + r51023;
        return r51024;
}

double f_od(double x) {
        double r51025 = x;
        double r51026 = r51025 * r51025;
        double r51027 = r51025 + r51026;
        return r51027;
}

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 r51028, r51029, r51030;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r51028);
        mpfr_init(r51029);
        mpfr_init(r51030);
}

double f_im(double x) {
        mpfr_set_d(r51028, x, MPFR_RNDN);
        mpfr_mul(r51029, r51028, r51028, MPFR_RNDN);
        mpfr_add(r51030, r51028, r51029, MPFR_RNDN);
        return mpfr_get_d(r51030, MPFR_RNDN);
}

static mpfr_t r51031, r51032, r51033;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51031);
        mpfr_init(r51032);
        mpfr_init(r51033);
}

double f_fm(double x) {
        mpfr_set_d(r51031, x, MPFR_RNDN);
        mpfr_mul(r51032, r51031, r51031, MPFR_RNDN);
        mpfr_add(r51033, r51031, r51032, MPFR_RNDN);
        return mpfr_get_d(r51033, MPFR_RNDN);
}

static mpfr_t r51034, r51035, r51036;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r51034);
        mpfr_init(r51035);
        mpfr_init(r51036);
}

double f_dm(double x) {
        mpfr_set_d(r51034, x, MPFR_RNDN);
        mpfr_mul(r51035, r51034, r51034, MPFR_RNDN);
        mpfr_add(r51036, r51034, r51035, MPFR_RNDN);
        return mpfr_get_d(r51036, MPFR_RNDN);
}