#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 r20239 = 3;
        float r20240 = 2;
        float r20241 = r;
        float r20242 = r20241 * r20241;
        float r20243 = r20240 / r20242;
        float r20244 = r20239 + r20243;
        float r20245 = 0.125;
        float r20246 = v;
        float r20247 = r20240 * r20246;
        float r20248 = r20239 - r20247;
        float r20249 = r20245 * r20248;
        float r20250 = w;
        float r20251 = r20250 * r20250;
        float r20252 = r20251 * r20241;
        float r20253 = r20252 * r20241;
        float r20254 = r20249 * r20253;
        float r20255 = 1;
        float r20256 = r20255 - r20246;
        float r20257 = r20254 / r20256;
        float r20258 = r20244 - r20257;
        float r20259 = 4.5;
        float r20260 = r20258 - r20259;
        return r20260;
}

double f_id(double v, double w, double r) {
        double r20261 = 3;
        double r20262 = 2;
        double r20263 = r;
        double r20264 = r20263 * r20263;
        double r20265 = r20262 / r20264;
        double r20266 = r20261 + r20265;
        double r20267 = 0.125;
        double r20268 = v;
        double r20269 = r20262 * r20268;
        double r20270 = r20261 - r20269;
        double r20271 = r20267 * r20270;
        double r20272 = w;
        double r20273 = r20272 * r20272;
        double r20274 = r20273 * r20263;
        double r20275 = r20274 * r20263;
        double r20276 = r20271 * r20275;
        double r20277 = 1;
        double r20278 = r20277 - r20268;
        double r20279 = r20276 / r20278;
        double r20280 = r20266 - r20279;
        double r20281 = 4.5;
        double r20282 = r20280 - r20281;
        return r20282;
}


double f_of(float v, float w, float r) {
        float r20283 = v;
        float r20284 = -8.982586312012436e+76;
        bool r20285 = r20283 <= r20284;
        float r20286 = 3;
        float r20287 = 4.5;
        float r20288 = r20286 - r20287;
        float r20289 = 2;
        float r20290 = r;
        float r20291 = r20289 / r20290;
        float r20292 = r20291 / r20290;
        float r20293 = r20288 + r20292;
        float r20294 = 0.25;
        float r20295 = w;
        float r20296 = r20295 * r20290;
        float r20297 = r20296 * r20296;
        float r20298 = r20294 * r20297;
        float r20299 = r20293 - r20298;
        float r20300 = 0.125;
        float r20301 = r20300 / r20283;
        float r20302 = r20301 + r20300;
        float r20303 = r20302 * r20297;
        float r20304 = r20303 / r20283;
        float r20305 = r20299 + r20304;
        float r20306 = 2.9582720470002017e-15;
        bool r20307 = r20283 <= r20306;
        float r20308 = 1;
        float r20309 = r20283 * r20283;
        float r20310 = r20308 + r20309;
        float r20311 = r20283 + r20283;
        float r20312 = r20286 - r20311;
        float r20313 = r20310 * r20312;
        float r20314 = r20308 + r20283;
        float r20315 = r20314 * r20300;
        float r20316 = r20315 * r20297;
        float r20317 = r20309 * r20309;
        float r20318 = r20308 - r20317;
        float r20319 = r20316 / r20318;
        float r20320 = r20313 * r20319;
        float r20321 = r20293 - r20320;
        float r20322 = r20307 ? r20321 : r20305;
        float r20323 = r20285 ? r20305 : r20322;
        return r20323;
}

