#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 r27578 = a;
        float r27579 = r27578 * r27578;
        float r27580 = b;
        float r27581 = r27580 * r27580;
        float r27582 = r27579 + r27581;
        float r27583 = 2;
        float r27584 = pow(r27582, r27583);
        float r27585 = 4;
        float r27586 = 1;
        float r27587 = r27586 + r27578;
        float r27588 = r27579 * r27587;
        float r27589 = 3;
        float r27590 = r27589 * r27578;
        float r27591 = r27586 - r27590;
        float r27592 = r27581 * r27591;
        float r27593 = r27588 + r27592;
        float r27594 = r27585 * r27593;
        float r27595 = r27584 + r27594;
        float r27596 = r27595 - r27586;
        return r27596;
}

double f_id(double a, double b) {
        double r27597 = a;
        double r27598 = r27597 * r27597;
        double r27599 = b;
        double r27600 = r27599 * r27599;
        double r27601 = r27598 + r27600;
        double r27602 = 2;
        double r27603 = pow(r27601, r27602);
        double r27604 = 4;
        double r27605 = 1;
        double r27606 = r27605 + r27597;
        double r27607 = r27598 * r27606;
        double r27608 = 3;
        double r27609 = r27608 * r27597;
        double r27610 = r27605 - r27609;
        double r27611 = r27600 * r27610;
        double r27612 = r27607 + r27611;
        double r27613 = r27604 * r27612;
        double r27614 = r27603 + r27613;
        double r27615 = r27614 - r27605;
        return r27615;
}


double f_of(float a, float b) {
        float r27616 = a;
        float r27617 = r27616 * r27616;
        float r27618 = b;
        float r27619 = r27618 * r27618;
        float r27620 = r27617 + r27619;
        float r27621 = 2;
        float r27622 = pow(r27620, r27621);
        float r27623 = 4;
        float r27624 = 1;
        float r27625 = r27624 + r27616;
        float r27626 = r27617 * r27625;
        float r27627 = 3;
        float r27628 = r27627 * r27616;
        float r27629 = r27624 - r27628;
        float r27630 = r27619 * r27629;
        float r27631 = r27626 + r27630;
        float r27632 = r27623 * r27631;
        float r27633 = r27622 + r27632;
        float r27634 = r27633 - r27624;
        return r27634;
}

double f_od(double a, double b) {
        double r27635 = a;
        double r27636 = r27635 * r27635;
        double r27637 = b;
        double r27638 = r27637 * r27637;
        double r27639 = r27636 + r27638;
        double r27640 = 2;
        double r27641 = pow(r27639, r27640);
        double r27642 = 4;
        double r27643 = 1;
        double r27644 = r27643 + r27635;
        double r27645 = r27636 * r27644;
        double r27646 = 3;
        double r27647 = r27646 * r27635;
        double r27648 = r27643 - r27647;
        double r27649 = r27638 * r27648;
        double r27650 = r27645 + r27649;
        double r27651 = r27642 * r27650;
        double r27652 = r27641 + r27651;
        double r27653 = r27652 - r27643;
        return r27653;
}

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 r27654, r27655, r27656, r27657, r27658, r27659, r27660, r27661, r27662, r27663, r27664, r27665, r27666, r27667, r27668, r27669, r27670, r27671, r27672;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27654);
        mpfr_init(r27655);
        mpfr_init(r27656);
        mpfr_init(r27657);
        mpfr_init(r27658);
        mpfr_init_set_str(r27659, "2", 10, MPFR_RNDN);
        mpfr_init(r27660);
        mpfr_init_set_str(r27661, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27662, "1", 10, MPFR_RNDN);
        mpfr_init(r27663);
        mpfr_init(r27664);
        mpfr_init_set_str(r27665, "3", 10, MPFR_RNDN);
        mpfr_init(r27666);
        mpfr_init(r27667);
        mpfr_init(r27668);
        mpfr_init(r27669);
        mpfr_init(r27670);
        mpfr_init(r27671);
        mpfr_init(r27672);
}

double f_im(double a, double b) {
        mpfr_set_d(r27654, a, MPFR_RNDN);
        mpfr_mul(r27655, r27654, r27654, MPFR_RNDN);
        mpfr_set_d(r27656, b, MPFR_RNDN);
        mpfr_mul(r27657, r27656, r27656, MPFR_RNDN);
        mpfr_add(r27658, r27655, r27657, MPFR_RNDN);
        ;
        mpfr_pow(r27660, r27658, r27659, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27663, r27662, r27654, MPFR_RNDN);
        mpfr_mul(r27664, r27655, r27663, MPFR_RNDN);
        ;
        mpfr_mul(r27666, r27665, r27654, MPFR_RNDN);
        mpfr_sub(r27667, r27662, r27666, MPFR_RNDN);
        mpfr_mul(r27668, r27657, r27667, MPFR_RNDN);
        mpfr_add(r27669, r27664, r27668, MPFR_RNDN);
        mpfr_mul(r27670, r27661, r27669, MPFR_RNDN);
        mpfr_add(r27671, r27660, r27670, MPFR_RNDN);
        mpfr_sub(r27672, r27671, r27662, MPFR_RNDN);
        return mpfr_get_d(r27672, MPFR_RNDN);
}

