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

char *name = "Rosa's Benchmark";

double f_if(float x) {
        float r20971 = 0.954929658551372;
        float r20972 = x;
        float r20973 = r20971 * r20972;
        float r20974 = 0.12900613773279798;
        float r20975 = r20972 * r20972;
        float r20976 = r20975 * r20972;
        float r20977 = r20974 * r20976;
        float r20978 = r20973 - r20977;
        return r20978;
}

double f_id(double x) {
        double r20979 = 0.954929658551372;
        double r20980 = x;
        double r20981 = r20979 * r20980;
        double r20982 = 0.12900613773279798;
        double r20983 = r20980 * r20980;
        double r20984 = r20983 * r20980;
        double r20985 = r20982 * r20984;
        double r20986 = r20981 - r20985;
        return r20986;
}


double f_of(float x) {
        float r20987 = 0.954929658551372;
        float r20988 = x;
        float r20989 = 3;
        float r20990 = pow(r20988, r20989);
        float r20991 = 0.12900613773279798;
        float r20992 = -r20991;
        float r20993 = r20990 * r20992;
        float r20994 = fma(r20987, r20988, r20993);
        return r20994;
}

double f_od(double x) {
        double r20995 = 0.954929658551372;
        double r20996 = x;
        double r20997 = 3;
        double r20998 = pow(r20996, r20997);
        double r20999 = 0.12900613773279798;
        double r21000 = -r20999;
        double r21001 = r20998 * r21000;
        double r21002 = fma(r20995, r20996, r21001);
        return r21002;
}

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 r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r21003, "0.954929658551372", 10, MPFR_RNDN);
        mpfr_init(r21004);
        mpfr_init(r21005);
        mpfr_init_set_str(r21006, "0.12900613773279798", 10, MPFR_RNDN);
        mpfr_init(r21007);
        mpfr_init(r21008);
        mpfr_init(r21009);
        mpfr_init(r21010);
}

double f_im(double x) {
        ;
        mpfr_set_d(r21004, x, MPFR_RNDN);
        mpfr_mul(r21005, r21003, r21004, MPFR_RNDN);
        ;
        mpfr_mul(r21007, r21004, r21004, MPFR_RNDN);
        mpfr_mul(r21008, r21007, r21004, MPFR_RNDN);
        mpfr_mul(r21009, r21006, r21008, MPFR_RNDN);
        mpfr_sub(r21010, r21005, r21009, MPFR_RNDN);
        return mpfr_get_d(r21010, MPFR_RNDN);
}

static mpfr_t r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r21011, "0.954929658551372", 10, MPFR_RNDN);
        mpfr_init(r21012);
        mpfr_init_set_str(r21013, "3", 10, MPFR_RNDN);
        mpfr_init(r21014);
        mpfr_init_set_str(r21015, "0.12900613773279798", 10, MPFR_RNDN);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r21012, x, MPFR_RNDN);
        ;
        mpfr_pow(r21014, r21012, r21013, MPFR_RNDN);
        ;
        mpfr_neg(r21016, r21015, MPFR_RNDN);
        mpfr_mul(r21017, r21014, r21016, MPFR_RNDN);
        mpfr_fma(r21018, r21011, r21012, r21017, MPFR_RNDN);
        return mpfr_get_d(r21018, MPFR_RNDN);
}

static mpfr_t r21019, r21020, r21021, r21022, r21023, r21024, r21025, r21026;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r21019, "0.954929658551372", 10, MPFR_RNDN);
        mpfr_init(r21020);
        mpfr_init_set_str(r21021, "3", 10, MPFR_RNDN);
        mpfr_init(r21022);
        mpfr_init_set_str(r21023, "0.12900613773279798", 10, MPFR_RNDN);
        mpfr_init(r21024);
        mpfr_init(r21025);
        mpfr_init(r21026);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r21020, x, MPFR_RNDN);
        ;
        mpfr_pow(r21022, r21020, r21021, MPFR_RNDN);
        ;
        mpfr_neg(r21024, r21023, MPFR_RNDN);
        mpfr_mul(r21025, r21022, r21024, MPFR_RNDN);
        mpfr_fma(r21026, r21019, r21020, r21025, MPFR_RNDN);
        return mpfr_get_d(r21026, MPFR_RNDN);
}