#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 r27606 = a;
        float r27607 = r27606 * r27606;
        float r27608 = b;
        float r27609 = r27608 * r27608;
        float r27610 = r27607 + r27609;
        float r27611 = 2;
        float r27612 = pow(r27610, r27611);
        float r27613 = 4;
        float r27614 = 1;
        float r27615 = r27614 - r27606;
        float r27616 = r27607 * r27615;
        float r27617 = 3;
        float r27618 = r27617 + r27606;
        float r27619 = r27609 * r27618;
        float r27620 = r27616 + r27619;
        float r27621 = r27613 * r27620;
        float r27622 = r27612 + r27621;
        float r27623 = r27622 - r27614;
        return r27623;
}

double f_id(double a, double b) {
        double r27624 = a;
        double r27625 = r27624 * r27624;
        double r27626 = b;
        double r27627 = r27626 * r27626;
        double r27628 = r27625 + r27627;
        double r27629 = 2;
        double r27630 = pow(r27628, r27629);
        double r27631 = 4;
        double r27632 = 1;
        double r27633 = r27632 - r27624;
        double r27634 = r27625 * r27633;
        double r27635 = 3;
        double r27636 = r27635 + r27624;
        double r27637 = r27627 * r27636;
        double r27638 = r27634 + r27637;
        double r27639 = r27631 * r27638;
        double r27640 = r27630 + r27639;
        double r27641 = r27640 - r27632;
        return r27641;
}


double f_of(float a, float b) {
        float r27642 = a;
        float r27643 = r27642 * r27642;
        float r27644 = b;
        float r27645 = r27644 * r27644;
        float r27646 = r27643 + r27645;
        float r27647 = 2;
        float r27648 = pow(r27646, r27647);
        float r27649 = 4;
        float r27650 = 1;
        float r27651 = r27650 - r27642;
        float r27652 = r27643 * r27651;
        float r27653 = 3;
        float r27654 = r27653 + r27642;
        float r27655 = r27645 * r27654;
        float r27656 = r27652 + r27655;
        float r27657 = r27649 * r27656;
        float r27658 = r27648 + r27657;
        float r27659 = r27658 - r27650;
        return r27659;
}

double f_od(double a, double b) {
        double r27660 = a;
        double r27661 = r27660 * r27660;
        double r27662 = b;
        double r27663 = r27662 * r27662;
        double r27664 = r27661 + r27663;
        double r27665 = 2;
        double r27666 = pow(r27664, r27665);
        double r27667 = 4;
        double r27668 = 1;
        double r27669 = r27668 - r27660;
        double r27670 = r27661 * r27669;
        double r27671 = 3;
        double r27672 = r27671 + r27660;
        double r27673 = r27663 * r27672;
        double r27674 = r27670 + r27673;
        double r27675 = r27667 * r27674;
        double r27676 = r27666 + r27675;
        double r27677 = r27676 - r27668;
        return r27677;
}

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 r27678, r27679, r27680, r27681, r27682, r27683, r27684, r27685, r27686, r27687, r27688, r27689, r27690, r27691, r27692, r27693, r27694, r27695;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27678);
        mpfr_init(r27679);
        mpfr_init(r27680);
        mpfr_init(r27681);
        mpfr_init(r27682);
        mpfr_init_set_str(r27683, "2", 10, MPFR_RNDN);
        mpfr_init(r27684);
        mpfr_init_set_str(r27685, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27686, "1", 10, MPFR_RNDN);
        mpfr_init(r27687);
        mpfr_init(r27688);
        mpfr_init_set_str(r27689, "3", 10, MPFR_RNDN);
        mpfr_init(r27690);
        mpfr_init(r27691);
        mpfr_init(r27692);
        mpfr_init(r27693);
        mpfr_init(r27694);
        mpfr_init(r27695);
}

double f_im(double a, double b) {
        mpfr_set_d(r27678, a, MPFR_RNDN);
        mpfr_mul(r27679, r27678, r27678, MPFR_RNDN);
        mpfr_set_d(r27680, b, MPFR_RNDN);
        mpfr_mul(r27681, r27680, r27680, MPFR_RNDN);
        mpfr_add(r27682, r27679, r27681, MPFR_RNDN);
        ;
        mpfr_pow(r27684, r27682, r27683, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r27687, r27686, r27678, MPFR_RNDN);
        mpfr_mul(r27688, r27679, r27687, MPFR_RNDN);
        ;
        mpfr_add(r27690, r27689, r27678, MPFR_RNDN);
        mpfr_mul(r27691, r27681, r27690, MPFR_RNDN);
        mpfr_add(r27692, r27688, r27691, MPFR_RNDN);
        mpfr_mul(r27693, r27685, r27692, MPFR_RNDN);
        mpfr_add(r27694, r27684, r27693, MPFR_RNDN);
        mpfr_sub(r27695, r27694, r27686, MPFR_RNDN);
        return mpfr_get_d(r27695, MPFR_RNDN);
}

