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

char *name = "turbine2";

double f_if(float v, float w, float r) {
        float r20602 = 6;
        float r20603 = v;
        float r20604 = r20602 * r20603;
        float r20605 = 0.5;
        float r20606 = r20605 * r20603;
        float r20607 = w;
        float r20608 = r20607 * r20607;
        float r20609 = r;
        float r20610 = r20608 * r20609;
        float r20611 = r20610 * r20609;
        float r20612 = r20606 * r20611;
        float r20613 = 1;
        float r20614 = r20613 - r20603;
        float r20615 = r20612 / r20614;
        float r20616 = r20604 - r20615;
        float r20617 = 2.5;
        float r20618 = r20616 - r20617;
        return r20618;
}

double f_id(double v, double w, double r) {
        double r20619 = 6;
        double r20620 = v;
        double r20621 = r20619 * r20620;
        double r20622 = 0.5;
        double r20623 = r20622 * r20620;
        double r20624 = w;
        double r20625 = r20624 * r20624;
        double r20626 = r;
        double r20627 = r20625 * r20626;
        double r20628 = r20627 * r20626;
        double r20629 = r20623 * r20628;
        double r20630 = 1;
        double r20631 = r20630 - r20620;
        double r20632 = r20629 / r20631;
        double r20633 = r20621 - r20632;
        double r20634 = 2.5;
        double r20635 = r20633 - r20634;
        return r20635;
}


double f_of(float v, float w, float r) {
        float r20636 = v;
        float r20637 = exp(r20636);
        float r20638 = 6;
        float r20639 = pow(r20637, r20638);
        float r20640 = 0.5;
        float r20641 = w;
        float r20642 = r20640 * r20641;
        float r20643 = r20642 * r20641;
        float r20644 = 1;
        float r20645 = r20644 - r20636;
        float r20646 = r;
        float r20647 = r20646 * r20646;
        float r20648 = r20645 / r20647;
        float r20649 = r20643 / r20648;
        float r20650 = pow(r20637, r20649);
        float r20651 = r20639 / r20650;
        float r20652 = log(r20651);
        float r20653 = 2.5;
        float r20654 = r20652 - r20653;
        return r20654;
}

double f_od(double v, double w, double r) {
        double r20655 = v;
        double r20656 = exp(r20655);
        double r20657 = 6;
        double r20658 = pow(r20656, r20657);
        double r20659 = 0.5;
        double r20660 = w;
        double r20661 = r20659 * r20660;
        double r20662 = r20661 * r20660;
        double r20663 = 1;
        double r20664 = r20663 - r20655;
        double r20665 = r;
        double r20666 = r20665 * r20665;
        double r20667 = r20664 / r20666;
        double r20668 = r20662 / r20667;
        double r20669 = pow(r20656, r20668);
        double r20670 = r20658 / r20669;
        double r20671 = log(r20670);
        double r20672 = 2.5;
        double r20673 = r20671 - r20672;
        return r20673;
}

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 r20674, r20675, r20676, r20677, r20678, r20679, r20680, r20681, r20682, r20683, r20684, r20685, r20686, r20687, r20688, r20689, r20690;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20674, "6", 10, MPFR_RNDN);
        mpfr_init(r20675);
        mpfr_init(r20676);
        mpfr_init_set_str(r20677, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20678);
        mpfr_init(r20679);
        mpfr_init(r20680);
        mpfr_init(r20681);
        mpfr_init(r20682);
        mpfr_init(r20683);
        mpfr_init(r20684);
        mpfr_init_set_str(r20685, "1", 10, MPFR_RNDN);
        mpfr_init(r20686);
        mpfr_init(r20687);
        mpfr_init(r20688);
        mpfr_init_set_str(r20689, "2.5", 10, MPFR_RNDN);
        mpfr_init(r20690);
}

double f_im(double v, double w, double r) {
        ;
        mpfr_set_d(r20675, v, MPFR_RNDN);
        mpfr_mul(r20676, r20674, r20675, MPFR_RNDN);
        ;
        mpfr_mul(r20678, r20677, r20675, MPFR_RNDN);
        mpfr_set_d(r20679, w, MPFR_RNDN);
        mpfr_mul(r20680, r20679, r20679, MPFR_RNDN);
        mpfr_set_d(r20681, r, MPFR_RNDN);
        mpfr_mul(r20682, r20680, r20681, MPFR_RNDN);
        mpfr_mul(r20683, r20682, r20681, MPFR_RNDN);
        mpfr_mul(r20684, r20678, r20683, MPFR_RNDN);
        ;
        mpfr_sub(r20686, r20685, r20675, MPFR_RNDN);
        mpfr_div(r20687, r20684, r20686, MPFR_RNDN);
        mpfr_sub(r20688, r20676, r20687, MPFR_RNDN);
        ;
        mpfr_sub(r20690, r20688, r20689, MPFR_RNDN);
        return mpfr_get_d(r20690, MPFR_RNDN);
}

