#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 r9990 = 1.0f;
        float r9991 = eps;
        float r9992 = r9990 - r9991;
        float r9993 = r9990 + r9991;
        float r9994 = r9992 / r9993;
        float r9995 = log(r9994);
        return r9995;
}

double f_id(double eps) {
        double r9996 = 1.0;
        double r9997 = eps;
        double r9998 = r9996 - r9997;
        double r9999 = r9996 + r9997;
        double r10000 = r9998 / r9999;
        double r10001 = log(r10000);
        return r10001;
}


double f_of(float eps) {
        float r10002 = eps;
        float r10003 = r10002 * r10002;
        float r10004 = 0.6666666666666666f;
        float r10005 = 2.0f;
        float r10006 = fma(r10003, r10004, r10005);
        float r10007 = -r10002;
        float r10008 = 5.0f;
        float r10009 = pow(r10002, r10008);
        float r10010 = 0.4f;
        float r10011 = -r10010;
        float r10012 = r10009 * r10011;
        float r10013 = fma(r10006, r10007, r10012);
        return r10013;
}

double f_od(double eps) {
        double r10014 = eps;
        double r10015 = r10014 * r10014;
        double r10016 = 0.6666666666666666;
        double r10017 = 2.0;
        double r10018 = fma(r10015, r10016, r10017);
        double r10019 = -r10014;
        double r10020 = 5.0;
        double r10021 = pow(r10014, r10020);
        double r10022 = 0.4;
        double r10023 = -r10022;
        double r10024 = r10021 * r10023;
        double r10025 = fma(r10018, r10019, r10024);
        return r10025;
}

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 r10026, r10027, r10028, r10029, r10030, r10031;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10026, "1", 10, MPFR_RNDN);
        mpfr_init(r10027);
        mpfr_init(r10028);
        mpfr_init(r10029);
        mpfr_init(r10030);
        mpfr_init(r10031);
}

double f_im(double eps) {
        ;
        mpfr_set_d(r10027, eps, MPFR_RNDN);
        mpfr_sub(r10028, r10026, r10027, MPFR_RNDN);
        mpfr_add(r10029, r10026, r10027, MPFR_RNDN);
        mpfr_div(r10030, r10028, r10029, MPFR_RNDN);
        mpfr_log(r10031, r10030, MPFR_RNDN);
        return mpfr_get_d(r10031, MPFR_RNDN);
}

static mpfr_t r10032, r10033, r10034, r10035, r10036, r10037, r10038, r10039, r10040, r10041, r10042, r10043;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10032);
        mpfr_init(r10033);
        mpfr_init_set_str(r10034, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r10035, "2", 10, MPFR_RNDN);
        mpfr_init(r10036);
        mpfr_init(r10037);
        mpfr_init_set_str(r10038, "5", 10, MPFR_RNDN);
        mpfr_init(r10039);
        mpfr_init_set_str(r10040, "2/5", 10, MPFR_RNDN);
        mpfr_init(r10041);
        mpfr_init(r10042);
        mpfr_init(r10043);
}

double f_fm(double eps) {
        mpfr_set_d(r10032, eps, MPFR_RNDN);
        mpfr_mul(r10033, r10032, r10032, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r10036, r10033, r10034, r10035, MPFR_RNDN);
        mpfr_neg(r10037, r10032, MPFR_RNDN);
        ;
        mpfr_pow(r10039, r10032, r10038, MPFR_RNDN);
        ;
        mpfr_neg(r10041, r10040, MPFR_RNDN);
        mpfr_mul(r10042, r10039, r10041, MPFR_RNDN);
        mpfr_fma(r10043, r10036, r10037, r10042, MPFR_RNDN);
        return mpfr_get_d(r10043, MPFR_RNDN);
}

static mpfr_t r10044, r10045, r10046, r10047, r10048, r10049, r10050, r10051, r10052, r10053, r10054, r10055;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10044);
        mpfr_init(r10045);
        mpfr_init_set_str(r10046, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r10047, "2", 10, MPFR_RNDN);
        mpfr_init(r10048);
        mpfr_init(r10049);
        mpfr_init_set_str(r10050, "5", 10, MPFR_RNDN);
        mpfr_init(r10051);
        mpfr_init_set_str(r10052, "2/5", 10, MPFR_RNDN);
        mpfr_init(r10053);
        mpfr_init(r10054);
        mpfr_init(r10055);
}

double f_dm(double eps) {
        mpfr_set_d(r10044, eps, MPFR_RNDN);
        mpfr_mul(r10045, r10044, r10044, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r10048, r10045, r10046, r10047, MPFR_RNDN);
        mpfr_neg(r10049, r10044, MPFR_RNDN);
        ;
        mpfr_pow(r10051, r10044, r10050, MPFR_RNDN);
        ;
        mpfr_neg(r10053, r10052, MPFR_RNDN);
        mpfr_mul(r10054, r10051, r10053, MPFR_RNDN);
        mpfr_fma(r10055, r10048, r10049, r10054, MPFR_RNDN);
        return mpfr_get_d(r10055, MPFR_RNDN);
}