#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 r27280 = 1;
        float r27281 = 5;
        float r27282 = v;
        float r27283 = r27282 * r27282;
        float r27284 = r27281 * r27283;
        float r27285 = r27280 - r27284;
        float r27286 = r27283 - r27280;
        float r27287 = r27285 / r27286;
        float r27288 = acos(r27287);
        return r27288;
}

double f_id(double v) {
        double r27289 = 1;
        double r27290 = 5;
        double r27291 = v;
        double r27292 = r27291 * r27291;
        double r27293 = r27290 * r27292;
        double r27294 = r27289 - r27293;
        double r27295 = r27292 - r27289;
        double r27296 = r27294 / r27295;
        double r27297 = acos(r27296);
        return r27297;
}


double f_of(float v) {
        float r27298 = 1;
        float r27299 = 5;
        float r27300 = v;
        float r27301 = r27300 * r27300;
        float r27302 = r27299 * r27301;
        float r27303 = exp(r27302);
        float r27304 = log(r27303);
        float r27305 = r27298 - r27304;
        float r27306 = r27301 - r27298;
        float r27307 = r27305 / r27306;
        float r27308 = acos(r27307);
        return r27308;
}

double f_od(double v) {
        double r27309 = 1;
        double r27310 = 5;
        double r27311 = v;
        double r27312 = r27311 * r27311;
        double r27313 = r27310 * r27312;
        double r27314 = exp(r27313);
        double r27315 = log(r27314);
        double r27316 = r27309 - r27315;
        double r27317 = r27312 - r27309;
        double r27318 = r27316 / r27317;
        double r27319 = acos(r27318);
        return r27319;
}

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 r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27320, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27321, "5", 10, MPFR_RNDN);
        mpfr_init(r27322);
        mpfr_init(r27323);
        mpfr_init(r27324);
        mpfr_init(r27325);
        mpfr_init(r27326);
        mpfr_init(r27327);
        mpfr_init(r27328);
}

double f_im(double v) {
        ;
        ;
        mpfr_set_d(r27322, v, MPFR_RNDN);
        mpfr_mul(r27323, r27322, r27322, MPFR_RNDN);
        mpfr_mul(r27324, r27321, r27323, MPFR_RNDN);
        mpfr_sub(r27325, r27320, r27324, MPFR_RNDN);
        mpfr_sub(r27326, r27323, r27320, MPFR_RNDN);
        mpfr_div(r27327, r27325, r27326, MPFR_RNDN);
        mpfr_acos(r27328, r27327, MPFR_RNDN);
        return mpfr_get_d(r27328, MPFR_RNDN);
}

static mpfr_t r27329, r27330, r27331, r27332, r27333, r27334, r27335, r27336, r27337, r27338, r27339;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27329, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27330, "5", 10, MPFR_RNDN);
        mpfr_init(r27331);
        mpfr_init(r27332);
        mpfr_init(r27333);
        mpfr_init(r27334);
        mpfr_init(r27335);
        mpfr_init(r27336);
        mpfr_init(r27337);
        mpfr_init(r27338);
        mpfr_init(r27339);
}

double f_fm(double v) {
        ;
        ;
        mpfr_set_d(r27331, v, MPFR_RNDN);
        mpfr_mul(r27332, r27331, r27331, MPFR_RNDN);
        mpfr_mul(r27333, r27330, r27332, MPFR_RNDN);
        mpfr_exp(r27334, r27333, MPFR_RNDN);
        mpfr_log(r27335, r27334, MPFR_RNDN);
        mpfr_sub(r27336, r27329, r27335, MPFR_RNDN);
        mpfr_sub(r27337, r27332, r27329, MPFR_RNDN);
        mpfr_div(r27338, r27336, r27337, MPFR_RNDN);
        mpfr_acos(r27339, r27338, MPFR_RNDN);
        return mpfr_get_d(r27339, MPFR_RNDN);
}

static mpfr_t r27340, r27341, r27342, r27343, r27344, r27345, r27346, r27347, r27348, r27349, r27350;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27340, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27341, "5", 10, MPFR_RNDN);
        mpfr_init(r27342);
        mpfr_init(r27343);
        mpfr_init(r27344);
        mpfr_init(r27345);
        mpfr_init(r27346);
        mpfr_init(r27347);
        mpfr_init(r27348);
        mpfr_init(r27349);
        mpfr_init(r27350);
}

double f_dm(double v) {
        ;
        ;
        mpfr_set_d(r27342, v, MPFR_RNDN);
        mpfr_mul(r27343, r27342, r27342, MPFR_RNDN);
        mpfr_mul(r27344, r27341, r27343, MPFR_RNDN);
        mpfr_exp(r27345, r27344, MPFR_RNDN);
        mpfr_log(r27346, r27345, MPFR_RNDN);
        mpfr_sub(r27347, r27340, r27346, MPFR_RNDN);
        mpfr_sub(r27348, r27343, r27340, MPFR_RNDN);
        mpfr_div(r27349, r27347, r27348, MPFR_RNDN);
        mpfr_acos(r27350, r27349, MPFR_RNDN);
        return mpfr_get_d(r27350, MPFR_RNDN);
}