static mpfr_t r27673, r27674, r27675, r27676, r27677, r27678, r27679, r27680, r27681, r27682, r27683, r27684, r27685, r27686, r27687, r27688, r27689, r27690, r27691;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27673);
        mpfr_init(r27674);
        mpfr_init(r27675);
        mpfr_init(r27676);
        mpfr_init(r27677);
        mpfr_init_set_str(r27678, "2", 10, MPFR_RNDN);
        mpfr_init(r27679);
        mpfr_init_set_str(r27680, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27681, "1", 10, MPFR_RNDN);
        mpfr_init(r27682);
        mpfr_init(r27683);
        mpfr_init_set_str(r27684, "3", 10, MPFR_RNDN);
        mpfr_init(r27685);
        mpfr_init(r27686);
        mpfr_init(r27687);
        mpfr_init(r27688);
        mpfr_init(r27689);
        mpfr_init(r27690);
        mpfr_init(r27691);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27673, a, MPFR_RNDN);
        mpfr_mul(r27674, r27673, r27673, MPFR_RNDN);
        mpfr_set_d(r27675, b, MPFR_RNDN);
        mpfr_mul(r27676, r27675, r27675, MPFR_RNDN);
        mpfr_add(r27677, r27674, r27676, MPFR_RNDN);
        ;
        mpfr_pow(r27679, r27677, r27678, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27682, r27681, r27673, MPFR_RNDN);
        mpfr_mul(r27683, r27674, r27682, MPFR_RNDN);
        ;
        mpfr_mul(r27685, r27684, r27673, MPFR_RNDN);
        mpfr_sub(r27686, r27681, r27685, MPFR_RNDN);
        mpfr_mul(r27687, r27676, r27686, MPFR_RNDN);
        mpfr_add(r27688, r27683, r27687, MPFR_RNDN);
        mpfr_mul(r27689, r27680, r27688, MPFR_RNDN);
        mpfr_add(r27690, r27679, r27689, MPFR_RNDN);
        mpfr_sub(r27691, r27690, r27681, MPFR_RNDN);
        return mpfr_get_d(r27691, MPFR_RNDN);
}

static mpfr_t r27692, r27693, r27694, r27695, r27696, r27697, r27698, r27699, r27700, r27701, r27702, r27703, r27704, r27705, r27706, r27707, r27708, r27709, r27710;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27692);
        mpfr_init(r27693);
        mpfr_init(r27694);
        mpfr_init(r27695);
        mpfr_init(r27696);
        mpfr_init_set_str(r27697, "2", 10, MPFR_RNDN);
        mpfr_init(r27698);
        mpfr_init_set_str(r27699, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r27700, "1", 10, MPFR_RNDN);
        mpfr_init(r27701);
        mpfr_init(r27702);
        mpfr_init_set_str(r27703, "3", 10, MPFR_RNDN);
        mpfr_init(r27704);
        mpfr_init(r27705);
        mpfr_init(r27706);
        mpfr_init(r27707);
        mpfr_init(r27708);
        mpfr_init(r27709);
        mpfr_init(r27710);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27692, a, MPFR_RNDN);
        mpfr_mul(r27693, r27692, r27692, MPFR_RNDN);
        mpfr_set_d(r27694, b, MPFR_RNDN);
        mpfr_mul(r27695, r27694, r27694, MPFR_RNDN);
        mpfr_add(r27696, r27693, r27695, MPFR_RNDN);
        ;
        mpfr_pow(r27698, r27696, r27697, MPFR_RNDN);
        ;
        ;
        mpfr_add(r27701, r27700, r27692, MPFR_RNDN);
        mpfr_mul(r27702, r27693, r27701, MPFR_RNDN);
        ;
        mpfr_mul(r27704, r27703, r27692, MPFR_RNDN);
        mpfr_sub(r27705, r27700, r27704, MPFR_RNDN);
        mpfr_mul(r27706, r27695, r27705, MPFR_RNDN);
        mpfr_add(r27707, r27702, r27706, MPFR_RNDN);
        mpfr_mul(r27708, r27699, r27707, MPFR_RNDN);
        mpfr_add(r27709, r27698, r27708, MPFR_RNDN);
        mpfr_sub(r27710, r27709, r27700, MPFR_RNDN);
        return mpfr_get_d(r27710, MPFR_RNDN);
}