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

char *name = "VandenBroeck and Keller, Equation (20)";

double f_if(float f) {
        float r21540 = 1;
        float r21541 = atan2(1.0, 0.0);
        float r21542 = 4;
        float r21543 = r21541 / r21542;
        float r21544 = r21540 / r21543;
        float r21545 = f;
        float r21546 = r21543 * r21545;
        float r21547 = exp(r21546);
        float r21548 = -r21546;
        float r21549 = exp(r21548);
        float r21550 = r21547 + r21549;
        float r21551 = r21547 - r21549;
        float r21552 = r21550 / r21551;
        float r21553 = log(r21552);
        float r21554 = r21544 * r21553;
        float r21555 = -r21554;
        return r21555;
}

double f_id(double f) {
        double r21556 = 1;
        double r21557 = atan2(1.0, 0.0);
        double r21558 = 4;
        double r21559 = r21557 / r21558;
        double r21560 = r21556 / r21559;
        double r21561 = f;
        double r21562 = r21559 * r21561;
        double r21563 = exp(r21562);
        double r21564 = -r21562;
        double r21565 = exp(r21564);
        double r21566 = r21563 + r21565;
        double r21567 = r21563 - r21565;
        double r21568 = r21566 / r21567;
        double r21569 = log(r21568);
        double r21570 = r21560 * r21569;
        double r21571 = -r21570;
        return r21571;
}


double f_of(float f) {
        float r21572 = 4;
        float r21573 = atan2(1.0, 0.0);
        float r21574 = r21572 / r21573;
        float r21575 = log(r21574);
        float r21576 = r21575 / r21573;
        float r21577 = r21572 * r21576;
        float r21578 = 1/12;
        float r21579 = f;
        float r21580 = 2;
        float r21581 = pow(r21579, r21580);
        float r21582 = r21573 * r21581;
        float r21583 = r21578 * r21582;
        float r21584 = r21577 + r21583;
        float r21585 = 7/5760;
        float r21586 = exp(1.0);
        float r21587 = 3;
        float r21588 = pow(r21573, r21587);
        float r21589 = log(r21588);
        float r21590 = pow(r21586, r21589);
        float r21591 = pow(r21579, r21572);
        float r21592 = r21590 * r21591;
        float r21593 = r21585 * r21592;
        float r21594 = log(r21579);
        float r21595 = r21594 / r21573;
        float r21596 = r21572 * r21595;
        float r21597 = r21593 + r21596;
        float r21598 = r21584 - r21597;
        float r21599 = -r21598;
        return r21599;
}

double f_od(double f) {
        double r21600 = 4;
        double r21601 = atan2(1.0, 0.0);
        double r21602 = r21600 / r21601;
        double r21603 = log(r21602);
        double r21604 = r21603 / r21601;
        double r21605 = r21600 * r21604;
        double r21606 = 1/12;
        double r21607 = f;
        double r21608 = 2;
        double r21609 = pow(r21607, r21608);
        double r21610 = r21601 * r21609;
        double r21611 = r21606 * r21610;
        double r21612 = r21605 + r21611;
        double r21613 = 7/5760;
        double r21614 = exp(1.0);
        double r21615 = 3;
        double r21616 = pow(r21601, r21615);
        double r21617 = log(r21616);
        double r21618 = pow(r21614, r21617);
        double r21619 = pow(r21607, r21600);
        double r21620 = r21618 * r21619;
        double r21621 = r21613 * r21620;
        double r21622 = log(r21607);
        double r21623 = r21622 / r21601;
        double r21624 = r21600 * r21623;
        double r21625 = r21621 + r21624;
        double r21626 = r21612 - r21625;
        double r21627 = -r21626;
        return r21627;
}

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 r21628, r21629, r21630, r21631, r21632, r21633, r21634, r21635, r21636, r21637, r21638, r21639, r21640, r21641, r21642, r21643;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21628, "1", 10, MPFR_RNDN);
        mpfr_init(r21629);
        mpfr_init_set_str(r21630, "4", 10, MPFR_RNDN);
        mpfr_init(r21631);
        mpfr_init(r21632);
        mpfr_init(r21633);
        mpfr_init(r21634);
        mpfr_init(r21635);
        mpfr_init(r21636);
        mpfr_init(r21637);
        mpfr_init(r21638);
        mpfr_init(r21639);
        mpfr_init(r21640);
        mpfr_init(r21641);
        mpfr_init(r21642);
        mpfr_init(r21643);
}

