#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Falkner and Boettcher, Appendix B, 2";

double f_if(float v) {
        float r27367 = 2;
        float r27368 = sqrt(r27367);
        float r27369 = 4;
        float r27370 = r27368 / r27369;
        float r27371 = 1;
        float r27372 = 3;
        float r27373 = v;
        float r27374 = r27373 * r27373;
        float r27375 = r27372 * r27374;
        float r27376 = r27371 - r27375;
        float r27377 = sqrt(r27376);
        float r27378 = r27370 * r27377;
        float r27379 = r27371 - r27374;
        float r27380 = r27378 * r27379;
        return r27380;
}

double f_id(double v) {
        double r27381 = 2;
        double r27382 = sqrt(r27381);
        double r27383 = 4;
        double r27384 = r27382 / r27383;
        double r27385 = 1;
        double r27386 = 3;
        double r27387 = v;
        double r27388 = r27387 * r27387;
        double r27389 = r27386 * r27388;
        double r27390 = r27385 - r27389;
        double r27391 = sqrt(r27390);
        double r27392 = r27384 * r27391;
        double r27393 = r27385 - r27388;
        double r27394 = r27392 * r27393;
        return r27394;
}


double f_of(float v) {
        float r27395 = 2;
        float r27396 = sqrt(r27395);
        float r27397 = 4;
        float r27398 = r27396 / r27397;
        float r27399 = 1;
        float r27400 = 3;
        float r27401 = v;
        float r27402 = r27401 * r27401;
        float r27403 = r27400 * r27402;
        float r27404 = r27399 - r27403;
        float r27405 = sqrt(r27404);
        float r27406 = r27398 * r27405;
        float r27407 = r27399 - r27402;
        float r27408 = r27406 * r27407;
        return r27408;
}

double f_od(double v) {
        double r27409 = 2;
        double r27410 = sqrt(r27409);
        double r27411 = 4;
        double r27412 = r27410 / r27411;
        double r27413 = 1;
        double r27414 = 3;
        double r27415 = v;
        double r27416 = r27415 * r27415;
        double r27417 = r27414 * r27416;
        double r27418 = r27413 - r27417;
        double r27419 = sqrt(r27418);
        double r27420 = r27412 * r27419;
        double r27421 = r27413 - r27416;
        double r27422 = r27420 * r27421;
        return r27422;
}

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 r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430, r27431, r27432, r27433, r27434, r27435, r27436;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27423, "2", 10, MPFR_RNDN);
        mpfr_init(r27424);
        mpfr_init_set_str(r27425, "4", 10, MPFR_RNDN);
        mpfr_init(r27426);
        mpfr_init_set_str(r27427, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27428, "3", 10, MPFR_RNDN);
        mpfr_init(r27429);
        mpfr_init(r27430);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
        mpfr_init(r27434);
        mpfr_init(r27435);
        mpfr_init(r27436);
}

double f_im(double v) {
        ;
        mpfr_sqrt(r27424, r27423, MPFR_RNDN);
        ;
        mpfr_div(r27426, r27424, r27425, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r27429, v, MPFR_RNDN);
        mpfr_mul(r27430, r27429, r27429, MPFR_RNDN);
        mpfr_mul(r27431, r27428, r27430, MPFR_RNDN);
        mpfr_sub(r27432, r27427, r27431, MPFR_RNDN);
        mpfr_sqrt(r27433, r27432, MPFR_RNDN);
        mpfr_mul(r27434, r27426, r27433, MPFR_RNDN);
        mpfr_sub(r27435, r27427, r27430, MPFR_RNDN);
        mpfr_mul(r27436, r27434, r27435, MPFR_RNDN);
        return mpfr_get_d(r27436, MPFR_RNDN);
}

static mpfr_t r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27437, "2", 10, MPFR_RNDN);
        mpfr_init(r27438);
        mpfr_init_set_str(r27439, "4", 10, MPFR_RNDN);
        mpfr_init(r27440);
        mpfr_init_set_str(r27441, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27442, "3", 10, MPFR_RNDN);
        mpfr_init(r27443);
        mpfr_init(r27444);
        mpfr_init(r27445);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init(r27449);
        mpfr_init(r27450);
}

double f_fm(double v) {
        ;
        mpfr_sqrt(r27438, r27437, MPFR_RNDN);
        ;
        mpfr_div(r27440, r27438, r27439, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r27443, v, MPFR_RNDN);
        mpfr_mul(r27444, r27443, r27443, MPFR_RNDN);
        mpfr_mul(r27445, r27442, r27444, MPFR_RNDN);
        mpfr_sub(r27446, r27441, r27445, MPFR_RNDN);
        mpfr_sqrt(r27447, r27446, MPFR_RNDN);
        mpfr_mul(r27448, r27440, r27447, MPFR_RNDN);
        mpfr_sub(r27449, r27441, r27444, MPFR_RNDN);
        mpfr_mul(r27450, r27448, r27449, MPFR_RNDN);
        return mpfr_get_d(r27450, MPFR_RNDN);
}

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

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27451, "2", 10, MPFR_RNDN);
        mpfr_init(r27452);
        mpfr_init_set_str(r27453, "4", 10, MPFR_RNDN);
        mpfr_init(r27454);
        mpfr_init_set_str(r27455, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27456, "3", 10, MPFR_RNDN);
        mpfr_init(r27457);
        mpfr_init(r27458);
        mpfr_init(r27459);
        mpfr_init(r27460);
        mpfr_init(r27461);
        mpfr_init(r27462);
        mpfr_init(r27463);
        mpfr_init(r27464);
}

double f_dm(double v) {
        ;
        mpfr_sqrt(r27452, r27451, MPFR_RNDN);
        ;
        mpfr_div(r27454, r27452, r27453, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r27457, v, MPFR_RNDN);
        mpfr_mul(r27458, r27457, r27457, MPFR_RNDN);
        mpfr_mul(r27459, r27456, r27458, MPFR_RNDN);
        mpfr_sub(r27460, r27455, r27459, MPFR_RNDN);
        mpfr_sqrt(r27461, r27460, MPFR_RNDN);
        mpfr_mul(r27462, r27454, r27461, MPFR_RNDN);
        mpfr_sub(r27463, r27455, r27458, MPFR_RNDN);
        mpfr_mul(r27464, r27462, r27463, MPFR_RNDN);
        return mpfr_get_d(r27464, MPFR_RNDN);
}