#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 r19629 = 3;
        float r19630 = 2;
        float r19631 = r;
        float r19632 = r19631 * r19631;
        float r19633 = r19630 / r19632;
        float r19634 = r19629 + r19633;
        float r19635 = 0.125;
        float r19636 = v;
        float r19637 = r19630 * r19636;
        float r19638 = r19629 - r19637;
        float r19639 = r19635 * r19638;
        float r19640 = w;
        float r19641 = r19640 * r19640;
        float r19642 = r19641 * r19631;
        float r19643 = r19642 * r19631;
        float r19644 = r19639 * r19643;
        float r19645 = 1;
        float r19646 = r19645 - r19636;
        float r19647 = r19644 / r19646;
        float r19648 = r19634 - r19647;
        float r19649 = 4.5;
        float r19650 = r19648 - r19649;
        return r19650;
}

double f_id(double v, double w, double r) {
        double r19651 = 3;
        double r19652 = 2;
        double r19653 = r;
        double r19654 = r19653 * r19653;
        double r19655 = r19652 / r19654;
        double r19656 = r19651 + r19655;
        double r19657 = 0.125;
        double r19658 = v;
        double r19659 = r19652 * r19658;
        double r19660 = r19651 - r19659;
        double r19661 = r19657 * r19660;
        double r19662 = w;
        double r19663 = r19662 * r19662;
        double r19664 = r19663 * r19653;
        double r19665 = r19664 * r19653;
        double r19666 = r19661 * r19665;
        double r19667 = 1;
        double r19668 = r19667 - r19658;
        double r19669 = r19666 / r19668;
        double r19670 = r19656 - r19669;
        double r19671 = 4.5;
        double r19672 = r19670 - r19671;
        return r19672;
}


double f_of(float v, float w, float r) {
        float r19673 = 2;
        float r19674 = r;
        float r19675 = r19673 / r19674;
        float r19676 = r19675 / r19674;
        float r19677 = 3;
        float r19678 = 4.5;
        float r19679 = r19677 - r19678;
        float r19680 = r19676 + r19679;
        float r19681 = w;
        float r19682 = r19681 * r19674;
        float r19683 = fabs(r19682);
        float r19684 = sqrt(r19683);
        float r19685 = r19683 * r19684;
        float r19686 = r19685 * r19684;
        float r19687 = 1;
        float r19688 = v;
        float r19689 = r19687 - r19688;
        float r19690 = 0.375;
        float r19691 = 0.25;
        float r19692 = r19688 * r19691;
        float r19693 = r19690 - r19692;
        float r19694 = r19689 / r19693;
        float r19695 = r19686 / r19694;
        float r19696 = r19680 - r19695;
        return r19696;
}

double f_od(double v, double w, double r) {
        double r19697 = 2;
        double r19698 = r;
        double r19699 = r19697 / r19698;
        double r19700 = r19699 / r19698;
        double r19701 = 3;
        double r19702 = 4.5;
        double r19703 = r19701 - r19702;
        double r19704 = r19700 + r19703;
        double r19705 = w;
        double r19706 = r19705 * r19698;
        double r19707 = fabs(r19706);
        double r19708 = sqrt(r19707);
        double r19709 = r19707 * r19708;
        double r19710 = r19709 * r19708;
        double r19711 = 1;
        double r19712 = v;
        double r19713 = r19711 - r19712;
        double r19714 = 0.375;
        double r19715 = 0.25;
        double r19716 = r19712 * r19715;
        double r19717 = r19714 - r19716;
        double r19718 = r19713 / r19717;
        double r19719 = r19710 / r19718;
        double r19720 = r19704 - r19719;
        return r19720;
}

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 r19721, r19722, r19723, r19724, r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r19721, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19722, "2", 10, MPFR_RNDN);
        mpfr_init(r19723);
        mpfr_init(r19724);
        mpfr_init(r19725);
        mpfr_init(r19726);
        mpfr_init_set_str(r19727, "0.125", 10, MPFR_RNDN);
        mpfr_init(r19728);
        mpfr_init(r19729);
        mpfr_init(r19730);
        mpfr_init(r19731);
        mpfr_init(r19732);
        mpfr_init(r19733);
        mpfr_init(r19734);
        mpfr_init(r19735);
        mpfr_init(r19736);
        mpfr_init_set_str(r19737, "1", 10, MPFR_RNDN);
        mpfr_init(r19738);
        mpfr_init(r19739);
        mpfr_init(r19740);
        mpfr_init_set_str(r19741, "4.5", 10, MPFR_RNDN);
        mpfr_init(r19742);
}

