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

char *name = "Rosa's TurbineBenchmark";

double f_if(float v, float w, float r) {
        float r7853 = 3;
        float r7854 = 2;
        float r7855 = r;
        float r7856 = r7855 * r7855;
        float r7857 = r7854 / r7856;
        float r7858 = r7853 + r7857;
        float r7859 = 0.125;
        float r7860 = v;
        float r7861 = r7854 * r7860;
        float r7862 = r7853 - r7861;
        float r7863 = r7859 * r7862;
        float r7864 = w;
        float r7865 = r7864 * r7864;
        float r7866 = r7865 * r7855;
        float r7867 = r7866 * r7855;
        float r7868 = r7863 * r7867;
        float r7869 = 1;
        float r7870 = r7869 - r7860;
        float r7871 = r7868 / r7870;
        float r7872 = r7858 - r7871;
        float r7873 = 4.5;
        float r7874 = r7872 - r7873;
        return r7874;
}

double f_id(double v, double w, double r) {
        double r7875 = 3;
        double r7876 = 2;
        double r7877 = r;
        double r7878 = r7877 * r7877;
        double r7879 = r7876 / r7878;
        double r7880 = r7875 + r7879;
        double r7881 = 0.125;
        double r7882 = v;
        double r7883 = r7876 * r7882;
        double r7884 = r7875 - r7883;
        double r7885 = r7881 * r7884;
        double r7886 = w;
        double r7887 = r7886 * r7886;
        double r7888 = r7887 * r7877;
        double r7889 = r7888 * r7877;
        double r7890 = r7885 * r7889;
        double r7891 = 1;
        double r7892 = r7891 - r7882;
        double r7893 = r7890 / r7892;
        double r7894 = r7880 - r7893;
        double r7895 = 4.5;
        double r7896 = r7894 - r7895;
        return r7896;
}


double f_of(float v, float w, float r) {
        float r7897 = 3;
        float r7898 = 2;
        float r7899 = r;
        float r7900 = r7899 * r7899;
        float r7901 = r7898 / r7900;
        float r7902 = r7897 + r7901;
        float r7903 = 0.125;
        float r7904 = v;
        float r7905 = r7898 * r7904;
        float r7906 = r7897 - r7905;
        float r7907 = r7903 * r7906;
        float r7908 = 1;
        float r7909 = r7908 - r7904;
        float r7910 = cbrt(r7909);
        float r7911 = r7910 * r7910;
        float r7912 = r7907 / r7911;
        float r7913 = w;
        float r7914 = r7913 * r7899;
        float r7915 = fabs(r7914);
        float r7916 = r7915 * r7915;
        float r7917 = r7916 / r7910;
        float r7918 = r7912 * r7917;
        float r7919 = r7902 - r7918;
        float r7920 = 4.5;
        float r7921 = r7919 - r7920;
        return r7921;
}

double f_od(double v, double w, double r) {
        double r7922 = 3;
        double r7923 = 2;
        double r7924 = r;
        double r7925 = r7924 * r7924;
        double r7926 = r7923 / r7925;
        double r7927 = r7922 + r7926;
        double r7928 = 0.125;
        double r7929 = v;
        double r7930 = r7923 * r7929;
        double r7931 = r7922 - r7930;
        double r7932 = r7928 * r7931;
        double r7933 = 1;
        double r7934 = r7933 - r7929;
        double r7935 = cbrt(r7934);
        double r7936 = r7935 * r7935;
        double r7937 = r7932 / r7936;
        double r7938 = w;
        double r7939 = r7938 * r7924;
        double r7940 = fabs(r7939);
        double r7941 = r7940 * r7940;
        double r7942 = r7941 / r7935;
        double r7943 = r7937 * r7942;
        double r7944 = r7927 - r7943;
        double r7945 = 4.5;
        double r7946 = r7944 - r7945;
        return r7946;
}

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 r7947, r7948, r7949, r7950, r7951, r7952, r7953, r7954, r7955, r7956, r7957, r7958, r7959, r7960, r7961, r7962, r7963, r7964, r7965, r7966, r7967, r7968;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7947, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r7948, "2", 10, MPFR_RNDN);
        mpfr_init(r7949);
        mpfr_init(r7950);
        mpfr_init(r7951);
        mpfr_init(r7952);
        mpfr_init_set_str(r7953, "0.125", 10, MPFR_RNDN);
        mpfr_init(r7954);
        mpfr_init(r7955);
        mpfr_init(r7956);
        mpfr_init(r7957);
        mpfr_init(r7958);
        mpfr_init(r7959);
        mpfr_init(r7960);
        mpfr_init(r7961);
        mpfr_init(r7962);
        mpfr_init_set_str(r7963, "1", 10, MPFR_RNDN);
        mpfr_init(r7964);
        mpfr_init(r7965);
        mpfr_init(r7966);
        mpfr_init_set_str(r7967, "4.5", 10, MPFR_RNDN);
        mpfr_init(r7968);
}

