#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 r24747 = t1;
        float r24748 = -r24747;
        float r24749 = v;
        float r24750 = r24748 * r24749;
        float r24751 = u;
        float r24752 = r24747 + r24751;
        float r24753 = r24752 * r24752;
        float r24754 = r24750 / r24753;
        return r24754;
}

double f_id(double u, double v, double t1) {
        double r24755 = t1;
        double r24756 = -r24755;
        double r24757 = v;
        double r24758 = r24756 * r24757;
        double r24759 = u;
        double r24760 = r24755 + r24759;
        double r24761 = r24760 * r24760;
        double r24762 = r24758 / r24761;
        return r24762;
}


double f_of(float u, float v, float t1) {
        float r24763 = t1;
        float r24764 = -r24763;
        float r24765 = cbrt(r24764);
        float r24766 = r24765 * r24765;
        float r24767 = u;
        float r24768 = r24763 + r24767;
        float r24769 = cbrt(r24768);
        float r24770 = r24769 * r24769;
        float r24771 = r24766 / r24770;
        float r24772 = r24765 / r24769;
        float r24773 = v;
        float r24774 = r24773 / r24768;
        float r24775 = r24772 * r24774;
        float r24776 = r24771 * r24775;
        return r24776;
}

double f_od(double u, double v, double t1) {
        double r24777 = t1;
        double r24778 = -r24777;
        double r24779 = cbrt(r24778);
        double r24780 = r24779 * r24779;
        double r24781 = u;
        double r24782 = r24777 + r24781;
        double r24783 = cbrt(r24782);
        double r24784 = r24783 * r24783;
        double r24785 = r24780 / r24784;
        double r24786 = r24779 / r24783;
        double r24787 = v;
        double r24788 = r24787 / r24782;
        double r24789 = r24786 * r24788;
        double r24790 = r24785 * r24789;
        return r24790;
}

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 r24791, r24792, r24793, r24794, r24795, r24796, r24797, r24798;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r24791);
        mpfr_init(r24792);
        mpfr_init(r24793);
        mpfr_init(r24794);
        mpfr_init(r24795);
        mpfr_init(r24796);
        mpfr_init(r24797);
        mpfr_init(r24798);
}

double f_im(double u, double v, double t1) {
        mpfr_set_d(r24791, t1, MPFR_RNDN);
        mpfr_neg(r24792, r24791, MPFR_RNDN);
        mpfr_set_d(r24793, v, MPFR_RNDN);
        mpfr_mul(r24794, r24792, r24793, MPFR_RNDN);
        mpfr_set_d(r24795, u, MPFR_RNDN);
        mpfr_add(r24796, r24791, r24795, MPFR_RNDN);
        mpfr_mul(r24797, r24796, r24796, MPFR_RNDN);
        mpfr_div(r24798, r24794, r24797, MPFR_RNDN);
        return mpfr_get_d(r24798, MPFR_RNDN);
}

static mpfr_t r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r24799);
        mpfr_init(r24800);
        mpfr_init(r24801);
        mpfr_init(r24802);
        mpfr_init(r24803);
        mpfr_init(r24804);
        mpfr_init(r24805);
        mpfr_init(r24806);
        mpfr_init(r24807);
        mpfr_init(r24808);
        mpfr_init(r24809);
        mpfr_init(r24810);
        mpfr_init(r24811);
        mpfr_init(r24812);
}

double f_fm(double u, double v, double t1) {
        mpfr_set_d(r24799, t1, MPFR_RNDN);
        mpfr_neg(r24800, r24799, MPFR_RNDN);
        mpfr_cbrt(r24801, r24800, MPFR_RNDN);
        mpfr_mul(r24802, r24801, r24801, MPFR_RNDN);
        mpfr_set_d(r24803, u, MPFR_RNDN);
        mpfr_add(r24804, r24799, r24803, MPFR_RNDN);
        mpfr_cbrt(r24805, r24804, MPFR_RNDN);
        mpfr_mul(r24806, r24805, r24805, MPFR_RNDN);
        mpfr_div(r24807, r24802, r24806, MPFR_RNDN);
        mpfr_div(r24808, r24801, r24805, MPFR_RNDN);
        mpfr_set_d(r24809, v, MPFR_RNDN);
        mpfr_div(r24810, r24809, r24804, MPFR_RNDN);
        mpfr_mul(r24811, r24808, r24810, MPFR_RNDN);
        mpfr_mul(r24812, r24807, r24811, MPFR_RNDN);
        return mpfr_get_d(r24812, MPFR_RNDN);
}

static mpfr_t r24813, r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r24813);
        mpfr_init(r24814);
        mpfr_init(r24815);
        mpfr_init(r24816);
        mpfr_init(r24817);
        mpfr_init(r24818);
        mpfr_init(r24819);
        mpfr_init(r24820);
        mpfr_init(r24821);
        mpfr_init(r24822);
        mpfr_init(r24823);
        mpfr_init(r24824);
        mpfr_init(r24825);
        mpfr_init(r24826);
}

double f_dm(double u, double v, double t1) {
        mpfr_set_d(r24813, t1, MPFR_RNDN);
        mpfr_neg(r24814, r24813, MPFR_RNDN);
        mpfr_cbrt(r24815, r24814, MPFR_RNDN);
        mpfr_mul(r24816, r24815, r24815, MPFR_RNDN);
        mpfr_set_d(r24817, u, MPFR_RNDN);
        mpfr_add(r24818, r24813, r24817, MPFR_RNDN);
        mpfr_cbrt(r24819, r24818, MPFR_RNDN);
        mpfr_mul(r24820, r24819, r24819, MPFR_RNDN);
        mpfr_div(r24821, r24816, r24820, MPFR_RNDN);
        mpfr_div(r24822, r24815, r24819, MPFR_RNDN);
        mpfr_set_d(r24823, v, MPFR_RNDN);
        mpfr_div(r24824, r24823, r24818, MPFR_RNDN);
        mpfr_mul(r24825, r24822, r24824, MPFR_RNDN);
        mpfr_mul(r24826, r24821, r24825, MPFR_RNDN);
        return mpfr_get_d(r24826, MPFR_RNDN);
}