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

char *name = "Bouland and Aaronson, Equation (26)";

double f_if(float a, float b) {
        float r27700 = a;
        float r27701 = r27700 * r27700;
        float r27702 = b;
        float r27703 = r27702 * r27702;
        float r27704 = r27701 + r27703;
        float r27705 = 2;
        float r27706 = pow(r27704, r27705);
        float r27707 = 4;
        float r27708 = r27707 * r27703;
        float r27709 = r27706 + r27708;
        float r27710 = 1;
        float r27711 = r27709 - r27710;
        return r27711;
}

double f_id(double a, double b) {
        double r27712 = a;
        double r27713 = r27712 * r27712;
        double r27714 = b;
        double r27715 = r27714 * r27714;
        double r27716 = r27713 + r27715;
        double r27717 = 2;
        double r27718 = pow(r27716, r27717);
        double r27719 = 4;
        double r27720 = r27719 * r27715;
        double r27721 = r27718 + r27720;
        double r27722 = 1;
        double r27723 = r27721 - r27722;
        return r27723;
}


double f_of(float a, float b) {
        float r27724 = b;
        float r27725 = a;
        float r27726 = r27725 * r27725;
        float r27727 = fma(r27724, r27724, r27726);
        float r27728 = r27727 * r27727;
        float r27729 = 4;
        float r27730 = r27724 * r27724;
        float r27731 = r27729 * r27730;
        float r27732 = r27728 + r27731;
        float r27733 = 1;
        float r27734 = r27732 - r27733;
        return r27734;
}

double f_od(double a, double b) {
        double r27735 = b;
        double r27736 = a;
        double r27737 = r27736 * r27736;
        double r27738 = fma(r27735, r27735, r27737);
        double r27739 = r27738 * r27738;
        double r27740 = 4;
        double r27741 = r27735 * r27735;
        double r27742 = r27740 * r27741;
        double r27743 = r27739 + r27742;
        double r27744 = 1;
        double r27745 = r27743 - r27744;
        return r27745;
}

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 r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27746);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init_set_str(r27751, "2", 10, MPFR_RNDN);
        mpfr_init(r27752);
        mpfr_init_set_str(r27753, "4", 10, MPFR_RNDN);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init_set_str(r27756, "1", 10, MPFR_RNDN);
        mpfr_init(r27757);
}

double f_im(double a, double b) {
        mpfr_set_d(r27746, a, MPFR_RNDN);
        mpfr_mul(r27747, r27746, r27746, MPFR_RNDN);
        mpfr_set_d(r27748, b, MPFR_RNDN);
        mpfr_mul(r27749, r27748, r27748, MPFR_RNDN);
        mpfr_add(r27750, r27747, r27749, MPFR_RNDN);
        ;
        mpfr_pow(r27752, r27750, r27751, MPFR_RNDN);
        ;
        mpfr_mul(r27754, r27753, r27749, MPFR_RNDN);
        mpfr_add(r27755, r27752, r27754, MPFR_RNDN);
        ;
        mpfr_sub(r27757, r27755, r27756, MPFR_RNDN);
        return mpfr_get_d(r27757, MPFR_RNDN);
}

static mpfr_t r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27758);
        mpfr_init(r27759);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
        mpfr_init_set_str(r27763, "4", 10, MPFR_RNDN);
        mpfr_init(r27764);
        mpfr_init(r27765);
        mpfr_init(r27766);
        mpfr_init_set_str(r27767, "1", 10, MPFR_RNDN);
        mpfr_init(r27768);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27758, b, MPFR_RNDN);
        mpfr_set_d(r27759, a, MPFR_RNDN);
        mpfr_mul(r27760, r27759, r27759, MPFR_RNDN);
        mpfr_fma(r27761, r27758, r27758, r27760, MPFR_RNDN);
        mpfr_mul(r27762, r27761, r27761, MPFR_RNDN);
        ;
        mpfr_mul(r27764, r27758, r27758, MPFR_RNDN);
        mpfr_mul(r27765, r27763, r27764, MPFR_RNDN);
        mpfr_add(r27766, r27762, r27765, MPFR_RNDN);
        ;
        mpfr_sub(r27768, r27766, r27767, MPFR_RNDN);
        return mpfr_get_d(r27768, MPFR_RNDN);
}

static mpfr_t r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27769);
        mpfr_init(r27770);
        mpfr_init(r27771);
        mpfr_init(r27772);
        mpfr_init(r27773);
        mpfr_init_set_str(r27774, "4", 10, MPFR_RNDN);
        mpfr_init(r27775);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init_set_str(r27778, "1", 10, MPFR_RNDN);
        mpfr_init(r27779);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27769, b, MPFR_RNDN);
        mpfr_set_d(r27770, a, MPFR_RNDN);
        mpfr_mul(r27771, r27770, r27770, MPFR_RNDN);
        mpfr_fma(r27772, r27769, r27769, r27771, MPFR_RNDN);
        mpfr_mul(r27773, r27772, r27772, MPFR_RNDN);
        ;
        mpfr_mul(r27775, r27769, r27769, MPFR_RNDN);
        mpfr_mul(r27776, r27774, r27775, MPFR_RNDN);
        mpfr_add(r27777, r27773, r27776, MPFR_RNDN);
        ;
        mpfr_sub(r27779, r27777, r27778, MPFR_RNDN);
        return mpfr_get_d(r27779, MPFR_RNDN);
}