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

char *name = "Falkner and Boettcher, Appendix A";

double f_if(float a, float k, float m) {
        float r27649 = a;
        float r27650 = k;
        float r27651 = m;
        float r27652 = pow(r27650, r27651);
        float r27653 = r27649 * r27652;
        float r27654 = 1;
        float r27655 = 10;
        float r27656 = r27655 * r27650;
        float r27657 = r27654 + r27656;
        float r27658 = r27650 * r27650;
        float r27659 = r27657 + r27658;
        float r27660 = r27653 / r27659;
        return r27660;
}

double f_id(double a, double k, double m) {
        double r27661 = a;
        double r27662 = k;
        double r27663 = m;
        double r27664 = pow(r27662, r27663);
        double r27665 = r27661 * r27664;
        double r27666 = 1;
        double r27667 = 10;
        double r27668 = r27667 * r27662;
        double r27669 = r27666 + r27668;
        double r27670 = r27662 * r27662;
        double r27671 = r27669 + r27670;
        double r27672 = r27665 / r27671;
        return r27672;
}


double f_of(float a, float k, float m) {
        float r27673 = k;
        float r27674 = 3.9501795571041663e+151;
        bool r27675 = r27673 <= r27674;
        float r27676 = m;
        float r27677 = pow(r27673, r27676);
        float r27678 = a;
        float r27679 = r27677 * r27678;
        float r27680 = 10;
        float r27681 = r27680 + r27673;
        float r27682 = 1;
        float r27683 = fma(r27673, r27681, r27682);
        float r27684 = sqrt(r27683);
        float r27685 = r27679 / r27684;
        float r27686 = r27680 * r27673;
        float r27687 = r27682 + r27686;
        float r27688 = r27673 * r27673;
        float r27689 = r27687 + r27688;
        float r27690 = sqrt(r27689);
        float r27691 = r27685 / r27690;
        float r27692 = r27673 / r27678;
        float r27693 = -r27676;
        float r27694 = pow(r27673, r27693);
        float r27695 = r27673 + r27680;
        float r27696 = r27694 * r27695;
        float r27697 = r27694 / r27678;
        float r27698 = fma(r27692, r27696, r27697);
        float r27699 = r27682 / r27698;
        float r27700 = r27675 ? r27691 : r27699;
        return r27700;
}

double f_od(double a, double k, double m) {
        double r27701 = k;
        double r27702 = 3.9501795571041663e+151;
        bool r27703 = r27701 <= r27702;
        double r27704 = m;
        double r27705 = pow(r27701, r27704);
        double r27706 = a;
        double r27707 = r27705 * r27706;
        double r27708 = 10;
        double r27709 = r27708 + r27701;
        double r27710 = 1;
        double r27711 = fma(r27701, r27709, r27710);
        double r27712 = sqrt(r27711);
        double r27713 = r27707 / r27712;
        double r27714 = r27708 * r27701;
        double r27715 = r27710 + r27714;
        double r27716 = r27701 * r27701;
        double r27717 = r27715 + r27716;
        double r27718 = sqrt(r27717);
        double r27719 = r27713 / r27718;
        double r27720 = r27701 / r27706;
        double r27721 = -r27704;
        double r27722 = pow(r27701, r27721);
        double r27723 = r27701 + r27708;
        double r27724 = r27722 * r27723;
        double r27725 = r27722 / r27706;
        double r27726 = fma(r27720, r27724, r27725);
        double r27727 = r27710 / r27726;
        double r27728 = r27703 ? r27719 : r27727;
        return r27728;
}

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 r27729, r27730, r27731, r27732, r27733, r27734, r27735, r27736, r27737, r27738, r27739, r27740;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27729);
        mpfr_init(r27730);
        mpfr_init(r27731);
        mpfr_init(r27732);
        mpfr_init(r27733);
        mpfr_init_set_str(r27734, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27735, "10", 10, MPFR_RNDN);
        mpfr_init(r27736);
        mpfr_init(r27737);
        mpfr_init(r27738);
        mpfr_init(r27739);
        mpfr_init(r27740);
}

double f_im(double a, double k, double m) {
        mpfr_set_d(r27729, a, MPFR_RNDN);
        mpfr_set_d(r27730, k, MPFR_RNDN);
        mpfr_set_d(r27731, m, MPFR_RNDN);
        mpfr_pow(r27732, r27730, r27731, MPFR_RNDN);
        mpfr_mul(r27733, r27729, r27732, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r27736, r27735, r27730, MPFR_RNDN);
        mpfr_add(r27737, r27734, r27736, MPFR_RNDN);
        mpfr_mul(r27738, r27730, r27730, MPFR_RNDN);
        mpfr_add(r27739, r27737, r27738, MPFR_RNDN);
        mpfr_div(r27740, r27733, r27739, MPFR_RNDN);
        return mpfr_get_d(r27740, MPFR_RNDN);
}

static mpfr_t r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27741);
        mpfr_init_set_str(r27742, "3.9501795571041663e+151", 10, MPFR_RNDN);
        mpfr_init(r27743);
        mpfr_init(r27744);
        mpfr_init(r27745);
        mpfr_init(r27746);
        mpfr_init(r27747);
        mpfr_init_set_str(r27748, "10", 10, MPFR_RNDN);
        mpfr_init(r27749);
        mpfr_init_set_str(r27750, "1", 10, MPFR_RNDN);
        mpfr_init(r27751);
        mpfr_init(r27752);
        mpfr_init(r27753);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init(r27756);
        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);
        mpfr_init(r27768);
}

