#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 = 2; float r21396 = r; float r21397 = r21395 / r21396; float r21398 = r21397 / r21396; float r21399 = 3; float r21400 = r21398 + r21399; float r21401 = v; float r21402 = r21401 * r21395; float r21403 = r21399 - r21402; float r21404 = 1; float r21405 = r21404 - r21401; float r21406 = 0.125; float r21407 = r21405 / r21406; float r21408 = r21403 / r21407; float r21409 = sqrt(r21408); float r21410 = r21409 * r21409; float r21411 = w; float r21412 = r21411 * r21396; float r21413 = r21412 * r21412; float r21414 = 4.5; float r21415 = fma(r21410, r21413, r21414); float r21416 = r21400 - r21415; return r21416; } double f_od(double v, double w, double r) { double r21417 = 2; double r21418 = r; double r21419 = r21417 / r21418; double r21420 = r21419 / r21418; double r21421 = 3; double r21422 = r21420 + r21421; double r21423 = v; double r21424 = r21423 * r21417; double r21425 = r21421 - r21424; double r21426 = 1; double r21427 = r21426 - r21423; double r21428 = 0.125; double r21429 = r21427 / r21428; double r21430 = r21425 / r21429; double r21431 = sqrt(r21430); double r21432 = r21431 * r21431; double r21433 = w; double r21434 = r21433 * r21418; double r21435 = r21434 * r21434; double r21436 = 4.5; double r21437 = fma(r21432, r21435, r21436); double r21438 = r21422 - r21437; return r21438; } 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 r21439, r21440, r21441, r21442, r21443, r21444, r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459, r21460; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21439, "3", 10, MPFR_RNDN); mpfr_init_set_str(r21440, "2", 10, MPFR_RNDN); mpfr_init(r21441); mpfr_init(r21442); mpfr_init(r21443); mpfr_init(r21444); mpfr_init_set_str(r21445, "0.125", 10, MPFR_RNDN); mpfr_init(r21446); mpfr_init(r21447); mpfr_init(r21448); mpfr_init(r21449); mpfr_init(r21450); mpfr_init(r21451); mpfr_init(r21452); mpfr_init(r21453); mpfr_init(r21454); mpfr_init_set_str(r21455, "1", 10, MPFR_RNDN); mpfr_init(r21456); mpfr_init(r21457); mpfr_init(r21458); mpfr_init_set_str(r21459, "4.5", 10, MPFR_RNDN); mpfr_init(r21460); } double f_im(double v, double w, double r) { ; ; mpfr_set_d(r21441, r, MPFR_RNDN); mpfr_mul(r21442, r21441, r21441, MPFR_RNDN); mpfr_div(r21443, r21440, r21442, MPFR_RNDN); mpfr_add(r21444, r21439, r21443, MPFR_RNDN); ; mpfr_set_d(r21446, v, MPFR_RNDN); mpfr_mul(r21447, r21440, r21446, MPFR_RNDN); mpfr_sub(r21448, r21439, r21447, MPFR_RNDN); mpfr_mul(r21449, r21445, r21448, MPFR_RNDN); mpfr_set_d(r21450, w, MPFR_RNDN); mpfr_mul(r21451, r21450, r21450, MPFR_RNDN); mpfr_mul(r21452, r21451, r21441, MPFR_RNDN); mpfr_mul(r21453, r21452, r21441, MPFR_RNDN); mpfr_mul(r21454, r21449, r21453, MPFR_RNDN); ; mpfr_sub(r21456, r21455, r21446, MPFR_RNDN); mpfr_div(r21457, r21454, r21456, MPFR_RNDN); mpfr_sub(r21458, r21444, r21457, MPFR_RNDN); ; mpfr_sub(r21460, r21458, r21459, MPFR_RNDN); return mpfr_get_d(r21460, MPFR_RNDN); } static mpfr_t r21461, r21462, r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476, r21477, r21478, r21479, r21480, r21481, r21482; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21461, "2", 10, MPFR_RNDN); mpfr_init(r21462); mpfr_init(r21463); mpfr_init(r21464); mpfr_init_set_str(r21465, "3", 10, MPFR_RNDN); mpfr_init(r21466); mpfr_init(r21467); mpfr_init(r21468); mpfr_init(r21469); mpfr_init_set_str(r21470, "1", 10, MPFR_RNDN); mpfr_init(r21471); mpfr_init_set_str(r21472, "0.125", 10, MPFR_RNDN); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init(r21477); mpfr_init(r21478); mpfr_init(r21479); mpfr_init_set_str(r21480, "4.5", 10, MPFR_RNDN); mpfr_init(r21481); mpfr_init(r21482); } double f_fm(double v, double w, double r) { ; mpfr_set_d(r21462, r, MPFR_RNDN); mpfr_div(r21463, r21461, r21462, MPFR_RNDN); mpfr_div(r21464, r21463, r21462, MPFR_RNDN); ; mpfr_add(r21466, r21464, r21465, MPFR_RNDN); mpfr_set_d(r21467, v, MPFR_RNDN); mpfr_mul(r21468, r21467, r21461, MPFR_RNDN); mpfr_sub(r21469, r21465, r21468, MPFR_RNDN); ; mpfr_sub(r21471, r21470, r21467, MPFR_RNDN); ; mpfr_div(r21473, r21471, r21472, MPFR_RNDN); mpfr_div(r21474, r21469, r21473, MPFR_RNDN); mpfr_sqrt(r21475, r21474, MPFR_RNDN); mpfr_mul(r21476, r21475, r21475, MPFR_RNDN); mpfr_set_d(r21477, w, MPFR_RNDN); mpfr_mul(r21478, r21477, r21462, MPFR_RNDN); mpfr_mul(r21479, r21478, r21478, MPFR_RNDN); ; mpfr_fma(r21481, r21476, r21479, r21480, MPFR_RNDN); mpfr_sub(r21482, r21466, r21481, MPFR_RNDN); return mpfr_get_d(r21482, MPFR_RNDN); } static mpfr_t r21483, r21484, r21485, r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494, r21495, r21496, r21497, r21498, r21499, r21500, r21501, r21502, r21503, r21504; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21483, "2", 10, MPFR_RNDN); mpfr_init(r21484); mpfr_init(r21485); mpfr_init(r21486); mpfr_init_set_str(r21487, "3", 10, MPFR_RNDN); mpfr_init(r21488); mpfr_init(r21489); mpfr_init(r21490); mpfr_init(r21491); mpfr_init_set_str(r21492, "1", 10, MPFR_RNDN); mpfr_init(r21493); mpfr_init_set_str(r21494, "0.125", 10, MPFR_RNDN); mpfr_init(r21495); mpfr_init(r21496); mpfr_init(r21497); mpfr_init(r21498); mpfr_init(r21499); mpfr_init(r21500); mpfr_init(r21501); mpfr_init_set_str(r21502, "4.5", 10, MPFR_RNDN); mpfr_init(r21503); mpfr_init(r21504); } double f_dm(double v, double w, double r) { ; mpfr_set_d(r21484, r, MPFR_RNDN); mpfr_div(r21485, r21483, r21484, MPFR_RNDN); mpfr_div(r21486, r21485, r21484, MPFR_RNDN); ; mpfr_add(r21488, r21486, r21487, MPFR_RNDN); mpfr_set_d(r21489, v, MPFR_RNDN); mpfr_mul(r21490, r21489, r21483, MPFR_RNDN); mpfr_sub(r21491, r21487, r21490, MPFR_RNDN); ; mpfr_sub(r21493, r21492, r21489, MPFR_RNDN); ; mpfr_div(r21495, r21493, r21494, MPFR_RNDN); mpfr_div(r21496, r21491, r21495, MPFR_RNDN); mpfr_sqrt(r21497, r21496, MPFR_RNDN); mpfr_mul(r21498, r21497, r21497, MPFR_RNDN); mpfr_set_d(r21499, w, MPFR_RNDN); mpfr_mul(r21500, r21499, r21484, MPFR_RNDN); mpfr_mul(r21501, r21500, r21500, MPFR_RNDN); ; mpfr_fma(r21503, r21498, r21501, r21502, MPFR_RNDN); mpfr_sub(r21504, r21488, r21503, MPFR_RNDN); return mpfr_get_d(r21504, MPFR_RNDN); }