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

char *name = "Numeric.Integration.TanhSinh:nonNegative from integration-0.2.1";

double f_if(float x) {
        float r28064 = x;
        float r28065 = 1.0;
        float r28066 = r28065 - r28064;
        float r28067 = r28064 / r28066;
        return r28067;
}

double f_id(double x) {
        double r28068 = x;
        double r28069 = 1.0;
        double r28070 = r28069 - r28068;
        double r28071 = r28068 / r28070;
        return r28071;
}


double f_of(float x) {
        float r28072 = x;
        float r28073 = 1.0;
        float r28074 = r28073 - r28072;
        float r28075 = r28072 / r28074;
        return r28075;
}

double f_od(double x) {
        double r28076 = x;
        double r28077 = 1.0;
        double r28078 = r28077 - r28076;
        double r28079 = r28076 / r28078;
        return r28079;
}

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 r28080, r28081, r28082, r28083;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r28080);
        mpfr_init_set_str(r28081, "1.0", 10, MPFR_RNDN);
        mpfr_init(r28082);
        mpfr_init(r28083);
}

double f_im(double x) {
        mpfr_set_d(r28080, x, MPFR_RNDN);
        ;
        mpfr_sub(r28082, r28081, r28080, MPFR_RNDN);
        mpfr_div(r28083, r28080, r28082, MPFR_RNDN);
        return mpfr_get_d(r28083, MPFR_RNDN);
}

static mpfr_t r28084, r28085, r28086, r28087;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r28084);
        mpfr_init_set_str(r28085, "1.0", 10, MPFR_RNDN);
        mpfr_init(r28086);
        mpfr_init(r28087);
}

double f_fm(double x) {
        mpfr_set_d(r28084, x, MPFR_RNDN);
        ;
        mpfr_sub(r28086, r28085, r28084, MPFR_RNDN);
        mpfr_div(r28087, r28084, r28086, MPFR_RNDN);
        return mpfr_get_d(r28087, MPFR_RNDN);
}

static mpfr_t r28088, r28089, r28090, r28091;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r28088);
        mpfr_init_set_str(r28089, "1.0", 10, MPFR_RNDN);
        mpfr_init(r28090);
        mpfr_init(r28091);
}

double f_dm(double x) {
        mpfr_set_d(r28088, x, MPFR_RNDN);
        ;
        mpfr_sub(r28090, r28089, r28088, MPFR_RNDN);
        mpfr_div(r28091, r28088, r28090, MPFR_RNDN);
        return mpfr_get_d(r28091, MPFR_RNDN);
}