double f_od(double v, double w, double r) {
        double r20324 = v;
        double r20325 = -8.982586312012436e+76;
        bool r20326 = r20324 <= r20325;
        double r20327 = 3;
        double r20328 = 4.5;
        double r20329 = r20327 - r20328;
        double r20330 = 2;
        double r20331 = r;
        double r20332 = r20330 / r20331;
        double r20333 = r20332 / r20331;
        double r20334 = r20329 + r20333;
        double r20335 = 0.25;
        double r20336 = w;
        double r20337 = r20336 * r20331;
        double r20338 = r20337 * r20337;
        double r20339 = r20335 * r20338;
        double r20340 = r20334 - r20339;
        double r20341 = 0.125;
        double r20342 = r20341 / r20324;
        double r20343 = r20342 + r20341;
        double r20344 = r20343 * r20338;
        double r20345 = r20344 / r20324;
        double r20346 = r20340 + r20345;
        double r20347 = 2.9582720470002017e-15;
        bool r20348 = r20324 <= r20347;
        double r20349 = 1;
        double r20350 = r20324 * r20324;
        double r20351 = r20349 + r20350;
        double r20352 = r20324 + r20324;
        double r20353 = r20327 - r20352;
        double r20354 = r20351 * r20353;
        double r20355 = r20349 + r20324;
        double r20356 = r20355 * r20341;
        double r20357 = r20356 * r20338;
        double r20358 = r20350 * r20350;
        double r20359 = r20349 - r20358;
        double r20360 = r20357 / r20359;
        double r20361 = r20354 * r20360;
        double r20362 = r20334 - r20361;
        double r20363 = r20348 ? r20362 : r20346;
        double r20364 = r20326 ? r20346 : r20363;
        return r20364;
}

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 r20365, r20366, r20367, r20368, r20369, r20370, r20371, r20372, r20373, r20374, r20375, r20376, r20377, r20378, r20379, r20380, r20381, r20382, r20383, r20384, r20385, r20386;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20365, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r20366, "2", 10, MPFR_RNDN);
        mpfr_init(r20367);
        mpfr_init(r20368);
        mpfr_init(r20369);
        mpfr_init(r20370);
        mpfr_init_set_str(r20371, "0.125", 10, MPFR_RNDN);
        mpfr_init(r20372);
        mpfr_init(r20373);
        mpfr_init(r20374);
        mpfr_init(r20375);
        mpfr_init(r20376);
        mpfr_init(r20377);
        mpfr_init(r20378);
        mpfr_init(r20379);
        mpfr_init(r20380);
        mpfr_init_set_str(r20381, "1", 10, MPFR_RNDN);
        mpfr_init(r20382);
        mpfr_init(r20383);
        mpfr_init(r20384);
        mpfr_init_set_str(r20385, "4.5", 10, MPFR_RNDN);
        mpfr_init(r20386);
}

double f_im(double v, double w, double r) {
        ;
        ;
        mpfr_set_d(r20367, r, MPFR_RNDN);
        mpfr_mul(r20368, r20367, r20367, MPFR_RNDN);
        mpfr_div(r20369, r20366, r20368, MPFR_RNDN);
        mpfr_add(r20370, r20365, r20369, MPFR_RNDN);
        ;
        mpfr_set_d(r20372, v, MPFR_RNDN);
        mpfr_mul(r20373, r20366, r20372, MPFR_RNDN);
        mpfr_sub(r20374, r20365, r20373, MPFR_RNDN);
        mpfr_mul(r20375, r20371, r20374, MPFR_RNDN);
        mpfr_set_d(r20376, w, MPFR_RNDN);
        mpfr_mul(r20377, r20376, r20376, MPFR_RNDN);
        mpfr_mul(r20378, r20377, r20367, MPFR_RNDN);
        mpfr_mul(r20379, r20378, r20367, MPFR_RNDN);
        mpfr_mul(r20380, r20375, r20379, MPFR_RNDN);
        ;
        mpfr_sub(r20382, r20381, r20372, MPFR_RNDN);
        mpfr_div(r20383, r20380, r20382, MPFR_RNDN);
        mpfr_sub(r20384, r20370, r20383, MPFR_RNDN);
        ;
        mpfr_sub(r20386, r20384, r20385, MPFR_RNDN);
        return mpfr_get_d(r20386, MPFR_RNDN);
}

static mpfr_t r20387, r20388, r20389, r20390, r20391, r20392, r20393, r20394, r20395, r20396, r20397, r20398, r20399, r20400, r20401, r20402, r20403, r20404, r20405, r20406, r20407, r20408, r20409, r20410, r20411, r20412, r20413, r20414, r20415, r20416, r20417, r20418, r20419, r20420, r20421, r20422, r20423, r20424, r20425, r20426, r20427;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20387);
        mpfr_init_set_str(r20388, "-8.982586312012436e+76", 10, MPFR_RNDN);
        mpfr_init(r20389);
        mpfr_init_set_str(r20390, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r20391, "4.5", 10, MPFR_RNDN);
        mpfr_init(r20392);
        mpfr_init_set_str(r20393, "2", 10, MPFR_RNDN);
        mpfr_init(r20394);
        mpfr_init(r20395);
        mpfr_init(r20396);
        mpfr_init(r20397);
        mpfr_init_set_str(r20398, "0.25", 10, MPFR_RNDN);
        mpfr_init(r20399);
        mpfr_init(r20400);
        mpfr_init(r20401);
        mpfr_init(r20402);
        mpfr_init(r20403);
        mpfr_init_set_str(r20404, "0.125", 10, MPFR_RNDN);
        mpfr_init(r20405);
        mpfr_init(r20406);
        mpfr_init(r20407);
        mpfr_init(r20408);
        mpfr_init(r20409);
        mpfr_init_set_str(r20410, "2.9582720470002017e-15", 10, MPFR_RNDN);
        mpfr_init(r20411);
        mpfr_init_set_str(r20412, "1", 10, MPFR_RNDN);
        mpfr_init(r20413);
        mpfr_init(r20414);
        mpfr_init(r20415);
        mpfr_init(r20416);
        mpfr_init(r20417);
        mpfr_init(r20418);
        mpfr_init(r20419);
        mpfr_init(r20420);
        mpfr_init(r20421);
        mpfr_init(r20422);
        mpfr_init(r20423);
        mpfr_init(r20424);
        mpfr_init(r20425);
        mpfr_init(r20426);
        mpfr_init(r20427);
}