double f_fm(double a, double k, double m) {
        mpfr_set_d(r27741, k, MPFR_RNDN);
        ;
        mpfr_set_si(r27743, mpfr_cmp(r27741, r27742) <= 0, MPFR_RNDN);
        mpfr_set_d(r27744, m, MPFR_RNDN);
        mpfr_pow(r27745, r27741, r27744, MPFR_RNDN);
        mpfr_set_d(r27746, a, MPFR_RNDN);
        mpfr_mul(r27747, r27745, r27746, MPFR_RNDN);
        ;
        mpfr_add(r27749, r27748, r27741, MPFR_RNDN);
        ;
        mpfr_fma(r27751, r27741, r27749, r27750, MPFR_RNDN);
        mpfr_sqrt(r27752, r27751, MPFR_RNDN);
        mpfr_div(r27753, r27747, r27752, MPFR_RNDN);
        mpfr_mul(r27754, r27748, r27741, MPFR_RNDN);
        mpfr_add(r27755, r27750, r27754, MPFR_RNDN);
        mpfr_mul(r27756, r27741, r27741, MPFR_RNDN);
        mpfr_add(r27757, r27755, r27756, MPFR_RNDN);
        mpfr_sqrt(r27758, r27757, MPFR_RNDN);
        mpfr_div(r27759, r27753, r27758, MPFR_RNDN);
        mpfr_div(r27760, r27741, r27746, MPFR_RNDN);
        mpfr_neg(r27761, r27744, MPFR_RNDN);
        mpfr_pow(r27762, r27741, r27761, MPFR_RNDN);
        mpfr_add(r27763, r27741, r27748, MPFR_RNDN);
        mpfr_mul(r27764, r27762, r27763, MPFR_RNDN);
        mpfr_div(r27765, r27762, r27746, MPFR_RNDN);
        mpfr_fma(r27766, r27760, r27764, r27765, MPFR_RNDN);
        mpfr_div(r27767, r27750, r27766, MPFR_RNDN);
        if (mpfr_get_si(r27743, MPFR_RNDN)) { mpfr_set(r27768, r27759, MPFR_RNDN); } else { mpfr_set(r27768, r27767, MPFR_RNDN); };
        return mpfr_get_d(r27768, MPFR_RNDN);
}

static mpfr_t r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27769);
        mpfr_init_set_str(r27770, "3.9501795571041663e+151", 10, MPFR_RNDN);
        mpfr_init(r27771);
        mpfr_init(r27772);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init(r27775);
        mpfr_init_set_str(r27776, "10", 10, MPFR_RNDN);
        mpfr_init(r27777);
        mpfr_init_set_str(r27778, "1", 10, MPFR_RNDN);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init(r27782);
        mpfr_init(r27783);
        mpfr_init(r27784);
        mpfr_init(r27785);
        mpfr_init(r27786);
        mpfr_init(r27787);
        mpfr_init(r27788);
        mpfr_init(r27789);
        mpfr_init(r27790);
        mpfr_init(r27791);
        mpfr_init(r27792);
        mpfr_init(r27793);
        mpfr_init(r27794);
        mpfr_init(r27795);
        mpfr_init(r27796);
}

double f_dm(double a, double k, double m) {
        mpfr_set_d(r27769, k, MPFR_RNDN);
        ;
        mpfr_set_si(r27771, mpfr_cmp(r27769, r27770) <= 0, MPFR_RNDN);
        mpfr_set_d(r27772, m, MPFR_RNDN);
        mpfr_pow(r27773, r27769, r27772, MPFR_RNDN);
        mpfr_set_d(r27774, a, MPFR_RNDN);
        mpfr_mul(r27775, r27773, r27774, MPFR_RNDN);
        ;
        mpfr_add(r27777, r27776, r27769, MPFR_RNDN);
        ;
        mpfr_fma(r27779, r27769, r27777, r27778, MPFR_RNDN);
        mpfr_sqrt(r27780, r27779, MPFR_RNDN);
        mpfr_div(r27781, r27775, r27780, MPFR_RNDN);
        mpfr_mul(r27782, r27776, r27769, MPFR_RNDN);
        mpfr_add(r27783, r27778, r27782, MPFR_RNDN);
        mpfr_mul(r27784, r27769, r27769, MPFR_RNDN);
        mpfr_add(r27785, r27783, r27784, MPFR_RNDN);
        mpfr_sqrt(r27786, r27785, MPFR_RNDN);
        mpfr_div(r27787, r27781, r27786, MPFR_RNDN);
        mpfr_div(r27788, r27769, r27774, MPFR_RNDN);
        mpfr_neg(r27789, r27772, MPFR_RNDN);
        mpfr_pow(r27790, r27769, r27789, MPFR_RNDN);
        mpfr_add(r27791, r27769, r27776, MPFR_RNDN);
        mpfr_mul(r27792, r27790, r27791, MPFR_RNDN);
        mpfr_div(r27793, r27790, r27774, MPFR_RNDN);
        mpfr_fma(r27794, r27788, r27792, r27793, MPFR_RNDN);
        mpfr_div(r27795, r27778, r27794, MPFR_RNDN);
        if (mpfr_get_si(r27771, MPFR_RNDN)) { mpfr_set(r27796, r27787, MPFR_RNDN); } else { mpfr_set(r27796, r27795, MPFR_RNDN); };
        return mpfr_get_d(r27796, MPFR_RNDN);
}