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

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

double f_if(float a, float b) {
        float r27671 = a;
        float r27672 = r27671 * r27671;
        float r27673 = b;
        float r27674 = r27673 * r27673;
        float r27675 = r27672 + r27674;
        float r27676 = 2;
        float r27677 = pow(r27675, r27676);
        float r27678 = 4;
        float r27679 = 1;
        float r27680 = r27679 - r27671;
        float r27681 = r27672 * r27680;
        float r27682 = 3;
        float r27683 = r27682 + r27671;
        float r27684 = r27674 * r27683;
        float r27685 = r27681 + r27684;
        float r27686 = r27678 * r27685;
        float r27687 = r27677 + r27686;
        float r27688 = r27687 - r27679;
        return r27688;
}

double f_id(double a, double b) {
        double r27689 = a;
        double r27690 = r27689 * r27689;
        double r27691 = b;
        double r27692 = r27691 * r27691;
        double r27693 = r27690 + r27692;
        double r27694 = 2;
        double r27695 = pow(r27693, r27694);
        double r27696 = 4;
        double r27697 = 1;
        double r27698 = r27697 - r27689;
        double r27699 = r27690 * r27698;
        double r27700 = 3;
        double r27701 = r27700 + r27689;
        double r27702 = r27692 * r27701;
        double r27703 = r27699 + r27702;
        double r27704 = r27696 * r27703;
        double r27705 = r27695 + r27704;
        double r27706 = r27705 - r27697;
        return r27706;
}


double f_of(float a, float b) {
        float r27707 = 4;
        float r27708 = a;
        float r27709 = 1;
        float r27710 = r27709 - r27708;
        float r27711 = r27708 * r27710;
        float r27712 = 3;
        float r27713 = r27708 + r27712;
        float r27714 = b;
        float r27715 = r27714 * r27714;
        float r27716 = r27713 * r27715;
        float r27717 = fma(r27708, r27711, r27716);
        float r27718 = r27708 * r27714;
        float r27719 = r27718 * r27718;
        float r27720 = 2;
        float r27721 = pow(r27714, r27707);
        float r27722 = pow(r27708, r27707);
        float r27723 = r27721 + r27722;
        float r27724 = fma(r27719, r27720, r27723);
        float r27725 = fma(r27707, r27717, r27724);
        float r27726 = r27725 - r27709;
        return r27726;
}

double f_od(double a, double b) {
        double r27727 = 4;
        double r27728 = a;
        double r27729 = 1;
        double r27730 = r27729 - r27728;
        double r27731 = r27728 * r27730;
        double r27732 = 3;
        double r27733 = r27728 + r27732;
        double r27734 = b;
        double r27735 = r27734 * r27734;
        double r27736 = r27733 * r27735;
        double r27737 = fma(r27728, r27731, r27736);
        double r27738 = r27728 * r27734;
        double r27739 = r27738 * r27738;
        double r27740 = 2;
        double r27741 = pow(r27734, r27727);
        double r27742 = pow(r27728, r27727);
        double r27743 = r27741 + r27742;
        double r27744 = fma(r27739, r27740, r27743);
        double r27745 = fma(r27727, r27737, r27744);
        double r27746 = r27745 - r27729;
        return r27746;
}

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 r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init_set_str(r27752, "2", 10, MPFR_RNDN);
        mpfr_init(r27753);
        mpfr_init_set_str(r27754, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27755, "1", 10, MPFR_RNDN);
        mpfr_init(r27756);
        mpfr_init(r27757);
        mpfr_init_set_str(r27758, "3", 10, MPFR_RNDN);
        mpfr_init(r27759);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
        mpfr_init(r27763);
        mpfr_init(r27764);
}

double f_im(double a, double b) {
        mpfr_set_d(r27747, a, MPFR_RNDN);
        mpfr_mul(r27748, r27747, r27747, MPFR_RNDN);
        mpfr_set_d(r27749, b, MPFR_RNDN);
        mpfr_mul(r27750, r27749, r27749, MPFR_RNDN);
        mpfr_add(r27751, r27748, r27750, MPFR_RNDN);
        ;
        mpfr_pow(r27753, r27751, r27752, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r27756, r27755, r27747, MPFR_RNDN);
        mpfr_mul(r27757, r27748, r27756, MPFR_RNDN);
        ;
        mpfr_add(r27759, r27758, r27747, MPFR_RNDN);
        mpfr_mul(r27760, r27750, r27759, MPFR_RNDN);
        mpfr_add(r27761, r27757, r27760, MPFR_RNDN);
        mpfr_mul(r27762, r27754, r27761, MPFR_RNDN);
        mpfr_add(r27763, r27753, r27762, MPFR_RNDN);
        mpfr_sub(r27764, r27763, r27755, MPFR_RNDN);
        return mpfr_get_d(r27764, MPFR_RNDN);
}

