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

char *name = "Rosa's DopplerBench";

double f_if(float u, float v, float t1) {
        float r7694 = t1;
        float r7695 = -r7694;
        float r7696 = v;
        float r7697 = r7695 * r7696;
        float r7698 = u;
        float r7699 = r7694 + r7698;
        float r7700 = r7699 * r7699;
        float r7701 = r7697 / r7700;
        return r7701;
}

double f_id(double u, double v, double t1) {
        double r7702 = t1;
        double r7703 = -r7702;
        double r7704 = v;
        double r7705 = r7703 * r7704;
        double r7706 = u;
        double r7707 = r7702 + r7706;
        double r7708 = r7707 * r7707;
        double r7709 = r7705 / r7708;
        return r7709;
}


double f_of(float u, float v, float t1) {
        float r7710 = 1.0f;
        float r7711 = t1;
        float r7712 = u;
        float r7713 = r7711 + r7712;
        float r7714 = r7713 / r7711;
        float r7715 = r7710 / r7714;
        float r7716 = v;
        float r7717 = r7715 * r7716;
        float r7718 = -r7710;
        float r7719 = r7718 / r7713;
        float r7720 = r7717 * r7719;
        return r7720;
}

double f_od(double u, double v, double t1) {
        double r7721 = 1.0;
        double r7722 = t1;
        double r7723 = u;
        double r7724 = r7722 + r7723;
        double r7725 = r7724 / r7722;
        double r7726 = r7721 / r7725;
        double r7727 = v;
        double r7728 = r7726 * r7727;
        double r7729 = -r7721;
        double r7730 = r7729 / r7724;
        double r7731 = r7728 * r7730;
        return r7731;
}

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 r7732, r7733, r7734, r7735, r7736, r7737, r7738, r7739;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r7732);
        mpfr_init(r7733);
        mpfr_init(r7734);
        mpfr_init(r7735);
        mpfr_init(r7736);
        mpfr_init(r7737);
        mpfr_init(r7738);
        mpfr_init(r7739);
}

double f_im(double u, double v, double t1) {
        mpfr_set_d(r7732, t1, MPFR_RNDN);
        mpfr_neg(r7733, r7732, MPFR_RNDN);
        mpfr_set_d(r7734, v, MPFR_RNDN);
        mpfr_mul(r7735, r7733, r7734, MPFR_RNDN);
        mpfr_set_d(r7736, u, MPFR_RNDN);
        mpfr_add(r7737, r7732, r7736, MPFR_RNDN);
        mpfr_mul(r7738, r7737, r7737, MPFR_RNDN);
        mpfr_div(r7739, r7735, r7738, MPFR_RNDN);
        return mpfr_get_d(r7739, MPFR_RNDN);
}

static mpfr_t r7740, r7741, r7742, r7743, r7744, r7745, r7746, r7747, r7748, r7749, r7750;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7740, "1", 10, MPFR_RNDN);
        mpfr_init(r7741);
        mpfr_init(r7742);
        mpfr_init(r7743);
        mpfr_init(r7744);
        mpfr_init(r7745);
        mpfr_init(r7746);
        mpfr_init(r7747);
        mpfr_init(r7748);
        mpfr_init(r7749);
        mpfr_init(r7750);
}

double f_fm(double u, double v, double t1) {
        ;
        mpfr_set_d(r7741, t1, MPFR_RNDN);
        mpfr_set_d(r7742, u, MPFR_RNDN);
        mpfr_add(r7743, r7741, r7742, MPFR_RNDN);
        mpfr_div(r7744, r7743, r7741, MPFR_RNDN);
        mpfr_div(r7745, r7740, r7744, MPFR_RNDN);
        mpfr_set_d(r7746, v, MPFR_RNDN);
        mpfr_mul(r7747, r7745, r7746, MPFR_RNDN);
        mpfr_neg(r7748, r7740, MPFR_RNDN);
        mpfr_div(r7749, r7748, r7743, MPFR_RNDN);
        mpfr_mul(r7750, r7747, r7749, MPFR_RNDN);
        return mpfr_get_d(r7750, MPFR_RNDN);
}

static mpfr_t r7751, r7752, r7753, r7754, r7755, r7756, r7757, r7758, r7759, r7760, r7761;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7751, "1", 10, MPFR_RNDN);
        mpfr_init(r7752);
        mpfr_init(r7753);
        mpfr_init(r7754);
        mpfr_init(r7755);
        mpfr_init(r7756);
        mpfr_init(r7757);
        mpfr_init(r7758);
        mpfr_init(r7759);
        mpfr_init(r7760);
        mpfr_init(r7761);
}

double f_dm(double u, double v, double t1) {
        ;
        mpfr_set_d(r7752, t1, MPFR_RNDN);
        mpfr_set_d(r7753, u, MPFR_RNDN);
        mpfr_add(r7754, r7752, r7753, MPFR_RNDN);
        mpfr_div(r7755, r7754, r7752, MPFR_RNDN);
        mpfr_div(r7756, r7751, r7755, MPFR_RNDN);
        mpfr_set_d(r7757, v, MPFR_RNDN);
        mpfr_mul(r7758, r7756, r7757, MPFR_RNDN);
        mpfr_neg(r7759, r7751, MPFR_RNDN);
        mpfr_div(r7760, r7759, r7754, MPFR_RNDN);
        mpfr_mul(r7761, r7758, r7760, MPFR_RNDN);
        return mpfr_get_d(r7761, MPFR_RNDN);
}