double f_im(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r7949, r, MPFR_RNDN);
        mpfr_mul(r7950, r7949, r7949, MPFR_RNDN);
        mpfr_div(r7951, r7948, r7950, MPFR_RNDN);
        mpfr_add(r7952, r7947, r7951, MPFR_RNDN);
        ;
        mpfr_set_d(r7954, v, MPFR_RNDN);
        mpfr_mul(r7955, r7948, r7954, MPFR_RNDN);
        mpfr_sub(r7956, r7947, r7955, MPFR_RNDN);
        mpfr_mul(r7957, r7953, r7956, MPFR_RNDN);
        mpfr_set_d(r7958, w, MPFR_RNDN);
        mpfr_mul(r7959, r7958, r7958, MPFR_RNDN);
        mpfr_mul(r7960, r7959, r7949, MPFR_RNDN);
        mpfr_mul(r7961, r7960, r7949, MPFR_RNDN);
        mpfr_mul(r7962, r7957, r7961, MPFR_RNDN);
        ;
        mpfr_sub(r7964, r7963, r7954, MPFR_RNDN);
        mpfr_div(r7965, r7962, r7964, MPFR_RNDN);
        mpfr_sub(r7966, r7952, r7965, MPFR_RNDN);
        ;
        mpfr_sub(r7968, r7966, r7967, MPFR_RNDN);
        return mpfr_get_d(r7968, MPFR_RNDN);
}

static mpfr_t r7969, r7970, r7971, r7972, r7973, r7974, r7975, r7976, r7977, r7978, r7979, r7980, r7981, r7982, r7983, r7984, r7985, r7986, r7987, r7988, r7989, r7990, r7991, r7992, r7993;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7969, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r7970, "2", 10, MPFR_RNDN);
        mpfr_init(r7971);
        mpfr_init(r7972);
        mpfr_init(r7973);
        mpfr_init(r7974);
        mpfr_init_set_str(r7975, "0.125", 10, MPFR_RNDN);
        mpfr_init(r7976);
        mpfr_init(r7977);
        mpfr_init(r7978);
        mpfr_init(r7979);
        mpfr_init_set_str(r7980, "1", 10, MPFR_RNDN);
        mpfr_init(r7981);
        mpfr_init(r7982);
        mpfr_init(r7983);
        mpfr_init(r7984);
        mpfr_init(r7985);
        mpfr_init(r7986);
        mpfr_init(r7987);
        mpfr_init(r7988);
        mpfr_init(r7989);
        mpfr_init(r7990);
        mpfr_init(r7991);
        mpfr_init_set_str(r7992, "4.5", 10, MPFR_RNDN);
        mpfr_init(r7993);
}

