#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 r27196 = 1;
        float r27197 = 5;
        float r27198 = v;
        float r27199 = r27198 * r27198;
        float r27200 = r27197 * r27199;
        float r27201 = r27196 - r27200;
        float r27202 = r27199 - r27196;
        float r27203 = r27201 / r27202;
        float r27204 = acos(r27203);
        return r27204;
}

double f_id(double v) {
        double r27205 = 1;
        double r27206 = 5;
        double r27207 = v;
        double r27208 = r27207 * r27207;
        double r27209 = r27206 * r27208;
        double r27210 = r27205 - r27209;
        double r27211 = r27208 - r27205;
        double r27212 = r27210 / r27211;
        double r27213 = acos(r27212);
        return r27213;
}


double f_of(float v) {
        float r27214 = 1;
        float r27215 = 5;
        float r27216 = v;
        float r27217 = r27216 * r27216;
        float r27218 = r27215 * r27217;
        float r27219 = r27214 - r27218;
        float r27220 = cbrt(r27219);
        float r27221 = r27220 * r27220;
        float r27222 = r27217 - r27214;
        float r27223 = cbrt(r27222);
        float r27224 = r27223 * r27223;
        float r27225 = r27221 / r27224;
        float r27226 = r27220 / r27223;
        float r27227 = r27225 * r27226;
        float r27228 = acos(r27227);
        return r27228;
}

double f_od(double v) {
        double r27229 = 1;
        double r27230 = 5;
        double r27231 = v;
        double r27232 = r27231 * r27231;
        double r27233 = r27230 * r27232;
        double r27234 = r27229 - r27233;
        double r27235 = cbrt(r27234);
        double r27236 = r27235 * r27235;
        double r27237 = r27232 - r27229;
        double r27238 = cbrt(r27237);
        double r27239 = r27238 * r27238;
        double r27240 = r27236 / r27239;
        double r27241 = r27235 / r27238;
        double r27242 = r27240 * r27241;
        double r27243 = acos(r27242);
        return r27243;
}

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 r27244, r27245, r27246, r27247, r27248, r27249, r27250, r27251, r27252;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27244, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27245, "5", 10, MPFR_RNDN);
        mpfr_init(r27246);
        mpfr_init(r27247);
        mpfr_init(r27248);
        mpfr_init(r27249);
        mpfr_init(r27250);
        mpfr_init(r27251);
        mpfr_init(r27252);
}

double f_im(double v) {
        ;
        ;
        mpfr_set_d(r27246, v, MPFR_RNDN);
        mpfr_mul(r27247, r27246, r27246, MPFR_RNDN);
        mpfr_mul(r27248, r27245, r27247, MPFR_RNDN);
        mpfr_sub(r27249, r27244, r27248, MPFR_RNDN);
        mpfr_sub(r27250, r27247, r27244, MPFR_RNDN);
        mpfr_div(r27251, r27249, r27250, MPFR_RNDN);
        mpfr_acos(r27252, r27251, MPFR_RNDN);
        return mpfr_get_d(r27252, MPFR_RNDN);
}

static mpfr_t r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261, r27262, r27263, r27264, r27265, r27266, r27267;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27253, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27254, "5", 10, MPFR_RNDN);
        mpfr_init(r27255);
        mpfr_init(r27256);
        mpfr_init(r27257);
        mpfr_init(r27258);
        mpfr_init(r27259);
        mpfr_init(r27260);
        mpfr_init(r27261);
        mpfr_init(r27262);
        mpfr_init(r27263);
        mpfr_init(r27264);
        mpfr_init(r27265);
        mpfr_init(r27266);
        mpfr_init(r27267);
}

double f_fm(double v) {
        ;
        ;
        mpfr_set_d(r27255, v, MPFR_RNDN);
        mpfr_mul(r27256, r27255, r27255, MPFR_RNDN);
        mpfr_mul(r27257, r27254, r27256, MPFR_RNDN);
        mpfr_sub(r27258, r27253, r27257, MPFR_RNDN);
        mpfr_cbrt(r27259, r27258, MPFR_RNDN);
        mpfr_mul(r27260, r27259, r27259, MPFR_RNDN);
        mpfr_sub(r27261, r27256, r27253, MPFR_RNDN);
        mpfr_cbrt(r27262, r27261, MPFR_RNDN);
        mpfr_mul(r27263, r27262, r27262, MPFR_RNDN);
        mpfr_div(r27264, r27260, r27263, MPFR_RNDN);
        mpfr_div(r27265, r27259, r27262, MPFR_RNDN);
        mpfr_mul(r27266, r27264, r27265, MPFR_RNDN);
        mpfr_acos(r27267, r27266, MPFR_RNDN);
        return mpfr_get_d(r27267, MPFR_RNDN);
}

static mpfr_t r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279, r27280, r27281, r27282;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27268, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27269, "5", 10, MPFR_RNDN);
        mpfr_init(r27270);
        mpfr_init(r27271);
        mpfr_init(r27272);
        mpfr_init(r27273);
        mpfr_init(r27274);
        mpfr_init(r27275);
        mpfr_init(r27276);
        mpfr_init(r27277);
        mpfr_init(r27278);
        mpfr_init(r27279);
        mpfr_init(r27280);
        mpfr_init(r27281);
        mpfr_init(r27282);
}

double f_dm(double v) {
        ;
        ;
        mpfr_set_d(r27270, v, MPFR_RNDN);
        mpfr_mul(r27271, r27270, r27270, MPFR_RNDN);
        mpfr_mul(r27272, r27269, r27271, MPFR_RNDN);
        mpfr_sub(r27273, r27268, r27272, MPFR_RNDN);
        mpfr_cbrt(r27274, r27273, MPFR_RNDN);
        mpfr_mul(r27275, r27274, r27274, MPFR_RNDN);
        mpfr_sub(r27276, r27271, r27268, MPFR_RNDN);
        mpfr_cbrt(r27277, r27276, MPFR_RNDN);
        mpfr_mul(r27278, r27277, r27277, MPFR_RNDN);
        mpfr_div(r27279, r27275, r27278, MPFR_RNDN);
        mpfr_div(r27280, r27274, r27277, MPFR_RNDN);
        mpfr_mul(r27281, r27279, r27280, MPFR_RNDN);
        mpfr_acos(r27282, r27281, MPFR_RNDN);
        return mpfr_get_d(r27282, MPFR_RNDN);
}