#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 r27656 = a;
        float r27657 = r27656 * r27656;
        float r27658 = b;
        float r27659 = r27658 * r27658;
        float r27660 = r27657 + r27659;
        float r27661 = 2;
        float r27662 = pow(r27660, r27661);
        float r27663 = 4;
        float r27664 = r27663 * r27659;
        float r27665 = r27662 + r27664;
        float r27666 = 1;
        float r27667 = r27665 - r27666;
        return r27667;
}

double f_id(double a, double b) {
        double r27668 = a;
        double r27669 = r27668 * r27668;
        double r27670 = b;
        double r27671 = r27670 * r27670;
        double r27672 = r27669 + r27671;
        double r27673 = 2;
        double r27674 = pow(r27672, r27673);
        double r27675 = 4;
        double r27676 = r27675 * r27671;
        double r27677 = r27674 + r27676;
        double r27678 = 1;
        double r27679 = r27677 - r27678;
        return r27679;
}


double f_of(float a, float b) {
        float r27680 = b;
        float r27681 = -118.38282236026332;
        bool r27682 = r27680 <= r27681;
        float r27683 = 2.5969064672601682e-17;
        bool r27684 = r27680 <= r27683;
        float r27685 = !r27684;
        bool r27686 = r27682 || r27685;
        float r27687 = 4;
        float r27688 = r27687 * r27680;
        float r27689 = pow(r27680, r27687);
        float r27690 = 1;
        float r27691 = r27689 - r27690;
        float r27692 = fma(r27688, r27680, r27691);
        float r27693 = 2;
        float r27694 = a;
        float r27695 = r27693 * r27694;
        float r27696 = fma(r27695, r27694, r27687);
        float r27697 = r27680 * r27680;
        float r27698 = pow(r27694, r27687);
        float r27699 = r27698 - r27690;
        float r27700 = fma(r27696, r27697, r27699);
        float r27701 = r27686 ? r27692 : r27700;
        return r27701;
}

double f_od(double a, double b) {
        double r27702 = b;
        double r27703 = -118.38282236026332;
        bool r27704 = r27702 <= r27703;
        double r27705 = 2.5969064672601682e-17;
        bool r27706 = r27702 <= r27705;
        double r27707 = !r27706;
        bool r27708 = r27704 || r27707;
        double r27709 = 4;
        double r27710 = r27709 * r27702;
        double r27711 = pow(r27702, r27709);
        double r27712 = 1;
        double r27713 = r27711 - r27712;
        double r27714 = fma(r27710, r27702, r27713);
        double r27715 = 2;
        double r27716 = a;
        double r27717 = r27715 * r27716;
        double r27718 = fma(r27717, r27716, r27709);
        double r27719 = r27702 * r27702;
        double r27720 = pow(r27716, r27709);
        double r27721 = r27720 - r27712;
        double r27722 = fma(r27718, r27719, r27721);
        double r27723 = r27708 ? r27714 : r27722;
        return r27723;
}

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 r27724, r27725, r27726, r27727, r27728, r27729, r27730, r27731, r27732, r27733, r27734, r27735;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27724);
        mpfr_init(r27725);
        mpfr_init(r27726);
        mpfr_init(r27727);
        mpfr_init(r27728);
        mpfr_init_set_str(r27729, "2", 10, MPFR_RNDN);
        mpfr_init(r27730);
        mpfr_init_set_str(r27731, "4", 10, MPFR_RNDN);
        mpfr_init(r27732);
        mpfr_init(r27733);
        mpfr_init_set_str(r27734, "1", 10, MPFR_RNDN);
        mpfr_init(r27735);
}

double f_im(double a, double b) {
        mpfr_set_d(r27724, a, MPFR_RNDN);
        mpfr_mul(r27725, r27724, r27724, MPFR_RNDN);
        mpfr_set_d(r27726, b, MPFR_RNDN);
        mpfr_mul(r27727, r27726, r27726, MPFR_RNDN);
        mpfr_add(r27728, r27725, r27727, MPFR_RNDN);
        ;
        mpfr_pow(r27730, r27728, r27729, MPFR_RNDN);
        ;
        mpfr_mul(r27732, r27731, r27727, MPFR_RNDN);
        mpfr_add(r27733, r27730, r27732, MPFR_RNDN);
        ;
        mpfr_sub(r27735, r27733, r27734, MPFR_RNDN);
        return mpfr_get_d(r27735, MPFR_RNDN);
}

