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

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

double f_if(float a, float b) {
        float r27654 = a;
        float r27655 = r27654 * r27654;
        float r27656 = b;
        float r27657 = r27656 * r27656;
        float r27658 = r27655 + r27657;
        float r27659 = 2;
        float r27660 = pow(r27658, r27659);
        float r27661 = 4;
        float r27662 = 1;
        float r27663 = r27662 + r27654;
        float r27664 = r27655 * r27663;
        float r27665 = 3;
        float r27666 = r27665 * r27654;
        float r27667 = r27662 - r27666;
        float r27668 = r27657 * r27667;
        float r27669 = r27664 + r27668;
        float r27670 = r27661 * r27669;
        float r27671 = r27660 + r27670;
        float r27672 = r27671 - r27662;
        return r27672;
}

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


double f_of(float a, float b) {
        float r27692 = b;
        float r27693 = a;
        float r27694 = hypot(r27692, r27693);
        float r27695 = 3;
        float r27696 = 1;
        float r27697 = r27695 + r27696;
        float r27698 = pow(r27694, r27697);
        float r27699 = 4;
        float r27700 = r27693 * r27693;
        float r27701 = r27696 + r27693;
        float r27702 = r27700 * r27701;
        float r27703 = r27692 * r27692;
        float r27704 = r27695 * r27693;
        float r27705 = r27696 - r27704;
        float r27706 = r27703 * r27705;
        float r27707 = r27702 + r27706;
        float r27708 = r27699 * r27707;
        float r27709 = r27698 + r27708;
        float r27710 = r27709 - r27696;
        return r27710;
}

double f_od(double a, double b) {
        double r27711 = b;
        double r27712 = a;
        double r27713 = hypot(r27711, r27712);
        double r27714 = 3;
        double r27715 = 1;
        double r27716 = r27714 + r27715;
        double r27717 = pow(r27713, r27716);
        double r27718 = 4;
        double r27719 = r27712 * r27712;
        double r27720 = r27715 + r27712;
        double r27721 = r27719 * r27720;
        double r27722 = r27711 * r27711;
        double r27723 = r27714 * r27712;
        double r27724 = r27715 - r27723;
        double r27725 = r27722 * r27724;
        double r27726 = r27721 + r27725;
        double r27727 = r27718 * r27726;
        double r27728 = r27717 + r27727;
        double r27729 = r27728 - r27715;
        return r27729;
}

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 r27730, r27731, r27732, r27733, r27734, r27735, r27736, r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27730);
        mpfr_init(r27731);
        mpfr_init(r27732);
        mpfr_init(r27733);
        mpfr_init(r27734);
        mpfr_init_set_str(r27735, "2", 10, MPFR_RNDN);
        mpfr_init(r27736);
        mpfr_init_set_str(r27737, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27738, "1", 10, MPFR_RNDN);
        mpfr_init(r27739);
        mpfr_init(r27740);
        mpfr_init_set_str(r27741, "3", 10, MPFR_RNDN);
        mpfr_init(r27742);
        mpfr_init(r27743);
        mpfr_init(r27744);
        mpfr_init(r27745);
        mpfr_init(r27746);
        mpfr_init(r27747);
        mpfr_init(r27748);
}

double f_im(double a, double b) {
        mpfr_set_d(r27730, a, MPFR_RNDN);
        mpfr_mul(r27731, r27730, r27730, MPFR_RNDN);
        mpfr_set_d(r27732, b, MPFR_RNDN);
        mpfr_mul(r27733, r27732, r27732, MPFR_RNDN);
        mpfr_add(r27734, r27731, r27733, MPFR_RNDN);
        ;
        mpfr_pow(r27736, r27734, r27735, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27739, r27738, r27730, MPFR_RNDN);
        mpfr_mul(r27740, r27731, r27739, MPFR_RNDN);
        ;
        mpfr_mul(r27742, r27741, r27730, MPFR_RNDN);
        mpfr_sub(r27743, r27738, r27742, MPFR_RNDN);
        mpfr_mul(r27744, r27733, r27743, MPFR_RNDN);
        mpfr_add(r27745, r27740, r27744, MPFR_RNDN);
        mpfr_mul(r27746, r27737, r27745, MPFR_RNDN);
        mpfr_add(r27747, r27736, r27746, MPFR_RNDN);
        mpfr_sub(r27748, r27747, r27738, MPFR_RNDN);
        return mpfr_get_d(r27748, MPFR_RNDN);
}

static mpfr_t r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init_set_str(r27752, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r27753, "1", 10, MPFR_RNDN);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init_set_str(r27756, "4", 10, MPFR_RNDN);
        mpfr_init(r27757);
        mpfr_init(r27758);
        mpfr_init(r27759);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
        mpfr_init(r27763);
        mpfr_init(r27764);
        mpfr_init(r27765);
        mpfr_init(r27766);
        mpfr_init(r27767);
}

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

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

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27768);
        mpfr_init(r27769);
        mpfr_init(r27770);
        mpfr_init_set_str(r27771, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r27772, "1", 10, MPFR_RNDN);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init_set_str(r27775, "4", 10, MPFR_RNDN);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init(r27778);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init(r27782);
        mpfr_init(r27783);
        mpfr_init(r27784);
        mpfr_init(r27785);
        mpfr_init(r27786);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27768, b, MPFR_RNDN);
        mpfr_set_d(r27769, a, MPFR_RNDN);
        mpfr_hypot(r27770, r27768, r27769, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27773, r27771, r27772, MPFR_RNDN);
        mpfr_pow(r27774, r27770, r27773, MPFR_RNDN);
        ;
        mpfr_mul(r27776, r27769, r27769, MPFR_RNDN);
        mpfr_add(r27777, r27772, r27769, MPFR_RNDN);
        mpfr_mul(r27778, r27776, r27777, MPFR_RNDN);
        mpfr_mul(r27779, r27768, r27768, MPFR_RNDN);
        mpfr_mul(r27780, r27771, r27769, MPFR_RNDN);
        mpfr_sub(r27781, r27772, r27780, MPFR_RNDN);
        mpfr_mul(r27782, r27779, r27781, MPFR_RNDN);
        mpfr_add(r27783, r27778, r27782, MPFR_RNDN);
        mpfr_mul(r27784, r27775, r27783, MPFR_RNDN);
        mpfr_add(r27785, r27774, r27784, MPFR_RNDN);
        mpfr_sub(r27786, r27785, r27772, MPFR_RNDN);
        return mpfr_get_d(r27786, MPFR_RNDN);
}