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

char *name = "Diagrams.Trail:splitAtParam  from diagrams-lib-1.3.0.3, C";

double f_if(float x, float y) {
        float r33832 = x;
        float r33833 = y;
        float r33834 = r33832 - r33833;
        float r33835 = 1.0;
        float r33836 = r33835 - r33833;
        float r33837 = r33834 / r33836;
        return r33837;
}

double f_id(double x, double y) {
        double r33838 = x;
        double r33839 = y;
        double r33840 = r33838 - r33839;
        double r33841 = 1.0;
        double r33842 = r33841 - r33839;
        double r33843 = r33840 / r33842;
        return r33843;
}


double f_of(float x, float y) {
        float r33844 = x;
        float r33845 = y;
        float r33846 = r33844 - r33845;
        float r33847 = 1.0;
        float r33848 = r33847 - r33845;
        float r33849 = r33846 / r33848;
        return r33849;
}

double f_od(double x, double y) {
        double r33850 = x;
        double r33851 = y;
        double r33852 = r33850 - r33851;
        double r33853 = 1.0;
        double r33854 = r33853 - r33851;
        double r33855 = r33852 / r33854;
        return r33855;
}

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 r33856, r33857, r33858, r33859, r33860, r33861;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r33856);
        mpfr_init(r33857);
        mpfr_init(r33858);
        mpfr_init_set_str(r33859, "1.0", 10, MPFR_RNDN);
        mpfr_init(r33860);
        mpfr_init(r33861);
}

double f_im(double x, double y) {
        mpfr_set_d(r33856, x, MPFR_RNDN);
        mpfr_set_d(r33857, y, MPFR_RNDN);
        mpfr_sub(r33858, r33856, r33857, MPFR_RNDN);
        ;
        mpfr_sub(r33860, r33859, r33857, MPFR_RNDN);
        mpfr_div(r33861, r33858, r33860, MPFR_RNDN);
        return mpfr_get_d(r33861, MPFR_RNDN);
}

static mpfr_t r33862, r33863, r33864, r33865, r33866, r33867;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r33862);
        mpfr_init(r33863);
        mpfr_init(r33864);
        mpfr_init_set_str(r33865, "1.0", 10, MPFR_RNDN);
        mpfr_init(r33866);
        mpfr_init(r33867);
}

double f_fm(double x, double y) {
        mpfr_set_d(r33862, x, MPFR_RNDN);
        mpfr_set_d(r33863, y, MPFR_RNDN);
        mpfr_sub(r33864, r33862, r33863, MPFR_RNDN);
        ;
        mpfr_sub(r33866, r33865, r33863, MPFR_RNDN);
        mpfr_div(r33867, r33864, r33866, MPFR_RNDN);
        return mpfr_get_d(r33867, MPFR_RNDN);
}

static mpfr_t r33868, r33869, r33870, r33871, r33872, r33873;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r33868);
        mpfr_init(r33869);
        mpfr_init(r33870);
        mpfr_init_set_str(r33871, "1.0", 10, MPFR_RNDN);
        mpfr_init(r33872);
        mpfr_init(r33873);
}

double f_dm(double x, double y) {
        mpfr_set_d(r33868, x, MPFR_RNDN);
        mpfr_set_d(r33869, y, MPFR_RNDN);
        mpfr_sub(r33870, r33868, r33869, MPFR_RNDN);
        ;
        mpfr_sub(r33872, r33871, r33869, MPFR_RNDN);
        mpfr_div(r33873, r33870, r33872, MPFR_RNDN);
        return mpfr_get_d(r33873, MPFR_RNDN);
}