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