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

char *name = "logs (example 3.8)";

double f_if(float n) {
        float r26403 = n;
        float r26404 = 1;
        float r26405 = r26403 + r26404;
        float r26406 = log(r26405);
        float r26407 = r26405 * r26406;
        float r26408 = log(r26403);
        float r26409 = r26403 * r26408;
        float r26410 = r26407 - r26409;
        float r26411 = r26410 - r26404;
        return r26411;
}

double f_id(double n) {
        double r26412 = n;
        double r26413 = 1;
        double r26414 = r26412 + r26413;
        double r26415 = log(r26414);
        double r26416 = r26414 * r26415;
        double r26417 = log(r26412);
        double r26418 = r26412 * r26417;
        double r26419 = r26416 - r26418;
        double r26420 = r26419 - r26413;
        return r26420;
}


double f_of(float n) {
        float r26421 = n;
        float r26422 = log(r26421);
        float r26423 = 1/6;
        float r26424 = r26423 / r26421;
        float r26425 = r26424 / r26421;
        float r26426 = r26422 - r26425;
        float r26427 = 1/2;
        float r26428 = r26427 / r26421;
        float r26429 = r26426 + r26428;
        return r26429;
}

double f_od(double n) {
        double r26430 = n;
        double r26431 = log(r26430);
        double r26432 = 1/6;
        double r26433 = r26432 / r26430;
        double r26434 = r26433 / r26430;
        double r26435 = r26431 - r26434;
        double r26436 = 1/2;
        double r26437 = r26436 / r26430;
        double r26438 = r26435 + r26437;
        return r26438;
}

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 r26439, r26440, r26441, r26442, r26443, r26444, r26445, r26446, r26447;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26439);
        mpfr_init_set_str(r26440, "1", 10, MPFR_RNDN);
        mpfr_init(r26441);
        mpfr_init(r26442);
        mpfr_init(r26443);
        mpfr_init(r26444);
        mpfr_init(r26445);
        mpfr_init(r26446);
        mpfr_init(r26447);
}

double f_im(double n) {
        mpfr_set_d(r26439, n, MPFR_RNDN);
        ;
        mpfr_add(r26441, r26439, r26440, MPFR_RNDN);
        mpfr_log(r26442, r26441, MPFR_RNDN);
        mpfr_mul(r26443, r26441, r26442, MPFR_RNDN);
        mpfr_log(r26444, r26439, MPFR_RNDN);
        mpfr_mul(r26445, r26439, r26444, MPFR_RNDN);
        mpfr_sub(r26446, r26443, r26445, MPFR_RNDN);
        mpfr_sub(r26447, r26446, r26440, MPFR_RNDN);
        return mpfr_get_d(r26447, MPFR_RNDN);
}

static mpfr_t r26448, r26449, r26450, r26451, r26452, r26453, r26454, r26455, r26456;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26448);
        mpfr_init(r26449);
        mpfr_init_set_str(r26450, "1/6", 10, MPFR_RNDN);
        mpfr_init(r26451);
        mpfr_init(r26452);
        mpfr_init(r26453);
        mpfr_init_set_str(r26454, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26455);
        mpfr_init(r26456);
}

double f_fm(double n) {
        mpfr_set_d(r26448, n, MPFR_RNDN);
        mpfr_log(r26449, r26448, MPFR_RNDN);
        ;
        mpfr_div(r26451, r26450, r26448, MPFR_RNDN);
        mpfr_div(r26452, r26451, r26448, MPFR_RNDN);
        mpfr_sub(r26453, r26449, r26452, MPFR_RNDN);
        ;
        mpfr_div(r26455, r26454, r26448, MPFR_RNDN);
        mpfr_add(r26456, r26453, r26455, MPFR_RNDN);
        return mpfr_get_d(r26456, MPFR_RNDN);
}

static mpfr_t r26457, r26458, r26459, r26460, r26461, r26462, r26463, r26464, r26465;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26457);
        mpfr_init(r26458);
        mpfr_init_set_str(r26459, "1/6", 10, MPFR_RNDN);
        mpfr_init(r26460);
        mpfr_init(r26461);
        mpfr_init(r26462);
        mpfr_init_set_str(r26463, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26464);
        mpfr_init(r26465);
}

double f_dm(double n) {
        mpfr_set_d(r26457, n, MPFR_RNDN);
        mpfr_log(r26458, r26457, MPFR_RNDN);
        ;
        mpfr_div(r26460, r26459, r26457, MPFR_RNDN);
        mpfr_div(r26461, r26460, r26457, MPFR_RNDN);
        mpfr_sub(r26462, r26458, r26461, MPFR_RNDN);
        ;
        mpfr_div(r26464, r26463, r26457, MPFR_RNDN);
        mpfr_add(r26465, r26462, r26464, MPFR_RNDN);
        return mpfr_get_d(r26465, MPFR_RNDN);
}