#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 r27408 = 1;
        float r27409 = 5;
        float r27410 = v;
        float r27411 = r27410 * r27410;
        float r27412 = r27409 * r27411;
        float r27413 = r27408 - r27412;
        float r27414 = r27411 - r27408;
        float r27415 = r27413 / r27414;
        float r27416 = acos(r27415);
        return r27416;
}

double f_id(double v) {
        double r27417 = 1;
        double r27418 = 5;
        double r27419 = v;
        double r27420 = r27419 * r27419;
        double r27421 = r27418 * r27420;
        double r27422 = r27417 - r27421;
        double r27423 = r27420 - r27417;
        double r27424 = r27422 / r27423;
        double r27425 = acos(r27424);
        return r27425;
}


double f_of(float v) {
        float r27426 = 1;
        float r27427 = v;
        float r27428 = r27427 * r27427;
        float r27429 = 5;
        float r27430 = r27428 * r27429;
        float r27431 = r27426 - r27430;
        float r27432 = r27428 - r27426;
        float r27433 = r27431 / r27432;
        float r27434 = 3;
        float r27435 = pow(r27433, r27434);
        float r27436 = cbrt(r27435);
        float r27437 = acos(r27436);
        return r27437;
}

double f_od(double v) {
        double r27438 = 1;
        double r27439 = v;
        double r27440 = r27439 * r27439;
        double r27441 = 5;
        double r27442 = r27440 * r27441;
        double r27443 = r27438 - r27442;
        double r27444 = r27440 - r27438;
        double r27445 = r27443 / r27444;
        double r27446 = 3;
        double r27447 = pow(r27445, r27446);
        double r27448 = cbrt(r27447);
        double r27449 = acos(r27448);
        return r27449;
}

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 r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27450, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27451, "5", 10, MPFR_RNDN);
        mpfr_init(r27452);
        mpfr_init(r27453);
        mpfr_init(r27454);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init(r27457);
        mpfr_init(r27458);
}

double f_im(double v) {
        ;
        ;
        mpfr_set_d(r27452, v, MPFR_RNDN);
        mpfr_mul(r27453, r27452, r27452, MPFR_RNDN);
        mpfr_mul(r27454, r27451, r27453, MPFR_RNDN);
        mpfr_sub(r27455, r27450, r27454, MPFR_RNDN);
        mpfr_sub(r27456, r27453, r27450, MPFR_RNDN);
        mpfr_div(r27457, r27455, r27456, MPFR_RNDN);
        mpfr_acos(r27458, r27457, MPFR_RNDN);
        return mpfr_get_d(r27458, MPFR_RNDN);
}

static mpfr_t r27459, r27460, r27461, r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27459, "1", 10, MPFR_RNDN);
        mpfr_init(r27460);
        mpfr_init(r27461);
        mpfr_init_set_str(r27462, "5", 10, MPFR_RNDN);
        mpfr_init(r27463);
        mpfr_init(r27464);
        mpfr_init(r27465);
        mpfr_init(r27466);
        mpfr_init_set_str(r27467, "3", 10, MPFR_RNDN);
        mpfr_init(r27468);
        mpfr_init(r27469);
        mpfr_init(r27470);
}

double f_fm(double v) {
        ;
        mpfr_set_d(r27460, v, MPFR_RNDN);
        mpfr_mul(r27461, r27460, r27460, MPFR_RNDN);
        ;
        mpfr_mul(r27463, r27461, r27462, MPFR_RNDN);
        mpfr_sub(r27464, r27459, r27463, MPFR_RNDN);
        mpfr_sub(r27465, r27461, r27459, MPFR_RNDN);
        mpfr_div(r27466, r27464, r27465, MPFR_RNDN);
        ;
        mpfr_pow(r27468, r27466, r27467, MPFR_RNDN);
        mpfr_cbrt(r27469, r27468, MPFR_RNDN);
        mpfr_acos(r27470, r27469, MPFR_RNDN);
        return mpfr_get_d(r27470, MPFR_RNDN);
}

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

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

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