#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 r20895 = t1;
        float r20896 = -r20895;
        float r20897 = v;
        float r20898 = r20896 * r20897;
        float r20899 = u;
        float r20900 = r20895 + r20899;
        float r20901 = r20900 * r20900;
        float r20902 = r20898 / r20901;
        return r20902;
}

double f_id(double u, double v, double t1) {
        double r20903 = t1;
        double r20904 = -r20903;
        double r20905 = v;
        double r20906 = r20904 * r20905;
        double r20907 = u;
        double r20908 = r20903 + r20907;
        double r20909 = r20908 * r20908;
        double r20910 = r20906 / r20909;
        return r20910;
}


double f_of(float u, float v, float t1) {
        float r20911 = t1;
        float r20912 = -r20911;
        float r20913 = u;
        float r20914 = r20911 + r20913;
        float r20915 = r20912 / r20914;
        float r20916 = v;
        float r20917 = r20916 / r20914;
        float r20918 = r20915 * r20917;
        float r20919 = -7.159429525164284e-244;
        bool r20920 = r20918 <= r20919;
        float r20921 = r20915 * r20916;
        float r20922 = 1;
        float r20923 = r20922 / r20914;
        float r20924 = r20921 * r20923;
        float r20925 = -0.0;
        bool r20926 = r20918 <= r20925;
        float r20927 = r20912 * r20916;
        float r20928 = r20927 / r20914;
        float r20929 = r20928 / r20914;
        float r20930 = r20926 ? r20929 : r20918;
        float r20931 = r20920 ? r20924 : r20930;
        return r20931;
}

double f_od(double u, double v, double t1) {
        double r20932 = t1;
        double r20933 = -r20932;
        double r20934 = u;
        double r20935 = r20932 + r20934;
        double r20936 = r20933 / r20935;
        double r20937 = v;
        double r20938 = r20937 / r20935;
        double r20939 = r20936 * r20938;
        double r20940 = -7.159429525164284e-244;
        bool r20941 = r20939 <= r20940;
        double r20942 = r20936 * r20937;
        double r20943 = 1;
        double r20944 = r20943 / r20935;
        double r20945 = r20942 * r20944;
        double r20946 = -0.0;
        bool r20947 = r20939 <= r20946;
        double r20948 = r20933 * r20937;
        double r20949 = r20948 / r20935;
        double r20950 = r20949 / r20935;
        double r20951 = r20947 ? r20950 : r20939;
        double r20952 = r20941 ? r20945 : r20951;
        return r20952;
}

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 r20953, r20954, r20955, r20956, r20957, r20958, r20959, r20960;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20953);
        mpfr_init(r20954);
        mpfr_init(r20955);
        mpfr_init(r20956);
        mpfr_init(r20957);
        mpfr_init(r20958);
        mpfr_init(r20959);
        mpfr_init(r20960);
}

double f_im(double u, double v, double t1) {
        mpfr_set_d(r20953, t1, MPFR_RNDN);
        mpfr_neg(r20954, r20953, MPFR_RNDN);
        mpfr_set_d(r20955, v, MPFR_RNDN);
        mpfr_mul(r20956, r20954, r20955, MPFR_RNDN);
        mpfr_set_d(r20957, u, MPFR_RNDN);
        mpfr_add(r20958, r20953, r20957, MPFR_RNDN);
        mpfr_mul(r20959, r20958, r20958, MPFR_RNDN);
        mpfr_div(r20960, r20956, r20959, MPFR_RNDN);
        return mpfr_get_d(r20960, MPFR_RNDN);
}

static mpfr_t r20961, r20962, r20963, r20964, r20965, r20966, r20967, r20968, r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976, r20977, r20978, r20979, r20980, r20981;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20961);
        mpfr_init(r20962);
        mpfr_init(r20963);
        mpfr_init(r20964);
        mpfr_init(r20965);
        mpfr_init(r20966);
        mpfr_init(r20967);
        mpfr_init(r20968);
        mpfr_init_set_str(r20969, "-7.159429525164284e-244", 10, MPFR_RNDN);
        mpfr_init(r20970);
        mpfr_init(r20971);
        mpfr_init_set_str(r20972, "1", 10, MPFR_RNDN);
        mpfr_init(r20973);
        mpfr_init(r20974);
        mpfr_init_set_str(r20975, "-0.0", 10, MPFR_RNDN);
        mpfr_init(r20976);
        mpfr_init(r20977);
        mpfr_init(r20978);
        mpfr_init(r20979);
        mpfr_init(r20980);
        mpfr_init(r20981);
}