static mpfr_t r27736, r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27736);
        mpfr_init_set_str(r27737, "-118.38282236026332", 10, MPFR_RNDN);
        mpfr_init(r27738);
        mpfr_init_set_str(r27739, "2.5969064672601682e-17", 10, MPFR_RNDN);
        mpfr_init(r27740);
        mpfr_init(r27741);
        mpfr_init(r27742);
        mpfr_init_set_str(r27743, "4", 10, MPFR_RNDN);
        mpfr_init(r27744);
        mpfr_init(r27745);
        mpfr_init_set_str(r27746, "1", 10, MPFR_RNDN);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init_set_str(r27749, "2", 10, MPFR_RNDN);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init(r27752);
        mpfr_init(r27753);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init(r27756);
        mpfr_init(r27757);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27736, b, MPFR_RNDN);
        ;
        mpfr_set_si(r27738, mpfr_cmp(r27736, r27737) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27740, mpfr_cmp(r27736, r27739) <= 0, MPFR_RNDN);
        mpfr_set_si(r27741, !mpfr_get_si(r27740, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27742, mpfr_get_si(r27738, MPFR_RNDN) || mpfr_get_si(r27741, MPFR_RNDN), MPFR_RNDN);
        ;
        mpfr_mul(r27744, r27743, r27736, MPFR_RNDN);
        mpfr_pow(r27745, r27736, r27743, MPFR_RNDN);
        ;
        mpfr_sub(r27747, r27745, r27746, MPFR_RNDN);
        mpfr_fma(r27748, r27744, r27736, r27747, MPFR_RNDN);
        ;
        mpfr_set_d(r27750, a, MPFR_RNDN);
        mpfr_mul(r27751, r27749, r27750, MPFR_RNDN);
        mpfr_fma(r27752, r27751, r27750, r27743, MPFR_RNDN);
        mpfr_mul(r27753, r27736, r27736, MPFR_RNDN);
        mpfr_pow(r27754, r27750, r27743, MPFR_RNDN);
        mpfr_sub(r27755, r27754, r27746, MPFR_RNDN);
        mpfr_fma(r27756, r27752, r27753, r27755, MPFR_RNDN);
        if (mpfr_get_si(r27742, MPFR_RNDN)) { mpfr_set(r27757, r27748, MPFR_RNDN); } else { mpfr_set(r27757, r27756, MPFR_RNDN); };
        return mpfr_get_d(r27757, MPFR_RNDN);
}

static mpfr_t r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27758);
        mpfr_init_set_str(r27759, "-118.38282236026332", 10, MPFR_RNDN);
        mpfr_init(r27760);
        mpfr_init_set_str(r27761, "2.5969064672601682e-17", 10, MPFR_RNDN);
        mpfr_init(r27762);
        mpfr_init(r27763);
        mpfr_init(r27764);
        mpfr_init_set_str(r27765, "4", 10, MPFR_RNDN);
        mpfr_init(r27766);
        mpfr_init(r27767);
        mpfr_init_set_str(r27768, "1", 10, MPFR_RNDN);
        mpfr_init(r27769);
        mpfr_init(r27770);
        mpfr_init_set_str(r27771, "2", 10, MPFR_RNDN);
        mpfr_init(r27772);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init(r27775);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init(r27778);
        mpfr_init(r27779);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27758, b, MPFR_RNDN);
        ;
        mpfr_set_si(r27760, mpfr_cmp(r27758, r27759) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27762, mpfr_cmp(r27758, r27761) <= 0, MPFR_RNDN);
        mpfr_set_si(r27763, !mpfr_get_si(r27762, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27764, mpfr_get_si(r27760, MPFR_RNDN) || mpfr_get_si(r27763, MPFR_RNDN), MPFR_RNDN);
        ;
        mpfr_mul(r27766, r27765, r27758, MPFR_RNDN);
        mpfr_pow(r27767, r27758, r27765, MPFR_RNDN);
        ;
        mpfr_sub(r27769, r27767, r27768, MPFR_RNDN);
        mpfr_fma(r27770, r27766, r27758, r27769, MPFR_RNDN);
        ;
        mpfr_set_d(r27772, a, MPFR_RNDN);
        mpfr_mul(r27773, r27771, r27772, MPFR_RNDN);
        mpfr_fma(r27774, r27773, r27772, r27765, MPFR_RNDN);
        mpfr_mul(r27775, r27758, r27758, MPFR_RNDN);
        mpfr_pow(r27776, r27772, r27765, MPFR_RNDN);
        mpfr_sub(r27777, r27776, r27768, MPFR_RNDN);
        mpfr_fma(r27778, r27774, r27775, r27777, MPFR_RNDN);
        if (mpfr_get_si(r27764, MPFR_RNDN)) { mpfr_set(r27779, r27770, MPFR_RNDN); } else { mpfr_set(r27779, r27778, MPFR_RNDN); };
        return mpfr_get_d(r27779, MPFR_RNDN);
}