static mpfr_t r20691, r20692, r20693, r20694, r20695, r20696, r20697, r20698, r20699, r20700, r20701, r20702, r20703, r20704, r20705, r20706, r20707, r20708, r20709;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20691);
        mpfr_init(r20692);
        mpfr_init_set_str(r20693, "6", 10, MPFR_RNDN);
        mpfr_init(r20694);
        mpfr_init_set_str(r20695, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20696);
        mpfr_init(r20697);
        mpfr_init(r20698);
        mpfr_init_set_str(r20699, "1", 10, MPFR_RNDN);
        mpfr_init(r20700);
        mpfr_init(r20701);
        mpfr_init(r20702);
        mpfr_init(r20703);
        mpfr_init(r20704);
        mpfr_init(r20705);
        mpfr_init(r20706);
        mpfr_init(r20707);
        mpfr_init_set_str(r20708, "2.5", 10, MPFR_RNDN);
        mpfr_init(r20709);
}

double f_fm(double v, double w, double r) {
        mpfr_set_d(r20691, v, MPFR_RNDN);
        mpfr_exp(r20692, r20691, MPFR_RNDN);
        ;
        mpfr_pow(r20694, r20692, r20693, MPFR_RNDN);
        ;
        mpfr_set_d(r20696, w, MPFR_RNDN);
        mpfr_mul(r20697, r20695, r20696, MPFR_RNDN);
        mpfr_mul(r20698, r20697, r20696, MPFR_RNDN);
        ;
        mpfr_sub(r20700, r20699, r20691, MPFR_RNDN);
        mpfr_set_d(r20701, r, MPFR_RNDN);
        mpfr_mul(r20702, r20701, r20701, MPFR_RNDN);
        mpfr_div(r20703, r20700, r20702, MPFR_RNDN);
        mpfr_div(r20704, r20698, r20703, MPFR_RNDN);
        mpfr_pow(r20705, r20692, r20704, MPFR_RNDN);
        mpfr_div(r20706, r20694, r20705, MPFR_RNDN);
        mpfr_log(r20707, r20706, MPFR_RNDN);
        ;
        mpfr_sub(r20709, r20707, r20708, MPFR_RNDN);
        return mpfr_get_d(r20709, MPFR_RNDN);
}

static mpfr_t r20710, r20711, r20712, r20713, r20714, r20715, r20716, r20717, r20718, r20719, r20720, r20721, r20722, r20723, r20724, r20725, r20726, r20727, r20728;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20710);
        mpfr_init(r20711);
        mpfr_init_set_str(r20712, "6", 10, MPFR_RNDN);
        mpfr_init(r20713);
        mpfr_init_set_str(r20714, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20715);
        mpfr_init(r20716);
        mpfr_init(r20717);
        mpfr_init_set_str(r20718, "1", 10, MPFR_RNDN);
        mpfr_init(r20719);
        mpfr_init(r20720);
        mpfr_init(r20721);
        mpfr_init(r20722);
        mpfr_init(r20723);
        mpfr_init(r20724);
        mpfr_init(r20725);
        mpfr_init(r20726);
        mpfr_init_set_str(r20727, "2.5", 10, MPFR_RNDN);
        mpfr_init(r20728);
}

double f_dm(double v, double w, double r) {
        mpfr_set_d(r20710, v, MPFR_RNDN);
        mpfr_exp(r20711, r20710, MPFR_RNDN);
        ;
        mpfr_pow(r20713, r20711, r20712, MPFR_RNDN);
        ;
        mpfr_set_d(r20715, w, MPFR_RNDN);
        mpfr_mul(r20716, r20714, r20715, MPFR_RNDN);
        mpfr_mul(r20717, r20716, r20715, MPFR_RNDN);
        ;
        mpfr_sub(r20719, r20718, r20710, MPFR_RNDN);
        mpfr_set_d(r20720, r, MPFR_RNDN);
        mpfr_mul(r20721, r20720, r20720, MPFR_RNDN);
        mpfr_div(r20722, r20719, r20721, MPFR_RNDN);
        mpfr_div(r20723, r20717, r20722, MPFR_RNDN);
        mpfr_pow(r20724, r20711, r20723, MPFR_RNDN);
        mpfr_div(r20725, r20713, r20724, MPFR_RNDN);
        mpfr_log(r20726, r20725, MPFR_RNDN);
        ;
        mpfr_sub(r20728, r20726, r20727, MPFR_RNDN);
        return mpfr_get_d(r20728, MPFR_RNDN);
}