double f_im(double f) {
        ;
        mpfr_const_pi(r21629, MPFR_RNDN);
        ;
        mpfr_div(r21631, r21629, r21630, MPFR_RNDN);
        mpfr_div(r21632, r21628, r21631, MPFR_RNDN);
        mpfr_set_d(r21633, f, MPFR_RNDN);
        mpfr_mul(r21634, r21631, r21633, MPFR_RNDN);
        mpfr_exp(r21635, r21634, MPFR_RNDN);
        mpfr_neg(r21636, r21634, MPFR_RNDN);
        mpfr_exp(r21637, r21636, MPFR_RNDN);
        mpfr_add(r21638, r21635, r21637, MPFR_RNDN);
        mpfr_sub(r21639, r21635, r21637, MPFR_RNDN);
        mpfr_div(r21640, r21638, r21639, MPFR_RNDN);
        mpfr_log(r21641, r21640, MPFR_RNDN);
        mpfr_mul(r21642, r21632, r21641, MPFR_RNDN);
        mpfr_neg(r21643, r21642, MPFR_RNDN);
        return mpfr_get_d(r21643, MPFR_RNDN);
}

static mpfr_t r21644, r21645, r21646, r21647, r21648, r21649, r21650, r21651, r21652, r21653, r21654, r21655, r21656, r21657, r21658, r21659, r21660, r21661, r21662, r21663, r21664, r21665, r21666, r21667, r21668, r21669, r21670, r21671;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21644, "4", 10, MPFR_RNDN);
        mpfr_init(r21645);
        mpfr_init(r21646);
        mpfr_init(r21647);
        mpfr_init(r21648);
        mpfr_init(r21649);
        mpfr_init_set_str(r21650, "1/12", 10, MPFR_RNDN);
        mpfr_init(r21651);
        mpfr_init_set_str(r21652, "2", 10, MPFR_RNDN);
        mpfr_init(r21653);
        mpfr_init(r21654);
        mpfr_init(r21655);
        mpfr_init(r21656);
        mpfr_init_set_str(r21657, "7/5760", 10, MPFR_RNDN);
        mpfr_init(r21658);
        mpfr_init_set_str(r21659, "3", 10, MPFR_RNDN);
        mpfr_init(r21660);
        mpfr_init(r21661);
        mpfr_init(r21662);
        mpfr_init(r21663);
        mpfr_init(r21664);
        mpfr_init(r21665);
        mpfr_init(r21666);
        mpfr_init(r21667);
        mpfr_init(r21668);
        mpfr_init(r21669);
        mpfr_init(r21670);
        mpfr_init(r21671);
}

double f_fm(double f) {
        ;
        mpfr_const_pi(r21645, MPFR_RNDN);
        mpfr_div(r21646, r21644, r21645, MPFR_RNDN);
        mpfr_log(r21647, r21646, MPFR_RNDN);
        mpfr_div(r21648, r21647, r21645, MPFR_RNDN);
        mpfr_mul(r21649, r21644, r21648, MPFR_RNDN);
        ;
        mpfr_set_d(r21651, f, MPFR_RNDN);
        ;
        mpfr_pow(r21653, r21651, r21652, MPFR_RNDN);
        mpfr_mul(r21654, r21645, r21653, MPFR_RNDN);
        mpfr_mul(r21655, r21650, r21654, MPFR_RNDN);
        mpfr_add(r21656, r21649, r21655, MPFR_RNDN);
        ;
        mpfr_set_si(r21658, 1, MPFR_RNDN), mpfr_const_exp(r21658, r21658, MPFR_RNDN);
        ;
        mpfr_pow(r21660, r21645, r21659, MPFR_RNDN);
        mpfr_log(r21661, r21660, MPFR_RNDN);
        mpfr_pow(r21662, r21658, r21661, MPFR_RNDN);
        mpfr_pow(r21663, r21651, r21644, MPFR_RNDN);
        mpfr_mul(r21664, r21662, r21663, MPFR_RNDN);
        mpfr_mul(r21665, r21657, r21664, MPFR_RNDN);
        mpfr_log(r21666, r21651, MPFR_RNDN);
        mpfr_div(r21667, r21666, r21645, MPFR_RNDN);
        mpfr_mul(r21668, r21644, r21667, MPFR_RNDN);
        mpfr_add(r21669, r21665, r21668, MPFR_RNDN);
        mpfr_sub(r21670, r21656, r21669, MPFR_RNDN);
        mpfr_neg(r21671, r21670, MPFR_RNDN);
        return mpfr_get_d(r21671, MPFR_RNDN);
}

