#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 r17641 = a;
        float r17642 = r17641 * r17641;
        float r17643 = b;
        float r17644 = r17643 * r17643;
        float r17645 = r17642 + r17644;
        float r17646 = r17645 * r17645;
        float r17647 = 4.0f;
        float r17648 = 1.0f;
        float r17649 = r17648 + r17641;
        float r17650 = r17642 * r17649;
        float r17651 = 3.0f;
        float r17652 = r17651 * r17641;
        float r17653 = r17648 - r17652;
        float r17654 = r17644 * r17653;
        float r17655 = r17650 + r17654;
        float r17656 = r17647 * r17655;
        float r17657 = r17646 + r17656;
        float r17658 = r17657 - r17648;
        return r17658;
}

double f_id(double a, double b) {
        double r17659 = a;
        double r17660 = r17659 * r17659;
        double r17661 = b;
        double r17662 = r17661 * r17661;
        double r17663 = r17660 + r17662;
        double r17664 = r17663 * r17663;
        double r17665 = 4.0;
        double r17666 = 1.0;
        double r17667 = r17666 + r17659;
        double r17668 = r17660 * r17667;
        double r17669 = 3.0;
        double r17670 = r17669 * r17659;
        double r17671 = r17666 - r17670;
        double r17672 = r17662 * r17671;
        double r17673 = r17668 + r17672;
        double r17674 = r17665 * r17673;
        double r17675 = r17664 + r17674;
        double r17676 = r17675 - r17666;
        return r17676;
}


double f_of(float a, float b) {
        float r17677 = a;
        float r17678 = r17677 * r17677;
        float r17679 = b;
        float r17680 = r17679 * r17679;
        float r17681 = r17678 + r17680;
        float r17682 = sqrt(r17681);
        float r17683 = r17682 * r17682;
        float r17684 = r17683 * r17683;
        float r17685 = 4.0f;
        float r17686 = 1.0f;
        float r17687 = r17686 + r17677;
        float r17688 = r17678 * r17687;
        float r17689 = 3.0f;
        float r17690 = r17689 * r17677;
        float r17691 = r17686 - r17690;
        float r17692 = r17680 * r17691;
        float r17693 = r17688 + r17692;
        float r17694 = r17685 * r17693;
        float r17695 = r17684 + r17694;
        float r17696 = r17695 - r17686;
        return r17696;
}

double f_od(double a, double b) {
        double r17697 = a;
        double r17698 = r17697 * r17697;
        double r17699 = b;
        double r17700 = r17699 * r17699;
        double r17701 = r17698 + r17700;
        double r17702 = sqrt(r17701);
        double r17703 = r17702 * r17702;
        double r17704 = r17703 * r17703;
        double r17705 = 4.0;
        double r17706 = 1.0;
        double r17707 = r17706 + r17697;
        double r17708 = r17698 * r17707;
        double r17709 = 3.0;
        double r17710 = r17709 * r17697;
        double r17711 = r17706 - r17710;
        double r17712 = r17700 * r17711;
        double r17713 = r17708 + r17712;
        double r17714 = r17705 * r17713;
        double r17715 = r17704 + r17714;
        double r17716 = r17715 - r17706;
        return r17716;
}

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 r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732, r17733, r17734;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17717);
        mpfr_init(r17718);
        mpfr_init(r17719);
        mpfr_init(r17720);
        mpfr_init(r17721);
        mpfr_init(r17722);
        mpfr_init_set_str(r17723, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r17724, "1", 10, MPFR_RNDN);
        mpfr_init(r17725);
        mpfr_init(r17726);
        mpfr_init_set_str(r17727, "3", 10, MPFR_RNDN);
        mpfr_init(r17728);
        mpfr_init(r17729);
        mpfr_init(r17730);
        mpfr_init(r17731);
        mpfr_init(r17732);
        mpfr_init(r17733);
        mpfr_init(r17734);
}

double f_im(double a, double b) {
        mpfr_set_d(r17717, a, MPFR_RNDN);
        mpfr_sqr(r17718, r17717, MPFR_RNDN);
        mpfr_set_d(r17719, b, MPFR_RNDN);
        mpfr_sqr(r17720, r17719, MPFR_RNDN);
        mpfr_add(r17721, r17718, r17720, MPFR_RNDN);
        mpfr_sqr(r17722, r17721, MPFR_RNDN);
        ;
        ;
        mpfr_add(r17725, r17724, r17717, MPFR_RNDN);
        mpfr_mul(r17726, r17718, r17725, MPFR_RNDN);
        ;
        mpfr_mul(r17728, r17727, r17717, MPFR_RNDN);
        mpfr_sub(r17729, r17724, r17728, MPFR_RNDN);
        mpfr_mul(r17730, r17720, r17729, MPFR_RNDN);
        mpfr_add(r17731, r17726, r17730, MPFR_RNDN);
        mpfr_mul(r17732, r17723, r17731, MPFR_RNDN);
        mpfr_add(r17733, r17722, r17732, MPFR_RNDN);
        mpfr_sub(r17734, r17733, r17724, MPFR_RNDN);
        return mpfr_get_d(r17734, MPFR_RNDN);
}