double f_fm(double u, double v, double t1) {
        mpfr_set_d(r20961, t1, MPFR_RNDN);
        mpfr_neg(r20962, r20961, MPFR_RNDN);
        mpfr_set_d(r20963, u, MPFR_RNDN);
        mpfr_add(r20964, r20961, r20963, MPFR_RNDN);
        mpfr_div(r20965, r20962, r20964, MPFR_RNDN);
        mpfr_set_d(r20966, v, MPFR_RNDN);
        mpfr_div(r20967, r20966, r20964, MPFR_RNDN);
        mpfr_mul(r20968, r20965, r20967, MPFR_RNDN);
        ;
        mpfr_set_si(r20970, mpfr_cmp(r20968, r20969) <= 0, MPFR_RNDN);
        mpfr_mul(r20971, r20965, r20966, MPFR_RNDN);
        ;
        mpfr_div(r20973, r20972, r20964, MPFR_RNDN);
        mpfr_mul(r20974, r20971, r20973, MPFR_RNDN);
        ;
        mpfr_set_si(r20976, mpfr_cmp(r20968, r20975) <= 0, MPFR_RNDN);
        mpfr_mul(r20977, r20962, r20966, MPFR_RNDN);
        mpfr_div(r20978, r20977, r20964, MPFR_RNDN);
        mpfr_div(r20979, r20978, r20964, MPFR_RNDN);
        if (mpfr_get_si(r20976, MPFR_RNDN)) { mpfr_set(r20980, r20979, MPFR_RNDN); } else { mpfr_set(r20980, r20968, MPFR_RNDN); };
        if (mpfr_get_si(r20970, MPFR_RNDN)) { mpfr_set(r20981, r20974, MPFR_RNDN); } else { mpfr_set(r20981, r20980, MPFR_RNDN); };
        return mpfr_get_d(r20981, MPFR_RNDN);
}

static mpfr_t r20982, r20983, r20984, r20985, r20986, r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994, r20995, r20996, r20997, r20998, r20999, r21000, r21001, r21002;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20982);
        mpfr_init(r20983);
        mpfr_init(r20984);
        mpfr_init(r20985);
        mpfr_init(r20986);
        mpfr_init(r20987);
        mpfr_init(r20988);
        mpfr_init(r20989);
        mpfr_init_set_str(r20990, "-7.159429525164284e-244", 10, MPFR_RNDN);
        mpfr_init(r20991);
        mpfr_init(r20992);
        mpfr_init_set_str(r20993, "1", 10, MPFR_RNDN);
        mpfr_init(r20994);
        mpfr_init(r20995);
        mpfr_init_set_str(r20996, "-0.0", 10, MPFR_RNDN);
        mpfr_init(r20997);
        mpfr_init(r20998);
        mpfr_init(r20999);
        mpfr_init(r21000);
        mpfr_init(r21001);
        mpfr_init(r21002);
}

double f_dm(double u, double v, double t1) {
        mpfr_set_d(r20982, t1, MPFR_RNDN);
        mpfr_neg(r20983, r20982, MPFR_RNDN);
        mpfr_set_d(r20984, u, MPFR_RNDN);
        mpfr_add(r20985, r20982, r20984, MPFR_RNDN);
        mpfr_div(r20986, r20983, r20985, MPFR_RNDN);
        mpfr_set_d(r20987, v, MPFR_RNDN);
        mpfr_div(r20988, r20987, r20985, MPFR_RNDN);
        mpfr_mul(r20989, r20986, r20988, MPFR_RNDN);
        ;
        mpfr_set_si(r20991, mpfr_cmp(r20989, r20990) <= 0, MPFR_RNDN);
        mpfr_mul(r20992, r20986, r20987, MPFR_RNDN);
        ;
        mpfr_div(r20994, r20993, r20985, MPFR_RNDN);
        mpfr_mul(r20995, r20992, r20994, MPFR_RNDN);
        ;
        mpfr_set_si(r20997, mpfr_cmp(r20989, r20996) <= 0, MPFR_RNDN);
        mpfr_mul(r20998, r20983, r20987, MPFR_RNDN);
        mpfr_div(r20999, r20998, r20985, MPFR_RNDN);
        mpfr_div(r21000, r20999, r20985, MPFR_RNDN);
        if (mpfr_get_si(r20997, MPFR_RNDN)) { mpfr_set(r21001, r21000, MPFR_RNDN); } else { mpfr_set(r21001, r20989, MPFR_RNDN); };
        if (mpfr_get_si(r20991, MPFR_RNDN)) { mpfr_set(r21002, r20995, MPFR_RNDN); } else { mpfr_set(r21002, r21001, MPFR_RNDN); };
        return mpfr_get_d(r21002, MPFR_RNDN);
}