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

char *name = "Falkner and Boettcher, Appendix B, 1";

double f_if(float v) {
        float r25397 = 1;
        float r25398 = 5;
        float r25399 = v;
        float r25400 = r25399 * r25399;
        float r25401 = r25398 * r25400;
        float r25402 = r25397 - r25401;
        float r25403 = r25400 - r25397;
        float r25404 = r25402 / r25403;
        float r25405 = acos(r25404);
        return r25405;
}

double f_id(double v) {
        double r25406 = 1;
        double r25407 = 5;
        double r25408 = v;
        double r25409 = r25408 * r25408;
        double r25410 = r25407 * r25409;
        double r25411 = r25406 - r25410;
        double r25412 = r25409 - r25406;
        double r25413 = r25411 / r25412;
        double r25414 = acos(r25413);
        return r25414;
}


double f_of(float v) {
        float r25415 = exp(1.0);
        float r25416 = 1;
        float r25417 = 5;
        float r25418 = v;
        float r25419 = r25418 * r25418;
        float r25420 = r25417 * r25419;
        float r25421 = r25416 - r25420;
        float r25422 = r25419 - r25416;
        float r25423 = r25421 / r25422;
        float r25424 = acos(r25423);
        float r25425 = log(r25424);
        float r25426 = pow(r25415, r25425);
        return r25426;
}

double f_od(double v) {
        double r25427 = exp(1.0);
        double r25428 = 1;
        double r25429 = 5;
        double r25430 = v;
        double r25431 = r25430 * r25430;
        double r25432 = r25429 * r25431;
        double r25433 = r25428 - r25432;
        double r25434 = r25431 - r25428;
        double r25435 = r25433 / r25434;
        double r25436 = acos(r25435);
        double r25437 = log(r25436);
        double r25438 = pow(r25427, r25437);
        return r25438;
}

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 r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r25439, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r25440, "5", 10, MPFR_RNDN);
        mpfr_init(r25441);
        mpfr_init(r25442);
        mpfr_init(r25443);
        mpfr_init(r25444);
        mpfr_init(r25445);
        mpfr_init(r25446);
        mpfr_init(r25447);
}

double f_im(double v) {
        ;
        ;
        mpfr_set_d(r25441, v, MPFR_RNDN);
        mpfr_mul(r25442, r25441, r25441, MPFR_RNDN);
        mpfr_mul(r25443, r25440, r25442, MPFR_RNDN);
        mpfr_sub(r25444, r25439, r25443, MPFR_RNDN);
        mpfr_sub(r25445, r25442, r25439, MPFR_RNDN);
        mpfr_div(r25446, r25444, r25445, MPFR_RNDN);
        mpfr_acos(r25447, r25446, MPFR_RNDN);
        return mpfr_get_d(r25447, MPFR_RNDN);
}

static mpfr_t r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25448);
        mpfr_init_set_str(r25449, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r25450, "5", 10, MPFR_RNDN);
        mpfr_init(r25451);
        mpfr_init(r25452);
        mpfr_init(r25453);
        mpfr_init(r25454);
        mpfr_init(r25455);
        mpfr_init(r25456);
        mpfr_init(r25457);
        mpfr_init(r25458);
        mpfr_init(r25459);
}

double f_fm(double v) {
        mpfr_set_si(r25448, 1, MPFR_RNDN), mpfr_const_exp(r25448, r25448, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r25451, v, MPFR_RNDN);
        mpfr_mul(r25452, r25451, r25451, MPFR_RNDN);
        mpfr_mul(r25453, r25450, r25452, MPFR_RNDN);
        mpfr_sub(r25454, r25449, r25453, MPFR_RNDN);
        mpfr_sub(r25455, r25452, r25449, MPFR_RNDN);
        mpfr_div(r25456, r25454, r25455, MPFR_RNDN);
        mpfr_acos(r25457, r25456, MPFR_RNDN);
        mpfr_log(r25458, r25457, MPFR_RNDN);
        mpfr_pow(r25459, r25448, r25458, MPFR_RNDN);
        return mpfr_get_d(r25459, MPFR_RNDN);
}

static mpfr_t r25460, r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25460);
        mpfr_init_set_str(r25461, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r25462, "5", 10, MPFR_RNDN);
        mpfr_init(r25463);
        mpfr_init(r25464);
        mpfr_init(r25465);
        mpfr_init(r25466);
        mpfr_init(r25467);
        mpfr_init(r25468);
        mpfr_init(r25469);
        mpfr_init(r25470);
        mpfr_init(r25471);
}

double f_dm(double v) {
        mpfr_set_si(r25460, 1, MPFR_RNDN), mpfr_const_exp(r25460, r25460, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r25463, v, MPFR_RNDN);
        mpfr_mul(r25464, r25463, r25463, MPFR_RNDN);
        mpfr_mul(r25465, r25462, r25464, MPFR_RNDN);
        mpfr_sub(r25466, r25461, r25465, MPFR_RNDN);
        mpfr_sub(r25467, r25464, r25461, MPFR_RNDN);
        mpfr_div(r25468, r25466, r25467, MPFR_RNDN);
        mpfr_acos(r25469, r25468, MPFR_RNDN);
        mpfr_log(r25470, r25469, MPFR_RNDN);
        mpfr_pow(r25471, r25460, r25470, MPFR_RNDN);
        return mpfr_get_d(r25471, MPFR_RNDN);
}