#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Rosa's TurbineBenchmark";

double f_if(float v, float w, float r) {
        float r25648 = 3;
        float r25649 = 2;
        float r25650 = r;
        float r25651 = r25650 * r25650;
        float r25652 = r25649 / r25651;
        float r25653 = r25648 + r25652;
        float r25654 = 0.125;
        float r25655 = v;
        float r25656 = r25649 * r25655;
        float r25657 = r25648 - r25656;
        float r25658 = r25654 * r25657;
        float r25659 = w;
        float r25660 = r25659 * r25659;
        float r25661 = r25660 * r25650;
        float r25662 = r25661 * r25650;
        float r25663 = r25658 * r25662;
        float r25664 = 1;
        float r25665 = r25664 - r25655;
        float r25666 = r25663 / r25665;
        float r25667 = r25653 - r25666;
        float r25668 = 4.5;
        float r25669 = r25667 - r25668;
        return r25669;
}

double f_id(double v, double w, double r) {
        double r25670 = 3;
        double r25671 = 2;
        double r25672 = r;
        double r25673 = r25672 * r25672;
        double r25674 = r25671 / r25673;
        double r25675 = r25670 + r25674;
        double r25676 = 0.125;
        double r25677 = v;
        double r25678 = r25671 * r25677;
        double r25679 = r25670 - r25678;
        double r25680 = r25676 * r25679;
        double r25681 = w;
        double r25682 = r25681 * r25681;
        double r25683 = r25682 * r25672;
        double r25684 = r25683 * r25672;
        double r25685 = r25680 * r25684;
        double r25686 = 1;
        double r25687 = r25686 - r25677;
        double r25688 = r25685 / r25687;
        double r25689 = r25675 - r25688;
        double r25690 = 4.5;
        double r25691 = r25689 - r25690;
        return r25691;
}


double f_of(float v, float w, float r) {
        float r25692 = v;
        float r25693 = -2.526418831017353e+31;
        bool r25694 = r25692 <= r25693;
        float r25695 = 2;
        float r25696 = r;
        float r25697 = r25696 * r25696;
        float r25698 = r25695 / r25697;
        float r25699 = 3;
        float r25700 = 4.5;
        float r25701 = r25699 - r25700;
        float r25702 = r25698 + r25701;
        float r25703 = w;
        float r25704 = r25703 * r25696;
        float r25705 = r25704 * r25704;
        float r25706 = 0.25;
        float r25707 = r25705 * r25706;
        float r25708 = r25705 / r25692;
        float r25709 = 0.125;
        float r25710 = r25709 / r25692;
        float r25711 = r25709 + r25710;
        float r25712 = r25708 * r25711;
        float r25713 = r25707 - r25712;
        float r25714 = r25702 - r25713;
        float r25715 = 4.907722326036034e-19;
        bool r25716 = r25692 <= r25715;
        float r25717 = r25695 / r25696;
        float r25718 = r25717 / r25696;
        float r25719 = r25699 + r25718;
        float r25720 = r25696 * r25703;
        float r25721 = r25720 * r25720;
        float r25722 = r25709 * r25721;
        float r25723 = 1;
        float r25724 = r25692 * r25692;
        float r25725 = r25723 - r25724;
        float r25726 = r25692 + r25692;
        float r25727 = r25699 - r25726;
        float r25728 = r25725 / r25727;
        float r25729 = r25722 / r25728;
        float r25730 = r25723 + r25692;
        float r25731 = r25729 * r25730;
        float r25732 = r25719 - r25731;
        float r25733 = r25732 - r25700;
        float r25734 = r25716 ? r25733 : r25714;
        float r25735 = r25694 ? r25714 : r25734;
        return r25735;
}

