#include #include #include #include #include char *name = "Rosa's DopplerBench"; double f_if(float u, float v, float t1) { float r24388 = t1; float r24389 = -r24388; float r24390 = v; float r24391 = r24389 * r24390; float r24392 = u; float r24393 = r24388 + r24392; float r24394 = r24393 * r24393; float r24395 = r24391 / r24394; return r24395; } double f_id(double u, double v, double t1) { double r24396 = t1; double r24397 = -r24396; double r24398 = v; double r24399 = r24397 * r24398; double r24400 = u; double r24401 = r24396 + r24400; double r24402 = r24401 * r24401; double r24403 = r24399 / r24402; return r24403; } double f_of(float u, float v, float t1) { float r24404 = t1; float r24405 = -4.47775029934645e-194; bool r24406 = r24404 <= r24405; float r24407 = 8.543332634879257e-235; bool r24408 = r24404 <= r24407; float r24409 = !r24408; bool r24410 = r24406 || r24409; float r24411 = u; float r24412 = r24411 + r24404; float r24413 = r24404 / r24412; float r24414 = v; float r24415 = -r24414; float r24416 = r24415 / r24412; float r24417 = r24413 * r24416; float r24418 = -r24404; float r24419 = r24418 * r24414; float r24420 = r24419 / r24412; float r24421 = r24420 / r24412; float r24422 = r24410 ? r24417 : r24421; return r24422; } double f_od(double u, double v, double t1) { double r24423 = t1; double r24424 = -4.47775029934645e-194; bool r24425 = r24423 <= r24424; double r24426 = 8.543332634879257e-235; bool r24427 = r24423 <= r24426; double r24428 = !r24427; bool r24429 = r24425 || r24428; double r24430 = u; double r24431 = r24430 + r24423; double r24432 = r24423 / r24431; double r24433 = v; double r24434 = -r24433; double r24435 = r24434 / r24431; double r24436 = r24432 * r24435; double r24437 = -r24423; double r24438 = r24437 * r24433; double r24439 = r24438 / r24431; double r24440 = r24439 / r24431; double r24441 = r24429 ? r24436 : r24440; return r24441; } 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 r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); } double f_im(double u, double v, double t1) { mpfr_set_d(r24442, t1, MPFR_RNDN); mpfr_neg(r24443, r24442, MPFR_RNDN); mpfr_set_d(r24444, v, MPFR_RNDN); mpfr_mul(r24445, r24443, r24444, MPFR_RNDN); mpfr_set_d(r24446, u, MPFR_RNDN); mpfr_add(r24447, r24442, r24446, MPFR_RNDN); mpfr_mul(r24448, r24447, r24447, MPFR_RNDN); mpfr_div(r24449, r24445, r24448, MPFR_RNDN); return mpfr_get_d(r24449, MPFR_RNDN); } static mpfr_t r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24450); mpfr_init_set_str(r24451, "-4.47775029934645e-194", 10, MPFR_RNDN); mpfr_init(r24452); mpfr_init_set_str(r24453, "8.543332634879257e-235", 10, MPFR_RNDN); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); mpfr_init(r24464); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); } double f_fm(double u, double v, double t1) { mpfr_set_d(r24450, t1, MPFR_RNDN); ; mpfr_set_si(r24452, mpfr_cmp(r24450, r24451) <= 0, MPFR_RNDN); ; mpfr_set_si(r24454, mpfr_cmp(r24450, r24453) <= 0, MPFR_RNDN); mpfr_set_si(r24455, !mpfr_get_si(r24454, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r24456, mpfr_get_si(r24452, MPFR_RNDN) || mpfr_get_si(r24455, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r24457, u, MPFR_RNDN); mpfr_add(r24458, r24457, r24450, MPFR_RNDN); mpfr_div(r24459, r24450, r24458, MPFR_RNDN); mpfr_set_d(r24460, v, MPFR_RNDN); mpfr_neg(r24461, r24460, MPFR_RNDN); mpfr_div(r24462, r24461, r24458, MPFR_RNDN); mpfr_mul(r24463, r24459, r24462, MPFR_RNDN); mpfr_neg(r24464, r24450, MPFR_RNDN); mpfr_mul(r24465, r24464, r24460, MPFR_RNDN); mpfr_div(r24466, r24465, r24458, MPFR_RNDN); mpfr_div(r24467, r24466, r24458, MPFR_RNDN); if (mpfr_get_si(r24456, MPFR_RNDN)) { mpfr_set(r24468, r24463, MPFR_RNDN); } else { mpfr_set(r24468, r24467, MPFR_RNDN); }; return mpfr_get_d(r24468, MPFR_RNDN); } static mpfr_t r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24469); mpfr_init_set_str(r24470, "-4.47775029934645e-194", 10, MPFR_RNDN); mpfr_init(r24471); mpfr_init_set_str(r24472, "8.543332634879257e-235", 10, MPFR_RNDN); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); mpfr_init(r24481); mpfr_init(r24482); mpfr_init(r24483); mpfr_init(r24484); mpfr_init(r24485); mpfr_init(r24486); mpfr_init(r24487); } double f_dm(double u, double v, double t1) { mpfr_set_d(r24469, t1, MPFR_RNDN); ; mpfr_set_si(r24471, mpfr_cmp(r24469, r24470) <= 0, MPFR_RNDN); ; mpfr_set_si(r24473, mpfr_cmp(r24469, r24472) <= 0, MPFR_RNDN); mpfr_set_si(r24474, !mpfr_get_si(r24473, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r24475, mpfr_get_si(r24471, MPFR_RNDN) || mpfr_get_si(r24474, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r24476, u, MPFR_RNDN); mpfr_add(r24477, r24476, r24469, MPFR_RNDN); mpfr_div(r24478, r24469, r24477, MPFR_RNDN); mpfr_set_d(r24479, v, MPFR_RNDN); mpfr_neg(r24480, r24479, MPFR_RNDN); mpfr_div(r24481, r24480, r24477, MPFR_RNDN); mpfr_mul(r24482, r24478, r24481, MPFR_RNDN); mpfr_neg(r24483, r24469, MPFR_RNDN); mpfr_mul(r24484, r24483, r24479, MPFR_RNDN); mpfr_div(r24485, r24484, r24477, MPFR_RNDN); mpfr_div(r24486, r24485, r24477, MPFR_RNDN); if (mpfr_get_si(r24475, MPFR_RNDN)) { mpfr_set(r24487, r24482, MPFR_RNDN); } else { mpfr_set(r24487, r24486, MPFR_RNDN); }; return mpfr_get_d(r24487, MPFR_RNDN); }