double f_fm(double v, double w, double r) {
        mpfr_set_d(r20387, v, MPFR_RNDN);
        ;
        mpfr_set_si(r20389, mpfr_cmp(r20387, r20388) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r20392, r20390, r20391, MPFR_RNDN);
        ;
        mpfr_set_d(r20394, r, MPFR_RNDN);
        mpfr_div(r20395, r20393, r20394, MPFR_RNDN);
        mpfr_div(r20396, r20395, r20394, MPFR_RNDN);
        mpfr_add(r20397, r20392, r20396, MPFR_RNDN);
        ;
        mpfr_set_d(r20399, w, MPFR_RNDN);
        mpfr_mul(r20400, r20399, r20394, MPFR_RNDN);
        mpfr_mul(r20401, r20400, r20400, MPFR_RNDN);
        mpfr_mul(r20402, r20398, r20401, MPFR_RNDN);
        mpfr_sub(r20403, r20397, r20402, MPFR_RNDN);
        ;
        mpfr_div(r20405, r20404, r20387, MPFR_RNDN);
        mpfr_add(r20406, r20405, r20404, MPFR_RNDN);
        mpfr_mul(r20407, r20406, r20401, MPFR_RNDN);
        mpfr_div(r20408, r20407, r20387, MPFR_RNDN);
        mpfr_add(r20409, r20403, r20408, MPFR_RNDN);
        ;
        mpfr_set_si(r20411, mpfr_cmp(r20387, r20410) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r20413, r20387, r20387, MPFR_RNDN);
        mpfr_add(r20414, r20412, r20413, MPFR_RNDN);
        mpfr_add(r20415, r20387, r20387, MPFR_RNDN);
        mpfr_sub(r20416, r20390, r20415, MPFR_RNDN);
        mpfr_mul(r20417, r20414, r20416, MPFR_RNDN);
        mpfr_add(r20418, r20412, r20387, MPFR_RNDN);
        mpfr_mul(r20419, r20418, r20404, MPFR_RNDN);
        mpfr_mul(r20420, r20419, r20401, MPFR_RNDN);
        mpfr_mul(r20421, r20413, r20413, MPFR_RNDN);
        mpfr_sub(r20422, r20412, r20421, MPFR_RNDN);
        mpfr_div(r20423, r20420, r20422, MPFR_RNDN);
        mpfr_mul(r20424, r20417, r20423, MPFR_RNDN);
        mpfr_sub(r20425, r20397, r20424, MPFR_RNDN);
        if (mpfr_get_si(r20411, MPFR_RNDN)) { mpfr_set(r20426, r20425, MPFR_RNDN); } else { mpfr_set(r20426, r20409, MPFR_RNDN); };
        if (mpfr_get_si(r20389, MPFR_RNDN)) { mpfr_set(r20427, r20409, MPFR_RNDN); } else { mpfr_set(r20427, r20426, MPFR_RNDN); };
        return mpfr_get_d(r20427, MPFR_RNDN);
}