double f_od(double v, double w, double r) {
        double r25736 = v;
        double r25737 = -2.526418831017353e+31;
        bool r25738 = r25736 <= r25737;
        double r25739 = 2;
        double r25740 = r;
        double r25741 = r25740 * r25740;
        double r25742 = r25739 / r25741;
        double r25743 = 3;
        double r25744 = 4.5;
        double r25745 = r25743 - r25744;
        double r25746 = r25742 + r25745;
        double r25747 = w;
        double r25748 = r25747 * r25740;
        double r25749 = r25748 * r25748;
        double r25750 = 0.25;
        double r25751 = r25749 * r25750;
        double r25752 = r25749 / r25736;
        double r25753 = 0.125;
        double r25754 = r25753 / r25736;
        double r25755 = r25753 + r25754;
        double r25756 = r25752 * r25755;
        double r25757 = r25751 - r25756;
        double r25758 = r25746 - r25757;
        double r25759 = 4.907722326036034e-19;
        bool r25760 = r25736 <= r25759;
        double r25761 = r25739 / r25740;
        double r25762 = r25761 / r25740;
        double r25763 = r25743 + r25762;
        double r25764 = r25740 * r25747;
        double r25765 = r25764 * r25764;
        double r25766 = r25753 * r25765;
        double r25767 = 1;
        double r25768 = r25736 * r25736;
        double r25769 = r25767 - r25768;
        double r25770 = r25736 + r25736;
        double r25771 = r25743 - r25770;
        double r25772 = r25769 / r25771;
        double r25773 = r25766 / r25772;
        double r25774 = r25767 + r25736;
        double r25775 = r25773 * r25774;
        double r25776 = r25763 - r25775;
        double r25777 = r25776 - r25744;
        double r25778 = r25760 ? r25777 : r25758;
        double r25779 = r25738 ? r25758 : r25778;
        return r25779;
}

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 r25780, r25781, r25782, r25783, r25784, r25785, r25786, r25787, r25788, r25789, r25790, r25791, r25792, r25793, r25794, r25795, r25796, r25797, r25798, r25799, r25800, r25801;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25780, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r25781, "2", 10, MPFR_RNDN);
        mpfr_init(r25782);
        mpfr_init(r25783);
        mpfr_init(r25784);
        mpfr_init(r25785);
        mpfr_init_set_str(r25786, "0.125", 10, MPFR_RNDN);
        mpfr_init(r25787);
        mpfr_init(r25788);
        mpfr_init(r25789);
        mpfr_init(r25790);
        mpfr_init(r25791);
        mpfr_init(r25792);
        mpfr_init(r25793);
        mpfr_init(r25794);
        mpfr_init(r25795);
        mpfr_init_set_str(r25796, "1", 10, MPFR_RNDN);
        mpfr_init(r25797);
        mpfr_init(r25798);
        mpfr_init(r25799);
        mpfr_init_set_str(r25800, "4.5", 10, MPFR_RNDN);
        mpfr_init(r25801);
}

double f_im(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r25782, r, MPFR_RNDN);
        mpfr_mul(r25783, r25782, r25782, MPFR_RNDN);
        mpfr_div(r25784, r25781, r25783, MPFR_RNDN);
        mpfr_add(r25785, r25780, r25784, MPFR_RNDN);
        ;
        mpfr_set_d(r25787, v, MPFR_RNDN);
        mpfr_mul(r25788, r25781, r25787, MPFR_RNDN);
        mpfr_sub(r25789, r25780, r25788, MPFR_RNDN);
        mpfr_mul(r25790, r25786, r25789, MPFR_RNDN);
        mpfr_set_d(r25791, w, MPFR_RNDN);
        mpfr_mul(r25792, r25791, r25791, MPFR_RNDN);
        mpfr_mul(r25793, r25792, r25782, MPFR_RNDN);
        mpfr_mul(r25794, r25793, r25782, MPFR_RNDN);
        mpfr_mul(r25795, r25790, r25794, MPFR_RNDN);
        ;
        mpfr_sub(r25797, r25796, r25787, MPFR_RNDN);
        mpfr_div(r25798, r25795, r25797, MPFR_RNDN);
        mpfr_sub(r25799, r25785, r25798, MPFR_RNDN);
        ;
        mpfr_sub(r25801, r25799, r25800, MPFR_RNDN);
        return mpfr_get_d(r25801, MPFR_RNDN);
}

static mpfr_t r25802, r25803, r25804, r25805, r25806, r25807, r25808, r25809, r25810, r25811, r25812, r25813, r25814, r25815, r25816, r25817, r25818, r25819, r25820, r25821, r25822, r25823, r25824, r25825, r25826, r25827, r25828, r25829, r25830, r25831, r25832, r25833, r25834, r25835, r25836, r25837, r25838, r25839, r25840, r25841, r25842, r25843, r25844, r25845;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25802);
        mpfr_init_set_str(r25803, "-2.526418831017353e+31", 10, MPFR_RNDN);
        mpfr_init(r25804);
        mpfr_init_set_str(r25805, "2", 10, MPFR_RNDN);
        mpfr_init(r25806);
        mpfr_init(r25807);
        mpfr_init(r25808);
        mpfr_init_set_str(r25809, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r25810, "4.5", 10, MPFR_RNDN);
        mpfr_init(r25811);
        mpfr_init(r25812);
        mpfr_init(r25813);
        mpfr_init(r25814);
        mpfr_init(r25815);
        mpfr_init_set_str(r25816, "0.25", 10, MPFR_RNDN);
        mpfr_init(r25817);
        mpfr_init(r25818);
        mpfr_init_set_str(r25819, "0.125", 10, MPFR_RNDN);
        mpfr_init(r25820);
        mpfr_init(r25821);
        mpfr_init(r25822);
        mpfr_init(r25823);
        mpfr_init(r25824);
        mpfr_init_set_str(r25825, "4.907722326036034e-19", 10, MPFR_RNDN);
        mpfr_init(r25826);
        mpfr_init(r25827);
        mpfr_init(r25828);
        mpfr_init(r25829);
        mpfr_init(r25830);
        mpfr_init(r25831);
        mpfr_init(r25832);
        mpfr_init_set_str(r25833, "1", 10, MPFR_RNDN);
        mpfr_init(r25834);
        mpfr_init(r25835);
        mpfr_init(r25836);
        mpfr_init(r25837);
        mpfr_init(r25838);
        mpfr_init(r25839);
        mpfr_init(r25840);
        mpfr_init(r25841);
        mpfr_init(r25842);
        mpfr_init(r25843);
        mpfr_init(r25844);
        mpfr_init(r25845);
}