static mpfr_t r21672, r21673, r21674, r21675, r21676, r21677, r21678, r21679, r21680, r21681, r21682, r21683, r21684, r21685, r21686, r21687, r21688, r21689, r21690, r21691, r21692, r21693, r21694, r21695, r21696, r21697, r21698, r21699;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21672, "4", 10, MPFR_RNDN);
        mpfr_init(r21673);
        mpfr_init(r21674);
        mpfr_init(r21675);
        mpfr_init(r21676);
        mpfr_init(r21677);
        mpfr_init_set_str(r21678, "1/12", 10, MPFR_RNDN);
        mpfr_init(r21679);
        mpfr_init_set_str(r21680, "2", 10, MPFR_RNDN);
        mpfr_init(r21681);
        mpfr_init(r21682);
        mpfr_init(r21683);
        mpfr_init(r21684);
        mpfr_init_set_str(r21685, "7/5760", 10, MPFR_RNDN);
        mpfr_init(r21686);
        mpfr_init_set_str(r21687, "3", 10, MPFR_RNDN);
        mpfr_init(r21688);
        mpfr_init(r21689);
        mpfr_init(r21690);
        mpfr_init(r21691);
        mpfr_init(r21692);
        mpfr_init(r21693);
        mpfr_init(r21694);
        mpfr_init(r21695);
        mpfr_init(r21696);
        mpfr_init(r21697);
        mpfr_init(r21698);
        mpfr_init(r21699);
}

double f_dm(double f) {
        ;
        mpfr_const_pi(r21673, MPFR_RNDN);
        mpfr_div(r21674, r21672, r21673, MPFR_RNDN);
        mpfr_log(r21675, r21674, MPFR_RNDN);
        mpfr_div(r21676, r21675, r21673, MPFR_RNDN);
        mpfr_mul(r21677, r21672, r21676, MPFR_RNDN);
        ;
        mpfr_set_d(r21679, f, MPFR_RNDN);
        ;
        mpfr_pow(r21681, r21679, r21680, MPFR_RNDN);
        mpfr_mul(r21682, r21673, r21681, MPFR_RNDN);
        mpfr_mul(r21683, r21678, r21682, MPFR_RNDN);
        mpfr_add(r21684, r21677, r21683, MPFR_RNDN);
        ;
        mpfr_set_si(r21686, 1, MPFR_RNDN), mpfr_const_exp(r21686, r21686, MPFR_RNDN);
        ;
        mpfr_pow(r21688, r21673, r21687, MPFR_RNDN);
        mpfr_log(r21689, r21688, MPFR_RNDN);
        mpfr_pow(r21690, r21686, r21689, MPFR_RNDN);
        mpfr_pow(r21691, r21679, r21672, MPFR_RNDN);
        mpfr_mul(r21692, r21690, r21691, MPFR_RNDN);
        mpfr_mul(r21693, r21685, r21692, MPFR_RNDN);
        mpfr_log(r21694, r21679, MPFR_RNDN);
        mpfr_div(r21695, r21694, r21673, MPFR_RNDN);
        mpfr_mul(r21696, r21672, r21695, MPFR_RNDN);
        mpfr_add(r21697, r21693, r21696, MPFR_RNDN);
        mpfr_sub(r21698, r21684, r21697, MPFR_RNDN);
        mpfr_neg(r21699, r21698, MPFR_RNDN);
        return mpfr_get_d(r21699, MPFR_RNDN);
}