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

char *name = "Kahan's exp quotient";

double f_if(float x) {
        float r25769 = x;
        float r25770 = exp(r25769);
        float r25771 = 1;
        float r25772 = r25770 - r25771;
        float r25773 = r25772 / r25769;
        return r25773;
}

double f_id(double x) {
        double r25774 = x;
        double r25775 = exp(r25774);
        double r25776 = 1;
        double r25777 = r25775 - r25776;
        double r25778 = r25777 / r25774;
        return r25778;
}


double f_of(float x) {
        float r25779 = x;
        float r25780 = expm1(r25779);
        float r25781 = r25780 / r25779;
        float r25782 = expm1(r25781);
        float r25783 = log1p(r25782);
        return r25783;
}

double f_od(double x) {
        double r25784 = x;
        double r25785 = expm1(r25784);
        double r25786 = r25785 / r25784;
        double r25787 = expm1(r25786);
        double r25788 = log1p(r25787);
        return r25788;
}

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 r25789, r25790, r25791, r25792, r25793;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r25789);
        mpfr_init(r25790);
        mpfr_init_set_str(r25791, "1", 10, MPFR_RNDN);
        mpfr_init(r25792);
        mpfr_init(r25793);
}

double f_im(double x) {
        mpfr_set_d(r25789, x, MPFR_RNDN);
        mpfr_exp(r25790, r25789, MPFR_RNDN);
        ;
        mpfr_sub(r25792, r25790, r25791, MPFR_RNDN);
        mpfr_div(r25793, r25792, r25789, MPFR_RNDN);
        return mpfr_get_d(r25793, MPFR_RNDN);
}

static mpfr_t r25794, r25795, r25796, r25797, r25798;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r25794);
        mpfr_init(r25795);
        mpfr_init(r25796);
        mpfr_init(r25797);
        mpfr_init(r25798);
}

double f_fm(double x) {
        mpfr_set_d(r25794, x, MPFR_RNDN);
        mpfr_expm1(r25795, r25794, MPFR_RNDN);
        mpfr_div(r25796, r25795, r25794, MPFR_RNDN);
        mpfr_expm1(r25797, r25796, MPFR_RNDN);
        mpfr_log1p(r25798, r25797, MPFR_RNDN);
        return mpfr_get_d(r25798, MPFR_RNDN);
}

static mpfr_t r25799, r25800, r25801, r25802, r25803;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r25799);
        mpfr_init(r25800);
        mpfr_init(r25801);
        mpfr_init(r25802);
        mpfr_init(r25803);
}

double f_dm(double x) {
        mpfr_set_d(r25799, x, MPFR_RNDN);
        mpfr_expm1(r25800, r25799, MPFR_RNDN);
        mpfr_div(r25801, r25800, r25799, MPFR_RNDN);
        mpfr_expm1(r25802, r25801, MPFR_RNDN);
        mpfr_log1p(r25803, r25802, MPFR_RNDN);
        return mpfr_get_d(r25803, MPFR_RNDN);
}