double f_fm(double v, double w, double r) {
        mpfr_set_d(r25802, v, MPFR_RNDN);
        ;
        mpfr_set_si(r25804, mpfr_cmp(r25802, r25803) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r25806, r, MPFR_RNDN);
        mpfr_mul(r25807, r25806, r25806, MPFR_RNDN);
        mpfr_div(r25808, r25805, r25807, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25811, r25809, r25810, MPFR_RNDN);
        mpfr_add(r25812, r25808, r25811, MPFR_RNDN);
        mpfr_set_d(r25813, w, MPFR_RNDN);
        mpfr_mul(r25814, r25813, r25806, MPFR_RNDN);
        mpfr_mul(r25815, r25814, r25814, MPFR_RNDN);
        ;
        mpfr_mul(r25817, r25815, r25816, MPFR_RNDN);
        mpfr_div(r25818, r25815, r25802, MPFR_RNDN);
        ;
        mpfr_div(r25820, r25819, r25802, MPFR_RNDN);
        mpfr_add(r25821, r25819, r25820, MPFR_RNDN);
        mpfr_mul(r25822, r25818, r25821, MPFR_RNDN);
        mpfr_sub(r25823, r25817, r25822, MPFR_RNDN);
        mpfr_sub(r25824, r25812, r25823, MPFR_RNDN);
        ;
        mpfr_set_si(r25826, mpfr_cmp(r25802, r25825) <= 0, MPFR_RNDN);
        mpfr_div(r25827, r25805, r25806, MPFR_RNDN);
        mpfr_div(r25828, r25827, r25806, MPFR_RNDN);
        mpfr_add(r25829, r25809, r25828, MPFR_RNDN);
        mpfr_mul(r25830, r25806, r25813, MPFR_RNDN);
        mpfr_mul(r25831, r25830, r25830, MPFR_RNDN);
        mpfr_mul(r25832, r25819, r25831, MPFR_RNDN);
        ;
        mpfr_mul(r25834, r25802, r25802, MPFR_RNDN);
        mpfr_sub(r25835, r25833, r25834, MPFR_RNDN);
        mpfr_add(r25836, r25802, r25802, MPFR_RNDN);
        mpfr_sub(r25837, r25809, r25836, MPFR_RNDN);
        mpfr_div(r25838, r25835, r25837, MPFR_RNDN);
        mpfr_div(r25839, r25832, r25838, MPFR_RNDN);
        mpfr_add(r25840, r25833, r25802, MPFR_RNDN);
        mpfr_mul(r25841, r25839, r25840, MPFR_RNDN);
        mpfr_sub(r25842, r25829, r25841, MPFR_RNDN);
        mpfr_sub(r25843, r25842, r25810, MPFR_RNDN);
        if (mpfr_get_si(r25826, MPFR_RNDN)) { mpfr_set(r25844, r25843, MPFR_RNDN); } else { mpfr_set(r25844, r25824, MPFR_RNDN); };
        if (mpfr_get_si(r25804, MPFR_RNDN)) { mpfr_set(r25845, r25824, MPFR_RNDN); } else { mpfr_set(r25845, r25844, MPFR_RNDN); };
        return mpfr_get_d(r25845, MPFR_RNDN);
}