static mpfr_t r20428, r20429, r20430, r20431, r20432, r20433, r20434, r20435, r20436, r20437, r20438, r20439, r20440, r20441, r20442, r20443, r20444, r20445, r20446, r20447, r20448, r20449, r20450, r20451, r20452, r20453, r20454, r20455, r20456, r20457, r20458, r20459, r20460, r20461, r20462, r20463, r20464, r20465, r20466, r20467, r20468;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20428);
        mpfr_init_set_str(r20429, "-8.982586312012436e+76", 10, MPFR_RNDN);
        mpfr_init(r20430);
        mpfr_init_set_str(r20431, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r20432, "4.5", 10, MPFR_RNDN);
        mpfr_init(r20433);
        mpfr_init_set_str(r20434, "2", 10, MPFR_RNDN);
        mpfr_init(r20435);
        mpfr_init(r20436);
        mpfr_init(r20437);
        mpfr_init(r20438);
        mpfr_init_set_str(r20439, "0.25", 10, MPFR_RNDN);
        mpfr_init(r20440);
        mpfr_init(r20441);
        mpfr_init(r20442);
        mpfr_init(r20443);
        mpfr_init(r20444);
        mpfr_init_set_str(r20445, "0.125", 10, MPFR_RNDN);
        mpfr_init(r20446);
        mpfr_init(r20447);
        mpfr_init(r20448);
        mpfr_init(r20449);
        mpfr_init(r20450);
        mpfr_init_set_str(r20451, "2.9582720470002017e-15", 10, MPFR_RNDN);
        mpfr_init(r20452);
        mpfr_init_set_str(r20453, "1", 10, MPFR_RNDN);
        mpfr_init(r20454);
        mpfr_init(r20455);
        mpfr_init(r20456);
        mpfr_init(r20457);
        mpfr_init(r20458);
        mpfr_init(r20459);
        mpfr_init(r20460);
        mpfr_init(r20461);
        mpfr_init(r20462);
        mpfr_init(r20463);
        mpfr_init(r20464);
        mpfr_init(r20465);
        mpfr_init(r20466);
        mpfr_init(r20467);
        mpfr_init(r20468);
}

double f_dm(double v, double w, double r) {
        mpfr_set_d(r20428, v, MPFR_RNDN);
        ;
        mpfr_set_si(r20430, mpfr_cmp(r20428, r20429) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r20433, r20431, r20432, MPFR_RNDN);
        ;
        mpfr_set_d(r20435, r, MPFR_RNDN);
        mpfr_div(r20436, r20434, r20435, MPFR_RNDN);
        mpfr_div(r20437, r20436, r20435, MPFR_RNDN);
        mpfr_add(r20438, r20433, r20437, MPFR_RNDN);
        ;
        mpfr_set_d(r20440, w, MPFR_RNDN);
        mpfr_mul(r20441, r20440, r20435, MPFR_RNDN);
        mpfr_mul(r20442, r20441, r20441, MPFR_RNDN);
        mpfr_mul(r20443, r20439, r20442, MPFR_RNDN);
        mpfr_sub(r20444, r20438, r20443, MPFR_RNDN);
        ;
        mpfr_div(r20446, r20445, r20428, MPFR_RNDN);
        mpfr_add(r20447, r20446, r20445, MPFR_RNDN);
        mpfr_mul(r20448, r20447, r20442, MPFR_RNDN);
        mpfr_div(r20449, r20448, r20428, MPFR_RNDN);
        mpfr_add(r20450, r20444, r20449, MPFR_RNDN);
        ;
        mpfr_set_si(r20452, mpfr_cmp(r20428, r20451) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r20454, r20428, r20428, MPFR_RNDN);
        mpfr_add(r20455, r20453, r20454, MPFR_RNDN);
        mpfr_add(r20456, r20428, r20428, MPFR_RNDN);
        mpfr_sub(r20457, r20431, r20456, MPFR_RNDN);
        mpfr_mul(r20458, r20455, r20457, MPFR_RNDN);
        mpfr_add(r20459, r20453, r20428, MPFR_RNDN);
        mpfr_mul(r20460, r20459, r20445, MPFR_RNDN);
        mpfr_mul(r20461, r20460, r20442, MPFR_RNDN);
        mpfr_mul(r20462, r20454, r20454, MPFR_RNDN);
        mpfr_sub(r20463, r20453, r20462, MPFR_RNDN);
        mpfr_div(r20464, r20461, r20463, MPFR_RNDN);
        mpfr_mul(r20465, r20458, r20464, MPFR_RNDN);
        mpfr_sub(r20466, r20438, r20465, MPFR_RNDN);
        if (mpfr_get_si(r20452, MPFR_RNDN)) { mpfr_set(r20467, r20466, MPFR_RNDN); } else { mpfr_set(r20467, r20450, MPFR_RNDN); };
        if (mpfr_get_si(r20430, MPFR_RNDN)) { mpfr_set(r20468, r20450, MPFR_RNDN); } else { mpfr_set(r20468, r20467, MPFR_RNDN); };
        return mpfr_get_d(r20468, MPFR_RNDN);
}

