#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 r27756 = x;
        float r27757 = 1.0;
        float r27758 = r27757 - r27756;
        float r27759 = r27756 / r27758;
        return r27759;
}

double f_id(double x) {
        double r27760 = x;
        double r27761 = 1.0;
        double r27762 = r27761 - r27760;
        double r27763 = r27760 / r27762;
        return r27763;
}


double f_of(float x) {
        float r27764 = x;
        float r27765 = 1.0;
        float r27766 = r27765 - r27764;
        float r27767 = r27764 / r27766;
        return r27767;
}

double f_od(double x) {
        double r27768 = x;
        double r27769 = 1.0;
        double r27770 = r27769 - r27768;
        double r27771 = r27768 / r27770;
        return r27771;
}

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 r27772, r27773, r27774, r27775;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r27772);
        mpfr_init_set_str(r27773, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27774);
        mpfr_init(r27775);
}

double f_im(double x) {
        mpfr_set_d(r27772, x, MPFR_RNDN);
        ;
        mpfr_sub(r27774, r27773, r27772, MPFR_RNDN);
        mpfr_div(r27775, r27772, r27774, MPFR_RNDN);
        return mpfr_get_d(r27775, MPFR_RNDN);
}

static mpfr_t r27776, r27777, r27778, r27779;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27776);
        mpfr_init_set_str(r27777, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27778);
        mpfr_init(r27779);
}

double f_fm(double x) {
        mpfr_set_d(r27776, x, MPFR_RNDN);
        ;
        mpfr_sub(r27778, r27777, r27776, MPFR_RNDN);
        mpfr_div(r27779, r27776, r27778, MPFR_RNDN);
        return mpfr_get_d(r27779, MPFR_RNDN);
}

static mpfr_t r27780, r27781, r27782, r27783;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r27780);
        mpfr_init_set_str(r27781, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27782);
        mpfr_init(r27783);
}

double f_dm(double x) {
        mpfr_set_d(r27780, x, MPFR_RNDN);
        ;
        mpfr_sub(r27782, r27781, r27780, MPFR_RNDN);
        mpfr_div(r27783, r27780, r27782, MPFR_RNDN);
        return mpfr_get_d(r27783, MPFR_RNDN);
}