static mpfr_t r27696, r27697, r27698, r27699, r27700, r27701, r27702, r27703, r27704, r27705, r27706, r27707, r27708, r27709, r27710, r27711, r27712, r27713;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27696);
        mpfr_init(r27697);
        mpfr_init(r27698);
        mpfr_init(r27699);
        mpfr_init(r27700);
        mpfr_init_set_str(r27701, "2", 10, MPFR_RNDN);
        mpfr_init(r27702);
        mpfr_init_set_str(r27703, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27704, "1", 10, MPFR_RNDN);
        mpfr_init(r27705);
        mpfr_init(r27706);
        mpfr_init_set_str(r27707, "3", 10, MPFR_RNDN);
        mpfr_init(r27708);
        mpfr_init(r27709);
        mpfr_init(r27710);
        mpfr_init(r27711);
        mpfr_init(r27712);
        mpfr_init(r27713);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27696, a, MPFR_RNDN);
        mpfr_mul(r27697, r27696, r27696, MPFR_RNDN);
        mpfr_set_d(r27698, b, MPFR_RNDN);
        mpfr_mul(r27699, r27698, r27698, MPFR_RNDN);
        mpfr_add(r27700, r27697, r27699, MPFR_RNDN);
        ;
        mpfr_pow(r27702, r27700, r27701, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r27705, r27704, r27696, MPFR_RNDN);
        mpfr_mul(r27706, r27697, r27705, MPFR_RNDN);
        ;
        mpfr_add(r27708, r27707, r27696, MPFR_RNDN);
        mpfr_mul(r27709, r27699, r27708, MPFR_RNDN);
        mpfr_add(r27710, r27706, r27709, MPFR_RNDN);
        mpfr_mul(r27711, r27703, r27710, MPFR_RNDN);
        mpfr_add(r27712, r27702, r27711, MPFR_RNDN);
        mpfr_sub(r27713, r27712, r27704, MPFR_RNDN);
        return mpfr_get_d(r27713, MPFR_RNDN);
}

static mpfr_t r27714, r27715, r27716, r27717, r27718, r27719, r27720, r27721, r27722, r27723, r27724, r27725, r27726, r27727, r27728, r27729, r27730, r27731;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27714);
        mpfr_init(r27715);
        mpfr_init(r27716);
        mpfr_init(r27717);
        mpfr_init(r27718);
        mpfr_init_set_str(r27719, "2", 10, MPFR_RNDN);
        mpfr_init(r27720);
        mpfr_init_set_str(r27721, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27722, "1", 10, MPFR_RNDN);
        mpfr_init(r27723);
        mpfr_init(r27724);
        mpfr_init_set_str(r27725, "3", 10, MPFR_RNDN);
        mpfr_init(r27726);
        mpfr_init(r27727);
        mpfr_init(r27728);
        mpfr_init(r27729);
        mpfr_init(r27730);
        mpfr_init(r27731);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27714, a, MPFR_RNDN);
        mpfr_mul(r27715, r27714, r27714, MPFR_RNDN);
        mpfr_set_d(r27716, b, MPFR_RNDN);
        mpfr_mul(r27717, r27716, r27716, MPFR_RNDN);
        mpfr_add(r27718, r27715, r27717, MPFR_RNDN);
        ;
        mpfr_pow(r27720, r27718, r27719, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r27723, r27722, r27714, MPFR_RNDN);
        mpfr_mul(r27724, r27715, r27723, MPFR_RNDN);
        ;
        mpfr_add(r27726, r27725, r27714, MPFR_RNDN);
        mpfr_mul(r27727, r27717, r27726, MPFR_RNDN);
        mpfr_add(r27728, r27724, r27727, MPFR_RNDN);
        mpfr_mul(r27729, r27721, r27728, MPFR_RNDN);
        mpfr_add(r27730, r27720, r27729, MPFR_RNDN);
        mpfr_sub(r27731, r27730, r27722, MPFR_RNDN);
        return mpfr_get_d(r27731, MPFR_RNDN);
}