double f_im(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r19723, r, MPFR_RNDN);
        mpfr_mul(r19724, r19723, r19723, MPFR_RNDN);
        mpfr_div(r19725, r19722, r19724, MPFR_RNDN);
        mpfr_add(r19726, r19721, r19725, MPFR_RNDN);
        ;
        mpfr_set_d(r19728, v, MPFR_RNDN);
        mpfr_mul(r19729, r19722, r19728, MPFR_RNDN);
        mpfr_sub(r19730, r19721, r19729, MPFR_RNDN);
        mpfr_mul(r19731, r19727, r19730, MPFR_RNDN);
        mpfr_set_d(r19732, w, MPFR_RNDN);
        mpfr_mul(r19733, r19732, r19732, MPFR_RNDN);
        mpfr_mul(r19734, r19733, r19723, MPFR_RNDN);
        mpfr_mul(r19735, r19734, r19723, MPFR_RNDN);
        mpfr_mul(r19736, r19731, r19735, MPFR_RNDN);
        ;
        mpfr_sub(r19738, r19737, r19728, MPFR_RNDN);
        mpfr_div(r19739, r19736, r19738, MPFR_RNDN);
        mpfr_sub(r19740, r19726, r19739, MPFR_RNDN);
        ;
        mpfr_sub(r19742, r19740, r19741, MPFR_RNDN);
        return mpfr_get_d(r19742, MPFR_RNDN);
}

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, r19765, r19766;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r19743, "2", 10, MPFR_RNDN);
        mpfr_init(r19744);
        mpfr_init(r19745);
        mpfr_init(r19746);
        mpfr_init_set_str(r19747, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19748, "4.5", 10, MPFR_RNDN);
        mpfr_init(r19749);
        mpfr_init(r19750);
        mpfr_init(r19751);
        mpfr_init(r19752);
        mpfr_init(r19753);
        mpfr_init(r19754);
        mpfr_init(r19755);
        mpfr_init(r19756);
        mpfr_init_set_str(r19757, "1", 10, MPFR_RNDN);
        mpfr_init(r19758);
        mpfr_init(r19759);
        mpfr_init_set_str(r19760, "0.375", 10, MPFR_RNDN);
        mpfr_init_set_str(r19761, "0.25", 10, MPFR_RNDN);
        mpfr_init(r19762);
        mpfr_init(r19763);
        mpfr_init(r19764);
        mpfr_init(r19765);
        mpfr_init(r19766);
}

double f_fm(double v, double w, double r) {
        ;
        mpfr_set_d(r19744, r, MPFR_RNDN);
        mpfr_div(r19745, r19743, r19744, MPFR_RNDN);
        mpfr_div(r19746, r19745, r19744, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r19749, r19747, r19748, MPFR_RNDN);
        mpfr_add(r19750, r19746, r19749, MPFR_RNDN);
        mpfr_set_d(r19751, w, MPFR_RNDN);
        mpfr_mul(r19752, r19751, r19744, MPFR_RNDN);
        mpfr_abs(r19753, r19752, MPFR_RNDN);
        mpfr_sqrt(r19754, r19753, MPFR_RNDN);
        mpfr_mul(r19755, r19753, r19754, MPFR_RNDN);
        mpfr_mul(r19756, r19755, r19754, MPFR_RNDN);
        ;
        mpfr_set_d(r19758, v, MPFR_RNDN);
        mpfr_sub(r19759, r19757, r19758, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r19762, r19758, r19761, MPFR_RNDN);
        mpfr_sub(r19763, r19760, r19762, MPFR_RNDN);
        mpfr_div(r19764, r19759, r19763, MPFR_RNDN);
        mpfr_div(r19765, r19756, r19764, MPFR_RNDN);
        mpfr_sub(r19766, r19750, r19765, MPFR_RNDN);
        return mpfr_get_d(r19766, MPFR_RNDN);
}

static mpfr_t r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789, r19790;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r19767, "2", 10, MPFR_RNDN);
        mpfr_init(r19768);
        mpfr_init(r19769);
        mpfr_init(r19770);
        mpfr_init_set_str(r19771, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19772, "4.5", 10, MPFR_RNDN);
        mpfr_init(r19773);
        mpfr_init(r19774);
        mpfr_init(r19775);
        mpfr_init(r19776);
        mpfr_init(r19777);
        mpfr_init(r19778);
        mpfr_init(r19779);
        mpfr_init(r19780);
        mpfr_init_set_str(r19781, "1", 10, MPFR_RNDN);
        mpfr_init(r19782);
        mpfr_init(r19783);
        mpfr_init_set_str(r19784, "0.375", 10, MPFR_RNDN);
        mpfr_init_set_str(r19785, "0.25", 10, MPFR_RNDN);
        mpfr_init(r19786);
        mpfr_init(r19787);
        mpfr_init(r19788);
        mpfr_init(r19789);
        mpfr_init(r19790);
}

double f_dm(double v, double w, double r) {
        ;
        mpfr_set_d(r19768, r, MPFR_RNDN);
        mpfr_div(r19769, r19767, r19768, MPFR_RNDN);
        mpfr_div(r19770, r19769, r19768, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r19773, r19771, r19772, MPFR_RNDN);
        mpfr_add(r19774, r19770, r19773, MPFR_RNDN);
        mpfr_set_d(r19775, w, MPFR_RNDN);
        mpfr_mul(r19776, r19775, r19768, MPFR_RNDN);
        mpfr_abs(r19777, r19776, MPFR_RNDN);
        mpfr_sqrt(r19778, r19777, MPFR_RNDN);
        mpfr_mul(r19779, r19777, r19778, MPFR_RNDN);
        mpfr_mul(r19780, r19779, r19778, MPFR_RNDN);
        ;
        mpfr_set_d(r19782, v, MPFR_RNDN);
        mpfr_sub(r19783, r19781, r19782, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r19786, r19782, r19785, MPFR_RNDN);
        mpfr_sub(r19787, r19784, r19786, MPFR_RNDN);
        mpfr_div(r19788, r19783, r19787, MPFR_RNDN);
        mpfr_div(r19789, r19780, r19788, MPFR_RNDN);
        mpfr_sub(r19790, r19774, r19789, MPFR_RNDN);
        return mpfr_get_d(r19790, MPFR_RNDN);
}

