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

char *name = "NMSE problem 3.3.3";

double f_if(float x) {
        float r15974 = 1.0f;
        float r15975 = x;
        float r15976 = r15975 + r15974;
        float r15977 = r15974 / r15976;
        float r15978 = 2.0f;
        float r15979 = r15978 / r15975;
        float r15980 = r15977 - r15979;
        float r15981 = r15975 - r15974;
        float r15982 = r15974 / r15981;
        float r15983 = r15980 + r15982;
        return r15983;
}

double f_id(double x) {
        double r15984 = 1.0;
        double r15985 = x;
        double r15986 = r15985 + r15984;
        double r15987 = r15984 / r15986;
        double r15988 = 2.0;
        double r15989 = r15988 / r15985;
        double r15990 = r15987 - r15989;
        double r15991 = r15985 - r15984;
        double r15992 = r15984 / r15991;
        double r15993 = r15990 + r15992;
        return r15993;
}


double f_of(float x) {
        float r15994 = x;
        float r15995 = -2.633003449024554f;
        bool r15996 = r15994 <= r15995;
        float r15997 = 2.0f;
        float r15998 = 5.0f;
        float r15999 = pow(r15994, r15998);
        float r16000 = r15997 / r15999;
        float r16001 = 7.0f;
        float r16002 = pow(r15994, r16001);
        float r16003 = r15997 / r16002;
        float r16004 = r16000 + r16003;
        float r16005 = r15994 * (r15994 * r15994);
        float r16006 = r15997 / r16005;
        float r16007 = r16004 + r16006;
        float r16008 = 228205724867.64624f;
        bool r16009 = r15994 <= r16008;
        float r16010 = -r15997;
        float r16011 = 1.0f;
        float r16012 = fma(r15994, r15994, r16011);
        float r16013 = r16011 / r15994;
        float r16014 = fma(r16012, r15994, r16013);
        float r16015 = r16010 * r16014;
        float r16016 = r16009 ? r16015 : r16007;
        float r16017 = r15996 ? r16007 : r16016;
        return r16017;
}

double f_od(double x) {
        double r16018 = x;
        double r16019 = -2.633003449024554;
        bool r16020 = r16018 <= r16019;
        double r16021 = 2.0;
        double r16022 = 5.0;
        double r16023 = pow(r16018, r16022);
        double r16024 = r16021 / r16023;
        double r16025 = 7.0;
        double r16026 = pow(r16018, r16025);
        double r16027 = r16021 / r16026;
        double r16028 = r16024 + r16027;
        double r16029 = r16018 * (r16018 * r16018);
        double r16030 = r16021 / r16029;
        double r16031 = r16028 + r16030;
        double r16032 = 228205724867.64624;
        bool r16033 = r16018 <= r16032;
        double r16034 = -r16021;
        double r16035 = 1.0;
        double r16036 = fma(r16018, r16018, r16035);
        double r16037 = r16035 / r16018;
        double r16038 = fma(r16036, r16018, r16037);
        double r16039 = r16034 * r16038;
        double r16040 = r16033 ? r16039 : r16031;
        double r16041 = r16020 ? r16031 : r16040;
        return r16041;
}

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 r16042, r16043, r16044, r16045, r16046, r16047, r16048, r16049, r16050, r16051;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r16042, "1", 10, MPFR_RNDN);
        mpfr_init(r16043);
        mpfr_init(r16044);
        mpfr_init(r16045);
        mpfr_init_set_str(r16046, "2", 10, MPFR_RNDN);
        mpfr_init(r16047);
        mpfr_init(r16048);
        mpfr_init(r16049);
        mpfr_init(r16050);
        mpfr_init(r16051);
}

double f_im(double x) {
        ;
        mpfr_set_d(r16043, x, MPFR_RNDN);
        mpfr_add(r16044, r16043, r16042, MPFR_RNDN);
        mpfr_div(r16045, r16042, r16044, MPFR_RNDN);
        ;
        mpfr_div(r16047, r16046, r16043, MPFR_RNDN);
        mpfr_sub(r16048, r16045, r16047, MPFR_RNDN);
        mpfr_sub(r16049, r16043, r16042, MPFR_RNDN);
        mpfr_div(r16050, r16042, r16049, MPFR_RNDN);
        mpfr_add(r16051, r16048, r16050, MPFR_RNDN);
        return mpfr_get_d(r16051, MPFR_RNDN);
}

static mpfr_t r16052, r16053, r16054, r16055, r16056, r16057, r16058, r16059, r16060, r16061, r16062, r16063, r16064, r16065, r16066, r16067, r16068, r16069, r16070, r16071, r16072, r16073, r16074, r16075;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16052);
        mpfr_init_set_str(r16053, "-2.633003449024554", 10, MPFR_RNDN);
        mpfr_init(r16054);
        mpfr_init_set_str(r16055, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r16056, "5", 10, MPFR_RNDN);
        mpfr_init(r16057);
        mpfr_init(r16058);
        mpfr_init_set_str(r16059, "7", 10, MPFR_RNDN);
        mpfr_init(r16060);
        mpfr_init(r16061);
        mpfr_init(r16062);
        mpfr_init(r16063);
        mpfr_init(r16064);
        mpfr_init(r16065);
        mpfr_init_set_str(r16066, "228205724867.64624", 10, MPFR_RNDN);
        mpfr_init(r16067);
        mpfr_init(r16068);
        mpfr_init_set_str(r16069, "1", 10, MPFR_RNDN);
        mpfr_init(r16070);
        mpfr_init(r16071);
        mpfr_init(r16072);
        mpfr_init(r16073);
        mpfr_init(r16074);
        mpfr_init(r16075);
}

