#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r22401 = x; float r22402 = y; float r22403 = z; float r22404 = fma(r22401, r22402, r22403); float r22405 = 1; float r22406 = r22401 * r22402; float r22407 = r22406 + r22403; float r22408 = r22405 + r22407; float r22409 = r22404 - r22408; return r22409; } double f_id(double x, double y, double z) { double r22410 = x; double r22411 = y; double r22412 = z; double r22413 = fma(r22410, r22411, r22412); double r22414 = 1; double r22415 = r22410 * r22411; double r22416 = r22415 + r22412; double r22417 = r22414 + r22416; double r22418 = r22413 - r22417; return r22418; } double f_of(float x, float y, float z) { float r22419 = x; float r22420 = y; float r22421 = z; float r22422 = fma(r22419, r22420, r22421); float r22423 = r22422 - r22421; float r22424 = r22420 * r22419; float r22425 = r22423 - r22424; float r22426 = 1; float r22427 = r22425 - r22426; float r22428 = exp(r22427); float r22429 = log(r22428); float r22430 = exp(r22429); float r22431 = log(r22430); float r22432 = -1.0000000000293188; bool r22433 = r22431 <= r22432; float r22434 = r22422 - r22424; float r22435 = r22434 - r22421; float r22436 = r22435 - r22426; float r22437 = exp(r22436); float r22438 = log(r22437); float r22439 = -1.0; bool r22440 = r22431 <= r22439; float r22441 = r22440 ? r22431 : r22438; float r22442 = r22433 ? r22438 : r22441; return r22442; } double f_od(double x, double y, double z) { double r22443 = x; double r22444 = y; double r22445 = z; double r22446 = fma(r22443, r22444, r22445); double r22447 = r22446 - r22445; double r22448 = r22444 * r22443; double r22449 = r22447 - r22448; double r22450 = 1; double r22451 = r22449 - r22450; double r22452 = exp(r22451); double r22453 = log(r22452); double r22454 = exp(r22453); double r22455 = log(r22454); double r22456 = -1.0000000000293188; bool r22457 = r22455 <= r22456; double r22458 = r22446 - r22448; double r22459 = r22458 - r22445; double r22460 = r22459 - r22450; double r22461 = exp(r22460); double r22462 = log(r22461); double r22463 = -1.0; bool r22464 = r22455 <= r22463; double r22465 = r22464 ? r22455 : r22462; double r22466 = r22457 ? r22462 : r22465; return r22466; } 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 r22467, r22468, r22469, r22470, r22471, r22472, r22473, r22474, r22475; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r22467); mpfr_init(r22468); mpfr_init(r22469); mpfr_init(r22470); mpfr_init_set_str(r22471, "1", 10, MPFR_RNDN); mpfr_init(r22472); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); } double f_im(double x, double y, double z) { mpfr_set_d(r22467, x, MPFR_RNDN); mpfr_set_d(r22468, y, MPFR_RNDN); mpfr_set_d(r22469, z, MPFR_RNDN); mpfr_fma(r22470, r22467, r22468, r22469, MPFR_RNDN); ; mpfr_mul(r22472, r22467, r22468, MPFR_RNDN); mpfr_add(r22473, r22472, r22469, MPFR_RNDN); mpfr_add(r22474, r22471, r22473, MPFR_RNDN); mpfr_sub(r22475, r22470, r22474, MPFR_RNDN); return mpfr_get_d(r22475, MPFR_RNDN); } static mpfr_t r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r22476); mpfr_init(r22477); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init_set_str(r22483, "1", 10, MPFR_RNDN); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init_set_str(r22489, "-1.0000000000293188", 10, MPFR_RNDN); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); mpfr_init(r22494); mpfr_init(r22495); mpfr_init_set_str(r22496, "-1.0", 10, MPFR_RNDN); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); } double f_fm(double x, double y, double z) { mpfr_set_d(r22476, x, MPFR_RNDN); mpfr_set_d(r22477, y, MPFR_RNDN); mpfr_set_d(r22478, z, MPFR_RNDN); mpfr_fma(r22479, r22476, r22477, r22478, MPFR_RNDN); mpfr_sub(r22480, r22479, r22478, MPFR_RNDN); mpfr_mul(r22481, r22477, r22476, MPFR_RNDN); mpfr_sub(r22482, r22480, r22481, MPFR_RNDN); ; mpfr_sub(r22484, r22482, r22483, MPFR_RNDN); mpfr_exp(r22485, r22484, MPFR_RNDN); mpfr_log(r22486, r22485, MPFR_RNDN); mpfr_exp(r22487, r22486, MPFR_RNDN); mpfr_log(r22488, r22487, MPFR_RNDN); ; mpfr_set_si(r22490, mpfr_cmp(r22488, r22489) <= 0, MPFR_RNDN); mpfr_sub(r22491, r22479, r22481, MPFR_RNDN); mpfr_sub(r22492, r22491, r22478, MPFR_RNDN); mpfr_sub(r22493, r22492, r22483, MPFR_RNDN); mpfr_exp(r22494, r22493, MPFR_RNDN); mpfr_log(r22495, r22494, MPFR_RNDN); ; mpfr_set_si(r22497, mpfr_cmp(r22488, r22496) <= 0, MPFR_RNDN); if (mpfr_get_si(r22497, MPFR_RNDN)) { mpfr_set(r22498, r22488, MPFR_RNDN); } else { mpfr_set(r22498, r22495, MPFR_RNDN); }; if (mpfr_get_si(r22490, MPFR_RNDN)) { mpfr_set(r22499, r22495, MPFR_RNDN); } else { mpfr_set(r22499, r22498, MPFR_RNDN); }; return mpfr_get_d(r22499, MPFR_RNDN); } static mpfr_t r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519, r22520, r22521, r22522, r22523; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init_set_str(r22507, "1", 10, MPFR_RNDN); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init_set_str(r22513, "-1.0000000000293188", 10, MPFR_RNDN); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); mpfr_init_set_str(r22520, "-1.0", 10, MPFR_RNDN); mpfr_init(r22521); mpfr_init(r22522); mpfr_init(r22523); } double f_dm(double x, double y, double z) { mpfr_set_d(r22500, x, MPFR_RNDN); mpfr_set_d(r22501, y, MPFR_RNDN); mpfr_set_d(r22502, z, MPFR_RNDN); mpfr_fma(r22503, r22500, r22501, r22502, MPFR_RNDN); mpfr_sub(r22504, r22503, r22502, MPFR_RNDN); mpfr_mul(r22505, r22501, r22500, MPFR_RNDN); mpfr_sub(r22506, r22504, r22505, MPFR_RNDN); ; mpfr_sub(r22508, r22506, r22507, MPFR_RNDN); mpfr_exp(r22509, r22508, MPFR_RNDN); mpfr_log(r22510, r22509, MPFR_RNDN); mpfr_exp(r22511, r22510, MPFR_RNDN); mpfr_log(r22512, r22511, MPFR_RNDN); ; mpfr_set_si(r22514, mpfr_cmp(r22512, r22513) <= 0, MPFR_RNDN); mpfr_sub(r22515, r22503, r22505, MPFR_RNDN); mpfr_sub(r22516, r22515, r22502, MPFR_RNDN); mpfr_sub(r22517, r22516, r22507, MPFR_RNDN); mpfr_exp(r22518, r22517, MPFR_RNDN); mpfr_log(r22519, r22518, MPFR_RNDN); ; mpfr_set_si(r22521, mpfr_cmp(r22512, r22520) <= 0, MPFR_RNDN); if (mpfr_get_si(r22521, MPFR_RNDN)) { mpfr_set(r22522, r22512, MPFR_RNDN); } else { mpfr_set(r22522, r22519, MPFR_RNDN); }; if (mpfr_get_si(r22514, MPFR_RNDN)) { mpfr_set(r22523, r22519, MPFR_RNDN); } else { mpfr_set(r22523, r22522, MPFR_RNDN); }; return mpfr_get_d(r22523, MPFR_RNDN); }