#include #include #include #include #include char *name = "Rosa's TurbineBenchmark"; double f_if(float v, float w, float r) { float r19649 = 3; float r19650 = 2; float r19651 = r; float r19652 = r19651 * r19651; float r19653 = r19650 / r19652; float r19654 = r19649 + r19653; float r19655 = 0.125; float r19656 = v; float r19657 = r19650 * r19656; float r19658 = r19649 - r19657; float r19659 = r19655 * r19658; float r19660 = w; float r19661 = r19660 * r19660; float r19662 = r19661 * r19651; float r19663 = r19662 * r19651; float r19664 = r19659 * r19663; float r19665 = 1; float r19666 = r19665 - r19656; float r19667 = r19664 / r19666; float r19668 = r19654 - r19667; float r19669 = 4.5; float r19670 = r19668 - r19669; return r19670; } double f_id(double v, double w, double r) { double r19671 = 3; double r19672 = 2; double r19673 = r; double r19674 = r19673 * r19673; double r19675 = r19672 / r19674; double r19676 = r19671 + r19675; double r19677 = 0.125; double r19678 = v; double r19679 = r19672 * r19678; double r19680 = r19671 - r19679; double r19681 = r19677 * r19680; double r19682 = w; double r19683 = r19682 * r19682; double r19684 = r19683 * r19673; double r19685 = r19684 * r19673; double r19686 = r19681 * r19685; double r19687 = 1; double r19688 = r19687 - r19678; double r19689 = r19686 / r19688; double r19690 = r19676 - r19689; double r19691 = 4.5; double r19692 = r19690 - r19691; return r19692; } double f_of(float v, float w, float r) { float r19693 = 3; float r19694 = 2; float r19695 = r; float r19696 = r19695 * r19695; float r19697 = r19694 / r19696; float r19698 = r19693 + r19697; float r19699 = 0.125; float r19700 = v; float r19701 = r19694 * r19700; float r19702 = r19693 - r19701; float r19703 = r19699 * r19702; float r19704 = 1; float r19705 = r19704 - r19700; float r19706 = w; float r19707 = r19706 * r19695; float r19708 = fabs(r19707); float r19709 = cbrt(r19708); float r19710 = r19709 * r19709; float r19711 = r19708 * r19710; float r19712 = r19711 * r19709; float r19713 = r19705 / r19712; float r19714 = r19703 / r19713; float r19715 = r19698 - r19714; float r19716 = 4.5; float r19717 = r19715 - r19716; return r19717; } double f_od(double v, double w, double r) { double r19718 = 3; double r19719 = 2; double r19720 = r; double r19721 = r19720 * r19720; double r19722 = r19719 / r19721; double r19723 = r19718 + r19722; double r19724 = 0.125; double r19725 = v; double r19726 = r19719 * r19725; double r19727 = r19718 - r19726; double r19728 = r19724 * r19727; double r19729 = 1; double r19730 = r19729 - r19725; double r19731 = w; double r19732 = r19731 * r19720; double r19733 = fabs(r19732); double r19734 = cbrt(r19733); double r19735 = r19734 * r19734; double r19736 = r19733 * r19735; double r19737 = r19736 * r19734; double r19738 = r19730 / r19737; double r19739 = r19728 / r19738; double r19740 = r19723 - r19739; double r19741 = 4.5; double r19742 = r19740 - r19741; return r19742; } 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 r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r19743, "3", 10, MPFR_RNDN); mpfr_init_set_str(r19744, "2", 10, MPFR_RNDN); mpfr_init(r19745); mpfr_init(r19746); mpfr_init(r19747); mpfr_init(r19748); mpfr_init_set_str(r19749, "0.125", 10, MPFR_RNDN); mpfr_init(r19750); mpfr_init(r19751); mpfr_init(r19752); mpfr_init(r19753); mpfr_init(r19754); mpfr_init(r19755); mpfr_init(r19756); mpfr_init(r19757); mpfr_init(r19758); mpfr_init_set_str(r19759, "1", 10, MPFR_RNDN); mpfr_init(r19760); mpfr_init(r19761); mpfr_init(r19762); mpfr_init_set_str(r19763, "4.5", 10, MPFR_RNDN); mpfr_init(r19764); } double f_im(double v, double w, double r) { ; ; mpfr_set_d(r19745, r, MPFR_RNDN); mpfr_mul(r19746, r19745, r19745, MPFR_RNDN); mpfr_div(r19747, r19744, r19746, MPFR_RNDN); mpfr_add(r19748, r19743, r19747, MPFR_RNDN); ; mpfr_set_d(r19750, v, MPFR_RNDN); mpfr_mul(r19751, r19744, r19750, MPFR_RNDN); mpfr_sub(r19752, r19743, r19751, MPFR_RNDN); mpfr_mul(r19753, r19749, r19752, MPFR_RNDN); mpfr_set_d(r19754, w, MPFR_RNDN); mpfr_mul(r19755, r19754, r19754, MPFR_RNDN); mpfr_mul(r19756, r19755, r19745, MPFR_RNDN); mpfr_mul(r19757, r19756, r19745, MPFR_RNDN); mpfr_mul(r19758, r19753, r19757, MPFR_RNDN); ; mpfr_sub(r19760, r19759, r19750, MPFR_RNDN); mpfr_div(r19761, r19758, r19760, MPFR_RNDN); mpfr_sub(r19762, r19748, r19761, MPFR_RNDN); ; mpfr_sub(r19764, r19762, r19763, MPFR_RNDN); return mpfr_get_d(r19764, MPFR_RNDN); } static mpfr_t r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r19765, "3", 10, MPFR_RNDN); mpfr_init_set_str(r19766, "2", 10, MPFR_RNDN); mpfr_init(r19767); mpfr_init(r19768); mpfr_init(r19769); mpfr_init(r19770); mpfr_init_set_str(r19771, "0.125", 10, MPFR_RNDN); mpfr_init(r19772); mpfr_init(r19773); mpfr_init(r19774); mpfr_init(r19775); mpfr_init_set_str(r19776, "1", 10, MPFR_RNDN); mpfr_init(r19777); mpfr_init(r19778); mpfr_init(r19779); mpfr_init(r19780); mpfr_init(r19781); mpfr_init(r19782); mpfr_init(r19783); mpfr_init(r19784); mpfr_init(r19785); mpfr_init(r19786); mpfr_init(r19787); mpfr_init_set_str(r19788, "4.5", 10, MPFR_RNDN); mpfr_init(r19789); } double f_fm(double v, double w, double r) { ; ; mpfr_set_d(r19767, r, MPFR_RNDN); mpfr_mul(r19768, r19767, r19767, MPFR_RNDN); mpfr_div(r19769, r19766, r19768, MPFR_RNDN); mpfr_add(r19770, r19765, r19769, MPFR_RNDN); ; mpfr_set_d(r19772, v, MPFR_RNDN); mpfr_mul(r19773, r19766, r19772, MPFR_RNDN); mpfr_sub(r19774, r19765, r19773, MPFR_RNDN); mpfr_mul(r19775, r19771, r19774, MPFR_RNDN); ; mpfr_sub(r19777, r19776, r19772, MPFR_RNDN); mpfr_set_d(r19778, w, MPFR_RNDN); mpfr_mul(r19779, r19778, r19767, MPFR_RNDN); mpfr_abs(r19780, r19779, MPFR_RNDN); mpfr_cbrt(r19781, r19780, MPFR_RNDN); mpfr_mul(r19782, r19781, r19781, MPFR_RNDN); mpfr_mul(r19783, r19780, r19782, MPFR_RNDN); mpfr_mul(r19784, r19783, r19781, MPFR_RNDN); mpfr_div(r19785, r19777, r19784, MPFR_RNDN); mpfr_div(r19786, r19775, r19785, MPFR_RNDN); mpfr_sub(r19787, r19770, r19786, MPFR_RNDN); ; mpfr_sub(r19789, r19787, r19788, MPFR_RNDN); return mpfr_get_d(r19789, MPFR_RNDN); } static mpfr_t r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803, r19804, r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r19790, "3", 10, MPFR_RNDN); mpfr_init_set_str(r19791, "2", 10, MPFR_RNDN); mpfr_init(r19792); mpfr_init(r19793); mpfr_init(r19794); mpfr_init(r19795); mpfr_init_set_str(r19796, "0.125", 10, MPFR_RNDN); mpfr_init(r19797); mpfr_init(r19798); mpfr_init(r19799); mpfr_init(r19800); mpfr_init_set_str(r19801, "1", 10, MPFR_RNDN); mpfr_init(r19802); mpfr_init(r19803); mpfr_init(r19804); mpfr_init(r19805); mpfr_init(r19806); mpfr_init(r19807); mpfr_init(r19808); mpfr_init(r19809); mpfr_init(r19810); mpfr_init(r19811); mpfr_init(r19812); mpfr_init_set_str(r19813, "4.5", 10, MPFR_RNDN); mpfr_init(r19814); } double f_dm(double v, double w, double r) { ; ; mpfr_set_d(r19792, r, MPFR_RNDN); mpfr_mul(r19793, r19792, r19792, MPFR_RNDN); mpfr_div(r19794, r19791, r19793, MPFR_RNDN); mpfr_add(r19795, r19790, r19794, MPFR_RNDN); ; mpfr_set_d(r19797, v, MPFR_RNDN); mpfr_mul(r19798, r19791, r19797, MPFR_RNDN); mpfr_sub(r19799, r19790, r19798, MPFR_RNDN); mpfr_mul(r19800, r19796, r19799, MPFR_RNDN); ; mpfr_sub(r19802, r19801, r19797, MPFR_RNDN); mpfr_set_d(r19803, w, MPFR_RNDN); mpfr_mul(r19804, r19803, r19792, MPFR_RNDN); mpfr_abs(r19805, r19804, MPFR_RNDN); mpfr_cbrt(r19806, r19805, MPFR_RNDN); mpfr_mul(r19807, r19806, r19806, MPFR_RNDN); mpfr_mul(r19808, r19805, r19807, MPFR_RNDN); mpfr_mul(r19809, r19808, r19806, MPFR_RNDN); mpfr_div(r19810, r19802, r19809, MPFR_RNDN); mpfr_div(r19811, r19800, r19810, MPFR_RNDN); mpfr_sub(r19812, r19795, r19811, MPFR_RNDN); ; mpfr_sub(r19814, r19812, r19813, MPFR_RNDN); return mpfr_get_d(r19814, MPFR_RNDN); }