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

char *name = "Asymptote A";

double f_if(float x) {
        float r26685 = 1;
        float r26686 = x;
        float r26687 = r26686 + r26685;
        float r26688 = r26685 / r26687;
        float r26689 = r26686 - r26685;
        float r26690 = r26685 / r26689;
        float r26691 = r26688 - r26690;
        return r26691;
}

double f_id(double x) {
        double r26692 = 1;
        double r26693 = x;
        double r26694 = r26693 + r26692;
        double r26695 = r26692 / r26694;
        double r26696 = r26693 - r26692;
        double r26697 = r26692 / r26696;
        double r26698 = r26695 - r26697;
        return r26698;
}


double f_of(float x) {
        float r26699 = 1;
        float r26700 = r26699 + r26699;
        float r26701 = x;
        float r26702 = r26699 + r26701;
        float r26703 = r26700 / r26702;
        float r26704 = -r26703;
        float r26705 = r26701 - r26699;
        float r26706 = r26704 / r26705;
        return r26706;
}

double f_od(double x) {
        double r26707 = 1;
        double r26708 = r26707 + r26707;
        double r26709 = x;
        double r26710 = r26707 + r26709;
        double r26711 = r26708 / r26710;
        double r26712 = -r26711;
        double r26713 = r26709 - r26707;
        double r26714 = r26712 / r26713;
        return r26714;
}

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 r26715, r26716, r26717, r26718, r26719, r26720, r26721;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r26715, "1", 10, MPFR_RNDN);
        mpfr_init(r26716);
        mpfr_init(r26717);
        mpfr_init(r26718);
        mpfr_init(r26719);
        mpfr_init(r26720);
        mpfr_init(r26721);
}

double f_im(double x) {
        ;
        mpfr_set_d(r26716, x, MPFR_RNDN);
        mpfr_add(r26717, r26716, r26715, MPFR_RNDN);
        mpfr_div(r26718, r26715, r26717, MPFR_RNDN);
        mpfr_sub(r26719, r26716, r26715, MPFR_RNDN);
        mpfr_div(r26720, r26715, r26719, MPFR_RNDN);
        mpfr_sub(r26721, r26718, r26720, MPFR_RNDN);
        return mpfr_get_d(r26721, MPFR_RNDN);
}

static mpfr_t r26722, r26723, r26724, r26725, r26726, r26727, r26728, r26729;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r26722, "1", 10, MPFR_RNDN);
        mpfr_init(r26723);
        mpfr_init(r26724);
        mpfr_init(r26725);
        mpfr_init(r26726);
        mpfr_init(r26727);
        mpfr_init(r26728);
        mpfr_init(r26729);
}

double f_fm(double x) {
        ;
        mpfr_add(r26723, r26722, r26722, MPFR_RNDN);
        mpfr_set_d(r26724, x, MPFR_RNDN);
        mpfr_add(r26725, r26722, r26724, MPFR_RNDN);
        mpfr_div(r26726, r26723, r26725, MPFR_RNDN);
        mpfr_neg(r26727, r26726, MPFR_RNDN);
        mpfr_sub(r26728, r26724, r26722, MPFR_RNDN);
        mpfr_div(r26729, r26727, r26728, MPFR_RNDN);
        return mpfr_get_d(r26729, MPFR_RNDN);
}

static mpfr_t r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r26730, "1", 10, MPFR_RNDN);
        mpfr_init(r26731);
        mpfr_init(r26732);
        mpfr_init(r26733);
        mpfr_init(r26734);
        mpfr_init(r26735);
        mpfr_init(r26736);
        mpfr_init(r26737);
}

double f_dm(double x) {
        ;
        mpfr_add(r26731, r26730, r26730, MPFR_RNDN);
        mpfr_set_d(r26732, x, MPFR_RNDN);
        mpfr_add(r26733, r26730, r26732, MPFR_RNDN);
        mpfr_div(r26734, r26731, r26733, MPFR_RNDN);
        mpfr_neg(r26735, r26734, MPFR_RNDN);
        mpfr_sub(r26736, r26732, r26730, MPFR_RNDN);
        mpfr_div(r26737, r26735, r26736, MPFR_RNDN);
        return mpfr_get_d(r26737, MPFR_RNDN);
}