double f_fm(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r7971, r, MPFR_RNDN);
        mpfr_mul(r7972, r7971, r7971, MPFR_RNDN);
        mpfr_div(r7973, r7970, r7972, MPFR_RNDN);
        mpfr_add(r7974, r7969, r7973, MPFR_RNDN);
        ;
        mpfr_set_d(r7976, v, MPFR_RNDN);
        mpfr_mul(r7977, r7970, r7976, MPFR_RNDN);
        mpfr_sub(r7978, r7969, r7977, MPFR_RNDN);
        mpfr_mul(r7979, r7975, r7978, MPFR_RNDN);
        ;
        mpfr_sub(r7981, r7980, r7976, MPFR_RNDN);
        mpfr_cbrt(r7982, r7981, MPFR_RNDN);
        mpfr_mul(r7983, r7982, r7982, MPFR_RNDN);
        mpfr_div(r7984, r7979, r7983, MPFR_RNDN);
        mpfr_set_d(r7985, w, MPFR_RNDN);
        mpfr_mul(r7986, r7985, r7971, MPFR_RNDN);
        mpfr_abs(r7987, r7986, MPFR_RNDN);
        mpfr_mul(r7988, r7987, r7987, MPFR_RNDN);
        mpfr_div(r7989, r7988, r7982, MPFR_RNDN);
        mpfr_mul(r7990, r7984, r7989, MPFR_RNDN);
        mpfr_sub(r7991, r7974, r7990, MPFR_RNDN);
        ;
        mpfr_sub(r7993, r7991, r7992, MPFR_RNDN);
        return mpfr_get_d(r7993, MPFR_RNDN);
}

static mpfr_t r7994, r7995, r7996, r7997, r7998, r7999, r8000, r8001, r8002, r8003, r8004, r8005, r8006, r8007, r8008, r8009, r8010, r8011, r8012, r8013, r8014, r8015, r8016, r8017, r8018;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7994, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r7995, "2", 10, MPFR_RNDN);
        mpfr_init(r7996);
        mpfr_init(r7997);
        mpfr_init(r7998);
        mpfr_init(r7999);
        mpfr_init_set_str(r8000, "0.125", 10, MPFR_RNDN);
        mpfr_init(r8001);
        mpfr_init(r8002);
        mpfr_init(r8003);
        mpfr_init(r8004);
        mpfr_init_set_str(r8005, "1", 10, MPFR_RNDN);
        mpfr_init(r8006);
        mpfr_init(r8007);
        mpfr_init(r8008);
        mpfr_init(r8009);
        mpfr_init(r8010);
        mpfr_init(r8011);
        mpfr_init(r8012);
        mpfr_init(r8013);
        mpfr_init(r8014);
        mpfr_init(r8015);
        mpfr_init(r8016);
        mpfr_init_set_str(r8017, "4.5", 10, MPFR_RNDN);
        mpfr_init(r8018);
}

double f_dm(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r7996, r, MPFR_RNDN);
        mpfr_mul(r7997, r7996, r7996, MPFR_RNDN);
        mpfr_div(r7998, r7995, r7997, MPFR_RNDN);
        mpfr_add(r7999, r7994, r7998, MPFR_RNDN);
        ;
        mpfr_set_d(r8001, v, MPFR_RNDN);
        mpfr_mul(r8002, r7995, r8001, MPFR_RNDN);
        mpfr_sub(r8003, r7994, r8002, MPFR_RNDN);
        mpfr_mul(r8004, r8000, r8003, MPFR_RNDN);
        ;
        mpfr_sub(r8006, r8005, r8001, MPFR_RNDN);
        mpfr_cbrt(r8007, r8006, MPFR_RNDN);
        mpfr_mul(r8008, r8007, r8007, MPFR_RNDN);
        mpfr_div(r8009, r8004, r8008, MPFR_RNDN);
        mpfr_set_d(r8010, w, MPFR_RNDN);
        mpfr_mul(r8011, r8010, r7996, MPFR_RNDN);
        mpfr_abs(r8012, r8011, MPFR_RNDN);
        mpfr_mul(r8013, r8012, r8012, MPFR_RNDN);
        mpfr_div(r8014, r8013, r8007, MPFR_RNDN);
        mpfr_mul(r8015, r8009, r8014, MPFR_RNDN);
        mpfr_sub(r8016, r7999, r8015, MPFR_RNDN);
        ;
        mpfr_sub(r8018, r8016, r8017, MPFR_RNDN);
        return mpfr_get_d(r8018, MPFR_RNDN);
}

