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

char *name = "Asymptote B";

double f_if(float x) {
        float r10243 = 1.0f;
        float r10244 = x;
        float r10245 = r10244 - r10243;
        float r10246 = r10243 / r10245;
        float r10247 = r10244 + r10243;
        float r10248 = r10244 / r10247;
        float r10249 = r10246 + r10248;
        return r10249;
}

double f_id(double x) {
        double r10250 = 1.0;
        double r10251 = x;
        double r10252 = r10251 - r10250;
        double r10253 = r10250 / r10252;
        double r10254 = r10251 + r10250;
        double r10255 = r10251 / r10254;
        double r10256 = r10253 + r10255;
        return r10256;
}


double f_of(float x) {
        float r10257 = 1.0f;
        float r10258 = x;
        float r10259 = r10258 - r10257;
        float r10260 = r10257 / r10259;
        float r10261 = r10258 + r10257;
        float r10262 = r10258 / r10261;
        float r10263 = r10260 + r10262;
        return r10263;
}

double f_od(double x) {
        double r10264 = 1.0;
        double r10265 = x;
        double r10266 = r10265 - r10264;
        double r10267 = r10264 / r10266;
        double r10268 = r10265 + r10264;
        double r10269 = r10265 / r10268;
        double r10270 = r10267 + r10269;
        return r10270;
}

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 r10271, r10272, r10273, r10274, r10275, r10276, r10277;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r10271, "1", 10, MPFR_RNDN);
        mpfr_init(r10272);
        mpfr_init(r10273);
        mpfr_init(r10274);
        mpfr_init(r10275);
        mpfr_init(r10276);
        mpfr_init(r10277);
}

double f_im(double x) {
        ;
        mpfr_set_d(r10272, x, MPFR_RNDN);
        mpfr_sub(r10273, r10272, r10271, MPFR_RNDN);
        mpfr_div(r10274, r10271, r10273, MPFR_RNDN);
        mpfr_add(r10275, r10272, r10271, MPFR_RNDN);
        mpfr_div(r10276, r10272, r10275, MPFR_RNDN);
        mpfr_add(r10277, r10274, r10276, MPFR_RNDN);
        return mpfr_get_d(r10277, MPFR_RNDN);
}

static mpfr_t r10278, r10279, r10280, r10281, r10282, r10283, r10284;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r10278, "1", 10, MPFR_RNDN);
        mpfr_init(r10279);
        mpfr_init(r10280);
        mpfr_init(r10281);
        mpfr_init(r10282);
        mpfr_init(r10283);
        mpfr_init(r10284);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10279, x, MPFR_RNDN);
        mpfr_sub(r10280, r10279, r10278, MPFR_RNDN);
        mpfr_div(r10281, r10278, r10280, MPFR_RNDN);
        mpfr_add(r10282, r10279, r10278, MPFR_RNDN);
        mpfr_div(r10283, r10279, r10282, MPFR_RNDN);
        mpfr_add(r10284, r10281, r10283, MPFR_RNDN);
        return mpfr_get_d(r10284, MPFR_RNDN);
}

static mpfr_t r10285, r10286, r10287, r10288, r10289, r10290, r10291;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r10285, "1", 10, MPFR_RNDN);
        mpfr_init(r10286);
        mpfr_init(r10287);
        mpfr_init(r10288);
        mpfr_init(r10289);
        mpfr_init(r10290);
        mpfr_init(r10291);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10286, x, MPFR_RNDN);
        mpfr_sub(r10287, r10286, r10285, MPFR_RNDN);
        mpfr_div(r10288, r10285, r10287, MPFR_RNDN);
        mpfr_add(r10289, r10286, r10285, MPFR_RNDN);
        mpfr_div(r10290, r10286, r10289, MPFR_RNDN);
        mpfr_add(r10291, r10288, r10290, MPFR_RNDN);
        return mpfr_get_d(r10291, MPFR_RNDN);
}