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

char *name = "Expanding a square";

double f_if(float x) {
        float r7786 = x;
        float r7787 = 1;
        float r7788 = r7786 + r7787;
        float r7789 = r7788 * r7788;
        float r7790 = r7789 - r7787;
        return r7790;
}

double f_id(double x) {
        double r7791 = x;
        double r7792 = 1;
        double r7793 = r7791 + r7792;
        double r7794 = r7793 * r7793;
        double r7795 = r7794 - r7792;
        return r7795;
}


double f_of(float x) {
        float r7796 = x;
        float r7797 = 1;
        float r7798 = r7796 + r7797;
        float r7799 = r7798 * r7798;
        float r7800 = r7799 - r7797;
        float r7801 = 0.06199215151062582;
        bool r7802 = r7800 <= r7801;
        float r7803 = r7796 + r7796;
        float r7804 = fma(r7796, r7796, r7803);
        float r7805 = expm1(r7804);
        float r7806 = log1p(r7805);
        float r7807 = r7802 ? r7806 : r7800;
        return r7807;
}

double f_od(double x) {
        double r7808 = x;
        double r7809 = 1;
        double r7810 = r7808 + r7809;
        double r7811 = r7810 * r7810;
        double r7812 = r7811 - r7809;
        double r7813 = 0.06199215151062582;
        bool r7814 = r7812 <= r7813;
        double r7815 = r7808 + r7808;
        double r7816 = fma(r7808, r7808, r7815);
        double r7817 = expm1(r7816);
        double r7818 = log1p(r7817);
        double r7819 = r7814 ? r7818 : r7812;
        return r7819;
}

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 r7820, r7821, r7822, r7823, r7824;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7820);
        mpfr_init_set_str(r7821, "1", 10, MPFR_RNDN);
        mpfr_init(r7822);
        mpfr_init(r7823);
        mpfr_init(r7824);
}

double f_im(double x) {
        mpfr_set_d(r7820, x, MPFR_RNDN);
        ;
        mpfr_add(r7822, r7820, r7821, MPFR_RNDN);
        mpfr_mul(r7823, r7822, r7822, MPFR_RNDN);
        mpfr_sub(r7824, r7823, r7821, MPFR_RNDN);
        return mpfr_get_d(r7824, MPFR_RNDN);
}

static mpfr_t r7825, r7826, r7827, r7828, r7829, r7830, r7831, r7832, r7833, r7834, r7835, r7836;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7825);
        mpfr_init_set_str(r7826, "1", 10, MPFR_RNDN);
        mpfr_init(r7827);
        mpfr_init(r7828);
        mpfr_init(r7829);
        mpfr_init_set_str(r7830, "0.06199215151062582", 10, MPFR_RNDN);
        mpfr_init(r7831);
        mpfr_init(r7832);
        mpfr_init(r7833);
        mpfr_init(r7834);
        mpfr_init(r7835);
        mpfr_init(r7836);
}

double f_fm(double x) {
        mpfr_set_d(r7825, x, MPFR_RNDN);
        ;
        mpfr_add(r7827, r7825, r7826, MPFR_RNDN);
        mpfr_mul(r7828, r7827, r7827, MPFR_RNDN);
        mpfr_sub(r7829, r7828, r7826, MPFR_RNDN);
        ;
        mpfr_set_si(r7831, mpfr_cmp(r7829, r7830) <= 0, MPFR_RNDN);
        mpfr_add(r7832, r7825, r7825, MPFR_RNDN);
        mpfr_fma(r7833, r7825, r7825, r7832, MPFR_RNDN);
        mpfr_expm1(r7834, r7833, MPFR_RNDN);
        mpfr_log1p(r7835, r7834, MPFR_RNDN);
        if (mpfr_get_si(r7831, MPFR_RNDN)) { mpfr_set(r7836, r7835, MPFR_RNDN); } else { mpfr_set(r7836, r7829, MPFR_RNDN); };
        return mpfr_get_d(r7836, MPFR_RNDN);
}

static mpfr_t r7837, r7838, r7839, r7840, r7841, r7842, r7843, r7844, r7845, r7846, r7847, r7848;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7837);
        mpfr_init_set_str(r7838, "1", 10, MPFR_RNDN);
        mpfr_init(r7839);
        mpfr_init(r7840);
        mpfr_init(r7841);
        mpfr_init_set_str(r7842, "0.06199215151062582", 10, MPFR_RNDN);
        mpfr_init(r7843);
        mpfr_init(r7844);
        mpfr_init(r7845);
        mpfr_init(r7846);
        mpfr_init(r7847);
        mpfr_init(r7848);
}

double f_dm(double x) {
        mpfr_set_d(r7837, x, MPFR_RNDN);
        ;
        mpfr_add(r7839, r7837, r7838, MPFR_RNDN);
        mpfr_mul(r7840, r7839, r7839, MPFR_RNDN);
        mpfr_sub(r7841, r7840, r7838, MPFR_RNDN);
        ;
        mpfr_set_si(r7843, mpfr_cmp(r7841, r7842) <= 0, MPFR_RNDN);
        mpfr_add(r7844, r7837, r7837, MPFR_RNDN);
        mpfr_fma(r7845, r7837, r7837, r7844, MPFR_RNDN);
        mpfr_expm1(r7846, r7845, MPFR_RNDN);
        mpfr_log1p(r7847, r7846, MPFR_RNDN);
        if (mpfr_get_si(r7843, MPFR_RNDN)) { mpfr_set(r7848, r7847, MPFR_RNDN); } else { mpfr_set(r7848, r7841, MPFR_RNDN); };
        return mpfr_get_d(r7848, MPFR_RNDN);
}