#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r25245 = x; float r25246 = y; float r25247 = z; float r25248 = fma(r25245, r25246, r25247); float r25249 = 1; float r25250 = r25245 * r25246; float r25251 = r25250 + r25247; float r25252 = r25249 + r25251; float r25253 = r25248 - r25252; return r25253; } double f_id(double x, double y, double z) { double r25254 = x; double r25255 = y; double r25256 = z; double r25257 = fma(r25254, r25255, r25256); double r25258 = 1; double r25259 = r25254 * r25255; double r25260 = r25259 + r25256; double r25261 = r25258 + r25260; double r25262 = r25257 - r25261; return r25262; } double f_of(float x, float y, float z) { float r25263 = x; float r25264 = y; float r25265 = z; float r25266 = fma(r25263, r25264, r25265); float r25267 = r25266 - r25265; float r25268 = r25264 * r25263; float r25269 = r25267 - r25268; float r25270 = 1; float r25271 = r25269 - r25270; float r25272 = exp(r25271); float r25273 = log(r25272); float r25274 = exp(r25273); float r25275 = log(r25274); float r25276 = -1.015625; bool r25277 = r25275 <= r25276; float r25278 = r25266 - r25268; float r25279 = r25278 - r25265; float r25280 = r25279 - r25270; float r25281 = exp(r25280); float r25282 = log(r25281); float r25283 = -0.9999999999999453; bool r25284 = r25275 <= r25283; float r25285 = r25284 ? r25275 : r25282; float r25286 = r25277 ? r25282 : r25285; return r25286; } double f_od(double x, double y, double z) { double r25287 = x; double r25288 = y; double r25289 = z; double r25290 = fma(r25287, r25288, r25289); double r25291 = r25290 - r25289; double r25292 = r25288 * r25287; double r25293 = r25291 - r25292; double r25294 = 1; double r25295 = r25293 - r25294; double r25296 = exp(r25295); double r25297 = log(r25296); double r25298 = exp(r25297); double r25299 = log(r25298); double r25300 = -1.015625; bool r25301 = r25299 <= r25300; double r25302 = r25290 - r25292; double r25303 = r25302 - r25289; double r25304 = r25303 - r25294; double r25305 = exp(r25304); double r25306 = log(r25305); double r25307 = -0.9999999999999453; bool r25308 = r25299 <= r25307; double r25309 = r25308 ? r25299 : r25306; double r25310 = r25301 ? r25306 : r25309; return r25310; } 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 r25311, r25312, r25313, r25314, r25315, r25316, r25317, r25318, r25319; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r25311); mpfr_init(r25312); mpfr_init(r25313); mpfr_init(r25314); mpfr_init_set_str(r25315, "1", 10, MPFR_RNDN); mpfr_init(r25316); mpfr_init(r25317); mpfr_init(r25318); mpfr_init(r25319); } double f_im(double x, double y, double z) { mpfr_set_d(r25311, x, MPFR_RNDN); mpfr_set_d(r25312, y, MPFR_RNDN); mpfr_set_d(r25313, z, MPFR_RNDN); mpfr_fma(r25314, r25311, r25312, r25313, MPFR_RNDN); ; mpfr_mul(r25316, r25311, r25312, MPFR_RNDN); mpfr_add(r25317, r25316, r25313, MPFR_RNDN); mpfr_add(r25318, r25315, r25317, MPFR_RNDN); mpfr_sub(r25319, r25314, r25318, MPFR_RNDN); return mpfr_get_d(r25319, MPFR_RNDN); } static mpfr_t r25320, r25321, r25322, r25323, r25324, r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340, r25341, r25342, r25343; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r25320); mpfr_init(r25321); mpfr_init(r25322); mpfr_init(r25323); mpfr_init(r25324); mpfr_init(r25325); mpfr_init(r25326); mpfr_init_set_str(r25327, "1", 10, MPFR_RNDN); mpfr_init(r25328); mpfr_init(r25329); mpfr_init(r25330); mpfr_init(r25331); mpfr_init(r25332); mpfr_init_set_str(r25333, "-1.015625", 10, MPFR_RNDN); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init(r25337); mpfr_init(r25338); mpfr_init(r25339); mpfr_init_set_str(r25340, "-0.9999999999999453", 10, MPFR_RNDN); mpfr_init(r25341); mpfr_init(r25342); mpfr_init(r25343); } double f_fm(double x, double y, double z) { mpfr_set_d(r25320, x, MPFR_RNDN); mpfr_set_d(r25321, y, MPFR_RNDN); mpfr_set_d(r25322, z, MPFR_RNDN); mpfr_fma(r25323, r25320, r25321, r25322, MPFR_RNDN); mpfr_sub(r25324, r25323, r25322, MPFR_RNDN); mpfr_mul(r25325, r25321, r25320, MPFR_RNDN); mpfr_sub(r25326, r25324, r25325, MPFR_RNDN); ; mpfr_sub(r25328, r25326, r25327, MPFR_RNDN); mpfr_exp(r25329, r25328, MPFR_RNDN); mpfr_log(r25330, r25329, MPFR_RNDN); mpfr_exp(r25331, r25330, MPFR_RNDN); mpfr_log(r25332, r25331, MPFR_RNDN); ; mpfr_set_si(r25334, mpfr_cmp(r25332, r25333) <= 0, MPFR_RNDN); mpfr_sub(r25335, r25323, r25325, MPFR_RNDN); mpfr_sub(r25336, r25335, r25322, MPFR_RNDN); mpfr_sub(r25337, r25336, r25327, MPFR_RNDN); mpfr_exp(r25338, r25337, MPFR_RNDN); mpfr_log(r25339, r25338, MPFR_RNDN); ; mpfr_set_si(r25341, mpfr_cmp(r25332, r25340) <= 0, MPFR_RNDN); if (mpfr_get_si(r25341, MPFR_RNDN)) { mpfr_set(r25342, r25332, MPFR_RNDN); } else { mpfr_set(r25342, r25339, MPFR_RNDN); }; if (mpfr_get_si(r25334, MPFR_RNDN)) { mpfr_set(r25343, r25339, MPFR_RNDN); } else { mpfr_set(r25343, r25342, MPFR_RNDN); }; return mpfr_get_d(r25343, MPFR_RNDN); } static mpfr_t r25344, r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364, r25365, r25366, r25367; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r25344); mpfr_init(r25345); mpfr_init(r25346); mpfr_init(r25347); mpfr_init(r25348); mpfr_init(r25349); mpfr_init(r25350); mpfr_init_set_str(r25351, "1", 10, MPFR_RNDN); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init(r25356); mpfr_init_set_str(r25357, "-1.015625", 10, MPFR_RNDN); mpfr_init(r25358); mpfr_init(r25359); mpfr_init(r25360); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init_set_str(r25364, "-0.9999999999999453", 10, MPFR_RNDN); mpfr_init(r25365); mpfr_init(r25366); mpfr_init(r25367); } double f_dm(double x, double y, double z) { mpfr_set_d(r25344, x, MPFR_RNDN); mpfr_set_d(r25345, y, MPFR_RNDN); mpfr_set_d(r25346, z, MPFR_RNDN); mpfr_fma(r25347, r25344, r25345, r25346, MPFR_RNDN); mpfr_sub(r25348, r25347, r25346, MPFR_RNDN); mpfr_mul(r25349, r25345, r25344, MPFR_RNDN); mpfr_sub(r25350, r25348, r25349, MPFR_RNDN); ; mpfr_sub(r25352, r25350, r25351, MPFR_RNDN); mpfr_exp(r25353, r25352, MPFR_RNDN); mpfr_log(r25354, r25353, MPFR_RNDN); mpfr_exp(r25355, r25354, MPFR_RNDN); mpfr_log(r25356, r25355, MPFR_RNDN); ; mpfr_set_si(r25358, mpfr_cmp(r25356, r25357) <= 0, MPFR_RNDN); mpfr_sub(r25359, r25347, r25349, MPFR_RNDN); mpfr_sub(r25360, r25359, r25346, MPFR_RNDN); mpfr_sub(r25361, r25360, r25351, MPFR_RNDN); mpfr_exp(r25362, r25361, MPFR_RNDN); mpfr_log(r25363, r25362, MPFR_RNDN); ; mpfr_set_si(r25365, mpfr_cmp(r25356, r25364) <= 0, MPFR_RNDN); if (mpfr_get_si(r25365, MPFR_RNDN)) { mpfr_set(r25366, r25356, MPFR_RNDN); } else { mpfr_set(r25366, r25363, MPFR_RNDN); }; if (mpfr_get_si(r25358, MPFR_RNDN)) { mpfr_set(r25367, r25363, MPFR_RNDN); } else { mpfr_set(r25367, r25366, MPFR_RNDN); }; return mpfr_get_d(r25367, MPFR_RNDN); }