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

char *name = "logq (problem 3.4.3)";

double f_if(float eps) {
        float r10009 = 1.0f;
        float r10010 = eps;
        float r10011 = r10009 - r10010;
        float r10012 = r10009 + r10010;
        float r10013 = r10011 / r10012;
        float r10014 = log(r10013);
        return r10014;
}

double f_id(double eps) {
        double r10015 = 1.0;
        double r10016 = eps;
        double r10017 = r10015 - r10016;
        double r10018 = r10015 + r10016;
        double r10019 = r10017 / r10018;
        double r10020 = log(r10019);
        return r10020;
}


double f_of(float eps) {
        float r10021 = eps;
        float r10022 = r10021 * r10021;
        float r10023 = 0.6666666666666666f;
        float r10024 = 2.0f;
        float r10025 = fma(r10022, r10023, r10024);
        float r10026 = -r10021;
        float r10027 = 5.0f;
        float r10028 = pow(r10021, r10027);
        float r10029 = 0.4f;
        float r10030 = -r10029;
        float r10031 = r10028 * r10030;
        float r10032 = fma(r10025, r10026, r10031);
        return r10032;
}

double f_od(double eps) {
        double r10033 = eps;
        double r10034 = r10033 * r10033;
        double r10035 = 0.6666666666666666;
        double r10036 = 2.0;
        double r10037 = fma(r10034, r10035, r10036);
        double r10038 = -r10033;
        double r10039 = 5.0;
        double r10040 = pow(r10033, r10039);
        double r10041 = 0.4;
        double r10042 = -r10041;
        double r10043 = r10040 * r10042;
        double r10044 = fma(r10037, r10038, r10043);
        return r10044;
}

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 r10045, r10046, r10047, r10048, r10049, r10050;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10045, "1", 10, MPFR_RNDN);
        mpfr_init(r10046);
        mpfr_init(r10047);
        mpfr_init(r10048);
        mpfr_init(r10049);
        mpfr_init(r10050);
}

double f_im(double eps) {
        ;
        mpfr_set_d(r10046, eps, MPFR_RNDN);
        mpfr_sub(r10047, r10045, r10046, MPFR_RNDN);
        mpfr_add(r10048, r10045, r10046, MPFR_RNDN);
        mpfr_div(r10049, r10047, r10048, MPFR_RNDN);
        mpfr_log(r10050, r10049, MPFR_RNDN);
        return mpfr_get_d(r10050, MPFR_RNDN);
}

static mpfr_t r10051, r10052, r10053, r10054, r10055, r10056, r10057, r10058, r10059, r10060, r10061, r10062;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10051);
        mpfr_init(r10052);
        mpfr_init_set_str(r10053, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r10054, "2", 10, MPFR_RNDN);
        mpfr_init(r10055);
        mpfr_init(r10056);
        mpfr_init_set_str(r10057, "5", 10, MPFR_RNDN);
        mpfr_init(r10058);
        mpfr_init_set_str(r10059, "2/5", 10, MPFR_RNDN);
        mpfr_init(r10060);
        mpfr_init(r10061);
        mpfr_init(r10062);
}

double f_fm(double eps) {
        mpfr_set_d(r10051, eps, MPFR_RNDN);
        mpfr_mul(r10052, r10051, r10051, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r10055, r10052, r10053, r10054, MPFR_RNDN);
        mpfr_neg(r10056, r10051, MPFR_RNDN);
        ;
        mpfr_pow(r10058, r10051, r10057, MPFR_RNDN);
        ;
        mpfr_neg(r10060, r10059, MPFR_RNDN);
        mpfr_mul(r10061, r10058, r10060, MPFR_RNDN);
        mpfr_fma(r10062, r10055, r10056, r10061, MPFR_RNDN);
        return mpfr_get_d(r10062, MPFR_RNDN);
}

static mpfr_t r10063, r10064, r10065, r10066, r10067, r10068, r10069, r10070, r10071, r10072, r10073, r10074;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10063);
        mpfr_init(r10064);
        mpfr_init_set_str(r10065, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r10066, "2", 10, MPFR_RNDN);
        mpfr_init(r10067);
        mpfr_init(r10068);
        mpfr_init_set_str(r10069, "5", 10, MPFR_RNDN);
        mpfr_init(r10070);
        mpfr_init_set_str(r10071, "2/5", 10, MPFR_RNDN);
        mpfr_init(r10072);
        mpfr_init(r10073);
        mpfr_init(r10074);
}

double f_dm(double eps) {
        mpfr_set_d(r10063, eps, MPFR_RNDN);
        mpfr_mul(r10064, r10063, r10063, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r10067, r10064, r10065, r10066, MPFR_RNDN);
        mpfr_neg(r10068, r10063, MPFR_RNDN);
        ;
        mpfr_pow(r10070, r10063, r10069, MPFR_RNDN);
        ;
        mpfr_neg(r10072, r10071, MPFR_RNDN);
        mpfr_mul(r10073, r10070, r10072, MPFR_RNDN);
        mpfr_fma(r10074, r10067, r10068, r10073, MPFR_RNDN);
        return mpfr_get_d(r10074, MPFR_RNDN);
}