#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 r27086 = 1;
        float r27087 = x;
        float r27088 = r27087 - r27086;
        float r27089 = r27086 / r27088;
        float r27090 = r27087 + r27086;
        float r27091 = r27087 / r27090;
        float r27092 = r27089 + r27091;
        return r27092;
}

double f_id(double x) {
        double r27093 = 1;
        double r27094 = x;
        double r27095 = r27094 - r27093;
        double r27096 = r27093 / r27095;
        double r27097 = r27094 + r27093;
        double r27098 = r27094 / r27097;
        double r27099 = r27096 + r27098;
        return r27099;
}


double f_of(float x) {
        float r27100 = 1;
        float r27101 = x;
        float r27102 = r27101 - r27100;
        float r27103 = r27100 / r27102;
        float r27104 = cbrt(r27103);
        float r27105 = r27104 * r27104;
        float r27106 = r27101 + r27100;
        float r27107 = r27101 / r27106;
        float r27108 = fma(r27105, r27104, r27107);
        return r27108;
}

double f_od(double x) {
        double r27109 = 1;
        double r27110 = x;
        double r27111 = r27110 - r27109;
        double r27112 = r27109 / r27111;
        double r27113 = cbrt(r27112);
        double r27114 = r27113 * r27113;
        double r27115 = r27110 + r27109;
        double r27116 = r27110 / r27115;
        double r27117 = fma(r27114, r27113, r27116);
        return r27117;
}

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 r27118, r27119, r27120, r27121, r27122, r27123, r27124;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27118, "1", 10, MPFR_RNDN);
        mpfr_init(r27119);
        mpfr_init(r27120);
        mpfr_init(r27121);
        mpfr_init(r27122);
        mpfr_init(r27123);
        mpfr_init(r27124);
}

double f_im(double x) {
        ;
        mpfr_set_d(r27119, x, MPFR_RNDN);
        mpfr_sub(r27120, r27119, r27118, MPFR_RNDN);
        mpfr_div(r27121, r27118, r27120, MPFR_RNDN);
        mpfr_add(r27122, r27119, r27118, MPFR_RNDN);
        mpfr_div(r27123, r27119, r27122, MPFR_RNDN);
        mpfr_add(r27124, r27121, r27123, MPFR_RNDN);
        return mpfr_get_d(r27124, MPFR_RNDN);
}

static mpfr_t r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27125, "1", 10, MPFR_RNDN);
        mpfr_init(r27126);
        mpfr_init(r27127);
        mpfr_init(r27128);
        mpfr_init(r27129);
        mpfr_init(r27130);
        mpfr_init(r27131);
        mpfr_init(r27132);
        mpfr_init(r27133);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27126, x, MPFR_RNDN);
        mpfr_sub(r27127, r27126, r27125, MPFR_RNDN);
        mpfr_div(r27128, r27125, r27127, MPFR_RNDN);
        mpfr_cbrt(r27129, r27128, MPFR_RNDN);
        mpfr_mul(r27130, r27129, r27129, MPFR_RNDN);
        mpfr_add(r27131, r27126, r27125, MPFR_RNDN);
        mpfr_div(r27132, r27126, r27131, MPFR_RNDN);
        mpfr_fma(r27133, r27130, r27129, r27132, MPFR_RNDN);
        return mpfr_get_d(r27133, MPFR_RNDN);
}

static mpfr_t r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27134, "1", 10, MPFR_RNDN);
        mpfr_init(r27135);
        mpfr_init(r27136);
        mpfr_init(r27137);
        mpfr_init(r27138);
        mpfr_init(r27139);
        mpfr_init(r27140);
        mpfr_init(r27141);
        mpfr_init(r27142);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27135, x, MPFR_RNDN);
        mpfr_sub(r27136, r27135, r27134, MPFR_RNDN);
        mpfr_div(r27137, r27134, r27136, MPFR_RNDN);
        mpfr_cbrt(r27138, r27137, MPFR_RNDN);
        mpfr_mul(r27139, r27138, r27138, MPFR_RNDN);
        mpfr_add(r27140, r27135, r27134, MPFR_RNDN);
        mpfr_div(r27141, r27135, r27140, MPFR_RNDN);
        mpfr_fma(r27142, r27139, r27138, r27141, MPFR_RNDN);
        return mpfr_get_d(r27142, MPFR_RNDN);
}