#include #include #include #include #include char *name = "Rosa's TurbineBenchmark"; double f_if(float v, float w, float r) { float r20996 = 3; float r20997 = 2; float r20998 = r; float r20999 = r20998 * r20998; float r21000 = r20997 / r20999; float r21001 = r20996 + r21000; float r21002 = 0.125; float r21003 = v; float r21004 = r20997 * r21003; float r21005 = r20996 - r21004; float r21006 = r21002 * r21005; float r21007 = w; float r21008 = r21007 * r21007; float r21009 = r21008 * r20998; float r21010 = r21009 * r20998; float r21011 = r21006 * r21010; float r21012 = 1; float r21013 = r21012 - r21003; float r21014 = r21011 / r21013; float r21015 = r21001 - r21014; float r21016 = 4.5; float r21017 = r21015 - r21016; return r21017; } double f_id(double v, double w, double r) { double r21018 = 3; double r21019 = 2; double r21020 = r; double r21021 = r21020 * r21020; double r21022 = r21019 / r21021; double r21023 = r21018 + r21022; double r21024 = 0.125; double r21025 = v; double r21026 = r21019 * r21025; double r21027 = r21018 - r21026; double r21028 = r21024 * r21027; double r21029 = w; double r21030 = r21029 * r21029; double r21031 = r21030 * r21020; double r21032 = r21031 * r21020; double r21033 = r21028 * r21032; double r21034 = 1; double r21035 = r21034 - r21025; double r21036 = r21033 / r21035; double r21037 = r21023 - r21036; double r21038 = 4.5; double r21039 = r21037 - r21038; return r21039; } double f_of(float v, float w, float r) { float r21040 = 2; float r21041 = r; float r21042 = r21041 * r21041; float r21043 = r21040 / r21042; float r21044 = 3; float r21045 = v; float r21046 = r21045 * r21040; float r21047 = r21044 - r21046; float r21048 = 1; float r21049 = r21048 - r21045; float r21050 = 0.125; float r21051 = r21049 / r21050; float r21052 = r21047 / r21051; float r21053 = w; float r21054 = r21053 * r21041; float r21055 = r21054 * r21054; float r21056 = 4.5; float r21057 = fma(r21052, r21055, r21056); float r21058 = r21044 - r21057; float r21059 = r21043 + r21058; return r21059; } double f_od(double v, double w, double r) { double r21060 = 2; double r21061 = r; double r21062 = r21061 * r21061; double r21063 = r21060 / r21062; double r21064 = 3; double r21065 = v; double r21066 = r21065 * r21060; double r21067 = r21064 - r21066; double r21068 = 1; double r21069 = r21068 - r21065; double r21070 = 0.125; double r21071 = r21069 / r21070; double r21072 = r21067 / r21071; double r21073 = w; double r21074 = r21073 * r21061; double r21075 = r21074 * r21074; double r21076 = 4.5; double r21077 = fma(r21072, r21075, r21076); double r21078 = r21064 - r21077; double r21079 = r21063 + r21078; return r21079; } 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 r21080, r21081, r21082, r21083, r21084, r21085, r21086, r21087, r21088, r21089, r21090, r21091, r21092, r21093, r21094, r21095, r21096, r21097, r21098, r21099, r21100, r21101; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21080, "3", 10, MPFR_RNDN); mpfr_init_set_str(r21081, "2", 10, MPFR_RNDN); mpfr_init(r21082); mpfr_init(r21083); mpfr_init(r21084); mpfr_init(r21085); mpfr_init_set_str(r21086, "0.125", 10, MPFR_RNDN); mpfr_init(r21087); mpfr_init(r21088); mpfr_init(r21089); mpfr_init(r21090); mpfr_init(r21091); mpfr_init(r21092); mpfr_init(r21093); mpfr_init(r21094); mpfr_init(r21095); mpfr_init_set_str(r21096, "1", 10, MPFR_RNDN); mpfr_init(r21097); mpfr_init(r21098); mpfr_init(r21099); mpfr_init_set_str(r21100, "4.5", 10, MPFR_RNDN); mpfr_init(r21101); } double f_im(double v, double w, double r) { ; ; mpfr_set_d(r21082, r, MPFR_RNDN); mpfr_mul(r21083, r21082, r21082, MPFR_RNDN); mpfr_div(r21084, r21081, r21083, MPFR_RNDN); mpfr_add(r21085, r21080, r21084, MPFR_RNDN); ; mpfr_set_d(r21087, v, MPFR_RNDN); mpfr_mul(r21088, r21081, r21087, MPFR_RNDN); mpfr_sub(r21089, r21080, r21088, MPFR_RNDN); mpfr_mul(r21090, r21086, r21089, MPFR_RNDN); mpfr_set_d(r21091, w, MPFR_RNDN); mpfr_mul(r21092, r21091, r21091, MPFR_RNDN); mpfr_mul(r21093, r21092, r21082, MPFR_RNDN); mpfr_mul(r21094, r21093, r21082, MPFR_RNDN); mpfr_mul(r21095, r21090, r21094, MPFR_RNDN); ; mpfr_sub(r21097, r21096, r21087, MPFR_RNDN); mpfr_div(r21098, r21095, r21097, MPFR_RNDN); mpfr_sub(r21099, r21085, r21098, MPFR_RNDN); ; mpfr_sub(r21101, r21099, r21100, MPFR_RNDN); return mpfr_get_d(r21101, MPFR_RNDN); } static mpfr_t r21102, r21103, r21104, r21105, r21106, r21107, r21108, r21109, r21110, r21111, r21112, r21113, r21114, r21115, r21116, r21117, r21118, r21119, r21120, r21121; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21102, "2", 10, MPFR_RNDN); mpfr_init(r21103); mpfr_init(r21104); mpfr_init(r21105); mpfr_init_set_str(r21106, "3", 10, MPFR_RNDN); mpfr_init(r21107); mpfr_init(r21108); mpfr_init(r21109); mpfr_init_set_str(r21110, "1", 10, MPFR_RNDN); mpfr_init(r21111); mpfr_init_set_str(r21112, "0.125", 10, MPFR_RNDN); mpfr_init(r21113); mpfr_init(r21114); mpfr_init(r21115); mpfr_init(r21116); mpfr_init(r21117); mpfr_init_set_str(r21118, "4.5", 10, MPFR_RNDN); mpfr_init(r21119); mpfr_init(r21120); mpfr_init(r21121); } double f_fm(double v, double w, double r) { ; mpfr_set_d(r21103, r, MPFR_RNDN); mpfr_mul(r21104, r21103, r21103, MPFR_RNDN); mpfr_div(r21105, r21102, r21104, MPFR_RNDN); ; mpfr_set_d(r21107, v, MPFR_RNDN); mpfr_mul(r21108, r21107, r21102, MPFR_RNDN); mpfr_sub(r21109, r21106, r21108, MPFR_RNDN); ; mpfr_sub(r21111, r21110, r21107, MPFR_RNDN); ; mpfr_div(r21113, r21111, r21112, MPFR_RNDN); mpfr_div(r21114, r21109, r21113, MPFR_RNDN); mpfr_set_d(r21115, w, MPFR_RNDN); mpfr_mul(r21116, r21115, r21103, MPFR_RNDN); mpfr_mul(r21117, r21116, r21116, MPFR_RNDN); ; mpfr_fma(r21119, r21114, r21117, r21118, MPFR_RNDN); mpfr_sub(r21120, r21106, r21119, MPFR_RNDN); mpfr_add(r21121, r21105, r21120, MPFR_RNDN); return mpfr_get_d(r21121, MPFR_RNDN); } static mpfr_t r21122, r21123, r21124, r21125, r21126, r21127, r21128, r21129, r21130, r21131, r21132, r21133, r21134, r21135, r21136, r21137, r21138, r21139, r21140, r21141; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21122, "2", 10, MPFR_RNDN); mpfr_init(r21123); mpfr_init(r21124); mpfr_init(r21125); mpfr_init_set_str(r21126, "3", 10, MPFR_RNDN); mpfr_init(r21127); mpfr_init(r21128); mpfr_init(r21129); mpfr_init_set_str(r21130, "1", 10, MPFR_RNDN); mpfr_init(r21131); mpfr_init_set_str(r21132, "0.125", 10, MPFR_RNDN); mpfr_init(r21133); mpfr_init(r21134); mpfr_init(r21135); mpfr_init(r21136); mpfr_init(r21137); mpfr_init_set_str(r21138, "4.5", 10, MPFR_RNDN); mpfr_init(r21139); mpfr_init(r21140); mpfr_init(r21141); } double f_dm(double v, double w, double r) { ; mpfr_set_d(r21123, r, MPFR_RNDN); mpfr_mul(r21124, r21123, r21123, MPFR_RNDN); mpfr_div(r21125, r21122, r21124, MPFR_RNDN); ; mpfr_set_d(r21127, v, MPFR_RNDN); mpfr_mul(r21128, r21127, r21122, MPFR_RNDN); mpfr_sub(r21129, r21126, r21128, MPFR_RNDN); ; mpfr_sub(r21131, r21130, r21127, MPFR_RNDN); ; mpfr_div(r21133, r21131, r21132, MPFR_RNDN); mpfr_div(r21134, r21129, r21133, MPFR_RNDN); mpfr_set_d(r21135, w, MPFR_RNDN); mpfr_mul(r21136, r21135, r21123, MPFR_RNDN); mpfr_mul(r21137, r21136, r21136, MPFR_RNDN); ; mpfr_fma(r21139, r21134, r21137, r21138, MPFR_RNDN); mpfr_sub(r21140, r21126, r21139, MPFR_RNDN); mpfr_add(r21141, r21125, r21140, MPFR_RNDN); return mpfr_get_d(r21141, MPFR_RNDN); }