double f_fm(double x) {
        mpfr_set_d(r16052, x, MPFR_RNDN);
        ;
        mpfr_set_si(r16054, mpfr_cmp(r16052, r16053) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r16057, r16052, r16056, MPFR_RNDN);
        mpfr_div(r16058, r16055, r16057, MPFR_RNDN);
        ;
        mpfr_pow(r16060, r16052, r16059, MPFR_RNDN);
        mpfr_div(r16061, r16055, r16060, MPFR_RNDN);
        mpfr_add(r16062, r16058, r16061, MPFR_RNDN);
        mpfr_mul(r16063, r16052, r16052, MPFR_RNDN); mpfr_mul(r16063, r16063, r16052, MPFR_RNDN);
        mpfr_div(r16064, r16055, r16063, MPFR_RNDN);
        mpfr_add(r16065, r16062, r16064, MPFR_RNDN);
        ;
        mpfr_set_si(r16067, mpfr_cmp(r16052, r16066) <= 0, MPFR_RNDN);
        mpfr_neg(r16068, r16055, MPFR_RNDN);
        ;
        mpfr_fma(r16070, r16052, r16052, r16069, MPFR_RNDN);
        mpfr_div(r16071, r16069, r16052, MPFR_RNDN);
        mpfr_fma(r16072, r16070, r16052, r16071, MPFR_RNDN);
        mpfr_mul(r16073, r16068, r16072, MPFR_RNDN);
        if (mpfr_get_si(r16067, MPFR_RNDN)) { mpfr_set(r16074, r16073, MPFR_RNDN); } else { mpfr_set(r16074, r16065, MPFR_RNDN); };
        if (mpfr_get_si(r16054, MPFR_RNDN)) { mpfr_set(r16075, r16065, MPFR_RNDN); } else { mpfr_set(r16075, r16074, MPFR_RNDN); };
        return mpfr_get_d(r16075, MPFR_RNDN);
}

static mpfr_t r16076, r16077, r16078, r16079, r16080, r16081, r16082, r16083, r16084, r16085, r16086, r16087, r16088, r16089, r16090, r16091, r16092, r16093, r16094, r16095, r16096, r16097, r16098, r16099;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16076);
        mpfr_init_set_str(r16077, "-2.633003449024554", 10, MPFR_RNDN);
        mpfr_init(r16078);
        mpfr_init_set_str(r16079, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r16080, "5", 10, MPFR_RNDN);
        mpfr_init(r16081);
        mpfr_init(r16082);
        mpfr_init_set_str(r16083, "7", 10, MPFR_RNDN);
        mpfr_init(r16084);
        mpfr_init(r16085);
        mpfr_init(r16086);
        mpfr_init(r16087);
        mpfr_init(r16088);
        mpfr_init(r16089);
        mpfr_init_set_str(r16090, "228205724867.64624", 10, MPFR_RNDN);
        mpfr_init(r16091);
        mpfr_init(r16092);
        mpfr_init_set_str(r16093, "1", 10, MPFR_RNDN);
        mpfr_init(r16094);
        mpfr_init(r16095);
        mpfr_init(r16096);
        mpfr_init(r16097);
        mpfr_init(r16098);
        mpfr_init(r16099);
}

double f_dm(double x) {
        mpfr_set_d(r16076, x, MPFR_RNDN);
        ;
        mpfr_set_si(r16078, mpfr_cmp(r16076, r16077) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r16081, r16076, r16080, MPFR_RNDN);
        mpfr_div(r16082, r16079, r16081, MPFR_RNDN);
        ;
        mpfr_pow(r16084, r16076, r16083, MPFR_RNDN);
        mpfr_div(r16085, r16079, r16084, MPFR_RNDN);
        mpfr_add(r16086, r16082, r16085, MPFR_RNDN);
        mpfr_mul(r16087, r16076, r16076, MPFR_RNDN); mpfr_mul(r16087, r16087, r16076, MPFR_RNDN);
        mpfr_div(r16088, r16079, r16087, MPFR_RNDN);
        mpfr_add(r16089, r16086, r16088, MPFR_RNDN);
        ;
        mpfr_set_si(r16091, mpfr_cmp(r16076, r16090) <= 0, MPFR_RNDN);
        mpfr_neg(r16092, r16079, MPFR_RNDN);
        ;
        mpfr_fma(r16094, r16076, r16076, r16093, MPFR_RNDN);
        mpfr_div(r16095, r16093, r16076, MPFR_RNDN);
        mpfr_fma(r16096, r16094, r16076, r16095, MPFR_RNDN);
        mpfr_mul(r16097, r16092, r16096, MPFR_RNDN);
        if (mpfr_get_si(r16091, MPFR_RNDN)) { mpfr_set(r16098, r16097, MPFR_RNDN); } else { mpfr_set(r16098, r16089, MPFR_RNDN); };
        if (mpfr_get_si(r16078, MPFR_RNDN)) { mpfr_set(r16099, r16089, MPFR_RNDN); } else { mpfr_set(r16099, r16098, MPFR_RNDN); };
        return mpfr_get_d(r16099, MPFR_RNDN);
}