static mpfr_t r25846, r25847, r25848, r25849, r25850, r25851, r25852, r25853, r25854, r25855, r25856, r25857, r25858, r25859, r25860, r25861, r25862, r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871, r25872, r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25846);
        mpfr_init_set_str(r25847, "-2.526418831017353e+31", 10, MPFR_RNDN);
        mpfr_init(r25848);
        mpfr_init_set_str(r25849, "2", 10, MPFR_RNDN);
        mpfr_init(r25850);
        mpfr_init(r25851);
        mpfr_init(r25852);
        mpfr_init_set_str(r25853, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r25854, "4.5", 10, MPFR_RNDN);
        mpfr_init(r25855);
        mpfr_init(r25856);
        mpfr_init(r25857);
        mpfr_init(r25858);
        mpfr_init(r25859);
        mpfr_init_set_str(r25860, "0.25", 10, MPFR_RNDN);
        mpfr_init(r25861);
        mpfr_init(r25862);
        mpfr_init_set_str(r25863, "0.125", 10, MPFR_RNDN);
        mpfr_init(r25864);
        mpfr_init(r25865);
        mpfr_init(r25866);
        mpfr_init(r25867);
        mpfr_init(r25868);
        mpfr_init_set_str(r25869, "4.907722326036034e-19", 10, MPFR_RNDN);
        mpfr_init(r25870);
        mpfr_init(r25871);
        mpfr_init(r25872);
        mpfr_init(r25873);
        mpfr_init(r25874);
        mpfr_init(r25875);
        mpfr_init(r25876);
        mpfr_init_set_str(r25877, "1", 10, MPFR_RNDN);
        mpfr_init(r25878);
        mpfr_init(r25879);
        mpfr_init(r25880);
        mpfr_init(r25881);
        mpfr_init(r25882);
        mpfr_init(r25883);
        mpfr_init(r25884);
        mpfr_init(r25885);
        mpfr_init(r25886);
        mpfr_init(r25887);
        mpfr_init(r25888);
        mpfr_init(r25889);
}

double f_dm(double v, double w, double r) {
        mpfr_set_d(r25846, v, MPFR_RNDN);
        ;
        mpfr_set_si(r25848, mpfr_cmp(r25846, r25847) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r25850, r, MPFR_RNDN);
        mpfr_mul(r25851, r25850, r25850, MPFR_RNDN);
        mpfr_div(r25852, r25849, r25851, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r25855, r25853, r25854, MPFR_RNDN);
        mpfr_add(r25856, r25852, r25855, MPFR_RNDN);
        mpfr_set_d(r25857, w, MPFR_RNDN);
        mpfr_mul(r25858, r25857, r25850, MPFR_RNDN);
        mpfr_mul(r25859, r25858, r25858, MPFR_RNDN);
        ;
        mpfr_mul(r25861, r25859, r25860, MPFR_RNDN);
        mpfr_div(r25862, r25859, r25846, MPFR_RNDN);
        ;
        mpfr_div(r25864, r25863, r25846, MPFR_RNDN);
        mpfr_add(r25865, r25863, r25864, MPFR_RNDN);
        mpfr_mul(r25866, r25862, r25865, MPFR_RNDN);
        mpfr_sub(r25867, r25861, r25866, MPFR_RNDN);
        mpfr_sub(r25868, r25856, r25867, MPFR_RNDN);
        ;
        mpfr_set_si(r25870, mpfr_cmp(r25846, r25869) <= 0, MPFR_RNDN);
        mpfr_div(r25871, r25849, r25850, MPFR_RNDN);
        mpfr_div(r25872, r25871, r25850, MPFR_RNDN);
        mpfr_add(r25873, r25853, r25872, MPFR_RNDN);
        mpfr_mul(r25874, r25850, r25857, MPFR_RNDN);
        mpfr_mul(r25875, r25874, r25874, MPFR_RNDN);
        mpfr_mul(r25876, r25863, r25875, MPFR_RNDN);
        ;
        mpfr_mul(r25878, r25846, r25846, MPFR_RNDN);
        mpfr_sub(r25879, r25877, r25878, MPFR_RNDN);
        mpfr_add(r25880, r25846, r25846, MPFR_RNDN);
        mpfr_sub(r25881, r25853, r25880, MPFR_RNDN);
        mpfr_div(r25882, r25879, r25881, MPFR_RNDN);
        mpfr_div(r25883, r25876, r25882, MPFR_RNDN);
        mpfr_add(r25884, r25877, r25846, MPFR_RNDN);
        mpfr_mul(r25885, r25883, r25884, MPFR_RNDN);
        mpfr_sub(r25886, r25873, r25885, MPFR_RNDN);
        mpfr_sub(r25887, r25886, r25854, MPFR_RNDN);
        if (mpfr_get_si(r25870, MPFR_RNDN)) { mpfr_set(r25888, r25887, MPFR_RNDN); } else { mpfr_set(r25888, r25868, MPFR_RNDN); };
        if (mpfr_get_si(r25848, MPFR_RNDN)) { mpfr_set(r25889, r25868, MPFR_RNDN); } else { mpfr_set(r25889, r25888, MPFR_RNDN); };
        return mpfr_get_d(r25889, MPFR_RNDN);
}