#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 r26401 = n;
        float r26402 = 1;
        float r26403 = r26401 + r26402;
        float r26404 = log(r26403);
        float r26405 = r26403 * r26404;
        float r26406 = log(r26401);
        float r26407 = r26401 * r26406;
        float r26408 = r26405 - r26407;
        float r26409 = r26408 - r26402;
        return r26409;
}

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


double f_of(float n) {
        float r26419 = n;
        float r26420 = log(r26419);
        float r26421 = 1/2;
        float r26422 = r26421 / r26419;
        float r26423 = r26420 + r26422;
        float r26424 = 1;
        float r26425 = r26424 - r26424;
        float r26426 = r26423 + r26425;
        return r26426;
}

double f_od(double n) {
        double r26427 = n;
        double r26428 = log(r26427);
        double r26429 = 1/2;
        double r26430 = r26429 / r26427;
        double r26431 = r26428 + r26430;
        double r26432 = 1;
        double r26433 = r26432 - r26432;
        double r26434 = r26431 + r26433;
        return r26434;
}

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 r26435, r26436, r26437, r26438, r26439, r26440, r26441, r26442, r26443;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26435);
        mpfr_init_set_str(r26436, "1", 10, MPFR_RNDN);
        mpfr_init(r26437);
        mpfr_init(r26438);
        mpfr_init(r26439);
        mpfr_init(r26440);
        mpfr_init(r26441);
        mpfr_init(r26442);
        mpfr_init(r26443);
}

double f_im(double n) {
        mpfr_set_d(r26435, n, MPFR_RNDN);
        ;
        mpfr_add(r26437, r26435, r26436, MPFR_RNDN);
        mpfr_log(r26438, r26437, MPFR_RNDN);
        mpfr_mul(r26439, r26437, r26438, MPFR_RNDN);
        mpfr_log(r26440, r26435, MPFR_RNDN);
        mpfr_mul(r26441, r26435, r26440, MPFR_RNDN);
        mpfr_sub(r26442, r26439, r26441, MPFR_RNDN);
        mpfr_sub(r26443, r26442, r26436, MPFR_RNDN);
        return mpfr_get_d(r26443, MPFR_RNDN);
}

static mpfr_t r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26444);
        mpfr_init(r26445);
        mpfr_init_set_str(r26446, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26447);
        mpfr_init(r26448);
        mpfr_init_set_str(r26449, "1", 10, MPFR_RNDN);
        mpfr_init(r26450);
        mpfr_init(r26451);
}

double f_fm(double n) {
        mpfr_set_d(r26444, n, MPFR_RNDN);
        mpfr_log(r26445, r26444, MPFR_RNDN);
        ;
        mpfr_div(r26447, r26446, r26444, MPFR_RNDN);
        mpfr_add(r26448, r26445, r26447, MPFR_RNDN);
        ;
        mpfr_sub(r26450, r26449, r26449, MPFR_RNDN);
        mpfr_add(r26451, r26448, r26450, MPFR_RNDN);
        return mpfr_get_d(r26451, MPFR_RNDN);
}

static mpfr_t r26452, r26453, r26454, r26455, r26456, r26457, r26458, r26459;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26452);
        mpfr_init(r26453);
        mpfr_init_set_str(r26454, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26455);
        mpfr_init(r26456);
        mpfr_init_set_str(r26457, "1", 10, MPFR_RNDN);
        mpfr_init(r26458);
        mpfr_init(r26459);
}

double f_dm(double n) {
        mpfr_set_d(r26452, n, MPFR_RNDN);
        mpfr_log(r26453, r26452, MPFR_RNDN);
        ;
        mpfr_div(r26455, r26454, r26452, MPFR_RNDN);
        mpfr_add(r26456, r26453, r26455, MPFR_RNDN);
        ;
        mpfr_sub(r26458, r26457, r26457, MPFR_RNDN);
        mpfr_add(r26459, r26456, r26458, MPFR_RNDN);
        return mpfr_get_d(r26459, MPFR_RNDN);
}