#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r25292 = x; float r25293 = y; float r25294 = z; float r25295 = fma(r25292, r25293, r25294); float r25296 = 1; float r25297 = r25292 * r25293; float r25298 = r25297 + r25294; float r25299 = r25296 + r25298; float r25300 = r25295 - r25299; return r25300; } double f_id(double x, double y, double z) { double r25301 = x; double r25302 = y; double r25303 = z; double r25304 = fma(r25301, r25302, r25303); double r25305 = 1; double r25306 = r25301 * r25302; double r25307 = r25306 + r25303; double r25308 = r25305 + r25307; double r25309 = r25304 - r25308; return r25309; } double f_of(float x, float y, float z) { float r25310 = x; float r25311 = y; float r25312 = z; float r25313 = fma(r25310, r25311, r25312); float r25314 = r25311 * r25310; float r25315 = r25313 - r25314; float r25316 = r25315 - r25312; float r25317 = -2.919718334020837e-16; bool r25318 = r25316 <= r25317; float r25319 = 0.0; bool r25320 = r25316 <= r25319; float r25321 = !r25320; bool r25322 = r25318 || r25321; float r25323 = r25313 - r25312; float r25324 = r25323 - r25314; float r25325 = 1; float r25326 = r25324 - r25325; float r25327 = exp(r25326); float r25328 = log(r25327); float r25329 = exp(r25328); float r25330 = log(r25329); float r25331 = r25316 - r25325; float r25332 = exp(r25331); float r25333 = log(r25332); float r25334 = r25322 ? r25330 : r25333; return r25334; } double f_od(double x, double y, double z) { double r25335 = x; double r25336 = y; double r25337 = z; double r25338 = fma(r25335, r25336, r25337); double r25339 = r25336 * r25335; double r25340 = r25338 - r25339; double r25341 = r25340 - r25337; double r25342 = -2.919718334020837e-16; bool r25343 = r25341 <= r25342; double r25344 = 0.0; bool r25345 = r25341 <= r25344; double r25346 = !r25345; bool r25347 = r25343 || r25346; double r25348 = r25338 - r25337; double r25349 = r25348 - r25339; double r25350 = 1; double r25351 = r25349 - r25350; double r25352 = exp(r25351); double r25353 = log(r25352); double r25354 = exp(r25353); double r25355 = log(r25354); double r25356 = r25341 - r25350; double r25357 = exp(r25356); double r25358 = log(r25357); double r25359 = r25347 ? r25355 : r25358; return r25359; } 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 r25360, r25361, r25362, r25363, r25364, r25365, r25366, r25367, r25368; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r25360); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init_set_str(r25364, "1", 10, MPFR_RNDN); mpfr_init(r25365); mpfr_init(r25366); mpfr_init(r25367); mpfr_init(r25368); } double f_im(double x, double y, double z) { mpfr_set_d(r25360, x, MPFR_RNDN); mpfr_set_d(r25361, y, MPFR_RNDN); mpfr_set_d(r25362, z, MPFR_RNDN); mpfr_fma(r25363, r25360, r25361, r25362, MPFR_RNDN); ; mpfr_mul(r25365, r25360, r25361, MPFR_RNDN); mpfr_add(r25366, r25365, r25362, MPFR_RNDN); mpfr_add(r25367, r25364, r25366, MPFR_RNDN); mpfr_sub(r25368, r25363, r25367, MPFR_RNDN); return mpfr_get_d(r25368, MPFR_RNDN); } static mpfr_t r25369, r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391, r25392, r25393; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r25369); mpfr_init(r25370); mpfr_init(r25371); mpfr_init(r25372); mpfr_init(r25373); mpfr_init(r25374); mpfr_init(r25375); mpfr_init_set_str(r25376, "-2.919718334020837e-16", 10, MPFR_RNDN); mpfr_init(r25377); mpfr_init_set_str(r25378, "0.0", 10, MPFR_RNDN); mpfr_init(r25379); mpfr_init(r25380); mpfr_init(r25381); mpfr_init(r25382); mpfr_init(r25383); mpfr_init_set_str(r25384, "1", 10, MPFR_RNDN); mpfr_init(r25385); mpfr_init(r25386); mpfr_init(r25387); mpfr_init(r25388); mpfr_init(r25389); mpfr_init(r25390); mpfr_init(r25391); mpfr_init(r25392); mpfr_init(r25393); } double f_fm(double x, double y, double z) { mpfr_set_d(r25369, x, MPFR_RNDN); mpfr_set_d(r25370, y, MPFR_RNDN); mpfr_set_d(r25371, z, MPFR_RNDN); mpfr_fma(r25372, r25369, r25370, r25371, MPFR_RNDN); mpfr_mul(r25373, r25370, r25369, MPFR_RNDN); mpfr_sub(r25374, r25372, r25373, MPFR_RNDN); mpfr_sub(r25375, r25374, r25371, MPFR_RNDN); ; mpfr_set_si(r25377, mpfr_cmp(r25375, r25376) <= 0, MPFR_RNDN); ; mpfr_set_si(r25379, mpfr_cmp(r25375, r25378) <= 0, MPFR_RNDN); mpfr_set_si(r25380, !mpfr_get_si(r25379, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25381, mpfr_get_si(r25377, MPFR_RNDN) || mpfr_get_si(r25380, MPFR_RNDN), MPFR_RNDN); mpfr_sub(r25382, r25372, r25371, MPFR_RNDN); mpfr_sub(r25383, r25382, r25373, MPFR_RNDN); ; mpfr_sub(r25385, r25383, r25384, MPFR_RNDN); mpfr_exp(r25386, r25385, MPFR_RNDN); mpfr_log(r25387, r25386, MPFR_RNDN); mpfr_exp(r25388, r25387, MPFR_RNDN); mpfr_log(r25389, r25388, MPFR_RNDN); mpfr_sub(r25390, r25375, r25384, MPFR_RNDN); mpfr_exp(r25391, r25390, MPFR_RNDN); mpfr_log(r25392, r25391, MPFR_RNDN); if (mpfr_get_si(r25381, MPFR_RNDN)) { mpfr_set(r25393, r25389, MPFR_RNDN); } else { mpfr_set(r25393, r25392, MPFR_RNDN); }; return mpfr_get_d(r25393, MPFR_RNDN); } static mpfr_t r25394, r25395, r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r25394); mpfr_init(r25395); mpfr_init(r25396); mpfr_init(r25397); mpfr_init(r25398); mpfr_init(r25399); mpfr_init(r25400); mpfr_init_set_str(r25401, "-2.919718334020837e-16", 10, MPFR_RNDN); mpfr_init(r25402); mpfr_init_set_str(r25403, "0.0", 10, MPFR_RNDN); mpfr_init(r25404); mpfr_init(r25405); mpfr_init(r25406); mpfr_init(r25407); mpfr_init(r25408); mpfr_init_set_str(r25409, "1", 10, MPFR_RNDN); mpfr_init(r25410); mpfr_init(r25411); mpfr_init(r25412); mpfr_init(r25413); mpfr_init(r25414); mpfr_init(r25415); mpfr_init(r25416); mpfr_init(r25417); mpfr_init(r25418); } double f_dm(double x, double y, double z) { mpfr_set_d(r25394, x, MPFR_RNDN); mpfr_set_d(r25395, y, MPFR_RNDN); mpfr_set_d(r25396, z, MPFR_RNDN); mpfr_fma(r25397, r25394, r25395, r25396, MPFR_RNDN); mpfr_mul(r25398, r25395, r25394, MPFR_RNDN); mpfr_sub(r25399, r25397, r25398, MPFR_RNDN); mpfr_sub(r25400, r25399, r25396, MPFR_RNDN); ; mpfr_set_si(r25402, mpfr_cmp(r25400, r25401) <= 0, MPFR_RNDN); ; mpfr_set_si(r25404, mpfr_cmp(r25400, r25403) <= 0, MPFR_RNDN); mpfr_set_si(r25405, !mpfr_get_si(r25404, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25406, mpfr_get_si(r25402, MPFR_RNDN) || mpfr_get_si(r25405, MPFR_RNDN), MPFR_RNDN); mpfr_sub(r25407, r25397, r25396, MPFR_RNDN); mpfr_sub(r25408, r25407, r25398, MPFR_RNDN); ; mpfr_sub(r25410, r25408, r25409, MPFR_RNDN); mpfr_exp(r25411, r25410, MPFR_RNDN); mpfr_log(r25412, r25411, MPFR_RNDN); mpfr_exp(r25413, r25412, MPFR_RNDN); mpfr_log(r25414, r25413, MPFR_RNDN); mpfr_sub(r25415, r25400, r25409, MPFR_RNDN); mpfr_exp(r25416, r25415, MPFR_RNDN); mpfr_log(r25417, r25416, MPFR_RNDN); if (mpfr_get_si(r25406, MPFR_RNDN)) { mpfr_set(r25418, r25414, MPFR_RNDN); } else { mpfr_set(r25418, r25417, MPFR_RNDN); }; return mpfr_get_d(r25418, MPFR_RNDN); }