#include #include #include #include #include char *name = "Rosa's TurbineBenchmark"; double f_if(float v, float w, float r) { float r21351 = 3; float r21352 = 2; float r21353 = r; float r21354 = r21353 * r21353; float r21355 = r21352 / r21354; float r21356 = r21351 + r21355; float r21357 = 0.125; float r21358 = v; float r21359 = r21352 * r21358; float r21360 = r21351 - r21359; float r21361 = r21357 * r21360; float r21362 = w; float r21363 = r21362 * r21362; float r21364 = r21363 * r21353; float r21365 = r21364 * r21353; float r21366 = r21361 * r21365; float r21367 = 1; float r21368 = r21367 - r21358; float r21369 = r21366 / r21368; float r21370 = r21356 - r21369; float r21371 = 4.5; float r21372 = r21370 - r21371; return r21372; } double f_id(double v, double w, double r) { double r21373 = 3; double r21374 = 2; double r21375 = r; double r21376 = r21375 * r21375; double r21377 = r21374 / r21376; double r21378 = r21373 + r21377; double r21379 = 0.125; double r21380 = v; double r21381 = r21374 * r21380; double r21382 = r21373 - r21381; double r21383 = r21379 * r21382; double r21384 = w; double r21385 = r21384 * r21384; double r21386 = r21385 * r21375; double r21387 = r21386 * r21375; double r21388 = r21383 * r21387; double r21389 = 1; double r21390 = r21389 - r21380; double r21391 = r21388 / r21390; double r21392 = r21378 - r21391; double r21393 = 4.5; double r21394 = r21392 - r21393; return r21394; } double f_of(float v, float w, float r) { float r21395 = 1; float r21396 = r; float r21397 = 2; float r21398 = r21397 / r21396; float r21399 = r21396 / r21398; float r21400 = sqrt(r21399); float r21401 = r21395 / r21400; float r21402 = r21401 / r21400; float r21403 = 3; float r21404 = r21402 + r21403; float r21405 = v; float r21406 = r21405 * r21397; float r21407 = r21403 - r21406; float r21408 = r21395 - r21405; float r21409 = 0.125; float r21410 = r21408 / r21409; float r21411 = r21407 / r21410; float r21412 = w; float r21413 = r21412 * r21396; float r21414 = r21413 * r21413; float r21415 = 4.5; float r21416 = fma(r21411, r21414, r21415); float r21417 = r21404 - r21416; return r21417; } double f_od(double v, double w, double r) { double r21418 = 1; double r21419 = r; double r21420 = 2; double r21421 = r21420 / r21419; double r21422 = r21419 / r21421; double r21423 = sqrt(r21422); double r21424 = r21418 / r21423; double r21425 = r21424 / r21423; double r21426 = 3; double r21427 = r21425 + r21426; double r21428 = v; double r21429 = r21428 * r21420; double r21430 = r21426 - r21429; double r21431 = r21418 - r21428; double r21432 = 0.125; double r21433 = r21431 / r21432; double r21434 = r21430 / r21433; double r21435 = w; double r21436 = r21435 * r21419; double r21437 = r21436 * r21436; double r21438 = 4.5; double r21439 = fma(r21434, r21437, r21438); double r21440 = r21427 - r21439; return r21440; } 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 r21441, r21442, r21443, r21444, r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459, r21460, r21461, r21462; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21441, "3", 10, MPFR_RNDN); mpfr_init_set_str(r21442, "2", 10, MPFR_RNDN); mpfr_init(r21443); mpfr_init(r21444); mpfr_init(r21445); mpfr_init(r21446); mpfr_init_set_str(r21447, "0.125", 10, MPFR_RNDN); mpfr_init(r21448); mpfr_init(r21449); mpfr_init(r21450); mpfr_init(r21451); mpfr_init(r21452); mpfr_init(r21453); mpfr_init(r21454); mpfr_init(r21455); mpfr_init(r21456); mpfr_init_set_str(r21457, "1", 10, MPFR_RNDN); mpfr_init(r21458); mpfr_init(r21459); mpfr_init(r21460); mpfr_init_set_str(r21461, "4.5", 10, MPFR_RNDN); mpfr_init(r21462); } double f_im(double v, double w, double r) { ; ; mpfr_set_d(r21443, r, MPFR_RNDN); mpfr_mul(r21444, r21443, r21443, MPFR_RNDN); mpfr_div(r21445, r21442, r21444, MPFR_RNDN); mpfr_add(r21446, r21441, r21445, MPFR_RNDN); ; mpfr_set_d(r21448, v, MPFR_RNDN); mpfr_mul(r21449, r21442, r21448, MPFR_RNDN); mpfr_sub(r21450, r21441, r21449, MPFR_RNDN); mpfr_mul(r21451, r21447, r21450, MPFR_RNDN); mpfr_set_d(r21452, w, MPFR_RNDN); mpfr_mul(r21453, r21452, r21452, MPFR_RNDN); mpfr_mul(r21454, r21453, r21443, MPFR_RNDN); mpfr_mul(r21455, r21454, r21443, MPFR_RNDN); mpfr_mul(r21456, r21451, r21455, MPFR_RNDN); ; mpfr_sub(r21458, r21457, r21448, MPFR_RNDN); mpfr_div(r21459, r21456, r21458, MPFR_RNDN); mpfr_sub(r21460, r21446, r21459, MPFR_RNDN); ; mpfr_sub(r21462, r21460, r21461, MPFR_RNDN); return mpfr_get_d(r21462, MPFR_RNDN); } static mpfr_t r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476, r21477, r21478, r21479, r21480, r21481, r21482, r21483, r21484, r21485; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21463, "1", 10, MPFR_RNDN); mpfr_init(r21464); mpfr_init_set_str(r21465, "2", 10, MPFR_RNDN); mpfr_init(r21466); mpfr_init(r21467); mpfr_init(r21468); mpfr_init(r21469); mpfr_init(r21470); mpfr_init_set_str(r21471, "3", 10, MPFR_RNDN); mpfr_init(r21472); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init_set_str(r21477, "0.125", 10, MPFR_RNDN); mpfr_init(r21478); mpfr_init(r21479); mpfr_init(r21480); mpfr_init(r21481); mpfr_init(r21482); mpfr_init_set_str(r21483, "4.5", 10, MPFR_RNDN); mpfr_init(r21484); mpfr_init(r21485); } double f_fm(double v, double w, double r) { ; mpfr_set_d(r21464, r, MPFR_RNDN); ; mpfr_div(r21466, r21465, r21464, MPFR_RNDN); mpfr_div(r21467, r21464, r21466, MPFR_RNDN); mpfr_sqrt(r21468, r21467, MPFR_RNDN); mpfr_div(r21469, r21463, r21468, MPFR_RNDN); mpfr_div(r21470, r21469, r21468, MPFR_RNDN); ; mpfr_add(r21472, r21470, r21471, MPFR_RNDN); mpfr_set_d(r21473, v, MPFR_RNDN); mpfr_mul(r21474, r21473, r21465, MPFR_RNDN); mpfr_sub(r21475, r21471, r21474, MPFR_RNDN); mpfr_sub(r21476, r21463, r21473, MPFR_RNDN); ; mpfr_div(r21478, r21476, r21477, MPFR_RNDN); mpfr_div(r21479, r21475, r21478, MPFR_RNDN); mpfr_set_d(r21480, w, MPFR_RNDN); mpfr_mul(r21481, r21480, r21464, MPFR_RNDN); mpfr_mul(r21482, r21481, r21481, MPFR_RNDN); ; mpfr_fma(r21484, r21479, r21482, r21483, MPFR_RNDN); mpfr_sub(r21485, r21472, r21484, MPFR_RNDN); return mpfr_get_d(r21485, MPFR_RNDN); } static mpfr_t r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494, r21495, r21496, r21497, r21498, r21499, r21500, r21501, r21502, r21503, r21504, r21505, r21506, r21507, r21508; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21486, "1", 10, MPFR_RNDN); mpfr_init(r21487); mpfr_init_set_str(r21488, "2", 10, MPFR_RNDN); mpfr_init(r21489); mpfr_init(r21490); mpfr_init(r21491); mpfr_init(r21492); mpfr_init(r21493); mpfr_init_set_str(r21494, "3", 10, MPFR_RNDN); mpfr_init(r21495); mpfr_init(r21496); mpfr_init(r21497); mpfr_init(r21498); mpfr_init(r21499); mpfr_init_set_str(r21500, "0.125", 10, MPFR_RNDN); mpfr_init(r21501); mpfr_init(r21502); mpfr_init(r21503); mpfr_init(r21504); mpfr_init(r21505); mpfr_init_set_str(r21506, "4.5", 10, MPFR_RNDN); mpfr_init(r21507); mpfr_init(r21508); } double f_dm(double v, double w, double r) { ; mpfr_set_d(r21487, r, MPFR_RNDN); ; mpfr_div(r21489, r21488, r21487, MPFR_RNDN); mpfr_div(r21490, r21487, r21489, MPFR_RNDN); mpfr_sqrt(r21491, r21490, MPFR_RNDN); mpfr_div(r21492, r21486, r21491, MPFR_RNDN); mpfr_div(r21493, r21492, r21491, MPFR_RNDN); ; mpfr_add(r21495, r21493, r21494, MPFR_RNDN); mpfr_set_d(r21496, v, MPFR_RNDN); mpfr_mul(r21497, r21496, r21488, MPFR_RNDN); mpfr_sub(r21498, r21494, r21497, MPFR_RNDN); mpfr_sub(r21499, r21486, r21496, MPFR_RNDN); ; mpfr_div(r21501, r21499, r21500, MPFR_RNDN); mpfr_div(r21502, r21498, r21501, MPFR_RNDN); mpfr_set_d(r21503, w, MPFR_RNDN); mpfr_mul(r21504, r21503, r21487, MPFR_RNDN); mpfr_mul(r21505, r21504, r21504, MPFR_RNDN); ; mpfr_fma(r21507, r21502, r21505, r21506, MPFR_RNDN); mpfr_sub(r21508, r21495, r21507, MPFR_RNDN); return mpfr_get_d(r21508, MPFR_RNDN); }