#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 r10664 = 1;
        float r10665 = x;
        float r10666 = r10665 + r10664;
        float r10667 = r10664 / r10666;
        float r10668 = r10665 - r10664;
        float r10669 = r10664 / r10668;
        float r10670 = r10667 - r10669;
        return r10670;
}

double f_id(double x) {
        double r10671 = 1;
        double r10672 = x;
        double r10673 = r10672 + r10671;
        double r10674 = r10671 / r10673;
        double r10675 = r10672 - r10671;
        double r10676 = r10671 / r10675;
        double r10677 = r10674 - r10676;
        return r10677;
}


double f_of(float x) {
        float r10678 = 1;
        float r10679 = -r10678;
        float r10680 = r10679 - r10678;
        float r10681 = x;
        float r10682 = r10681 + r10678;
        float r10683 = r10680 / r10682;
        float r10684 = r10681 - r10678;
        float r10685 = r10683 / r10684;
        return r10685;
}

double f_od(double x) {
        double r10686 = 1;
        double r10687 = -r10686;
        double r10688 = r10687 - r10686;
        double r10689 = x;
        double r10690 = r10689 + r10686;
        double r10691 = r10688 / r10690;
        double r10692 = r10689 - r10686;
        double r10693 = r10691 / r10692;
        return r10693;
}

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 r10694, r10695, r10696, r10697, r10698, r10699, r10700;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r10694, "1", 10, MPFR_RNDN);
        mpfr_init(r10695);
        mpfr_init(r10696);
        mpfr_init(r10697);
        mpfr_init(r10698);
        mpfr_init(r10699);
        mpfr_init(r10700);
}

double f_im(double x) {
        ;
        mpfr_set_d(r10695, x, MPFR_RNDN);
        mpfr_add(r10696, r10695, r10694, MPFR_RNDN);
        mpfr_div(r10697, r10694, r10696, MPFR_RNDN);
        mpfr_sub(r10698, r10695, r10694, MPFR_RNDN);
        mpfr_div(r10699, r10694, r10698, MPFR_RNDN);
        mpfr_sub(r10700, r10697, r10699, MPFR_RNDN);
        return mpfr_get_d(r10700, MPFR_RNDN);
}

static mpfr_t r10701, r10702, r10703, r10704, r10705, r10706, r10707, r10708;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r10701, "1", 10, MPFR_RNDN);
        mpfr_init(r10702);
        mpfr_init(r10703);
        mpfr_init(r10704);
        mpfr_init(r10705);
        mpfr_init(r10706);
        mpfr_init(r10707);
        mpfr_init(r10708);
}

double f_fm(double x) {
        ;
        mpfr_neg(r10702, r10701, MPFR_RNDN);
        mpfr_sub(r10703, r10702, r10701, MPFR_RNDN);
        mpfr_set_d(r10704, x, MPFR_RNDN);
        mpfr_add(r10705, r10704, r10701, MPFR_RNDN);
        mpfr_div(r10706, r10703, r10705, MPFR_RNDN);
        mpfr_sub(r10707, r10704, r10701, MPFR_RNDN);
        mpfr_div(r10708, r10706, r10707, MPFR_RNDN);
        return mpfr_get_d(r10708, MPFR_RNDN);
}

static mpfr_t r10709, r10710, r10711, r10712, r10713, r10714, r10715, r10716;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r10709, "1", 10, MPFR_RNDN);
        mpfr_init(r10710);
        mpfr_init(r10711);
        mpfr_init(r10712);
        mpfr_init(r10713);
        mpfr_init(r10714);
        mpfr_init(r10715);
        mpfr_init(r10716);
}

double f_dm(double x) {
        ;
        mpfr_neg(r10710, r10709, MPFR_RNDN);
        mpfr_sub(r10711, r10710, r10709, MPFR_RNDN);
        mpfr_set_d(r10712, x, MPFR_RNDN);
        mpfr_add(r10713, r10712, r10709, MPFR_RNDN);
        mpfr_div(r10714, r10711, r10713, MPFR_RNDN);
        mpfr_sub(r10715, r10712, r10709, MPFR_RNDN);
        mpfr_div(r10716, r10714, r10715, MPFR_RNDN);
        return mpfr_get_d(r10716, MPFR_RNDN);
}