#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 r26516 = 1;
        float r26517 = eps;
        float r26518 = r26516 - r26517;
        float r26519 = r26516 + r26517;
        float r26520 = r26518 / r26519;
        float r26521 = log(r26520);
        return r26521;
}

double f_id(double eps) {
        double r26522 = 1;
        double r26523 = eps;
        double r26524 = r26522 - r26523;
        double r26525 = r26522 + r26523;
        double r26526 = r26524 / r26525;
        double r26527 = log(r26526);
        return r26527;
}


double f_of(float eps) {
        float r26528 = -2;
        float r26529 = eps;
        float r26530 = r26528 * r26529;
        float r26531 = 4;
        float r26532 = pow(r26529, r26531);
        float r26533 = r26530 * r26532;
        float r26534 = 1/5;
        float r26535 = 1/3;
        float r26536 = r26535 * r26528;
        float r26537 = 3;
        float r26538 = pow(r26529, r26537);
        float r26539 = fma(r26536, r26538, r26530);
        float r26540 = fma(r26533, r26534, r26539);
        return r26540;
}

double f_od(double eps) {
        double r26541 = -2;
        double r26542 = eps;
        double r26543 = r26541 * r26542;
        double r26544 = 4;
        double r26545 = pow(r26542, r26544);
        double r26546 = r26543 * r26545;
        double r26547 = 1/5;
        double r26548 = 1/3;
        double r26549 = r26548 * r26541;
        double r26550 = 3;
        double r26551 = pow(r26542, r26550);
        double r26552 = fma(r26549, r26551, r26543);
        double r26553 = fma(r26546, r26547, r26552);
        return r26553;
}

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 r26554, r26555, r26556, r26557, r26558, r26559;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26554, "1", 10, MPFR_RNDN);
        mpfr_init(r26555);
        mpfr_init(r26556);
        mpfr_init(r26557);
        mpfr_init(r26558);
        mpfr_init(r26559);
}

double f_im(double eps) {
        ;
        mpfr_set_d(r26555, eps, MPFR_RNDN);
        mpfr_sub(r26556, r26554, r26555, MPFR_RNDN);
        mpfr_add(r26557, r26554, r26555, MPFR_RNDN);
        mpfr_div(r26558, r26556, r26557, MPFR_RNDN);
        mpfr_log(r26559, r26558, MPFR_RNDN);
        return mpfr_get_d(r26559, MPFR_RNDN);
}

static mpfr_t r26560, r26561, r26562, r26563, r26564, r26565, r26566, r26567, r26568, r26569, r26570, r26571, r26572;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26560, "-2", 10, MPFR_RNDN);
        mpfr_init(r26561);
        mpfr_init(r26562);
        mpfr_init_set_str(r26563, "4", 10, MPFR_RNDN);
        mpfr_init(r26564);
        mpfr_init(r26565);
        mpfr_init_set_str(r26566, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r26567, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26568);
        mpfr_init_set_str(r26569, "3", 10, MPFR_RNDN);
        mpfr_init(r26570);
        mpfr_init(r26571);
        mpfr_init(r26572);
}

double f_fm(double eps) {
        ;
        mpfr_set_d(r26561, eps, MPFR_RNDN);
        mpfr_mul(r26562, r26560, r26561, MPFR_RNDN);
        ;
        mpfr_pow(r26564, r26561, r26563, MPFR_RNDN);
        mpfr_mul(r26565, r26562, r26564, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26568, r26567, r26560, MPFR_RNDN);
        ;
        mpfr_pow(r26570, r26561, r26569, MPFR_RNDN);
        mpfr_fma(r26571, r26568, r26570, r26562, MPFR_RNDN);
        mpfr_fma(r26572, r26565, r26566, r26571, MPFR_RNDN);
        return mpfr_get_d(r26572, MPFR_RNDN);
}

static mpfr_t r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584, r26585;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26573, "-2", 10, MPFR_RNDN);
        mpfr_init(r26574);
        mpfr_init(r26575);
        mpfr_init_set_str(r26576, "4", 10, MPFR_RNDN);
        mpfr_init(r26577);
        mpfr_init(r26578);
        mpfr_init_set_str(r26579, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r26580, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26581);
        mpfr_init_set_str(r26582, "3", 10, MPFR_RNDN);
        mpfr_init(r26583);
        mpfr_init(r26584);
        mpfr_init(r26585);
}

double f_dm(double eps) {
        ;
        mpfr_set_d(r26574, eps, MPFR_RNDN);
        mpfr_mul(r26575, r26573, r26574, MPFR_RNDN);
        ;
        mpfr_pow(r26577, r26574, r26576, MPFR_RNDN);
        mpfr_mul(r26578, r26575, r26577, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26581, r26580, r26573, MPFR_RNDN);
        ;
        mpfr_pow(r26583, r26574, r26582, MPFR_RNDN);
        mpfr_fma(r26584, r26581, r26583, r26575, MPFR_RNDN);
        mpfr_fma(r26585, r26578, r26579, r26584, MPFR_RNDN);
        return mpfr_get_d(r26585, MPFR_RNDN);
}