#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 r27393 = 1;
        float r27394 = 5;
        float r27395 = v;
        float r27396 = r27395 * r27395;
        float r27397 = r27394 * r27396;
        float r27398 = r27393 - r27397;
        float r27399 = r27396 - r27393;
        float r27400 = r27398 / r27399;
        float r27401 = acos(r27400);
        return r27401;
}

double f_id(double v) {
        double r27402 = 1;
        double r27403 = 5;
        double r27404 = v;
        double r27405 = r27404 * r27404;
        double r27406 = r27403 * r27405;
        double r27407 = r27402 - r27406;
        double r27408 = r27405 - r27402;
        double r27409 = r27407 / r27408;
        double r27410 = acos(r27409);
        return r27410;
}


double f_of(float v) {
        float r27411 = 1;
        float r27412 = 5;
        float r27413 = v;
        float r27414 = r27413 * r27413;
        float r27415 = r27412 * r27414;
        float r27416 = r27411 - r27415;
        float r27417 = r27414 - r27411;
        float r27418 = r27416 / r27417;
        float r27419 = acos(r27418);
        float r27420 = cbrt(r27419);
        float r27421 = cbrt(r27420);
        float r27422 = r27421 * r27421;
        float r27423 = r27420 * r27422;
        float r27424 = r27423 * r27421;
        float r27425 = r27422 * r27421;
        float r27426 = r27424 * r27425;
        return r27426;
}

double f_od(double v) {
        double r27427 = 1;
        double r27428 = 5;
        double r27429 = v;
        double r27430 = r27429 * r27429;
        double r27431 = r27428 * r27430;
        double r27432 = r27427 - r27431;
        double r27433 = r27430 - r27427;
        double r27434 = r27432 / r27433;
        double r27435 = acos(r27434);
        double r27436 = cbrt(r27435);
        double r27437 = cbrt(r27436);
        double r27438 = r27437 * r27437;
        double r27439 = r27436 * r27438;
        double r27440 = r27439 * r27437;
        double r27441 = r27438 * r27437;
        double r27442 = r27440 * r27441;
        return r27442;
}

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 r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27443, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27444, "5", 10, MPFR_RNDN);
        mpfr_init(r27445);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init(r27449);
        mpfr_init(r27450);
        mpfr_init(r27451);
}

double f_im(double v) {
        ;
        ;
        mpfr_set_d(r27445, v, MPFR_RNDN);
        mpfr_mul(r27446, r27445, r27445, MPFR_RNDN);
        mpfr_mul(r27447, r27444, r27446, MPFR_RNDN);
        mpfr_sub(r27448, r27443, r27447, MPFR_RNDN);
        mpfr_sub(r27449, r27446, r27443, MPFR_RNDN);
        mpfr_div(r27450, r27448, r27449, MPFR_RNDN);
        mpfr_acos(r27451, r27450, MPFR_RNDN);
        return mpfr_get_d(r27451, MPFR_RNDN);
}

static mpfr_t r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461, r27462, r27463, r27464, r27465, r27466, r27467;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27452, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27453, "5", 10, MPFR_RNDN);
        mpfr_init(r27454);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init(r27457);
        mpfr_init(r27458);
        mpfr_init(r27459);
        mpfr_init(r27460);
        mpfr_init(r27461);
        mpfr_init(r27462);
        mpfr_init(r27463);
        mpfr_init(r27464);
        mpfr_init(r27465);
        mpfr_init(r27466);
        mpfr_init(r27467);
}

double f_fm(double v) {
        ;
        ;
        mpfr_set_d(r27454, v, MPFR_RNDN);
        mpfr_mul(r27455, r27454, r27454, MPFR_RNDN);
        mpfr_mul(r27456, r27453, r27455, MPFR_RNDN);
        mpfr_sub(r27457, r27452, r27456, MPFR_RNDN);
        mpfr_sub(r27458, r27455, r27452, MPFR_RNDN);
        mpfr_div(r27459, r27457, r27458, MPFR_RNDN);
        mpfr_acos(r27460, r27459, MPFR_RNDN);
        mpfr_cbrt(r27461, r27460, MPFR_RNDN);
        mpfr_cbrt(r27462, r27461, MPFR_RNDN);
        mpfr_mul(r27463, r27462, r27462, MPFR_RNDN);
        mpfr_mul(r27464, r27461, r27463, MPFR_RNDN);
        mpfr_mul(r27465, r27464, r27462, MPFR_RNDN);
        mpfr_mul(r27466, r27463, r27462, MPFR_RNDN);
        mpfr_mul(r27467, r27465, r27466, MPFR_RNDN);
        return mpfr_get_d(r27467, MPFR_RNDN);
}

static mpfr_t r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27468, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27469, "5", 10, MPFR_RNDN);
        mpfr_init(r27470);
        mpfr_init(r27471);
        mpfr_init(r27472);
        mpfr_init(r27473);
        mpfr_init(r27474);
        mpfr_init(r27475);
        mpfr_init(r27476);
        mpfr_init(r27477);
        mpfr_init(r27478);
        mpfr_init(r27479);
        mpfr_init(r27480);
        mpfr_init(r27481);
        mpfr_init(r27482);
        mpfr_init(r27483);
}

double f_dm(double v) {
        ;
        ;
        mpfr_set_d(r27470, v, MPFR_RNDN);
        mpfr_mul(r27471, r27470, r27470, MPFR_RNDN);
        mpfr_mul(r27472, r27469, r27471, MPFR_RNDN);
        mpfr_sub(r27473, r27468, r27472, MPFR_RNDN);
        mpfr_sub(r27474, r27471, r27468, MPFR_RNDN);
        mpfr_div(r27475, r27473, r27474, MPFR_RNDN);
        mpfr_acos(r27476, r27475, MPFR_RNDN);
        mpfr_cbrt(r27477, r27476, MPFR_RNDN);
        mpfr_cbrt(r27478, r27477, MPFR_RNDN);
        mpfr_mul(r27479, r27478, r27478, MPFR_RNDN);
        mpfr_mul(r27480, r27477, r27479, MPFR_RNDN);
        mpfr_mul(r27481, r27480, r27478, MPFR_RNDN);
        mpfr_mul(r27482, r27479, r27478, MPFR_RNDN);
        mpfr_mul(r27483, r27481, r27482, MPFR_RNDN);
        return mpfr_get_d(r27483, MPFR_RNDN);
}