static mpfr_t r17735, r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743, r17744, r17745, r17746, r17747, r17748, r17749, r17750, r17751, r17752, r17753, r17754;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17735);
        mpfr_init(r17736);
        mpfr_init(r17737);
        mpfr_init(r17738);
        mpfr_init(r17739);
        mpfr_init(r17740);
        mpfr_init(r17741);
        mpfr_init(r17742);
        mpfr_init_set_str(r17743, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r17744, "1", 10, MPFR_RNDN);
        mpfr_init(r17745);
        mpfr_init(r17746);
        mpfr_init_set_str(r17747, "3", 10, MPFR_RNDN);
        mpfr_init(r17748);
        mpfr_init(r17749);
        mpfr_init(r17750);
        mpfr_init(r17751);
        mpfr_init(r17752);
        mpfr_init(r17753);
        mpfr_init(r17754);
}

double f_fm(double a, double b) {
        mpfr_set_d(r17735, a, MPFR_RNDN);
        mpfr_sqr(r17736, r17735, MPFR_RNDN);
        mpfr_set_d(r17737, b, MPFR_RNDN);
        mpfr_sqr(r17738, r17737, MPFR_RNDN);
        mpfr_add(r17739, r17736, r17738, MPFR_RNDN);
        mpfr_sqrt(r17740, r17739, MPFR_RNDN);
        mpfr_mul(r17741, r17740, r17740, MPFR_RNDN);
        mpfr_sqr(r17742, r17741, MPFR_RNDN);
        ;
        ;
        mpfr_add(r17745, r17744, r17735, MPFR_RNDN);
        mpfr_mul(r17746, r17736, r17745, MPFR_RNDN);
        ;
        mpfr_mul(r17748, r17747, r17735, MPFR_RNDN);
        mpfr_sub(r17749, r17744, r17748, MPFR_RNDN);
        mpfr_mul(r17750, r17738, r17749, MPFR_RNDN);
        mpfr_add(r17751, r17746, r17750, MPFR_RNDN);
        mpfr_mul(r17752, r17743, r17751, MPFR_RNDN);
        mpfr_add(r17753, r17742, r17752, MPFR_RNDN);
        mpfr_sub(r17754, r17753, r17744, MPFR_RNDN);
        return mpfr_get_d(r17754, MPFR_RNDN);
}

static mpfr_t r17755, r17756, r17757, r17758, r17759, r17760, r17761, r17762, r17763, r17764, r17765, r17766, r17767, r17768, r17769, r17770, r17771, r17772, r17773, r17774;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17755);
        mpfr_init(r17756);
        mpfr_init(r17757);
        mpfr_init(r17758);
        mpfr_init(r17759);
        mpfr_init(r17760);
        mpfr_init(r17761);
        mpfr_init(r17762);
        mpfr_init_set_str(r17763, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r17764, "1", 10, MPFR_RNDN);
        mpfr_init(r17765);
        mpfr_init(r17766);
        mpfr_init_set_str(r17767, "3", 10, MPFR_RNDN);
        mpfr_init(r17768);
        mpfr_init(r17769);
        mpfr_init(r17770);
        mpfr_init(r17771);
        mpfr_init(r17772);
        mpfr_init(r17773);
        mpfr_init(r17774);
}

double f_dm(double a, double b) {
        mpfr_set_d(r17755, a, MPFR_RNDN);
        mpfr_sqr(r17756, r17755, MPFR_RNDN);
        mpfr_set_d(r17757, b, MPFR_RNDN);
        mpfr_sqr(r17758, r17757, MPFR_RNDN);
        mpfr_add(r17759, r17756, r17758, MPFR_RNDN);
        mpfr_sqrt(r17760, r17759, MPFR_RNDN);
        mpfr_mul(r17761, r17760, r17760, MPFR_RNDN);
        mpfr_sqr(r17762, r17761, MPFR_RNDN);
        ;
        ;
        mpfr_add(r17765, r17764, r17755, MPFR_RNDN);
        mpfr_mul(r17766, r17756, r17765, MPFR_RNDN);
        ;
        mpfr_mul(r17768, r17767, r17755, MPFR_RNDN);
        mpfr_sub(r17769, r17764, r17768, MPFR_RNDN);
        mpfr_mul(r17770, r17758, r17769, MPFR_RNDN);
        mpfr_add(r17771, r17766, r17770, MPFR_RNDN);
        mpfr_mul(r17772, r17763, r17771, MPFR_RNDN);
        mpfr_add(r17773, r17762, r17772, MPFR_RNDN);
        mpfr_sub(r17774, r17773, r17764, MPFR_RNDN);
        return mpfr_get_d(r17774, MPFR_RNDN);
}