static mpfr_t r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27765, "4", 10, MPFR_RNDN);
        mpfr_init(r27766);
        mpfr_init_set_str(r27767, "1", 10, MPFR_RNDN);
        mpfr_init(r27768);
        mpfr_init(r27769);
        mpfr_init_set_str(r27770, "3", 10, MPFR_RNDN);
        mpfr_init(r27771);
        mpfr_init(r27772);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init(r27775);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init_set_str(r27778, "2", 10, MPFR_RNDN);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init(r27782);
        mpfr_init(r27783);
        mpfr_init(r27784);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r27766, a, MPFR_RNDN);
        ;
        mpfr_sub(r27768, r27767, r27766, MPFR_RNDN);
        mpfr_mul(r27769, r27766, r27768, MPFR_RNDN);
        ;
        mpfr_add(r27771, r27766, r27770, MPFR_RNDN);
        mpfr_set_d(r27772, b, MPFR_RNDN);
        mpfr_mul(r27773, r27772, r27772, MPFR_RNDN);
        mpfr_mul(r27774, r27771, r27773, MPFR_RNDN);
        mpfr_fma(r27775, r27766, r27769, r27774, MPFR_RNDN);
        mpfr_mul(r27776, r27766, r27772, MPFR_RNDN);
        mpfr_mul(r27777, r27776, r27776, MPFR_RNDN);
        ;
        mpfr_pow(r27779, r27772, r27765, MPFR_RNDN);
        mpfr_pow(r27780, r27766, r27765, MPFR_RNDN);
        mpfr_add(r27781, r27779, r27780, MPFR_RNDN);
        mpfr_fma(r27782, r27777, r27778, r27781, MPFR_RNDN);
        mpfr_fma(r27783, r27765, r27775, r27782, MPFR_RNDN);
        mpfr_sub(r27784, r27783, r27767, MPFR_RNDN);
        return mpfr_get_d(r27784, MPFR_RNDN);
}

static mpfr_t r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802, r27803, r27804;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27785, "4", 10, MPFR_RNDN);
        mpfr_init(r27786);
        mpfr_init_set_str(r27787, "1", 10, MPFR_RNDN);
        mpfr_init(r27788);
        mpfr_init(r27789);
        mpfr_init_set_str(r27790, "3", 10, MPFR_RNDN);
        mpfr_init(r27791);
        mpfr_init(r27792);
        mpfr_init(r27793);
        mpfr_init(r27794);
        mpfr_init(r27795);
        mpfr_init(r27796);
        mpfr_init(r27797);
        mpfr_init_set_str(r27798, "2", 10, MPFR_RNDN);
        mpfr_init(r27799);
        mpfr_init(r27800);
        mpfr_init(r27801);
        mpfr_init(r27802);
        mpfr_init(r27803);
        mpfr_init(r27804);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r27786, a, MPFR_RNDN);
        ;
        mpfr_sub(r27788, r27787, r27786, MPFR_RNDN);
        mpfr_mul(r27789, r27786, r27788, MPFR_RNDN);
        ;
        mpfr_add(r27791, r27786, r27790, MPFR_RNDN);
        mpfr_set_d(r27792, b, MPFR_RNDN);
        mpfr_mul(r27793, r27792, r27792, MPFR_RNDN);
        mpfr_mul(r27794, r27791, r27793, MPFR_RNDN);
        mpfr_fma(r27795, r27786, r27789, r27794, MPFR_RNDN);
        mpfr_mul(r27796, r27786, r27792, MPFR_RNDN);
        mpfr_mul(r27797, r27796, r27796, MPFR_RNDN);
        ;
        mpfr_pow(r27799, r27792, r27785, MPFR_RNDN);
        mpfr_pow(r27800, r27786, r27785, MPFR_RNDN);
        mpfr_add(r27801, r27799, r27800, MPFR_RNDN);
        mpfr_fma(r27802, r27797, r27798, r27801, MPFR_RNDN);
        mpfr_fma(r27803, r27785, r27795, r27802, MPFR_RNDN);
        mpfr_sub(r27804, r27803, r27787, MPFR_RNDN);
        return mpfr_get_d(r27804, MPFR_RNDN);
}