#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 = 5;
        float r27215 = -r27214;
        float r27216 = v;
        float r27217 = r27216 * r27216;
        float r27218 = 1;
        float r27219 = fma(r27215, r27217, r27218);
        float r27220 = r27217 - r27218;
        float r27221 = r27219 / r27220;
        float r27222 = acos(r27221);
        float r27223 = 3;
        float r27224 = pow(r27222, r27223);
        float r27225 = cbrt(r27224);
        float r27226 = log1p(r27225);
        float r27227 = expm1(r27226);
        return r27227;
}

double f_od(double v) {
        double r27228 = 5;
        double r27229 = -r27228;
        double r27230 = v;
        double r27231 = r27230 * r27230;
        double r27232 = 1;
        double r27233 = fma(r27229, r27231, r27232);
        double r27234 = r27231 - r27232;
        double r27235 = r27233 / r27234;
        double r27236 = acos(r27235);
        double r27237 = 3;
        double r27238 = pow(r27236, r27237);
        double r27239 = cbrt(r27238);
        double r27240 = log1p(r27239);
        double r27241 = expm1(r27240);
        return r27241;
}

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

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

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

static mpfr_t r27251, r27252, r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261, r27262, r27263, r27264;

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

double f_fm(double v) {
        ;
        mpfr_neg(r27252, r27251, MPFR_RNDN);
        mpfr_set_d(r27253, v, MPFR_RNDN);
        mpfr_mul(r27254, r27253, r27253, MPFR_RNDN);
        ;
        mpfr_fma(r27256, r27252, r27254, r27255, MPFR_RNDN);
        mpfr_sub(r27257, r27254, r27255, MPFR_RNDN);
        mpfr_div(r27258, r27256, r27257, MPFR_RNDN);
        mpfr_acos(r27259, r27258, MPFR_RNDN);
        ;
        mpfr_pow(r27261, r27259, r27260, MPFR_RNDN);
        mpfr_cbrt(r27262, r27261, MPFR_RNDN);
        mpfr_log1p(r27263, r27262, MPFR_RNDN);
        mpfr_expm1(r27264, r27263, MPFR_RNDN);
        return mpfr_get_d(r27264, MPFR_RNDN);
}

static mpfr_t r27265, r27266, r27267, r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27265, "5", 10, MPFR_RNDN);
        mpfr_init(r27266);
        mpfr_init(r27267);
        mpfr_init(r27268);
        mpfr_init_set_str(r27269, "1", 10, MPFR_RNDN);
        mpfr_init(r27270);
        mpfr_init(r27271);
        mpfr_init(r27272);
        mpfr_init(r27273);
        mpfr_init_set_str(r27274, "3", 10, MPFR_RNDN);
        mpfr_init(r27275);
        mpfr_init(r27276);
        mpfr_init(r27277);
        mpfr_init(r27278);
}

double f_dm(double v) {
        ;
        mpfr_neg(r27266, r27265, MPFR_RNDN);
        mpfr_set_d(r27267, v, MPFR_RNDN);
        mpfr_mul(r27268, r27267, r27267, MPFR_RNDN);
        ;
        mpfr_fma(r27270, r27266, r27268, r27269, MPFR_RNDN);
        mpfr_sub(r27271, r27268, r27269, MPFR_RNDN);
        mpfr_div(r27272, r27270, r27271, MPFR_RNDN);
        mpfr_acos(r27273, r27272, MPFR_RNDN);
        ;
        mpfr_pow(r27275, r27273, r27274, MPFR_RNDN);
        mpfr_cbrt(r27276, r27275, MPFR_RNDN);
        mpfr_log1p(r27277, r27276, MPFR_RNDN);
        mpfr_expm1(r27278, r27277, MPFR_RNDN);
        return mpfr_get_d(r27278, MPFR_RNDN);
}