#include #include #include #include #include char *name = "turbine3"; double f_if(float v, float w, float r) { float r20444 = 3; float r20445 = 2; float r20446 = r; float r20447 = r20446 * r20446; float r20448 = r20445 / r20447; float r20449 = r20444 - r20448; float r20450 = 0.125; float r20451 = 1; float r20452 = v; float r20453 = r20445 * r20452; float r20454 = r20451 + r20453; float r20455 = r20450 * r20454; float r20456 = w; float r20457 = r20456 * r20456; float r20458 = r20457 * r20446; float r20459 = r20458 * r20446; float r20460 = r20455 * r20459; float r20461 = r20451 - r20452; float r20462 = r20460 / r20461; float r20463 = r20449 - r20462; float r20464 = 0.5; float r20465 = r20463 - r20464; return r20465; } double f_id(double v, double w, double r) { double r20466 = 3; double r20467 = 2; double r20468 = r; double r20469 = r20468 * r20468; double r20470 = r20467 / r20469; double r20471 = r20466 - r20470; double r20472 = 0.125; double r20473 = 1; double r20474 = v; double r20475 = r20467 * r20474; double r20476 = r20473 + r20475; double r20477 = r20472 * r20476; double r20478 = w; double r20479 = r20478 * r20478; double r20480 = r20479 * r20468; double r20481 = r20480 * r20468; double r20482 = r20477 * r20481; double r20483 = r20473 - r20474; double r20484 = r20482 / r20483; double r20485 = r20471 - r20484; double r20486 = 0.5; double r20487 = r20485 - r20486; return r20487; } double f_of(float v, float w, float r) { float r20488 = 3; float r20489 = 2; float r20490 = r; float r20491 = r20490 * r20490; float r20492 = r20489 / r20491; float r20493 = r20488 - r20492; float r20494 = 0.125; float r20495 = 1; float r20496 = v; float r20497 = r20489 * r20496; float r20498 = r20495 + r20497; float r20499 = r20494 * r20498; float r20500 = r20499 / r20495; float r20501 = w; float r20502 = r20501 * r20501; float r20503 = r20502 * r20490; float r20504 = r20503 * r20490; float r20505 = r20495 - r20496; float r20506 = r20504 / r20505; float r20507 = r20500 * r20506; float r20508 = r20493 - r20507; float r20509 = 0.5; float r20510 = r20508 - r20509; return r20510; } double f_od(double v, double w, double r) { double r20511 = 3; double r20512 = 2; double r20513 = r; double r20514 = r20513 * r20513; double r20515 = r20512 / r20514; double r20516 = r20511 - r20515; double r20517 = 0.125; double r20518 = 1; double r20519 = v; double r20520 = r20512 * r20519; double r20521 = r20518 + r20520; double r20522 = r20517 * r20521; double r20523 = r20522 / r20518; double r20524 = w; double r20525 = r20524 * r20524; double r20526 = r20525 * r20513; double r20527 = r20526 * r20513; double r20528 = r20518 - r20519; double r20529 = r20527 / r20528; double r20530 = r20523 * r20529; double r20531 = r20516 - r20530; double r20532 = 0.5; double r20533 = r20531 - r20532; return r20533; } 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 r20534, r20535, r20536, r20537, r20538, r20539, r20540, r20541, r20542, r20543, r20544, r20545, r20546, r20547, r20548, r20549, r20550, r20551, r20552, r20553, r20554, r20555; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20534, "3", 10, MPFR_RNDN); mpfr_init_set_str(r20535, "2", 10, MPFR_RNDN); mpfr_init(r20536); mpfr_init(r20537); mpfr_init(r20538); mpfr_init(r20539); mpfr_init_set_str(r20540, "0.125", 10, MPFR_RNDN); mpfr_init_set_str(r20541, "1", 10, MPFR_RNDN); mpfr_init(r20542); mpfr_init(r20543); mpfr_init(r20544); mpfr_init(r20545); mpfr_init(r20546); mpfr_init(r20547); mpfr_init(r20548); mpfr_init(r20549); mpfr_init(r20550); mpfr_init(r20551); mpfr_init(r20552); mpfr_init(r20553); mpfr_init_set_str(r20554, "0.5", 10, MPFR_RNDN); mpfr_init(r20555); } double f_im(double v, double w, double r) { ; ; mpfr_set_d(r20536, r, MPFR_RNDN); mpfr_mul(r20537, r20536, r20536, MPFR_RNDN); mpfr_div(r20538, r20535, r20537, MPFR_RNDN); mpfr_sub(r20539, r20534, r20538, MPFR_RNDN); ; ; mpfr_set_d(r20542, v, MPFR_RNDN); mpfr_mul(r20543, r20535, r20542, MPFR_RNDN); mpfr_add(r20544, r20541, r20543, MPFR_RNDN); mpfr_mul(r20545, r20540, r20544, MPFR_RNDN); mpfr_set_d(r20546, w, MPFR_RNDN); mpfr_mul(r20547, r20546, r20546, MPFR_RNDN); mpfr_mul(r20548, r20547, r20536, MPFR_RNDN); mpfr_mul(r20549, r20548, r20536, MPFR_RNDN); mpfr_mul(r20550, r20545, r20549, MPFR_RNDN); mpfr_sub(r20551, r20541, r20542, MPFR_RNDN); mpfr_div(r20552, r20550, r20551, MPFR_RNDN); mpfr_sub(r20553, r20539, r20552, MPFR_RNDN); ; mpfr_sub(r20555, r20553, r20554, MPFR_RNDN); return mpfr_get_d(r20555, MPFR_RNDN); } static mpfr_t r20556, r20557, r20558, r20559, r20560, r20561, r20562, r20563, r20564, r20565, r20566, r20567, r20568, r20569, r20570, r20571, r20572, r20573, r20574, r20575, r20576, r20577, r20578; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20556, "3", 10, MPFR_RNDN); mpfr_init_set_str(r20557, "2", 10, MPFR_RNDN); mpfr_init(r20558); mpfr_init(r20559); mpfr_init(r20560); mpfr_init(r20561); mpfr_init_set_str(r20562, "0.125", 10, MPFR_RNDN); mpfr_init_set_str(r20563, "1", 10, MPFR_RNDN); mpfr_init(r20564); mpfr_init(r20565); mpfr_init(r20566); mpfr_init(r20567); mpfr_init(r20568); mpfr_init(r20569); mpfr_init(r20570); mpfr_init(r20571); mpfr_init(r20572); mpfr_init(r20573); mpfr_init(r20574); mpfr_init(r20575); mpfr_init(r20576); mpfr_init_set_str(r20577, "0.5", 10, MPFR_RNDN); mpfr_init(r20578); } double f_fm(double v, double w, double r) { ; ; mpfr_set_d(r20558, r, MPFR_RNDN); mpfr_mul(r20559, r20558, r20558, MPFR_RNDN); mpfr_div(r20560, r20557, r20559, MPFR_RNDN); mpfr_sub(r20561, r20556, r20560, MPFR_RNDN); ; ; mpfr_set_d(r20564, v, MPFR_RNDN); mpfr_mul(r20565, r20557, r20564, MPFR_RNDN); mpfr_add(r20566, r20563, r20565, MPFR_RNDN); mpfr_mul(r20567, r20562, r20566, MPFR_RNDN); mpfr_div(r20568, r20567, r20563, MPFR_RNDN); mpfr_set_d(r20569, w, MPFR_RNDN); mpfr_mul(r20570, r20569, r20569, MPFR_RNDN); mpfr_mul(r20571, r20570, r20558, MPFR_RNDN); mpfr_mul(r20572, r20571, r20558, MPFR_RNDN); mpfr_sub(r20573, r20563, r20564, MPFR_RNDN); mpfr_div(r20574, r20572, r20573, MPFR_RNDN); mpfr_mul(r20575, r20568, r20574, MPFR_RNDN); mpfr_sub(r20576, r20561, r20575, MPFR_RNDN); ; mpfr_sub(r20578, r20576, r20577, MPFR_RNDN); return mpfr_get_d(r20578, MPFR_RNDN); } static mpfr_t r20579, r20580, r20581, r20582, r20583, r20584, r20585, r20586, r20587, r20588, r20589, r20590, r20591, r20592, r20593, r20594, r20595, r20596, r20597, r20598, r20599, r20600, r20601; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20579, "3", 10, MPFR_RNDN); mpfr_init_set_str(r20580, "2", 10, MPFR_RNDN); mpfr_init(r20581); mpfr_init(r20582); mpfr_init(r20583); mpfr_init(r20584); mpfr_init_set_str(r20585, "0.125", 10, MPFR_RNDN); mpfr_init_set_str(r20586, "1", 10, MPFR_RNDN); mpfr_init(r20587); mpfr_init(r20588); mpfr_init(r20589); mpfr_init(r20590); mpfr_init(r20591); mpfr_init(r20592); mpfr_init(r20593); mpfr_init(r20594); mpfr_init(r20595); mpfr_init(r20596); mpfr_init(r20597); mpfr_init(r20598); mpfr_init(r20599); mpfr_init_set_str(r20600, "0.5", 10, MPFR_RNDN); mpfr_init(r20601); } double f_dm(double v, double w, double r) { ; ; mpfr_set_d(r20581, r, MPFR_RNDN); mpfr_mul(r20582, r20581, r20581, MPFR_RNDN); mpfr_div(r20583, r20580, r20582, MPFR_RNDN); mpfr_sub(r20584, r20579, r20583, MPFR_RNDN); ; ; mpfr_set_d(r20587, v, MPFR_RNDN); mpfr_mul(r20588, r20580, r20587, MPFR_RNDN); mpfr_add(r20589, r20586, r20588, MPFR_RNDN); mpfr_mul(r20590, r20585, r20589, MPFR_RNDN); mpfr_div(r20591, r20590, r20586, MPFR_RNDN); mpfr_set_d(r20592, w, MPFR_RNDN); mpfr_mul(r20593, r20592, r20592, MPFR_RNDN); mpfr_mul(r20594, r20593, r20581, MPFR_RNDN); mpfr_mul(r20595, r20594, r20581, MPFR_RNDN); mpfr_sub(r20596, r20586, r20587, MPFR_RNDN); mpfr_div(r20597, r20595, r20596, MPFR_RNDN); mpfr_mul(r20598, r20591, r20597, MPFR_RNDN); mpfr_sub(r20599, r20584, r20598, MPFR_RNDN); ; mpfr_sub(r20601, r20599, r20600, MPFR_RNDN); return mpfr_get_d(r20601, MPFR_RNDN); }